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 matrixRun the code above in your browser using DataLab