Re: [Rd] Using response variable in interaction as explanatory variable in glm crashes R

2017-10-09 Thread Martin Maechler
> Jan van der Laan > on Fri, 6 Oct 2017 12:13:39 +0200 writes: > It is actually model.matrix that crashes, not glm. Same > crash occurs with e.g. lm. > model.matrix(dob_mon ~ dob_day*dob_mon, data = tab) > also crashes R. Yes, segmentation fault. It only happens wh

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

2017-10-09 Thread peter dalgaard
AFAIR, it is a little more subtle than that. If you have replication weights, then the estimates are right, it is "just" that the SE from summary.lm() are wrong. Somehow, the text should reflect this. It is of some importance when you put glm() into the mix, because you can in fact get correct