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