Re: [R] problem with svyglm

2012-11-28 Thread Pablo Menese
I achive something diferent, I replicated the t value, the std. error and the hypothesis test but differents betas. But, you are right, the thing is, I detach the dataset, but even with it, I couldn't. I going to describe all because perhaps I omitted something important. I have this vector for th

Re: [R] problem with svyglm

2012-11-27 Thread David Winsemius
On Nov 27, 2012, at 2:31 PM, Pablo Menese wrote: Sorry, it send it alone... When I use it: logit <- glm(bach ~ egp4 + programa, weight=wst7, family=quasibinomial(link"logit")) I reach the same betas that in STATA, but the hypothesis test, the t value, and the std. error is different.

Re: [R] problem with svyglm

2012-11-27 Thread Pablo Menese
Sorry, it send it alone... When I use it: logit <- glm(bach ~ egp4 + programa, weight=wst7, family=quasibinomial(link"logit")) I reach the same betas that in STATA, but the hypothesis test, the t value, and the std. error is different. I think that the solution can't be so far from this...

Re: [R] problem with svyglm

2012-11-27 Thread Pablo Menese
I colud not, even without attach the dataset. The thing is, when I use this On Fri, Nov 23, 2012 at 5:56 PM, David Winsemius wrote: > > On Nov 23, 2012, at 12:08 PM, Pablo Menese wrote: > > I have this problem. >> >> test <- svydesign(id=~1,weights=~peso) >> >> logit <- svyglm(bach ~ job2 +

Re: [R] problem with svyglm

2012-11-23 Thread David Winsemius
On Nov 23, 2012, at 12:08 PM, Pablo Menese wrote: I have this problem. test <- svydesign(id=~1,weights=~peso) logit <- svyglm(bach ~ job2 + mujer + egp4 + programa + delay + mdeo + str + evprivate, family=binomial,design=test) then appear: Error in svyglm.survey.design(bach ~ job2 + muje

Re: [R] problem with svyglm()

2009-04-06 Thread Thomas Lumley
On Thu, 2 Apr 2009, Federico Calboli wrote: Hello, I'm trying to use the function svyglm in the library survey. I create a data survey object: data_svy<- svydesign(id=~PSU, strata=~sample_domain, weights=~sample_weight, data=data, nest=TRUE) and I try to use svyglm() with little success: R<