vignettes/MSFA.R

In rdevito/MSFA: Maximum Likelihood Estimation and Bayesian analysis of MSFA models

## ----help, echo = TRUE, results = TRUE, tidy = TRUE----------------------
library(MSFA)
data(data_immune)
help(data_immune)

## ---- starting values, messages = FALSE----------------------------------
start_value <- start_msfa(X_s = data_immune, k = 3, j_s = c(3, 4))

## ---- get estimate, results = FALSE--------------------------------------
mle <-  ecm_msfa(data_immune, start_value, trace = FALSE)

## ---- heatmap, fig.show = 'hold', fig.width = 7.5, fig.height = 6.5, message = FALSE----
library(gplots)
heatmap.2(mle$Phi,dendrogram='none', Rowv=FALSE, Colv=FALSE,trace='none', density.info="none", col=heat.colors(256))