楼主: zdh1999
4373 1

[程序分享] 广义相加模型考虑混杂因素程序书写 [推广有奖]

  • 0关注
  • 0粉丝

等待验证会员

小学生

14%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
30 点
帖子
2
精华
0
在线时间
5 小时
注册时间
2022-2-7
最后登录
2024-3-5

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币

summary(model<-gam(all~lag1pm2.5+ns(day,4*5)+ns(temp,6)+ns(rh,3)+as.factor(dow)+as.factor(vocation),data=sz,family=quasipoisson))

其中混杂因素温度ns(temp,6)和湿度ns(rh,3)中6和3是什么意思,这两个数值是随便界定的嘛

那我还要考虑气压和风速的影响,后面的程序又该如何书写??


二维码

扫码加我 拉你入群

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

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

关键词:Summary model Summa Mary mode

沙发
owenqi 在职认证  学生认证  发表于 2022-2-11 08:05:59 |只看作者 |坛友微信交流群
首先建议你把你要做的东西搞清楚,然后再去套用代码,不然你连你做的什么东西都不知道,那修改起来更加是没有头绪。
这里你的gam里面用的ns是splines包里的自然样条函数(natrual cubic splines)其中里面一个比较重要的参数是自由度(degree of freedom,df)根据是否包括截距,自由度和间隔点(knots)之间有如下关系 knots = df - 1 - intercept,其中默认截距是0,也就是说 ,默认情况下,knots = df -1。
这个数字是怎么界定的当然要看你实际情况,取大取小对样条的拟合程度和计算的复杂程度都会有影响。具体情况请自己参考样条相关的资料。

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-28 16:11