R中mgcv包的gam模型拟合自由度的选择

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

初次使用R，用mgcv包的gam模型来拟合，想自定义自由度来拟合，不知道如何实现？恳请高手帮助，谢谢啦！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：模型拟合 自由度 GAM 如何实现 自定义 自由度 模型 如何

可以自己设定k的值

1. gam(y~s(x,k=3),data=data)

复制代码

具体可参看帮助

1. ?formula.gam

复制代码

请问k值设定有什么标准么？怎么确定k值取多少合适呢？

可以用残差PACF绝对值之和最小来选择。具体就是
mod<-gam(.......)
n<-pacf(mod$residuals,lag.max=30)
m<-sum(abs(n))
测试不同的K值下的m，m有个最小值，对应的K就是最佳自由度

朱耀辉 发表于 2014-9-26 20:15
可以用残差PACF绝对值之和最小来选择。具体就是
mod

请教一下，按照你的提示，测试了一下，但是无法计算PACF之和,R软件提示如下
> m<-sum(abs(n))
Error in abs(n) : non-numeric argument to mathematical function
麻烦再给点建议，不知道问题在哪？？

错误在n，你的n不是数值型的，你在看看

朱耀辉 发表于 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.

朱耀辉 发表于 2015-8-8 15:32
错误在n，你的n不是数值型的，你在看看

class(n)
[1] "acf
as.numeric(n)
Error: (list) object cannot be coerced to type 'double'

请问 n也不能转数值型的，怎么办~~~   多谢

Alwaysonline21 发表于 2015-10-13 22:06
您好！是这样的，我现在在做关于气象因素与某种疾病的关系的毕业论文，现在遇到一个自由度的问题，就是我 ...

你好，你做出来了吗，我现在也在研究类似的课题，能交流一下吗，df怎么控制