作者eileen0516 (...)
看板Statistics
标题[程式] SPSS存活分析问题
时间Wed Apr 30 18:22:24 2014
[软体程式类别]:
SPSS
[程式问题]:
存活分析
[软体熟悉度]:
低(1~3个月)
[问题叙述]:
跑存活分析时，有些变项的平均数及中位数无法显示，而报表出现警告字样
未计算统计量，因为所有观察值都已审查过
不晓得这句话表示什麽意思，是我哪个环节出问题
麻烦解惑 谢谢
[程式范例]:
无
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 36.237.229.51
※ 文章网址: http://webptt.com/cn.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