作者windguide (夜很深)
看板Statistics
标题[问题] 用R分析2^4 design data
时间Wed Apr 28 12:08:39 2010
如果是跟统计软体有关请重发文章
如果跟论文有关也烦请您重发文章
文章类别是为了帮助大家搜寻资料与解答，造成不便之处请见谅
其实我遇到的问题是这样的
我用R分析2^4的Factorial Design，在做ANOVA分析的时候没有问题
问题是出在做回归分析的时候：
我记得回归分析的系数就代表了那个因子效应大小的一半
可是我做出来的回归分析系数跟我实际计算的因子效应并没有关系
例如：
回归分析： x_1 的系数 = A因子效应大小的一半 = 8.2526
实际用手算： A因子大小 = 3.0188
请问为什麽会出现这样的问题，是我的Model Fit错误吗?
麻烦各位帮我看一下，谢谢
以下是CODE
y =
c(7.037,6.376,14.707,15.219,11.635,12.089,17.273,17.815,10.403,10.151,
4.368,4.098,9.360,9.253,13.440,12.923,8.561,8.951,16.867,17.052,
13.876,13.658,19.824,19.639,11.846,12.337,6.125,5.904,11.190,10.935,
15.653,15.053)
trtA <- gl(2,2,32)
trtB <- gl(2,4,32)
trtC <- gl(2,8,32)
trtD <- gl(2,16,32)
out = aov(y~trtA*trtB*trtC*trtD)
summary(out)
out1 = lm(y~trtA*trtB*trtC*trtD)
summary(out1)
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.75.121
1F：推 kazuki:在two-level design factorial effects跟 regression 04/28 21:03
2F：→ kazuki:系数的估值会呈比例，也就是你说得1/2还是两倍。 04/28 21:04
3F：→ kazuki:在但mixed-level experiments里就不一定了。以上是我猜的^^ 04/28 21:05