From a1baafe832473a15c0d4601ecbbaa1d62053b668 Mon Sep 17 00:00:00 2001
From: Layla Riemann <layla.riemann@ptb.de>
Date: Mon, 2 Aug 2021 09:39:30 +0000
Subject: [PATCH] help file for R-code
---
get_lme_sd_v1.r | 18 ++++++++++++++++++
1 file changed, 18 insertions(+)
create mode 100644 get_lme_sd_v1.r
diff --git a/get_lme_sd_v1.r b/get_lme_sd_v1.r
new file mode 100644
index 0000000..a482cd3
--- /dev/null
+++ b/get_lme_sd_v1.r
@@ -0,0 +1,18 @@
+#
+# get_lme_sd_v1.r
+#
+
+get.sd.lme <- function(mod) {
+ sds <- exp(attr(mod$apVar, "Pars"))
+ names(sds) <- c(gsub("reStruct\\.", "", names(sds)[1:2]), "Residual")
+ sds
+}
+
+get.sd.lmer <- function(mod) {
+ vc <- VarCorr(mod)
+ sds <- unlist( lapply(vc, function(v) attr(v, "stddev")) )
+ names(sds) <- names(vc)
+ #NB: SDs are in decreasing order in nested models so ...
+ c( rev(sds), Residual=attr(vc, "sc"))
+}
+
--
GitLab