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

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

[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.

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 the

[R] Package survey: Compute standard deviations from complex survey designs

2012-10-16 Thread Sebastian Weirich
Hello, svyvar from the survey package computes variances (with standard errors) from survey design objects. Is there any way to compute standard deviations and their standard errors in a similar manner? Thanks a lot, Sebastian __

Re: [R] Package survey: Compute standard deviations from complex survey designs

2012-10-16 Thread Thomas Lumley
On Tue, Oct 16, 2012 at 10:40 PM, Sebastian Weirich sebastian.weir...@iqb.hu-berlin.de wrote: Hello, svyvar from the survey package computes variances (with standard errors) from survey design objects. Is there any way to compute standard deviations and their standard errors in a similar

Re: [R] package survey

2010-12-20 Thread Martin Maechler
] package survey On Dec 18, 2010, at 8:11 PM, Joel Schwartz wrote: and does anyone know if it is possible to find the codes for functions in survey package? Yes, you can find the code by doing the following: 1) Go to the CRAN R package list

[R] package survey

2010-12-18 Thread andrija djurovic
Hi R users, could someone help me to find out which formulas, for standard error calculation, are used in following example: a=data.frame(weights=rep(c(10,1),c(4,1)),fpc=rep(41,5),uk=rep(1,5)) srs-svydesign(id=~1, weights=~weights, data=a) srs1-svydesign(id=~1, weights=~weights,fpc=~fpc,

Re: [R] package survey

2010-12-18 Thread Joel Schwartz
but ,if there is, I haven't learned it yet. Joel -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of andrija djurovic Sent: Saturday, December 18, 2010 4:23 PM To: r-help@r-project.org Subject: [R] package survey Hi R users, could someone

Re: [R] package survey

2010-12-18 Thread David Winsemius
- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of andrija djurovic Sent: Saturday, December 18, 2010 4:23 PM To: r-help@r-project.org Subject: [R] package survey Hi R users, could someone help me to find out which formulas, for standard error calculation

Re: [R] package survey

2010-12-18 Thread Joel Schwartz
-Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Saturday, December 18, 2010 5:54 PM To: Joel Schwartz Cc: r-help@r-project.org Subject: Re: [R] package survey On Dec 18, 2010, at 8:11 PM, Joel Schwartz wrote: and does anyone know

[R] Package survey

2008-10-22 Thread Rita Sousa
Hi, I’m using the svyby for total statistics, for example: svyby(~p_igov,~div_a,desenho_nps,svytotal,drop.empty.groups=TRUE,vartype =c(se,var,cvpct)) In the numerical variable p_igov (and others) I have many non responses but if I maintain the NA it doesn’t work.

Re: [R] Package survey

2008-10-22 Thread Thomas Lumley
I think this has already been fixed in version 3.9-1. Also, you probably want to add na.rm=TRUE to the call svyby(~p_igov,~div_a,desenho_nps,svytotal,drop.empty.groups=TRUE,vartype =c(se,var,cvpct), na.rm=TRUE) otherwise the estimated total will be NA whenever any observation is NA.

Re: [R] Package Survey

2007-09-24 Thread Thomas Lumley
On Sun, 23 Sep 2007, Rita Cristina Pinto Sousa wrote: I?m using the package survey to obtain the statistics, fundamentally the variance estimates. Can you explain why do I obtain the same result with the replicate weights (as.svrepdesign function), for a stratified sample, and without the

[R] Package Survey

2007-09-20 Thread Rita Sousa
Hello, How I use the function as.svrepdesign without memory.size problems? desenho_npc_JK - as.svrepdesign(desenho_npc,type=JKn) Error: cannot allocate vector of size 161.3 Mb In addition: Warning messages: 1: Reached total allocation of 1022Mb: see help(memory.size) 2: Reached total

Re: [R] Package Survey

2007-09-20 Thread Thomas Lumley
On Thu, 20 Sep 2007, Rita Sousa wrote: How I use the function as.svrepdesign without memory.size problems? desenho_npc_JK - as.svrepdesign(desenho_npc,type=JKn) Error: cannot allocate vector of size 161.3 Mb There is currently no easy way to affect the amount of memory that this uses,