使用“退款”包装拟合功能性广义添加剂模型

时间:2015-12-28 02:47:24

标签: r gam

我想拟合一个函数广义加法模型(带有标量结果和一个函数预测器),其中结果在积分上回归:

\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)

0 个答案:

没有答案