Иногда кажется, что это происходит, когда я пытаюсь построить графики.
Error in (Midpoint.NO.TEL1.1[k, 2] - 250):(Midpoint.NO.TEL1.1[k, 2] + : NA/NaN argument
Как так получается, что иногда программа может строить графики с использованием точно такого же кода, а иногда нет?Как я могу решить это?
attach(ChromoSizes)
NO.TEL1.X1<-c(1:ChrSizes[1])
NO.TEL1.Y1<-rep(0,ChrSizes[1])
attach(FullMap.TEL1)
results <- matrix(0, nrow=1182, ncol=501)
for (k in 1:1182) {
results[k,1:501] <- NO.TEL1.Y1[(Midpoint.NO.TEL1.1[k,2]-250) : (Midpoint.NO.TEL1.1[k,2]+250)]
}
pileup <- colMeans(results)
plot(pileup, type = "l",col="dark blue",lwd=2,xlab ="Distance relative to DSB (bp)",ylab = "Nucleosome Occupacy")