# Divide units into clusters
clusters = rep(1:5, 10)
# Default probability 0.5, default ICC 0
draw_binary_icc(clusters = clusters)
# Specify probability or ICC
corr_draw = draw_binary_icc(prob = 0.5, clusters = clusters, ICC = 0.5)
# Verify ICC of data.
summary(lm(corr_draw ~ as.factor(clusters)))$r.squared
Run the code above in your browser using DataLab