| 查看: 2400 | 回复: 0 | |||
[求助]
sas或R如何做空气污染物和入院人数的时间序列分析（广义相加模型GAM)，有程序和数据
|
|
本科毕业论文想做一个关于PM10和每日入院人数的时间序列分析，用的是GAM模型，但在选取模型样条函数的自由度上无从下手，还有控制时间的变量如何加入模型语句中？以及控制星期几效应的哑元变量（dow)如何纳入模型？ 下面为sas初步的程序语句，其中输入的变量分别为：时间（1~730），控制星期几效应的哑元变量dow（1~7），总入院人数n，心血管病入院人数n1，呼吸系统病人数n2，pm10污染物浓度，so2浓度，no2浓度，温度，湿度。 data a; input time dow n n1 n2 pm10 so2 no2 t rh; logn1=log(n1); cards; 1 4 25 7 18 .106 .90 .72 3.5 33.00 2 5 33 14 19 .078 .18 .68 4.0 38.00 3 6 25 9 16 .126 .98 .80 4.8 50.00 4 7 34 8 26 .134 .70 .68 3.4 78.00 5 1 33 10 23 .142 .84 .68 3.5 76.00 6 2 20 9 11 .110 .64 .51 1.0 96.00 7 3 19 5 14 .116 .52 .51 2.6 88.00 8 4 27 11 16 .116 .78 .59 4.5 79.00 . . . 729 4 1 1 0 .162 .74 .54 3.6 32.00 730 5 0 0 0 .102 .54 .50 1.0 30.00 ; proc gam; model logn1=param(pm10)spline(time)spline(t)spline(rh)/dist=poisson; score data=a; run; 因数据较多，没有全部贴上来 |
回复此楼» 猜你喜欢
聘U V热熔胶研究人员
已经有10人回复
-
职称评审没过，求安慰
已经有7人回复
-
求助文献
已经有3人回复
-
垃圾破二本职称评审标准
已经有10人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复