Bic:

Usage Arguments Examples

View source: R/Bic.R

Usage

1
Bic(f, Y = NULL, nbr_par = NULL, Xquali = NULL, Xquanti = NULL, K = length(unique(Y)))

Arguments

f
Y
nbr_par
Xquali
Xquanti
K

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (f, Y = NULL, nbr_par = NULL, Xquali = NULL, Xquanti = NULL, 
    K = length(unique(Y))) 
{
    nbr_par <- 0
    if (!is.null(Xquali)) {
        p <- ncol(Xquali)
        for (j in 1:p) {
            nbr_par <- nbr_par + length(unique(Xquali[, j])) - 
                1
        }
        nbr_par <- nbr_par * K
    }
    if (!is.null(Xquanti)) {
        p <- ncol(Xquanti)
        nbr_par <- K * (p * (p + 1)/2 + p + 1) - 1 + nbr_par
    }
    return(-2 * sum(log(f)) + nbr_par * log(length(Y)))
  }

SaraTouzani/EM documentation built on Dec. 16, 2019, 12:31 a.m.