plotConvergence.nempi: Plot convergence of EM

Description Usage Arguments Value Author(s) Examples

View source: R/nempi_main.r

Description

Produces different convergence plots based on a nempi object

Usage

1
2
## S3 method for class 'nempi'
plotConvergence(x, type = "b", ...)

Arguments

x

nempi object

type

see ?plot.default

...

additional parameters for plot

Value

plot

Author(s)

Martin Pirkl

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
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)

cbg-ethz/nempi documentation built on May 4, 2021, 9:44 a.m.