X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=test%2Fgenerate_test_data%2FEMGLLF.R;h=f1eeffa0605bef7afd30eca136155fe82dbcf1c0;hp=41eb7cf6291f7760a603dfff72968968ad1d85a0;hb=923a335ee5493de5ace95d946b7bf93fbedd0377;hpb=f72448157e2ca194959adcffced29448fbb1bb74 diff --git a/test/generate_test_data/EMGLLF.R b/test/generate_test_data/EMGLLF.R index 41eb7cf..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]) + 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)