| 查看: 2402 | 回复: 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; 因数据较多，没有全部贴上来 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有27人回复
-
所感
已经有3人回复
-
要不要辞职读博？
已经有7人回复
-
不自信的我
已经有11人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复