})
logGam <- logGam - max(logGam) #adjust without changing proportions
- gam[i, ] <- exp(logGam)
- norm_fact <- sum(gam[i, ])
- gam[i, ] <- gam[i, ] / norm_fact
+ gam <- exp(logGam)
+ print(gam)
+ norm_fact <- sum(gam)
sumLogLLH <- sumLogLLH + log(norm_fact) - log((2 * base::pi)^(m/2))
}
llhLambda <- c(sumLogLLH/n, (dimension + m + 1) * k - 1)