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<-svyglm(data_svy[,4]~(data_svy[,iCol]==listModality[[iVar]] 
[iMod]),design=data_svy, family=binomial(link="logit")

Error in svyglm.survey.design(data_svy[, 4] ~ (data_svy[, iCol] ==  
listModality[[iVar]][iMod]),  :
         all variables must be in design= argument

Puzzlingly the following code works:

R<-glm(data[,4]~(data[,iCol]==listModality[[iVar]][iMod]),  
family=binomial(link="logit")

Any help would be gratefully accepted.

Best,

Federico

--
Federico C. F. Calboli
Department of Epidemiology and Public Health
Imperial College, St. Mary's Campus
Norfolk Place, London W2 1PG

Tel +44 (0)20 75941602   Fax +44 (0)20 75943193

f.calboli [.a.t] imperial.ac.uk
f.calboli [.a.t] gmail.com




	[[alternative HTML version deleted]]


