#object;回归方程之间的比较
#writer: mike
#time:2020,11,14
a <- c(1,2,4,5,67,8)
b <- c(3,4,6,556,86,234)
d <- c(111,32,123,543,64,65)
#构造数据集
data <- data.frame(a,b,d)
#建立用于比较的模型
model1 <- lm(d~a,data=data)
model2 <- lm(d~a+b,data=data)
#比较模型,注意这两个模型是有嵌套关系的
res <- anova(model1,model2)
print(res)
#Analysis of Variance Table
#Model 1: d ~ a
#Model 2: d ~ a + b
#Res.Df RSS Df Sum of Sq F Pr(>F)
#1 4 176231
#2 3 38667 1 137563 10.673 0.04689 *
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#注意这里的结果的计算方法,(两个方程残差的平方和的差值除以自由度之差) / (最长的方程的残差的平方和除以自由度)