R/mgsem_output_proc_casewise_likelihood.R

In alexanderrobitzsch/sirt: Supplementary Item Response Theory Models

## File Name: mgsem_output_proc_casewise_likelihood.R
## File Version: 0.08


mgsem_output_proc_casewise_likelihood <- function(data_proc, implied, estimator="ML")
{
    case_ll <- NULL
    if (estimator=='ML' & ( ! is.null(data_proc) ) ){
        requireNamespace('mvtnorm')
        N <- data_proc$N
        G <- data_proc$G
        data <- data_proc$data
        idgroup <- data_proc$idgroup
        case_ll <- rep(NA,N)
        for (gg in 1L:G){
            implied_gg <- implied[[gg]]
            ind_gg <- which(idgroup==gg)
            Mu <- as.vector(implied_gg$Mu[,1])
            Sigma <- implied_gg$Sigma
            y <- mvtnorm::dmvnorm(x=data[ind_gg,], mean=Mu, sigma=Sigma, log=TRUE)
            case_ll[ind_gg] <- y
        }
    }
    #-- output
    return(case_ll)
}