code for issue 222

example code for #222
florianhartig · Jun 3, 2024 · e269198 · e269198
1 parent 5533367
commit e269198
Show file tree

Hide file tree

Showing 2 changed files with 3,027 additions and 0 deletions.
diff --git a/Code/DHARMaIssues/222.R b/Code/DHARMaIssues/222.R
@@ -0,0 +1,21 @@
+dat = read.table("222.txt", header = T, sep = "|")
+head(dat)
+library(lme4)
+mod1=glmer(RT~DISP+(DISP|subj),family=Gamma(link="inverse"),data = dat)
+
+simulationOutput <- simulateResiduals(fittedModel = mod1, plot = T)
+# produces infinity error 
+
+x = predict(mod1,  type = "response")
+min(x)
+max(x)
+
+# unconditional simulations, default in DHARMa
+simulate(mod1)
+simulate(mod1, re.form = ~0)
+
+# conditional simulations 
+simulate(mod1, re.form = NULL)
+
+simulationOutput <- simulateResiduals(fittedModel = mod1, plot = T, re.form = NULL)
+