作者firstwind (初聆)
看板Statistics
標題[問題] 雙因子變異數分析變方不同質
時間Tue Jul 20 16:05:03 2021
如果是跟統計軟體有關請重發文章，使用程式做為分類。
統計軟體，如SPSS, AMOS, SAS, R, STATA, Eviews，請都使用程式做為分類
請詳述問題內容，以利板友幫忙解答，過短文章依板規處置，請注意。
為避免版面混亂，請勿手動置底問題，善用E做檔案編輯
第一次發文，請多包涵
我想要問關於two-way ANOVA的問題
https://imgur.com/bF3Nq9X
這是我的資料，light跟fertilizers分別是兩個處理因子
在看boxplot的時候就覺得很不妙了.......
https://imgur.com/KVPGwnG
作檢查的時候，果然發現變異數不同質(符合常態分佈
想請問遇到這種狀況，要怎麼解決比較好？
1.換一種統計方法？
但統計修得不夠好，也不知道應該換哪種統計方法才好(用的統計軟體是R
2.資料轉換？
問題同上......不知道該怎麼轉換比較合適
有看到使用welch的方法，但似乎是在回答單因子變異數分析時變方不同質的狀況，我不確定雙因子會不會有所不同
先感謝願意善心解答的大家
我會將各位寫在我論文的誌謝中的QwQ
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 59.120.16.47 (臺灣)
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1626768305.A.52F.html
1F：推 andrew43: car package Anova() 裡的 white.adjuat 了解一下 07/20 17:21
安裝完package以後，寫出來多加了white.adjust這個部分，做出來的結果解讀方式是不是就跟一般的two way Anova一樣呢？
(以我的結果是兩處理交感不顯著，需要再回歸到個別單因子分析這樣
https://imgur.com/goSaL9I
2F：→ Pieteacher: 做 gls 07/20 18:11
做gls的時候是否也要先大概知道變數之間的關係呢？
這是直接做的結果(code:mod<-lm(content~light*fertilizers, data=rutin)
https://imgur.com/q2GHtmg
下面是用a大的方法，先得知兩處理交感不顯著後的結果
https://imgur.com/xsJZOcx
※ 編輯: firstwind (140.112.74.142 臺灣), 07/21/2021 10:42:35
3F：推 andrew43: gls要用package，不是內建的lm 07/21 13:56
4F：→ andrew43: Anova建議用type III 07/21 13:58
5F：→ andrew43: 不過沒交互作用就無所謂 07/21 13:59
package是nlme嗎？裝了以後說沒有gls這個函數QQ
我一開始其實不是mod<-lm，而是mod<-aov
不確定Anova對變數的對比項之類的有沒有要求，所以不用lm，回來用aov寫了兩種
https://imgur.com/NkI9hsA
https://imgur.com/makpjAO
兩種看起來結果都是交感不顯著，但還是想知道應該用哪一個才是對的......
或者應該要怎麼寫比較好
不好意思問題有點多，再次感謝幫忙
※ 編輯: firstwind (140.112.74.142 臺灣), 07/21/2021 15:17:03
6F：→ andrew43: 留意指定contrast會影響有交互作用的主效應意義。若丟 07/21 15:31
7F：→ andrew43: 掉交互作用則沒差。 07/21 15:31
8F：→ andrew43: nlme的gls可以用，且要指定合適的weights 07/21 15:34
非常感謝解答！gls我跟它沒有很熟，也真的叫不出這個函數QQ，先把它放在一邊好了
我的實驗是平衡設計，那就先不指定contrast了
會先用上面的第一個方法，也會在跟老師確認
真的真的非常感謝！
※ 編輯: firstwind (140.112.74.142 臺灣), 07/21/2021 16:27:37
9F：→ andrew43: 留意你指定contrast是外部的變數不是d裡面的。題外話。 07/21 19:03