require("VGAMdata")
fit <- vglm(educ ~ sm.bs(age) * sex + ethnic,
cumulative(parallel = TRUE), xs.nz[1: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