X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=pkg%2FR%2FEMGLLF.R;h=bf4476b63d839fc1f260ff7b005309f32e80d4af;hb=43f9cdc196c73bf1074924d4921445d63e89a2f4;hp=03f0a75b47246a9a2d33be68082031267fe4a68c;hpb=ea5860f1b4fc91f06e371a0b26915198474a849d;p=valse.git diff --git a/pkg/R/EMGLLF.R b/pkg/R/EMGLLF.R index 03f0a75..bf4476b 100644 --- a/pkg/R/EMGLLF.R +++ b/pkg/R/EMGLLF.R @@ -179,7 +179,7 @@ EMGLLF <- function(phiInit, rhoInit, piInit, gamInit, mini, maxi, gamma, lambda, sumPen <- sum(pi^gamma * b) last_llh <- llh - llh <- -sumLogLLH/n + lambda * sumPen + llh <- -sumLogLLH/n #+ lambda * sumPen dist <- ifelse(ite == 1, llh, (llh - last_llh)/(1 + abs(llh))) Dist1 <- max((abs(phi - Phi))/(1 + abs(phi))) Dist2 <- max((abs(rho - Rho))/(1 + abs(rho)))