Re: [R] Package survey: singularities in linear regression models

2013-05-05 Thread Thomas Lumley
Ok, that's more helpful. The problem is with replicate-weight designs, and it's because svyglm() uses the fitted coefficients from the point estimate as starting values for fitting the replicates. And even if that is changed, the computation of the replicate variance doesn't like all the replicat

Re: [R] Package survey: singularities in linear regression models

2013-05-03 Thread Sebastian Weirich
Well, I have uploaded the data in the public folder of my dropbox. Due to data confidentiality, I haved to change the labels. To load the data: con <- url( "http://dl.dropboxusercontent.com/u/101865137/datEx.rda"; ) print(load(con)) # The replicate weights were created according to the jackknife

Re: [R] Package survey: singularities in linear regression models

2013-05-02 Thread Thomas Lumley
On Fri, May 3, 2013 at 2:27 AM, Sebastian Weirich < sebastian.weir...@iqb.hu-berlin.de> wrote: > Hello, > > I want to specify a linear regression model in which the metric outcome is > predicted by two factors and their interaction. glm() computes effects for > each factor level and the levels of

[R] Package survey: singularities in linear regression models

2013-05-02 Thread Sebastian Weirich
Hello, I want to specify a linear regression model in which the metric outcome is predicted by two factors and their interaction. glm() computes effects for each factor level and the levels of the interaction. In the case of singularities glm() displays "NA" for the corresponding coefficients