plotConvergence.nempi {nempi} | R Documentation |
Produces different convergence plots based on a nempi object
## S3 method for class 'nempi' plotConvergence(x, type = "b", ...)
x |
nempi object |
type |
see ?plot.default |
... |
additional parameters for plot |
plot
Martin Pirkl
D <- matrix(rnorm(1000*100), 1000, 100) colnames(D) <- sample(seq_len(5), 100, replace = TRUE) Gamma <- matrix(sample(c(0,1), 5*100, replace = TRUE, p = c(0.9, 0.1)), 5, 100) Gamma <- apply(Gamma, 2, function(x) return(x/sum(x))) Gamma[is.na(Gamma)] <- 0 rownames(Gamma) <- seq_len(5) result <- nempi(D, Gamma = Gamma) par(mfrow=c(2,3)) plotConvergence(result)