2013-04-06 51 views
0

代碼是這裏下面的代碼。迴歸線不通過點的雲

altura <- read.table("altura.txt", header=T, quote="\"") 
altura <- cbind(altura, altura$Esposa/altura$X.Marido, altura$X.Marido/altura$Esposa) 
is.data.frame(altura) 
names(altura) <- c("marido","esposa","r1","r2") 
with(altura,plot(marido~esposa)) 
g1 <- lm(altura$esposa~altura$marido) 
summary(g1) 
abline(g1$coefficients) 
abline(0,1,lty=5) 
with(altura,plot(esposa~marido)) 
g2 <- lm(altura$marido~altura$esposa) 
summary(g2) 
abline(g2$coefficients) 
abline(0,1,lty=5) 
cor(altura$marido,altura$esposa) 

簡單迴歸線不通過點雲。 abline使用匯總函數中的正確截距。這不是第一次發生。如你所見,在這兩個圖表中,我就是這個問題。一條線通過點,另一條線在下面。

+0

任何人都可以幫忙嗎?我認爲計算錯誤的係數是可能的?這是正常的嗎? – 2013-04-08 23:50:13

+0

你可以發佈你的數據還是圖表的圖片? – 2013-04-12 18:10:31

回答

3

我想我看到:您G2型號: g2 <- lm(altura$marido~altura$esposa應該with(altura,plot(marido~esposa))去你with(altura,plot(esposa~marido))

有它如

set.seed(1021) 

x <- rnorm(100) 
y <- 3*x + rnorm(100) 
m1 <- lm(y~x) 
plot(y~x) 
abline(m1$coefficients) 
m2 <- lm(x~y) 
abline(m2$coefficients, col = 'red') 

您正在密謀要將黑的副紅線反之亦然。