【统计建模与R软件】第二章 R软件的使用 上机练习5 R语言作图 线性回归 多图绘制 线条格式 箭头 文字 图例

上机练习及解析

1）数据的读取

数据：

weight height
4.17 1.75
5.58 1.8
5.18 1.65
6.11 1.9
4.5 1.74
4.61 1.91
5.17 1.75
4.53 1.8
5.33 1.65
5.14 1.9
4.81 1.74
4.17 1.91
4.41 1.75
3.59 1.8
5.87 1.65
3.83 1.9
6.03 1.74
4.89 1.91
4.32 1.75
4.69 1.8

代码：

# T1.1
d = read.table("HW.txt",header=T)

2）绘图它们的回归曲线，并预测在weight 等于6.85时身高等于多少？

# T2.1 画图
model = lm(d$height~d$weight)
plot(d$weight,d$height,xlab="weight",ylab="height")
abline(model) # 加上回归的线

# T2.2 预测
# 踩坑警告：变量必须从data.frame里面分离出来！
# 否则predict结果会有length(你的data.frame)那么多个数！

# 1.变量分离
x = d$weight
y = d$height
# 2.重新建模、预测
model = lm(y~x)
predict(model,newdata=data.frame(x=6.85))
# 这样最终就只有一个结果

结果：