作者anovachen (囧)
看板Statistics
標題[程式] [SAS] 計算標準誤出現不一致的結果
時間Thu Apr 25 23:51:20 2013
[軟體程式類別]:
SAS
[程式問題]:
用proc means stderr算standard error
和用proc glm中lsmeans /stderr求得的standard erro不相等。
[軟體熟悉度]:
高(1年以上)
[問題敘述]:
一群病人存活時間的資料，
依照年齡分組，
跑ANOVA。
順便要算各分組之間的standard deviation。
本來想用lsmeans /stderr算得的標準誤乘上√(該組樣本數)
但後來改用proc means算，同時確認stderr有沒有一樣，
卻意外發現不相同。
[程式範例]:
proc sort data=brca; by age_group; run;
proc means data=brca n mean std stderr;
var time;
by age_group;
run;
proc glm data=brca;
class age_group;
model time=age_group;
lsmeans age_group/ stderr; run;
兩者求得之標準誤會不一樣。為何會這樣？
(ps.我知道存活時間的資料應該改用存活分析，
但這是某課程的題目，所以就暫且按照題目要求用ANOVA)
The SAS System
MEANS 程序
age_group=1
分析變數 : time Time (months)
N 平均值 標準差 標準誤差
324 48.8735597 28.6992884 1.5944049
age_group=2
分析變數 : time Time (months)
N 平均值 標準差 標準誤差
596 44.5850671 28.4323597 1.1646348
The SAS System
GLM 程序
最小平方平均值
age_group time LSMEAN 標準誤差 Pr > |t|
1 48.8735597 1.5848091 <.0001
2 44.5850671 1.1684936 <.0001
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.67.123
話說...剛才發現多寫這一行：
means age_group;
就能算standard deviation了XXD
這邊算出的標準差倒是跟proc means算出來的一樣。
※ 編輯: anovachen 來自: 140.116.67.123 (04/26 00:12)