less than full rank contrast methods
On Tue, Dec 7, 2010 at 7:54 AM, Gabor Grothendieck
<ggrothendieck at gmail.com> wrote:
'On Tue, Dec 7, 2010 at 5:19 AM, mxkuhn <mxkuhn at 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