require("VGAMdata")
fit <- vglm(educ ~ sm.bs(age) * sex + ethnicity,
cumulative(parallel = TRUE), head(xs.nz, 200))
is.parallel(fit)
is.parallel(fit, type = "lm") # For each column of the LM matrix
Run the code above in your browser using DataLab