我正在尝试计算落在 1000 个窗口内的条目,问题是我正在使用 for 循环,这使得需要执行的操作数量非常大(我对 R 相当陌生)并且我得到了的边界错误。我知道必须有更好的方法来做到这一点。
文件(警告文件超过 100mb):bamDF.txt
采用:
dget(file="bamDF.txt")
脚本:
attach(bamDF)
out <- matrix(0,1,ceiling((max(pos, na.rm=TRUE)-min(pos, na.rm=TRUE))/interval))
interval <- 1000
for(q in 1:nrow(bamDF)){
for(z in 1:ceiling((max(pos, na.rm=TRUE)-min(pos, na.rm=TRUE))/interval)){
if(min(pos, na.rm=TRUE)+interval*(z-1)<pos[q]&&pos[q]<(min(pos, na.rm=TRUE)+interval*(z))){
out[z,] <- out[z,]+1;
}
}
}
detach(bamDF)