# install.packages('pacman')
library(pacman)
p_load('betareg', 'tidyverse', 'tibble', 'datapasta')
df<-tibble::tribble(
~target, ~value, ~pid, ~subtreatment, ~age, ~gender,
"male", 65, 2L, "EI (corr. shown)", 3, 1,
"male", 75, 2L, "EI (corr. shown)", 1, 0,
"female", 70, 1L, "EI (corr. shown)", 3, 0,
"female", 70, 1L, "EI (corr. shown)", 3, 0,
"male", 25, 2L, "EI (corr. shown)", 2, 0,
"female", 55, 1L, "EI (corr. shown)", 3, 0,
"female", 70, 1L, "EI (corr. shown)", 4, 1,
"female", 60, 1L, "EI (corr. shown)", 4, 1,
"male", 50, 2L, "EI (corr. shown)", 2, 1,
"female", 100, 1L, "EI (corr. shown)", 1, 0,
"female", 60, 1L, "EI (corr. shown)", 2, 1,
"female", 90, 1L, "EI (corr. shown)", 2, 0,
"female", 70, 1L, "EI (corr. shown)", 4, 1,
"female", 80, 1L, "EI (corr. shown)", 4, 1,
"female", 80, 1L, "EI (corr. shown)", 5, 0,
"female", 70, 1L, "EI (corr. shown)", 2, 1,
"male", 51, 2L, "EI (corr. shown)", 2, 0,
"male", 50, 2L, "EI (corr. shown)", 2, 1,
"male", 83, 2L, "EI (corr. shown)", 2, 0,
"male", 80, 2L, "EI (corr. shown)", 2, 1,
"female", 80, 1L, "EI (corr. shown)", 1, 0,
"female", 65, 1L, "EI (corr. shown)", 4, 0,
"female", 70, 1L, "EI (corr. shown)", 2, 1,
"male", 45, 2L, "EI (corr. shown)", 4, 0,
"female", 60, 1L, "EI (corr. shown)", 3, 0,
"female", 80, 1L, "EI (corr. shown)", 1, 0,
"female", 85, 1L, "EI (corr. shown)", 4, 0,
"male", 77, 2L, "EI (corr. shown)", 2, 1,
"female", 25, 1L, "EI (corr. shown)", 2, 0,
"female", 90, 1L, "EI (corr. shown)", 2, 0,
"male", 70, 2L, "EI (corr. shown)", 1, 0,
"female", 45, 1L, "EI (corr. shown)", 4, 0,
"female", 90, 1L, "EI (corr. shown)", 3, 1,
"female", 75, 1L, "EI (corr. shown)", 2, 1,
"male", 32, 2L, "EI (corr. shown)", 5, 1,
"male", 65, 2L, "EI (corr. shown)", 2, 0,
"female", 90, 1L, "EI (corr. shown)", 3, 0,
"female", 66, 1L, "EI (corr. shown)", 2, 0,
"female", 80, 1L, "EI (corr. shown)", 3, 1,
"male", 80, 2L, "EI (corr. shown)", 4, 1,
"male", 30, 2L, "EI (corr. shown)", 3, 0,
"female", 85, 1L, "EI (corr. shown)", 1, 0,
"male", 76, 2L, "EI (corr. shown)", 2, 0,
"male", 50, 2L, "EI (corr. shown)", 2, 1,
"female", 80, 1L, "EI (corr. shown)", 3, 1,
"male", 50, 2L, "EI (corr. shown)", 2, 0,
"female", 65, 1L, "EI (corr. shown)", 2, 0,
"male", 48, 2L, "EI (corr. shown)", 4, 0,
"male", 45, 2L, "EI (corr. shown)", 5, 0,
"male", 40, 2L, "EI (corr. shown)", 4, 1,
"male", 72, 2L, "EI (corr. not shown)", 4, 0,
"male", 85, 2L, "EI (corr. not shown)", 3, 1,
"male", 70, 2L, "EI (corr. not shown)", 2, 0,
"male", 70, 2L, "EI (corr. not shown)", 3, 1,
"female", 89, 1L, "EI (corr. not shown)", 2, 0,
"male", 90, 2L, "EI (corr. not shown)", 1, 1,
"female", 45, 1L, "EI (corr. not shown)", 0, 0,
"male", 90, 2L, "EI (corr. not shown)", 2, 0,
"female", 90, 1L, "EI (corr. not shown)", 1, 1,
"female", 70, 1L, "EI (corr. not shown)", 1, 1,
"female", 70, 1L, "EI (corr. not shown)", 2, 0,
"female", 30, 1L, "EI (corr. not shown)", 2, 1,
"male", 80, 2L, "EI (corr. not shown)", 4, 1,
"male", 90, 2L, "EI (corr. not shown)", 1, 1,
"female", 80, 1L, "EI (corr. not shown)", 3, 1,
"female", 80, 1L, "EI (corr. not shown)", 2, 0,
"female", 62, 1L, "EI (corr. not shown)", 2, 1,
"female", 70, 1L, "EI (corr. not shown)", 2, 1,
"female", 80, 1L, "EI (corr. not shown)", 3, 1,
"female", 55, 1L, "EI (corr. not shown)", 2, 1,
"female", 75, 1L, "EI (corr. not shown)", 2, 0,
"female", 70, 1L, "EI (corr. not shown)", 1, 1,
"female", 70, 1L, "EI (corr. not shown)", 2, 0,
"male", 85, 2L, "EI (corr. not shown)", 3, 1,
"male", 95, 2L, "EI (corr. not shown)", 5, 1,
"male", 40, 2L, "EI (corr. not shown)", 2, 0,
"female", 85, 1L, "EI (corr. not shown)", 3, 0,
"male", 80, 2L, "EI (corr. not shown)", 4, 0,
"female", 50, 1L, "EI (corr. not shown)", 3, 1,
"female", 57, 1L, "EI (corr. not shown)", 6, 1,
"female", 10, 1L, "EI (corr. not shown)", 2, 1,
"female", 60, 1L, "EI (corr. not shown)", 2, 1,
"female", 70, 1L, "EI (corr. not shown)", 4, 0,
"female", 60, 1L, "EI (corr. not shown)", 3, 1,
"female", 25, 1L, "EI (corr. not shown)", 2, 0,
"female", 90, 1L, "EI (corr. not shown)", 5, 1,
"male", 70, 2L, "EI (corr. not shown)", 3, 0,
"male", 90, 2L, "EI (corr. not shown)", 3, 0,
"female", 80, 1L, "EI (corr. not shown)", 3, 1,
"male", 65, 2L, "EI (corr. not shown)", 3, 0,
"male", 50, 2L, "EI (corr. not shown)", 3, 1,
"male", 40, 2L, "EI (corr. not shown)", 3, 0,
"male", 70, 2L, "EI (corr. not shown)", 4, 1,
"male", 52, 2L, "EI (corr. not shown)", 2, 1,
"female", 65, 1L, "EI (corr. not shown)", 3, 0,
"male", 80, 2L, "EI (corr. not shown)", 2, 0,
"male", 40, 2L, "EI (corr. not shown)", 2, 1,
"female", 15, 1L, "EI (corr. not shown)", 2, 1,
"female", 95, 1L, "EI (corr. not shown)", 5, 1,
"male", 90, 2L, "EI (corr. not shown)", 4, 1,
"male", 90, 2L, "FIC", 1, 1,
"male", 97, 2L, "FIC", 2, 1,
"male", 30, 2L, "FIC", 3, 1,
"male", 90, 2L, "FIC", 2, 1,
"male", 75, 2L, "FIC", 2, 0,
"female", 90, 1L, "FIC", 2, 1,
"male", 60, 2L, "FIC", 5, 0,
"male", 70, 2L, "FIC", 2, 1,
"male", 70, 2L, "FIC", 2, 1,
"male", 90, 2L, "FIC", 4, 0,
"male", 90, 2L, "FIC", 3, 0,
"female", 90, 1L, "FIC", 3, 0,
"male", 100, 2L, "FIC", 2, 0,
"male", 81, 2L, "FIC", 2, 0,
"female", 60, 1L, "FIC", 1, 0,
"male", 50, 2L, "FIC", 3, 0,
"female", 0, 1L, "FIC", 3, 0,
"female", 75, 1L, "FIC", 3, 0,
"female", 75, 1L, "FIC", 1, 0,
"male", 60, 2L, "FIC", 3, 0,
"female", 90, 1L, "FIC", 2, 0,
"female", 50, 1L, "FIC", 2, 0,
"female", 50, 1L, "FIC", 4, 0,
"male", 50, 2L, "FIC", 2, 1,
"male", 90, 2L, "FIC", 2, 1,
"female", 71, 1L, "FIC", 3, 0,
"male", 99, 2L, "FIC", 5, 0,
"female", 70, 1L, "FIC", 4, 1,
"male", 100, 2L, "FIC", 3, 0,
"female", 100, 1L, "FIC", 3, 0,
"male", 75, 2L, "FIC", 3, 0,
"male", 55, 2L, "FIC", 2, 0,
"male", 91, 2L, "FIC", 2, 1,
"male", 50, 2L, "FIC", 5, 1,
"female", 80, 1L, "FIC", 4, 0,
"male", 80, 2L, "FIC", 3, 0,
"male", 80, 2L, "FIC", 4, 0,
"female", 90, 1L, "FIC", 5, 0,
"female", 70, 1L, "FIC", 3, 1,
"female", 80, 1L, "FIC", 1, 1,
"male", 75, 2L, "FIC", 3, 1,
"female", 90, 1L, "FIC", 2, 0,
"male", 100, 2L, "FIC", 2, 0,
"female", 80, 1L, "FIC", 3, 1,
"male", 65, 2L, "FIC", 2, 0,
"male", 60, 2L, "FIC", 2, 1,
"female", 95, 1L, "FIC", 4, 1,
"male", 90, 2L, "FIC", 2, 1,
"male", 50, 2L, "FIC", NA, NA,
"male", 60, 2L, "FIC", 3, 1,
"male", 60, 2L, "FIN", 3, 0,
"female", 80, 1L, "FIN", 3, 1,
"female", 90, 1L, "FIN", 1, 0,
"female", 100, 1L, "FIN", 3, 0,
"female", 30, 1L, "FIN", 3, 1,
"female", 80, 1L, "FIN", 3, 0,
"male", 85, 2L, "FIN", 5, 0,
"male", 70, 2L, "FIN", 4, 1,
"male", 90, 2L, "FIN", 3, 0,
"female", 70, 1L, "FIN", 1, 0,
"female", 50, 1L, "FIN", 2, 0,
"male", 30, 2L, "FIN", 1, 1,
"male", 90, 2L, "FIN", 3, 1,
"male", 90, 2L, "FIN", 2, 1,
"male", 55, 2L, "FIN", 3, 0,
"female", 90, 1L, "FIN", 2, 1,
"female", 70, 1L, "FIN", 2, 0,
"female", 65, 1L, "FIN", 1, 1,
"male", 70, 2L, "FIN", 1, 0,
"male", 100, 2L, "FIN", 3, 0,
"male", 100, 2L, "FIN", 3, 0,
"female", 70, 1L, "FIN", 4, 1,
"female", 60, 1L, "FIN", 1, 1,
"male", 75, 2L, "FIN", 3, 1,
"male", 80, 2L, "FIN", 2, 1,
"female", 90, 1L, "FIN", 1, 0,
"female", 37, 1L, "FIN", 2, 1,
"male", 70, 2L, "FIN", 3, 0,
"female", 80, 1L, "FIN", 4, 0,
"male", 70, 2L, "FIN", 2, 0,
"male", 95, 2L, "FIN", 5, 1,
"male", 60, 2L, "FIN", 3, 1,
"female", 30, 1L, "FIN", 2, 1,
"male", 80, 2L, "FIN", 3, 1,
"male", 40, 2L, "FIN", 2, 1,
"male", 68, 2L, "FIN", 2, 0,
"male", 80, 2L, "FIN", 3, 1,
"male", 60, 2L, "FIN", 3, 0,
"female", 85, 1L, "FIN", 2, 1,
"male", 53, 2L, "FIN", 2, 1,
"female", 50, 1L, "FIN", 3, 1,
"female", 55, 1L, "FIN", 3, 0,
"male", 70, 2L, "FIN", 3, 1,
"male", 60, 2L, "FIN", 3, 1,
"female", 80, 1L, "FIN", 3, 1,
"male", 70, 2L, "FIN", 3, 1,
"female", 90, 1L, "FIN", 3, 0,
"female", 100, 1L, "FIN", 2, 1,
"female", 60, 1L, "FIN", 3, 0,
"female", 50, 1L, "FIN", 2, 1
)temp<- df %>% mutate(value=(value+1)/102)
# m1 <- betareg(value~target*subtreatment +(1|pid),data=temp)
# m2 <- glmmTMB(value~target*subtreatment +(1|pid), data=temp,family = beta_family("logit"),
# REML = T)
# summary(m1)
# summary(m2)