| 查看: 2414 | 回复: 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; 因数据较多，没有全部贴上来 |
回复此楼» 猜你喜欢
081700，311，求调剂
已经有15人回复
-
一志愿北京化工085600 310分求调剂
已经有18人回复
-
085600材料与化工301分求调剂院校
已经有4人回复
-
材料与化工371求调剂
已经有14人回复
-
336材料与化工085600求调剂
已经有7人回复
-
材料334求调剂
已经有18人回复
-
331求调剂
已经有8人回复
-
332求调剂
已经有17人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有4人回复
-
材料专硕322
已经有7人回复
10