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

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

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

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

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: Compute standard deviations from complex survey designs

On Tue, Oct 16, 2012 at 10:40 PM, Sebastian Weirich 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 manner? Usually you can do this s

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

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 __ R-help@r-project.o

Re: [R] package survey

l 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 if it is possible to find the >> codes for >> functions in

Re: [R] package survey

> -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 wro

Re: [R] package survey

a half-page of code. 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 help me t

Re: [R] package survey

sier or quicker way to do it from within R 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] package survey

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, dat

Re: [R] Package survey

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.

[R] Package survey

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. summary(base

Re: [R] Package Survey

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 repli

Re: [R] Package Survey

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, un

[R] Package Survey

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 tota