R做方差分析出现问题了！！！！！！！！！！！！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

> lm2.2=lm(log.均价~as.factor(环线)+as.factor(物业))
> Anova(lm2.2,type="3")
Anova Table (Type III tests)

Response: log.均价
                 Sum Sq Df  F value    Pr(>F)   
(Intercept)     163.568  1 8516.692 < 2.2e-16 ***
as.factor(环线)   0.708  4    9.210 0.0002176 ***
as.factor(物业)   0.627  2   16.311  6.29e-05 ***
Residuals         0.384 20                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

（上面做双因素不带交互作用方差分析没有问题，下面做交互为什么就出现问题了呢？求解）
> lm2.2=lm(log.均价~as.factor(环线)* as.factor(物业))
> Anova(lm2.2,type="3")
错误于linear.hypothesis.lm(mod, hyp.matrix, summary.model = sumry,  :
  One or more terms aliased in model

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：方差分析 HYPOTHESIS Intercept Residuals Residual 方差分析

lm2.2=lm(log.均价~as.factor(环线)* as.factor(物业))

应该改为：lm2.2=lm(log.均价~as.factor(环线)＋as.factor(物业)＋as.factor(环线)* as.factor(物业) )

可见问题不在R，而是你没理解好*的含义。

另外，做方差分析，建议使用aov函数

呵呵，谢谢LS的。明天我试试。我是按照教材上讲的写的，可能是印错了吧。

还是不行啊？？

> lm2.1=lm(log.均价~as.factor(环线)+as.factor(位置))
> Anova(lm2.1,type="3")
Anova Table (Type III tests)

Response: log.均价
                 Sum Sq Df   F value    Pr(>F)   
(Intercept)     248.805  1 7440.2662 < 2.2e-16 ***
as.factor(环线)   0.667  4    4.9857  0.005518 **
as.factor(位置)   0.308  1    9.2223  0.006270 **
Residuals         0.702 21                       
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> summary(lm2.1)

Call:
lm(formula = log.均价 ~ as.factor(环线) + as.factor(位置))

Residuals:
      Min        1Q    Median        3Q       Max
-0.334729 -0.110345 -0.009307  0.093236  0.295793

Coefficients:
                     Estimate Std. Error t value Pr(>|t|)   
(Intercept)          8.032209   0.093120  86.257  < 2e-16 ***
as.factor(环线)三环 -0.003559   0.102048  -0.035  0.97251   
as.factor(环线)四环 -0.273823   0.118245  -2.316  0.03077 * 
as.factor(环线)五环 -0.240490   0.118857  -2.023  0.05595 . 
as.factor(环线)一环  0.220214   0.119351   1.845  0.07917 . 
as.factor(位置)县城 -0.253537   0.083488  -3.037  0.00627 **
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.1829 on 21 degrees of freedom
Multiple R-squared: 0.7096,     Adjusted R-squared: 0.6405
F-statistic: 10.26 on 5 and 21 DF,  p-value: 4.384e-05

> lm2.2=lm(log.均价~as.factor(环线)+as.factor(物业)+as.factor(环线)*as.factor(物业))
> Anova(lm2.2,type="3")
错误于linear.hypothesis.lm(mod, hyp.matrix, summary.model = sumry,  :
  One or more terms aliased in model.

为什么还是出错啊？

好像是我数据本身的问题。重新换了两个因素就可以了。

> lm2.2=lm(log.均价~as.factor(环线)+as.factor(装修)+as.factor(环线)*as.factor(装修))
> Anova(lm2.2,type="3")
Anova Table (Type III tests)

Response: log.均价
                                 Sum Sq Df   F value Pr(>F)   
(Intercept)                     124.765  1 2173.0384 <2e-16 ***
as.factor(环线)                   0.297  4    1.2936 0.3117   
as.factor(装修)                   0.004  1    0.0643 0.8029   
as.factor(环线):as.factor(装修)   0.035  4    0.1504 0.9602   
Residuals                         0.976 17                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

错误于linear.hypothesis.lm(mod, hyp.matrix, summary.model = sumry,  :
  One or more terms aliased in model （这个是什么含义啊？？？？？？？？？？？）