| 查看: 2406 | 回复: 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; 因数据较多，没有全部贴上来 |
回复此楼» 猜你喜欢
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有6人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复