作者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/cn.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