2014-02-14 42 views
0

我试图用一个连续的预测变量在一个二元结果变量的散点图上拟合黄土线。在二进制数据中拟合R中的黄土线

这里是我使用的代码:

lw1<-loess(y~x, data=df) 
    plot(y~x, data=df, pch=19, cex=0.1) 
    lines(df$x, lw1$fitted, col='blue') 

,这是图我得到: ! https://bitbucket.org/heatherjbaldwin/akos_open/src/ec2a78d093e6cdf988434c03c1b7c7df145892ba/loessgraph1.png?at=master

我也试着订购x变量:

j<-order(df$x) 
    lines(df$x[j], lw1$fitted, col='blue') 

并获得该图: ! https://bitbucket.org/heatherjbaldwin/akos_open/src/ec2a78d093e6cdf988434c03c1b7c7df145892ba/loessgraph2%28ordered_x%29.png?at=master

下面是数据:https://bitbucket.org/heatherjbaldwin/akos_open/src/ec2a78d093e6cdf988434c03c1b7c7df145892ba/loesscurvedata.txt?at=master

任何帮助深表感谢。

回答

1

使用ggplot2loess平滑我得到这个:

ggplot(data=dat,aes(x,y)) + 
    geom_line() + 
    geom_smooth(method='loess') 

enter image description here

但我认为你是正在寻找分类这里。