2015-07-09 56 views
0

在我的研究领域中,已经确定两个变量-log10(x)和年。在R中通过abline(lm(y〜x))添加回归线会产生奇数结果-log10

我中的R由下面的散点图,以代号:

plot(-log10(LDR2EUR$V5),LDR2EUR$V6,ylab="r2 to rs13169313", xlab="log10(association p-value)",col=ifelse(LDR2EUR$V6==1,'purple',LDR2EUR$V7), pch=20) 

,然后我试图通过添加回归线:

abline(lm(LDR2EUR$V6~-log10(LDR2EUR$V5))) 

enter image description here

然而,线不适合像最合适的线应该的数据。

我想知道如果可怜的线路适合可能与-log10?

由于

cor(LDR2$V6,-log10(LDR2$V5)) 

回报0.9776906这在我看来,一个水平线不应导致,而是类似于直线y = X的直线。

任何指导将不胜感激。

回答

3

这是一个公式的问题,无关LOG10反倒是因为"-" -signs如何式表达式解释:

lm(LDR2EUR$V6~-log10(LDR2EUR$V5) 

....不退步V6反对消极V5而是去除V5从考虑。 (你最终绘制在V6的平均线,而不是尝试。

abline(lm(LDR2EUR$V6 ~ I(-log10(LDR2EUR$V5))) 

这有可能是你真正想要的(如果理论支持通过(0,0)的一行:

abline(lm(LDR2EUR$V6 ~ I(-log10(LDR2EUR$V5) -1)) 

哪也不会从任何东西中减去1,而是移除y截距并强制穿过原点。

+0

是不是必须是'我(-log10(LDR2EUR $ V5))-1' '-1'去除拦截? –

相关问题