Построение графика перпендикуляра смещает в графике регрессии наименьших квадратов в R

Question

Построение графика перпендикуляра смещает в графике регрессии наименьших квадратов в R

Я интересуюсь созданием графика с линией регрессии наименьших квадратов и линейными сегментами, подключающими точки данных к линии регрессии, как проиллюстрировано здесь в графических названных перпендикулярных смещениях: http://mathworld.wolfram.com/LeastSquaresFitting.html
_{(от MathWorld - веб-ресурс вольфрама: 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")

Я должен использовать некоторый для цикла и сегментов http://www.iiap.res.in/astrostat/School07/R/html/graphics/html/segments.html, чтобы сделать перпендикулярные смещения? Существует ли более эффективный путь? Обеспечьте пример, если это возможно.

16

r statistics plot linear-regression least-squares

задан Glorfindel 23 July 2019 в 01:22

1 ответ

Другие вопросы по тегам:

r statistics plot linear-regression least-squares

Похожие вопросы:

score 19 · Accepted Answer

Сначала необходимо определить координаты основания перпендикулярных сегментов, затем вызвать функцию сегментов , которая может принимать векторы координаты в качестве входных данных (нет необходимости в цикле).

perp.segment.coord <- function(x0, y0, lm.mod){
 #finds endpoint for a perpendicular segment from the point (x0,y0) to the line
 # defined by lm.mod as y=a+b*x
  a <- coef(lm.mod)[1]  #intercept
  b <- coef(lm.mod)[2]  #slope
  x1 <- (x0+b*y0-a*b)/(1+b^2)
  y1 <- a + b*x1
  list(x0=x0, y0=y0, x1=x1, y1=y1)
}

Теперь просто вызовите сегменты:

ss <- perp.segment.coord(temperature, diseasesev, severity.lm)
do.call(segments, ss)
#which is the same as:
segments(x0=ss$x0, x1=ss$x1, y0=ss$y0, y1=ss$y1)

Обратите внимание, что результаты не будут выглядеть перпендикулярно, если вы не убедитесь, что единицы x и y вашего графика имеют одинаковую видимую длину (изометрические масштабы). Вы можете сделать это, используя pty = "s" , чтобы получить квадратный график и установить xlim и ylim в тот же диапазон.