update to print.summary

R/summary.galamm.R

@@ -147,7 +147,10 @@ print.summary.galamm <- function(
   cat("\n")
   lme4::.prt.aictab(x$AICtab)
   cat("\n")
-  lme4::.prt.resids(residuals(x) / sigma(x), digits = digits)
+  if(length(x$model$family) == 1) {
+    lme4::.prt.resids(residuals(x) / sigma(x), digits = digits)
+  }
+
   if (exists("Lambda", x)) {
     cat("Lambda:\n")
     x$Lambda[x$Lambda == 0] <- NA

dev-scripts/lifespan_demo.R

@@ -26,8 +26,11 @@ mod <- galamm(
   load.var = "test",
   lambda = lmat,
   factor = c("epmem_ability", "wmem_ability", "execfun_ability"),
+  start = list(
+    theta = mod_init$parameters$parameter_estimates[mod_init$parameters$theta_inds]
+  ),
   control = galamm_control(
-    optim_control = list(REPORT = 2, factr = 1e9, trace = 3, maxit = 10))
+    optim_control = list(REPORT = 2, factr = 1e8, trace = 3, maxit = 10))
 )
 beepr::beep()
 summary(mod)

news.md

@@ -1,5 +1,8 @@
 # galamm (development versions)
 
+- The print.summary.galamm() method does no longer print residual percentiles 
+  with mixed response models. The version implemented until now had a bug, and 
+  it is not clear how to present this information in a useful way.
 - Fixed bug causing galamm() to fail with mixed response types, when the first 
   argument to "family" was not "gaussian".
 - BREAKING CHANGE: argument "factor", "factor_interaction" and "lambda" to 

tests/testthat/_snaps/galamm-mixed-resp.md

@@ -24,10 +24,6 @@
            AIC      BIC   logLik deviance df.resid 
        12356.9  12388.3  -6173.4  31426.4     3995 
       
-      Scaled residuals: 
-         Min     1Q Median     3Q    Max 
-      -15.23  -0.40   0.13   0.68  16.03 
-      
       Random effects:
        Groups Name        Variance Std.Dev.
        id     (Intercept) 0        0