朱耀辉 发表于 2015-8-8 15:32 
错误在n,你的n不是数值型的,你在看看
您好!是这样的,我现在在做关于气象因素与某种疾病的关系的毕业论文,现在遇到一个自由度的问题,就是我想研究两个气象因素(共有四个气象因素)与疾病之间的关系,先要做单个气象因素与疾病的关系,就需要调整其他两个非研究的气象因素(4-2),我是想用NS()函数来调整其他两个气象因素,在调整的过程中需要同时确定这两个非研究气象因素的自由度,因此想请教你具体这两个调整的气象因素的自由度该如何确定,请您指导一下!在此非常感谢!!!具体的formula如下:
model1 <- gam(cases~ave_temp+ns(
ave_pres,df)+ns(
precipitation,df)+ns(day1,2*7)+as.factor(holiday)+as.factor(DOW), data=cs1 , fmily=poisson(link="log"))
需要确定的变量的df为ave_pres,precipitation.