From 9d1bc874713e9713e1173db080762e840726511d Mon Sep 17 00:00:00 2001 From: Steve Phelps Date: Sat, 9 Dec 2023 17:53:20 +0000 Subject: [PATCH] residuals analysis --- jupyter-book/Anova-r.py | 36 +++++++++++++++++++++++++++++++----- 1 file changed, 31 insertions(+), 5 deletions(-) diff --git a/jupyter-book/Anova-r.py b/jupyter-book/Anova-r.py index 2a44c4b..a73a956 100644 --- a/jupyter-book/Anova-r.py +++ b/jupyter-book/Anova-r.py @@ -106,14 +106,40 @@ results.clean$Cooperation_frequency[results.clean$Cooperation_frequency == 1] <- 1 - epsilon # %% -model <- glmmTMB(Cooperation_frequency ~ +results.clean$Num_cooperates = round(results.clean$Cooperation_frequency * 6) + +# %% +6 - results.clean$Num_cooperates + +# %% +model <- glmmTMB(cbind(Num_cooperates, 6 - Num_cooperates) ~ + Participant_group + Partner_condition + t + Model + Temperature + + Partner_condition:Model + Participant_group:Model + + (1|Participant_id), + data = results.clean, + family = betabinomial) +summary(model) + +# %% +model <- glmmTMB(cbind(Num_cooperates, 6 - Num_cooperates) ~ Participant_group + Partner_condition + t + Model + Temperature + Partner_condition:Model + Participant_group:Model + Participant_labels_reversed:Participant_label + Participant_labels_reversed:Model + - Participant_label + Participant_chain_of_thought + Participant_pronoun + Participant_defect_first + Participant_labels_reversed + - (1 | Participant_id), + Participant_label + Participant_chain_of_thought + Participant_pronoun + Participant_defect_first + Participant_labels_reversed, data = results.clean, - family = beta_family(link = "logit")) + family = betabinomial) summary(model) # %% -simulationOutput <- simulateResiduals(fittedModel = model, plot = TRUE) +simulationOutput <- simulateResiduals(fittedModel = model, plot = TRUE, integerResponse=TRUE) + +# %% +options(repr.plot.width = 20, repr.plot.height = 10) + +# %% +testQuantiles(simulationOutput) + +# %% +hist(residuals(simulationOutput)) + +# %% +plotResiduals(simulationOutput, results.clean$Participant_group)