Nothing
options(digits = 4)
suppressWarnings(RNGversion("3.5.0"))
## package and data
library("betareg")
data("ReadingSkills", package = "betareg")
## beta regression mixture model
set.seed(1071)
rs_mix <- betamix(accuracy ~ iq, data = ReadingSkills, k = 3,
extra_components = extraComponent(type = "uniform",
coef = 0.99, delta = 0.01), nstart = 10)
## fitted model
print(rs_mix)
summary(rs_mix)
## further methods
table(clusters(rs_mix), ReadingSkills$dyslexia)
round(posterior(rs_mix), digits = 4)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.