I\'m running a logistic regression with a categorical predictor variable with levels 0-6. By default, R considers level 0 as the reference group.
How can I tell R to
The C
function has already been suggested, also look at contrasts
, relevel
, and reorder
, among others.
You can use relevel function: dataframe$x1 <- relevel(datafrmae$x1, "type here the reference category")
Use the C
function to define your contrasts in the dataframe.
If your dataframe is DF
and the factor variable is fct
, then
DF$fct <- C(DF$fct, contr.treatment, base=3)
(untested).
This is very easy with Frank Harrell's packages (which allows for many additional functions).
For example, for a dataframe named 'df'
library(Hmisc)
library(rms)
dd=datadist(df)
options(datadist='dd')
(m=lrm(y ~ catvar, data=df)) #uses the mode as the reference group
summary(m, catvar=3) #using level=3 as the reference group
summary(m, catvar=0) #using level=0 as the reference group
(And, some sample data would have been nice - you can always use dput
to output a small dataset)