Displaying 2 results from an estimated 2 matches for "rangeto".
Did you mean:
rangeno
2011 Nov 21
1
coverage plot
...ile="~/q20snpref/illusmp454merCbed",sep="\t",header=F)
colnames(data)<-c("Scaffold","sca_position","coverage")
depth<-mean(data[,"coverage'])
#depth now has the mean (overall)coverage
#set the bin-size
window<-1001
rangefrom<-0
rangeto<-length(data[,"sca_position"])
data.1kb<-runmean(data[,"coverage"],k=window)
png(file="cov_1k.png",width=400,height=1000)
plot(x=data[rangefrom:rangeto,"sca_position"],y=data.1kb[rangefrom:rangeto],pch=".",cex1,xlab="bp position",yla...
2011 Nov 26
1
plot xy data
...le="~/q20snpref/illusmp454merbed",sep="\t",header=F)
colnames(data)<-c("Scaffold","sca_position","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=...