Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # Generate data
- data1 <- data.frame(matrix(NA,200,2))
- data1[,1] <- c(1,5)
- data1[,2] <- 1
- data2 <- data.frame(matrix(NA,200,2))
- data2[,1] <- 3
- data2[,2] <- 2
- data3 <- data.frame(matrix(NA,200,2))
- data3[,1] <- c(2,4)
- data3[,2] <- 3
- data1 <- rbind(data1,data2,data3)
- names(data1) <- c("rating","bottle")
- data1$rating <- ordered(data1$rating)
- data1$bottle <- as.factor(data1$bottle)
- # Model data
- wineM <- brm(rating ~ bottle, data = data1,
- family = cumulative(threshold = "flexible"),
- chains = 2, cores = 2, iter=4000)
- # Plot Average
- conditional_effects(wineM, categorical = F)
- # Plot Categorical
- conditional_effects(wineM, categorical = T)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement