黄土线散​​点图,黄土不显示线在某一地区(scatter plot with loess line,

2019-10-16 16:44发布

我绘制这里染色体值超过长度

不加分的中间区域不包含数据,不应该得到一个黄土线。 如何修改我的代码,停止在这一地区的黄土线? 数据是连续的,但我可以添加一些行与一些特殊的价值,以纪念空白区域或标签中添加一列? 但如何在命令中使用呢?

我目前的命令:

library(IDPmisc)

# plot settings (edit here)
spanv<-0.05
pointcol1="#E69F00"
pointcol2="#56B4E9"
pointcol3="#009E73"
points=20
linecol="green"
xlabs=paste(onechr, " position", " (loess-span=", spanv, ")", sep="")

data1<-NaRV.omit(data[,c(2,7)]) # keep only x and y for the relevant data 
                                # and clean NA and Inf
ylabs='E / A - ratio'
p1<-ggplot(data1, aes(x=start, y=E.R)) +
ylim(0,5) +
geom_point(shape=points, col=pointcol1, na.rm=T) +
geom_hline(aes(yintercept=1, col=linecol)) +
geom_smooth(method="loess", span=spanv, fullrange=F, se=T, na.rm=T) +
xlab(xlabs) +
ylab(ylabs)

Answer 1:

我会做一两件事情:

  1. 执行loess()拟合外ggplot()分别预测为两个地区和各组预测添加到情节与自己geom_line()层。
  2. 与上述类似,但在这段时间内ggplot()操作的境界。 添加两层的情节,而不是一个,这两个使用geom_smooth()但重要的改变data提供给每个指只有一个或数据的其他部分的说法。

对于后者,或许是这样的:

....
geom_smooth(data = data[1:n, ], method="loess", span=spanv, fullrange=FALSE, 
            se=TRUE, na.rm=TRUE) +
geom_smooth(data = data[m:k, ], method="loess", span=spanv, fullrange=FALSE, 
            se=TRUE, na.rm=TRUE)
....

其中nmk指的是标记组1的端部和所述开始和组2的端部,并需要由您直接定义或提供的索引。



文章来源: scatter plot with loess line, loess to not show line in a given region