我想拟合一个功能广义加法模型(具有标量结果和功能预测器),其中结果在积分上回归:
\int F(X(t)) dt
默认情况下,R 中的“退款”包采用更灵活的形式
\int F(X(t), t) dt
有没有办法限制函数 F 不依赖于 t 变量?我通过将 t 方向的自由度指定为“pfr”函数中的一个来尝试此操作,但该函数忽略了此约束。这是显示我得到的错误的最小工作示例(取自“af”函数的文档):
library(refund)
data(DTI)
DTI1 <- DTI[DTI$visit==1 & DTI$case==1, ]
DTI2 <- DTI1[complete.cases(DTI1), ]
fit1 <- pfr(pasat ~ af(cca, k=c(8,1), m=list(c(2,2), c(0,0)), basistype='te'), data=DTI2)
vis.pfr(fit1)