## File Name: mgsem_proc_data.R
## File Version: 0.057
mgsem_proc_data <- function(data, group, weights)
{
N <- nrow(data)
if (is.null(group)){
group <- rep(1,N)
}
groups <- unique(group)
G <- length(groups)
if (is.null(weights)){
weights <- rep(1,N)
}
idgroup <- match(group, groups)
suffstat <- list()
for (gg in 1L:G){
ind_gg <- which(group==groups[gg])
dat_gg <- data[ ind_gg, ]
w <- weights[ind_gg]
res <- stats::cov.wt( x=dat_gg, w=w, method='ML')
s_gg <- list(M=res$center, S=res$cov, N=sum(w) )
suffstat[[gg]] <- s_gg
}
#-- output
res <- list(suffstat=suffstat, weights=weights, groups=groups,
G=G, N=N, idgroup=idgroup, data=data)
return(res)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.