Thank you! But i'm wondering:
if you run area <- factor(area, levels=c("C", "A","B","D","E","F") ) then you are transforming only 'area', aren't you? isn't it possible to transform the whole data like i did for agecat but now for area and area C as baseline, or are you doing so when you run > area <- factor(area, levels=c("C", "A","B","D","E","F") ) > attach(data) and then run the model with area as predictorvariable: > model <- glm(Y~ agecat+gender+area,...) My question is if i can run it as follows and still have a right solution : > data <-transform(data, area=(factor(area, levels=c("C", > "A","B","D","E","F") ) ) I'll be very grateful for any help you can provide! Kind regards, Louisa -- View this message in context: http://r.789695.n4.nabble.com/Inverse-Gaussian-Distribution-tp3172533p3173468.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.