[程式] SPSS存活分析問題

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