作者nckuwagtail (鶺鴒)
看板Statistics
標題[問題] 重複測量的迴歸分析
時間Tue Aug 13 16:11:54 2013
小妹只有學過基礎統計,目前手邊資料有點令我傷腦筋,
查了許多文獻資料,但還是沒有頭緒。
我有一批樹苗,
每顆樹苗每隔一段時間就記錄其存活狀況(存活/死亡)與當時的環境值,
一共有八次的存活及環境資料。
環境值有:光照等級-類別、地表落葉厚度-連續、是否有病蟲害-類別,
但如果樹苗死亡,就沒有測量環境資料了。
我想知道樹苗存活是否跟環境值有關。不知該用哪種分析?
似乎是該用logistic regression(?),
但因為有連續時間的重複調查,所以每次的環境值之間是有關連性的...,
然後我就卡住了。
我手邊的軟體是SAS,
看了許多proc phreg跟proc logistic的介紹,還是找不到合用的模組。
另一個類似的狀況,
同一批數據,每顆樹苗每隔一段時間就記錄其高度(連續)與當時的環境值,
我想知道樹高是否跟環境值有關。
好像要用proc glm? 但還是不知道該怎麼佈局。
希望可以有版友可以指點一下方向。
感謝。
--
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.67.82.165
1F:推 ching0629:Generalized estimating equation(GEE)應該可以滿足你 08/13 21:56
2F:→ ching0629:glm對你來講不夠,你用glm等同於用logistic regression 08/13 21:57
3F:→ ching0629:你的兩個狀況基本上是類似的 08/13 21:58
4F:→ nckuwagtail:剛剛試了一下,但遇到問題:如果苗木死掉的那次, 08/14 10:16
5F:→ nckuwagtail:就沒有測量環境值,這樣資料不知道會不會有問題? 08/14 10:16
6F:推 ax4869:成大先推~~ 08/14 18:52
7F:→ ching0629:也許要有點折衷,死亡紀錄必須記錄在前一次 08/16 23:54
8F:→ ching0629:而預測的目標變為"下次觀察死亡" 08/16 23:54
9F:推 anovachen:跟存活有關的問題,為何不用存活分析呢XD 08/17 00:09
10F:→ anovachen:Cox proportional hazard model或許能用? 08/17 00:10
11F:→ anovachen:問題二應該要用GEE才能解決了... 08/17 00:10
12F:→ anovachen:把環境值的Hazard ratio算出,HR顯著不為1就代表有影響 08/17 00:11
13F:→ anovachen:呃...如果環境值也隨時間變動,那到底要怎麼處理= = 08/17 00:13
15F:→ anovachen:Time-dependent covariates in the Cox proportional- 08/17 00:17
16F:→ anovachen:hazards regression model. Fisher LD, Lin DY. 1999 08/17 00:17
18F:→ anovachen:Your “Survival” Guide to Using Time‐Dependent 08/17 00:21
19F:→ anovachen:Covariates. Teresa M. Powell, Melissa E. Bagnell 08/17 00:21
20F:→ ching0629:A大,cox只能使用到baseline的環境監測值 08/17 18:22
21F:→ ching0629:明明就存在著各階段的環境值卻假設它不變並不恰當 08/17 18:23
22F:推 bmka:GEE不適用在survival data 09/02 12:38
23F:→ bmka:Cox PHM可處理time-dependent covariate 09/02 12:40
24F:→ bmka:觀察時間點少則應該考慮proportional odds model 09/02 12:41