search for: x1p

...ble by read.table ?!) Note, we have to agree if read.table will be called with (...., header = TRUE) or not. Here, mostly for teaching purposes, we have (almost) all our data files in the format nam.1 col.2 var.3 ... var.<p> <x11> <x12> <x13> ... <x1p> <x21> <x22> <x23> ... <x2p> .. .. .. .. .. .. .. .. .. .. .. .. I.e. we don't have `rownames' most of the time and don't want to type '1' '2' ... 'n' in front o...

...:X2C 1 1 0 0 0 0 0 2 1 0 0 1 0 0 3 1 1 0 0 0 0 4 1 1 0 0 1 0 5 1 0 1 0 0 0 6 1 0 1 0 0 1 (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C 1 1 0 0 0 0 0 2 1 1 0 0 0 0 3 1 0 1 0 0 0 4 1 1 0 1 0 0 5 1 0 0 0 1 0 6 1 1...

Hi, I recently ran into an inconsistency in the way model.matrix.default handles factor encoding for higher level interactions with categorical variables when the full hierarchy of effects is not present. Depending on which lower level interactions are specified, the factor encoding changes for a higher level interaction. Consider the following minimal reproducible example: -------------- >

...1 0 0 1 0 0 >> 3 1 1 0 0 0 0 >> 4 1 1 0 0 1 0 >> 5 1 0 1 0 0 0 >> 6 1 0 1 0 0 1 >> >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C >> 1 1 0 0 0 0 0 >> 2 1 1 0 0 0 0 >> 3 1 0 1 0 0 0 >> 4 1 1 0 1 0 0 >> 5 1 0 0...

...gt;> 3 1 1 0 0 0 0 >> >> 4 1 1 0 0 1 0 >> >> 5 1 0 1 0 0 0 >> >> 6 1 0 1 0 0 1 >> >> >> >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C >> >> 1 1 0 0 0 0 0 >> >> 2 1 1 0 0 0 0 >> >> 3 1 0 1 0 0 0 >> >> 4 1 1 0 1 0 0...

...0 >> >> >> 4 1 1 0 0 1 0 >> >> >> 5 1 0 1 0 0 0 >> >> >> 6 1 0 1 0 0 1 >> >> >> >> >> >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C >> >> >> 1 1 0 0 0 0 0 >> >> >> 2 1 1 0 0 0 0 >> >> >> 3 1 0 1 0 0 0 >> >> >> 4 1...

...0 0 > >> 3 1 1 0 0 0 0 > >> 4 1 1 0 0 1 0 > >> 5 1 0 1 0 0 0 > >> 6 1 0 1 0 0 1 > >> > >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C > >> 1 1 0 0 0 0 0 > >> 2 1 1 0 0 0 0 > >> 3 1 0 1 0 0 0 > >> 4 1 1 0 1 0 0 > >> 5...

...4 1 1 0 0 1 0 >> >> >> >> 5 1 0 1 0 0 0 >> >> >> >> 6 1 0 1 0 0 1 >> >> >> >> >> >> >> >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C >> >> >> >> 1 1 0 0 0 0 0 >> >> >> >> 2 1 1 0 0 0 0 >> >> >> >> 3 1 0 1 0 0 0 >> >...

...0 0 0 0 > >> >> 4 1 1 0 0 1 0 > >> >> 5 1 0 1 0 0 0 > >> >> 6 1 0 1 0 0 1 > >> >> > >> >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C > >> >> 1 1 0 0 0 0 0 > >> >> 2 1 1 0 0 0 0 > >> >> 3 1 0 1 0 0 0 > >> >> 4 1 1 0...

...t;> >> 4 1 1 0 0 1 0 > >> >> >> 5 1 0 1 0 0 0 > >> >> >> 6 1 0 1 0 0 1 > >> >> >> > >> >> >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C > >> >> >> 1 1 0 0 0 0 0 > >> >> >> 2 1 1 0 0 0 0 > >> >> >> 3 1 0 1 0 0 0 > >> >> &g...

...0 0 1 0 > >> >> >> >> 5 1 0 1 0 0 0 > >> >> >> >> 6 1 0 1 0 0 1 > >> >> >> >> > >> >> >> >> (Intercept) X1q X1p:X2B X1q:X2B X1p:X2C X1q:X2C > >> >> >> >> 1 1 0 0 0 0 0 > >> >> >> >> 2 1 1 0 0 0 0 > >> >> >> >> 3 1 0 1 0 0...