pharsing improved in lmm

aliaksah · aliaksah · commit 014f6eb1caf3 · 2025-12-18T17:36:18.000+01:00
diff --git a/R_script/reproducible_script.R b/R_script/reproducible_script.R
@@ -278,7 +278,7 @@ mixed.model.loglik.lme4 <- function (y, x, model, complex, mlpost_params)
     mm <- lmer(as.formula(paste0("y ~ 1 + (1 | dr)")), data = data, REML = FALSE)
   }
   
-  mloglik <- as.numeric(logLik(mm))  -  0.5*log(length(y)) * (dim(data)[2] - 2) #Laplace approximation for beta prior
+  mloglik <- as.numeric(logLik(mm))  -  0.5*log(length(y)) * (dim(data)[2] - 2) #Laplace approximation
   
   # logarithm of model prior
   if (length(mlpost_params$r) == 0)  mlpost_params$r <- 1/dim(x)[1]  # default value or parameter r

Original file line number	Diff line number	Diff line change
`@@ -278,7 +278,7 @@ mixed.model.loglik.lme4 <- function (y, x, model, complex, mlpost_params)`
`278`	`278`	`mm <- lmer(as.formula(paste0("y ~ 1 + (1 \| dr)")), data = data, REML = FALSE)`
`279`	`279`	`}`
`280`	`280`
`281`		`- mloglik <- as.numeric(logLik(mm)) - 0.5log(length(y)) (dim(data)[2] - 2) #Laplace approximation for beta prior`
	`281`	`+ mloglik <- as.numeric(logLik(mm)) - 0.5log(length(y)) (dim(data)[2] - 2) #Laplace approximation`
`282`	`282`
`283`	`283`	`# logarithm of model prior`
`284`	`284`	`if (length(mlpost_params$r) == 0) mlpost_params$r <- 1/dim(x)[1] # default value or parameter r`