我想拟合一个函数广义加法模型(带有标量结果和一个函数预测器),其中结果在积分上回归:
\int F(X(t)) dt
默认情况下,R中的“退款”套餐采用更灵活的格式
\int F(X(t), t) dt
有没有办法约束函数F不依赖于t变量?我通过在“pfr”函数中将t方向的自由度指定为1来尝试此操作,但该函数忽略了此约束。这是显示我得到的错误的最小工作示例(取自“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)