2014-10-05 135 views
2

我正試圖在SAS中進行一個簡單的單向隨機效應ANOVA。我想知道人口差異是否顯着不同於零。SAS中的單向隨機效應ANOVA:PROC GLM還是混合?

UCLA's idre site,他們聲明如下使用PROC MIXED:

proc mixed data = in.hsb12 covtest noclprint; 
    class school; 
    model mathach =/solution; 
    random intercept/subject = school; 
run; 

這對我來說很有意義給予我以前用PROC MIXED經驗。然而,在文本Murray Logan的生物統計設計和分析使用R的文章中,他表示,對於單向方差分析,固定和隨機效應沒有被區分並進行(在R中)「標準」單向方差分析儘管他正在測試方差,而不是手段。我發現,在SAS,他[R過程等同於使用任何以下內容:

  1. PROC ANOVA
  2. PROC GLM(同方差分析,但GLM到位方差分析)
  3. PROC GLM用隨機聲明

p從以上三種模式 - 值是相同的,而是來自加州大學洛杉磯分校使用的PROC MIXED型號而異。根據我的數據,這是p = 0.2508和p = 0.3138的差異。儘管在這種情況下結論不會改變,但我對這種差異感到不舒服。

任何人都可以提供建議哪一個更合適,也是爲什麼有這種差異?

+0

您在PROC ANOVA調用中使用了哪些選項? – Kostya 2014-10-15 23:10:56

+0

無,只是PROC ANOVA,其次是CLASS和MODEL語句。 – Meg 2014-10-16 00:20:36

回答

1

對於您的模型,PROC ANOVA和PROC MIXED之間的區別僅僅是由於數值噪聲(PROC MIXED的REML估計器)。但是,您問題中提到的p值對應於不同的測試。爲了在PROC MIXED中使用COVTEST的輸出獲得F值,需要考慮不相等的樣本大小重新計算MS_groups(手動方式如http://bio.classes.ucsc.edu/bio286/MIcksBookPDFs/QK08.PDF的p.231所述,或者僅使用PROC MIXED和相同的固定模型spec如PROC ANOVA中所述)。本文(http://isites.harvard.edu/fs/docs/icb.topic1140782.files/S98.pdf)提供了除了SAS手冊之外還使用了PROC MIXED的一些示例。

+0

非常有幫助,謝謝。我將不得不尋找一些時間來閱讀這些參考資料,以瞭解爲什麼PROC MIXED無法處理不平衡的數據,但PROC ANOVA和GLM都很好。 – Meg 2014-11-23 15:12:24