Displaying 1 result from an estimated 1 matches for "ma_resid".
Did you mean:
asresid
2008 Oct 21
0
Major Axis residuals
...:100 + 50 * runif(100)
z <- data.frame(x = x, y = y)
plot(z)
lm_xy <- lm(z$y ~ z$x)
abline(lm_xy, col = "blue")
z$lm_resid <- z$y - (z$x * coef(lm_xy)[2]) - coef(lm_xy)[1]
ma_xy <- line.cis(y = z$y, x = z$x, method = 2)
abline(ma_xy[1,1], ma_xy[2,1], col = "red")
z$ma_resid <- (z$y - (ma_xy[2,1] * z$x) - ma_xy[1,1])/
sqrt((ma_xy[2,1])^2 + 1)
plot(z$x, z$lm_resid)
abline(lm(z$lm_resid ~ z$x), col = "blue")
plot(z$x, z$ma_resid)
abline(lm(z$ma_resid ~ z$x), col = "red")
abline(h = 0, col = "gray")
______________________
C. E. Timothy...