| 1 | 1/1 | 返回列表 |
| 查看: 2413 | 回復(fù): 0 | ||
hi11銅蟲 (初入文壇)
|
[求助]
sas或R如何做空氣污染物和入院人數(shù)的時(shí)間序列分析(廣義相加模型GAM),有程序和數(shù)據(jù)
|
|
本科畢業(yè)論文想做一個(gè)關(guān)于PM10和每日入院人數(shù)的時(shí)間序列分析,用的是GAM模型,但在選取模型樣條函數(shù)的自由度上無從下手,還有控制時(shí)間的變量如何加入模型語句中?以及控制星期幾效應(yīng)的啞元變量(dow)如何納入模型? 下面為sas初步的程序語句,其中輸入的變量分別為:時(shí)間(1~730),控制星期幾效應(yīng)的啞元變量dow(1~7),總?cè)朐喝藬?shù)n,心血管病入院人數(shù)n1,呼吸系統(tǒng)病人數(shù)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; 因數(shù)據(jù)較多,沒有全部貼上來 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)279求調(diào)劑 +5 | 哈哈哈^_^ 2026-03-31 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +12 | Youth_ 2026-03-30 | 12/600 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +10 | Li李魚 2026-03-26 | 10/500 |
|
|
[考研] 一志愿鄭州大學(xué),080500學(xué)碩,總分317分求調(diào)劑 +10 | 舉個(gè)栗子oi 2026-03-24 | 11/550 |
|
|
[考研] 287求調(diào)劑 +14 | land xuxu 2026-03-26 | 14/700 |
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
[考研] 材料專碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
[考研] 0703化學(xué) +14 | 妮妮ninicgb 2026-03-27 | 14/700 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 311(085601)求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|