R/getCombiCoefs.R

In benjaminschlegel/brant: Test for Parallel Regression Assumption

getCombiCoefs <- function(model){
  classes = attr(model$terms,"dataClasses")
  factors = !(classes[2:length(classes)] %in% c("numeric", "logical"))
  f = i = var = 1
  result = data.frame(i=1:length(coef(model)),var=NA)
  for(factor in factors){
    if(factor){
      n = length(unlist(model$xlevels[f]))
      for(j in 1:(n-1)){
        result[i,"var"] = var
        i = i + 1
      }
      var = var + 1
      f = f + 1
    }else{
      result[i,"var"] = var
      var = var + 1
      i = i + 1
    }
  }
  return(result)
}