## Data set: pss (for description run ?pss)
barplot(table(pss$difference), xlab = "PSS index difference", ylab = "Frequency")
boxplot(pss$difference ~ pss$group, xlab = "Group", ylab = "PSS index difference")
## Fit a double model (mode = 1)
fit0 <- sdlrm(difference ~ group | group, data = pss, xi = 1)
## Constant dispersion test
disp_test(fit0)
Run the code above in your browser using DataLab