df =  read.csv("belief_mixed_effects.csv")
df =  read.csv("belief_mixed_effects.csv")
install.packages("tidyverse")
install.packages("glmmtmb")
install.packages("glmmTMB")
install.packages("easystats")
install.packages("statsr")
df$belief_movement = df$post_belief - df$prior_belief
m = lmer(belief_movement ~ post_stage +  (1|usertoken) ,df)
library(plyr)
library(statsr)
library(lme4)
library(tidyverse)
library(glmmTMB)
library(effects)
library(gtsummary)
theme_set(theme_sjplot())
library(plyr)
library(statsr)
library(lme4)
library(tidyverse)
library(glmmTMB)
library(effects)
library(gtsummary)
install.packages("effects")
install.packages("gtsummary")
install.packages("see")
library(plyr)
library(statsr)
library(lme4)
library(tidyverse)
library(glmmTMB)
library(effects)
library(gtsummary)
m = lmer(belief_movement ~ post_stage +  (1|usertoken) ,df)
m = lmer(belief_movement ~ post_stage +  (1|ID) ,df)
install.packages("effectsize")
library(easystats)
library(parameters)
library(see)
plot(parameters(model))
plot(parameters(m))
m = lmer(belief_movement ~ prior_uncertainty * post_stage +  (1|ID) ,df)
m1 = lmer(belief_movement ~ prior_uncertainty * post_stage +  (1|ID) ,df)
plot(parameters(m1))
df$uncertainty_movement = df$post_uncertainty - df$prior_uncertainty
m2 = lmer(uncertainty_movement ~   post_stage +  (1|ID) ,df)
plot(parameters(m2))
install.packages("scales")
install.packages("scales")
install.packages("scales")
install.packages("scales")
library(scales)
df$belief_movement_bnd = scales::rescale(df$belief_movement, to = c(0.001,0.999))
beta_m = glmmTMB(belief_movement_bnd~ post_stage * v +  (1|ID) ,df,family=list(family="beta", link="logit"))
library(glmmTMB)
beta_m = glmmTMB(belief_movement_bnd~ post_stage * v +  (1|ID) ,df,family=list(family="beta", link="logit"))
beta_m = glmmTMB(belief_movement_bnd~ post_stage +  (1|ID) ,df,family=list(family="beta", link="logit"))
plot(parameters(beta_m))
library(plyr)
library(statsr)
library(lme4)
library(tidyverse)
library(glmmTMB)
library(effects)
library(gtsummary)
library(easystats)
library(parameters)
library(see)
plot(parameters(beta_m))
df$uncertainty_movement_bnd = scales::rescale(df$uncertainty_movement, to = c(0.001,0.999))
beta_m = lmer(uncertainty_movement_bnd ~   post_stage +  (1|ID) ,df,family=list(family="beta", link="logit"))
beta_m = glmmTMB(uncertainty_movement_bnd ~   post_stage +  (1|ID) ,df,family=list(family="beta", link="logit"))
beta_m = glmmTMB(belief_movement_bnd~ post_stage +  (1|ID) ,df,family=list(family="beta", link="logit"))
beta_m1 = glmmTMB(uncertainty_movement_bnd ~   post_stage +  (1|ID) ,df,family=list(family="beta", link="logit"))
plot(parameters(beta_m1))
install.packages("sjPlot")
library(sjPlot)
library(sjlabelled)
library(sjmisc)
library(ggplot2)
library(sjPlot)
library(sjlabelled)
library(sjmisc)
library(ggplot2)
plot_model(beta_m1)
plot_model(beta_m,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
scale_y_continuous(breaks=seq(-.75,0.75,.25)) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Absolute Belief Difference", title = "") +
ylim(-0.25, 0.9)
plot_model(beta_m,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
scale_y_continuous(breaks=seq(-.75,0.75,.25)) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Absolute Belief Difference", title = "")
plot_model(beta_m,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Absolute Belief Difference", title = "")
plot_model(beta_m,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Absolute Belief Difference", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "belief movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "effect of disruption on belief movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m1,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption on belief movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m1,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption on uncertianty movement", title = "")  +
ylim(0.5,1.5)
beta_m2 = glmmTMB(belief_movement_bnd~ prior_uncertainty * post_stage +  (1|ID) ,df,family=list(family="beta", link="logit"))
beta_m3 = glmmTMB(uncertainty_movement_bnd ~ prior_uncertainty * post_stage +  (1|ID) ,df,family=list(family="beta", link="logit"))
plot_model(beta_m2,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption on uncertianty movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m2,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption and prior belief strength on uncertianty movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m2,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption and prior belief strength on belief movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m3,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption and prior belief strength on uncertianty movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m3,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption and prior belief strength on uncertianty movement", title = "")
plot_model(beta_m,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption on belief movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m1,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption on uncertianty movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m2,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption and prior belief strength on belief movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m2,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption and prior uncertainty on belief movement", title = "")  +
ylim(0.5,1.5)
plot_model(beta_m3,show.values = TRUE, vline.color = "grey", value.offset = .4, value.size = 3, type="est", show.intercept = TRUE ) +
theme(axis.text.y = element_text(size = 8),
plot.subtitle=element_text(size=11), plot.title = element_text(size = 1)) +
labs(subtitle = "Effect of disruption and prior uncertainty on uncertianty movement", title = "")
summary(beta_m1)
confint(beta_m1)
df %>%
ggplot(aes(x=belief_movement,fill=post_stage)) +
geom_density(alpha=0.5) +
theme(legend.position = "none") +
labs(x = " ", y = " ", subtitle = "belief_movement") +
facet_grid(~variable)
df %>%
ggplot(aes(x=belief_movement,fill=post_stage)) +
geom_density(alpha=0.3) +
theme(legend.position = "none") +
labs(x = " ", y = " ", subtitle = "belief_movement") +
facet_grid(~variable)
df %>%
ggplot(aes(x=belief_movement,fill=post_stage)) +
geom_density(alpha=0.3) +
labs(x = " ", y = " ", subtitle = "belief_movement") +
facet_grid(~variable)
df %>%
ggplot(aes(x=belief_movement,fill=post_stage)) +
geom_density(alpha=0.3) +
labs(x = " ", y = " ", subtitle = "belief movement (posterior - prior)") +
facet_grid(~variable)
