我想你可能会扭转T的顺序他在公式中的论点,因为y与x的常规图将从xyplot(y ~ x)
返回。但是,如果你真的想用一个二次多项式拟合的线性回归,然后用公式中的poly
功能和predict
功能得到拟合值:
xyplot(x~y,data=vals)
plot(vals$y, predict(lm(x~poly(y, 2), data=vals)))
,因为你是这不是一个精确配合绘制x对x的平方根。
resid( lm(x~poly(y, 2),data=vals))
#------------------------------
1 2 3 4 5 6 7
-0.57799840 -0.01438720 0.28444720 0.35772841 0.26036947 0.06297286 -0.14816950
8 9 10
-0.27107622 -0.18807651 0.23418989
另一种方式(它给出一个确切的FIT)是使用公式内I
功能:
plot(vals$y, fitted(lm(x~I(y^(1/2)), data=vals)))
points(vals, pch="x", col='red')
points(vals$y, vals$x, pch="x", col='red')

你需要使用绘图或曲线行帮助最合适? –
绘制线....出于某种原因,我尝试过的语法不起作用(我发誓我在问之前阅读文档!) – Rilcon42
您是使用'ggplot2'还是其他库? –