ここで垂直オフセットと呼ばれる図に示されているように、最小二乗回帰線と、データポイントを回帰線に接続する線分を使用してプロットを作成することに興味があり
ます。
- Wolfram Web リソース: wolfram.com )
ここでプロットと回帰直線を作成しました:
## Dataset from http://www.apsnet.org/education/advancedplantpath/topics/RModules/doc1/04_Linear_regression.html
## Disease severity as a function of temperature
# Response variable, disease severity
diseasesev<-c(1.9,3.1,3.3,4.8,5.3,6.1,6.4,7.6,9.8,12.4)
# Predictor variable, (Centigrade)
temperature<-c(2,1,5,5,20,20,23,10,30,25)
## For convenience, the data may be formatted into a dataframe
severity <- as.data.frame(cbind(diseasesev,temperature))
## Fit a linear model for the data and summarize the output from function lm()
severity.lm <- lm(diseasesev~temperature,data=severity)
# Take a look at the data
plot(
diseasesev~temperature,
data=severity,
xlab="Temperature",
ylab="% Disease Severity",
pch=16,
pty="s",
xlim=c(0,30),
ylim=c(0,30)
)
abline(severity.lm,lty=1)
title(main="Graph of % Disease Severity vs Temperature")
垂直オフセットを行うには、何らかの for ループとセグメントhttp://www.iiap.res.in/astrostat/School07/R/html/graphics/html/segments.htmlを使用する必要がありますか? より効率的な方法はありますか?可能であれば例を挙げてください。