##
sequence <- unlist(lapply(1:50000, function(x) {paste(sample(c("A","T","C","G"),prob=c(0.3,0.3,0.2,0.2),25,replace=TRUE), collapse="")}))
values <- runif(50000,min=-2,max=2)
plotGCbias(values, sequence)
Run the code above in your browser using DataLab