Mll {epiNEM} | R Documentation |
Computes marginal log-likelihood for model Phi given observed data matrix D1
Mll(Phi, D1, D0, ltype = "marginal", para = c(0.13, 0.05))
Phi |
model to be evaluated |
D1 |
observed data matrix |
D0 |
complementary D1 |
ltype |
likelihood type either "marginal" or "maximum" |
para |
false positive and false negative rates |
list with likelihood poster probability, egene positions
Phi <- matrix(sample(c(0,1), 9, replace = TRUE), 3, 3) data <- matrix(sample(c(0,1), 3*10, replace = TRUE), 10, 3) rownames(Phi) <- colnames(Phi) <- colnames(data) <- c("Ikk1", "Ikk2", "RelA") score <- Mll(Phi, D1 <- data, D0 <- 1 - data)