Question

我有两个名为ret和raw的数据框。我单独处理数据帧，以便每个数据帧输出一系列图。每个图表代表不同的染色体。两个数据帧都由相同的染色体分组（但每个数据帧的x轴和y轴测量不同的东西，因此数据不能合并在一起）。我希望输出显示来自数据框ret（称为ret$chr）的染色体图，它与数据框raw（称为raw$Chromosome的同一染色体相邻）。目前我的代码是：

# Plot the bins and sums for each chromosome. (png,pdf,ps,tiff)
ind = unique(as.character(ret$chr))
ind2 = unique(as.character(raw$Chromosome))
png(filename=paste("/Users/blabla/Desktop/SequencingScripts/Plots/temp","png"), width = 600, height = 2000)
par(mfrow=c(8,3), mar = rep(0.5, 4), oma = rep(0.5, 4))

for(jj in ind){
indic = which(ret$chr == jj)
indic2 = which(raw$Chromosome == jj)
plot(ret$retxvalue[indic],ret$retyvalue[indic],pch=19,xlab="Locus",ylab="Summed ZScore",type="s",main=jj,lwd=2, ylim=c(-1,2.5))
plot(raw$rawxvalue[indic2],raw$rawyvalue[indic2],pch=19,xlab="Locus",ylab="Ratio",type="s",main=jj,lwd=2, ylim=c(-1,2.5)xlim=c(0,20000))
}
dev.off()

但是我一直收到错误：

Error in plot.window(...) : need finite 'xlim' values In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf

Answer 1

是复制粘贴错误吗？（似乎你失去了昏迷）

 ylim=c(-1,2.5)xlim=c(0,20000)

在同一个png上绘制来自不同数据帧的相邻相关图

1 个答案: