我正在R中工作。我想灵活地指定一个函数q(X,Y),其中q相对于X严格递增(但不一定相对于Y)。 我想使用二维样条,但我不知道该怎么做(如何在二维中施加单调性约束)。 到目前为止,我的方法是选择 Y 的多个值,对于 Y 的每个值,我计算相对于 X 的单调样条曲线。然后我进行插值以创建整个函数 q。但我想知道是否可以一步完成所有这些?
scam
软件包具有以下功能:
smooth.construct.tesmi1.smooth.spec
Tensor product smoothing constructor for a
bivariate function monotone increasing in the
first covariate
因此
scam(z ~ s(x, y, bs = "tesmi1"), data = ...)
应该拟合一个惩罚回归样条,该样条在
x
中单调递增,并且相对于 y
不受约束。