0
点赞
收藏
分享

微信扫一扫

【R】【lm()】初理解

史值拥 2022-04-06 阅读 14

文章目录

浅试lm()

a<-c(1,2,3,4,5)
b<-c(2,4,6,8,10)
mydata<-data.frame(a,b)

在这里插入图片描述

myfit<-lm(a~b,mydata)
#myfit<-lm(formula,data)
summary.lm(myfit)

在这里插入图片描述

警告

R语言错误的提示(中英文翻译)–bahui7562–CSDN

R语言填坑–丶L–CSDN

文档示例(暂未精读)

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)
summary(lm.D90)

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

在这里插入图片描述

举报

相关推荐

0 条评论