[問題] 用lsmeans做ANOVA post hoc分析

看板R_Language作者salutyz (寫給城市的詩)時間9年前 (2016/12/26 10:38)推噓0(0推 0噓 7→)

留言7則, 3人參與討論串1/1

[問題類型]: 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來) [軟體熟悉度]: 請把以下不需要的部份刪除新手(沒寫過程式，R 是我的第一次) [問題敘述]: 請簡略描述你所要做的事情，或是這個程式的目的有A,B,C三個factors的repeated anova (type III)，已經跑出anova的結果，想要繼續做post hoc 要用lsmean和constrast的方式（這是依據同領域朋友的建議，不過他是用SAS裡proc glm 去跑底下 lsmeans ( for simple main effect) 和contrast(for differences) 查了很多資料用下面的程式跑之後在ref.grid()和lsmeans()都出現Error message Warning message: In lsm.basis.aovlist(object, trms, xlev, grid, ...) : Some predictors are correlated with the intercept - results are biased. May help to re-fit with different contrasts, e.g. 'contr.sum' 不確定要怎樣修改才對..... orz [程式範例]: fit.1 <- aov(value~A*B*C+Error(Subject/(A*B*C)),data=work.1) (late.rg1 <- ref.grid(fit.1)) lsmeans(fit.1, "A") [環境敘述]: 請提供 sessionInfo() 的輸出結果，裡面含有所有你使用的作業系統、R 的版本和套件版本資訊，讓版友更容易找出錯誤 R version 3.1.3 (2015-03-09) Platform: x86_64-apple-darwin13.4.0 (64-bit) Running under: OS X 10.10.5 (Yosemite) locale: [1] C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] lsmeans_2.20-23 estimability_1.1-1 car_2.0-25 nlme_3.1-126 [關鍵字]: 選擇性，也許未來有用因為是超新手所以折騰了好幾個禮拜跪求各位高手幫忙感激不盡（跪） -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 111.240.212.180 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1482719915.A.EE7.html