作者eileen0516 (...)
看板Statistics
標題[程式] SPSS存活分析問題
時間Wed Apr 30 18:22:24 2014
[軟體程式類別]:
SPSS
[程式問題]:
存活分析
[軟體熟悉度]:
低(1~3個月)
[問題敘述]:
跑存活分析時，有些變項的平均數及中位數無法顯示，而報表出現警告字樣
未計算統計量，因為所有觀察值都已審查過
不曉得這句話表示什麼意思，是我哪個環節出問題
麻煩解惑 謝謝
[程式範例]:
無
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.237.229.51
※ 文章網址: http://webptt.com/m.aspx?n=bbs/Statistics/M.1398853346.A.707.html
1F：→ bmka:存活分析方法很多種,你先把問題講清楚 04/30 20:54
2F：→ bmka:其實我可以猜到你遇到什麼問題(其實根本不是問題) 04/30 20:55
3F：→ bmka:但是還是要勸你,跑分析之前至少要對統計方法有基本認識 04/30 20:56
4F：→ eileen0516:我是利用kaplan-meier分析在不同病理分期的存活率 05/01 18:47
5F：→ eileen0516:所以問題是出在哪裡，問題還有哪裡需要補充嗎?麻煩解答 05/01 18:49
6F：→ eileen0516:謝謝 05/01 18:49
7F：→ bmka:不是問題,你的data就是沒那麼多failure event. 05/01 19:16
8F：→ bmka:K-M curve是nonparametric estimator 05/01 19:48
9F：→ bmka:data到哪裡,就只能估計到哪裡 05/01 19:49
10F：→ eileen0516:那我大概瞭解謝謝你 05/01 22:23
11F：→ eileen0516:不好意思，如果我還是想知道平均數和中位數 05/01 23:31
12F：→ eileen0516:我需要調整哪個變項，還是用另外的方法，我這一塊真的 05/01 23:32
13F：→ eileen0516:不是很熟 05/01 23:32
14F：→ bmka:如果一定要,那麼你的選擇只有parametric modeling 05/01 23:48
15F：推 excimo:審查 = censor 意思就是說 你的資料中沒有"event"發生 05/02 15:18
16F：→ excimo:你的資料庫當中沒有人死亡 自然KM算不出來 05/02 15:19
17F：→ bmka:樓上別來亂, 原po的data不是沒有event,而是event數不夠多 05/02 20:02
18F：→ bmka:而且最大的觀察值也被censor掉,所以KM的值都維持在0.5以上 05/02 20:03
19F：→ bmka:以至於median是not estimable. 至於平均值不能估是survival 05/02 20:04
20F：→ bmka:data的觀察時間不夠長, KM這種nonparametric method無法在所 05/02 20:04
21F：→ bmka:有failure event可能出現的的範圍作估計,當然mean也是not 05/02 20:06
22F：→ bmka:estimable 05/02 20:06
23F：推 lin15:Median算不出來是沒問題但是SPSS算Mean只要有一個event就會 05/03 14:27
24F：→ lin15:計算了 但Estimation is limited to the largest survival 05/03 14:27
25F：→ lin15:time if it is censored. 所以原po應該真的沒有event或操作 05/03 14:28
26F：→ lin15:有錯誤吧... 05/03 14:28
27F：→ bmka:如果SPSS只要有event就會report mean,那就表示這個軟體很糟糕 05/03 19:43
28F：推 excimo:我剛試過了 要跑出原PO的錯誤訊息 真的是因為沒有Event 05/03 22:16
29F：→ excimo:只要有一個event 就可以算出Mean 05/03 22:16
30F：→ bmka:那真的很糟糕,這樣會誤導不懂統計的人 05/03 22:52
31F：→ bmka:換個軟體吧(笑) 05/03 22:53
32F：→ bmka:還是稍微解釋一下, 通常因為觀察時間不夠長,無母數方法不能 05/03 22:59
33F：→ bmka:用來估計mean (超過最大觀察值以後就不知道會發生什麼事) 05/03 23:00
34F：→ bmka:能做的就是用area under the KM curve來估計restricted mean 05/03 23:03
35F：→ bmka:假設有興趣的survival time r.v.是T,實際上的研究只追蹤了 05/03 23:08
36F：→ bmka:一年, 那麼除非所有的failure event 都會在一年內發生, 05/03 23:10
37F：→ bmka:否則T的期望值是無法用無母數方法來估計的 05/03 23:10
38F：→ bmka:但是我們可以估計 T*I(T<=1) 的期望值,也就是restricted mean 05/03 23:11
39F：→ bmka:(題外話)其實這幾年哈佛的魏立人老師一直在提倡用restricted 05/03 23:13
40F：→ bmka:mean來做clinical trial的endpoint 05/03 23:14
41F：→ bmka:而他的演講舉的例子也就是用其中一個subgroup沒有event的情況 05/03 23:17
42F：→ bmka:來說明為什麼restricted mean比較適合當作outcome來比較 05/03 23:17