X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=test%2Fgenerate_test_data%2FEMGLLF.R;h=f1eeffa0605bef7afd30eca136155fe82dbcf1c0;hp=fc3220f55ee7f1021cd1aa19cfb253dc07ecc9a8;hb=923a335ee5493de5ace95d946b7bf93fbedd0377;hpb=2652fb8797f3097361ba23c7f3790be9754db987 diff --git a/test/generate_test_data/EMGLLF.R b/test/generate_test_data/EMGLLF.R index fc3220f..f1eeffa 100644 --- a/test/generate_test_data/EMGLLF.R +++ b/test/generate_test_data/EMGLLF.R @@ -128,7 +128,7 @@ EMGLLF_R = function(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,lambda,X,Y,ta sumLLF1 = 0.0; for (r in 1:k) { - Gam[i,r] = pi[r] * exp(-0.5*sqNorm2[r]) #* det(rho[,,r]) #FIXME: still issues here ?!?! + Gam[i,r] = pi[r] * exp(-0.5*sqNorm2[r]) * det(rho[,,r]) #FIXME: still issues here ?!?! sumLLF1 = sumLLF1 + Gam[i,r] / (2*base::pi)^(m/2) } sumLogLLF2 = sumLogLLF2 + log(sumLLF1)