R lm摘要错误

Question

不知怎么的，我的summary()函数停止了对lm回归输出的工作。为了演示我得到的问题和错误消息，我从“？lm”运行了示例回归：

ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
lm.D90 <- lm(weight ~ group - 1) # omitting intercept

anova(lm.D9)
Analysis of Variance Table

Response: weight
          Df Sum Sq Mean Sq F value Pr(>F)
group      1 0.6882 0.68820  1.4191  0.249
Residuals 18 8.7292 0.48496               
summary(lm.D90)
Error in summary.lm(lm.D90) : 
  length of 'dimnames' [1] not equal to array extent
In addition: Warning message:
In cbind(est, se, tval, 2 * pt(abs(tval), rdf, lower.tail = FALSE)) :
  number of rows of result is not a multiple of vector length (arg 1)

opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0))
plot(lm.D9, las = 1)      # Residuals, Fitted, ...
par(opar)

对于其他lm命令，错误消息也会出现，我在过去获得了这些命令的摘要输出，但从一秒到另一秒，R现在只给我上面的错误消息。

即使我重新启动R studio或我的计算机，也会出现错误信息。有没有人有一个建议，为什么会出现这种情况？任何帮助都是非常感谢的。

Answer 1

它确实像这样在我的机器上运行：

ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
lm.D90 <- lm(weight ~ group - 1) # omitting intercept
anova(lm.D9) 
Analysis of Variance Table

Response: weight
          Df Sum Sq Mean Sq F value Pr(>F) group      1 0.6882 0.68820  1.4191  0.249 Residuals 18 8.7292 0.48496               
summary(lm.D90)

Call: lm(formula = weight ~ group - 1)

Residuals:
    Min      1Q  Median      3Q     Max 
-1.0710 -0.4938  0.0685  0.2462  1.3690 

Coefficients:
         Estimate Std. Error t value Pr(>|t|)     groupCtl   5.0320     0.2202   22.85 9.55e-15 *** groupTrt   4.6610     0.2202   21.16 3.62e-14 ***
--- Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.6964 on 18 degrees of freedom Multiple R-squared:  0.9818,   Adjusted R-squared:  0.9798  F-statistic: 485.1 on 2 and 18 DF,  p-value: < 2.2e-16