R/plotdiversity.R

Defines functions plotdiversity

Documented in plotdiversity

plotdiversity <-
function(X, sample.times, makeplot=TRUE, filter=FALSE, ...) {
  if (filter) {
    seqn <- sample.times
  } else {
    seqn <- 1:length(sample.times)
  }
  diversity <- numeric(length(sample.times))
  k <- 1
  for (i in seqn) {
    diversity[k] <- meansnps(X$obs.strain[[i]], X$obs.freq[[i]], X$libr, X$nuc, X$librstrains)
    k <- k+1
  }
  if (makeplot) {
    plot(sample.times, diversity, type="l", ...)
  }
  return(invisible(diversity))
}

Try the seedy package in your browser

Any scripts or data that you put into this service are public.

seedy documentation built on May 29, 2017, 10:58 a.m.