aboutsummaryrefslogtreecommitdiff
path: root/man/summary.saem.mmkin.Rd
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2022-12-15 14:50:28 +0100
committerJohannes Ranke <jranke@uni-bremen.de>2022-12-15 14:50:28 +0100
commita54bd290bc3884d0000c52c1b29bc557825d9eae (patch)
tree841dc41b03035677aac0b288bca4927a0bd18582 /man/summary.saem.mmkin.Rd
parent7c8b3db180be371cb03b9518f14117060fbc4239 (diff)
List random effects correlations in output if any
Update docs
Diffstat (limited to 'man/summary.saem.mmkin.Rd')
-rw-r--r--man/summary.saem.mmkin.Rd13
1 files changed, 12 insertions, 1 deletions
diff --git a/man/summary.saem.mmkin.Rd b/man/summary.saem.mmkin.Rd
index fb099899..0845d4d2 100644
--- a/man/summary.saem.mmkin.Rd
+++ b/man/summary.saem.mmkin.Rd
@@ -92,10 +92,21 @@ f_mmkin_dfop_sfo <- mmkin(list(dfop_sfo), ds_syn_dfop_sfo,
f_saem_dfop_sfo <- saem(f_mmkin_dfop_sfo)
print(f_saem_dfop_sfo)
illparms(f_saem_dfop_sfo)
-f_saem_dfop_sfo_2 <- update(f_saem_dfop_sfo, covariance.model = diag(c(0, 0, 1, 1, 1, 0)))
+f_saem_dfop_sfo_2 <- update(f_saem_dfop_sfo,
+ no_random_effect = c("parent_0", "log_k_m1"))
illparms(f_saem_dfop_sfo_2)
intervals(f_saem_dfop_sfo_2)
summary(f_saem_dfop_sfo_2, data = TRUE)
+# Add a correlation between random effects of g and k2
+cov_model_3 <- f_saem_dfop_sfo_2$so@model@covariance.model
+cov_model_3["log_k2", "g_qlogis"] <- 1
+cov_model_3["g_qlogis", "log_k2"] <- 1
+f_saem_dfop_sfo_3 <- update(f_saem_dfop_sfo,
+ covariance.model = cov_model_3)
+intervals(f_saem_dfop_sfo_3)
+# The correlation does not improve the fit judged by AIC and BIC, although
+# the likelihood is higher with the additional parameter
+anova(f_saem_dfop_sfo, f_saem_dfop_sfo_2, f_saem_dfop_sfo_3)
}
}

Contact - Imprint