search for: cov_10k

Displaying 1 result from an estimated 1 matches for "cov_10k".

Did you mean: cov210
2011 Nov 26
1
plot xy data
...osition","coverage") depth<-mean(data[,"coverage"]) #depth now has the mean (overall)coverage #set the bin-size window<-10001 rangefrom<-0 rangeto<-length(data[,"sca_position"]) data.10kb<-runmed(data[,"coverage"],k=window) png(file="cov_10k.png",width=1000,height=1000) plot(x=data.10kb[rangefrom:rangeto],y=data[rangefrom:rangeto,"sca_position"],pch=".",cex=1,xlab="depth",ylab="bp_position",type="p") dev.off() Best regards, Sutada