[Rd] Bug in model.matrix.default for higher-order interaction encoding when specific model terms are missing

2017-10-12 Thread Tyler
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

Re: [Rd] Discourage the weights= option of lm with summarized data

2017-10-12 Thread Arie ten Cate
OK. We have now three suggestions to repair the text: - remove the text - add "not" at the beginning of the text - add at the end of the text a warning; something like: "Note that in this case the standard estimates of the parameters are in general not correct, and hence also the t values an