R/loglikf.R

loglikf <-
function(y,w,psi) {
  t(log(t(sapply(y, dnorm, psi$mu,psi$sigma)) %*% psi$p)) %*% w   # %*% is matrix product
}

Try the fitTetra package in your browser

Any scripts or data that you put into this service are public.

fitTetra documentation built on May 2, 2019, 3:52 p.m.