On Tue, Dec 7, 2010 at 7:54 AM, Gabor Grothendieck <ggrothendi...@gmail.com> wrote: > 'On Tue, Dec 7, 2010 at 5:19 AM, mxkuhn <mxk...@gmail.com> wrote: >> Greg and Frank, >> >> Thanks for the replies. I didn't express myself very well; I'm not interest >> in the model fitting aspect. I'd just like to get the full set of dummy >> variables (optimally from model.matrix) >> > > Try this: > > levels(when$time) <- c("morning", "afternoon", "night") > levels(when$day) <- c("Mon", "Tue", "Wed", "Thu", "Fri", "Sat", "Sun") > > contr.dummy <- function(n, ...) diag(n) > mm <- suppressWarnings(model.matrix(~ day + time, when)) > mm[is.na(mm)] <- 1 > mm > > You might also want to set the levels of your factors first so that it > includes levels that are not in the data and so that the levels are > sorted in an order other than alphabetical: >
The levels(when$day) lines should be: levels(when$day) <- c("Monday", "Tuesday", "Wednesday", "Thursday", "Friday", "Saturday", "Sunday") -- Statistics & Software Consulting GKX Group, GKX Associates Inc. tel: 1-877-GKX-GROUP email: ggrothendieck at gmail.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.