| 查看: 2448 | 回复: 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; 因数据较多，没有全部贴上来 |
回复此楼» 猜你喜欢
售SCI文章，我:8.O.5.51O.5.4,科目齐全,可十急
已经有3人回复
-
O(∩_∩)O
已经有4人回复
-
T^T
已经有4人回复
-
面上项目没有好文章就没希望了吗？
已经有21人回复
-
[急招]吉林大学招收26年9月入学考核制博士生1名，截止今晚24点前
已经有4人回复
-
o>_<o
已经有3人回复
-
重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶
已经有7人回复
-
前几天时间戳更新了
已经有13人回复
-
AASSDDFF
已经有3人回复
-
→_→
已经有3人回复