dans les sous-populations avec des données pondérées à l'aide de trousse de sondage()

Question

ftp://cran.r-project.org/pub/R/web/packages/survey/vignettes/domain.pdf

L'ensemble de données complet est tch2012. Cependant, je ne m'intéresse qu'à la sous-population de tch2012 dans laquelle deux critères sont remplis: âge < = 5 et genre == "féminin". Et dans cette sous-population, je veux comparer ceux avec la maladie (maladie == "1") et sans la maladie (maladie == "0").

Voici le code que j'ai écrit:

bibliothèque (enquête)

tch2012.tsl.dsgn < - svydesign (id = ~ HOSP_KID, couches = ~ KID_STRATUM, poids = ~ DISCWT, data = tch2012 , nid = TRUE)

créer un pointeur vers l'âge pédiatrique sous-population des femmes 5 ans et moins

tch2012_f_age5.tsl.dsgn < - sous-ensemble (tch2012.tsl.dsgn, AGE < = 5 & sexe == "femme")

données pondérées du nombre total avec et sans la maladie chez les patients pédiatriques âgés de 5 ans et des femmes sous

svyby (~ nombre, ~ la maladie, la conception = tch2012_f_age5.tsl.dsgn , svytotal)

Cependant, je reçu le message d'erreur ci-dessous lorsque je courais le svyby()

erreur en somme (sapply (covmats, Ncol)): invalide 'type' (liste) des arguments

Depuis que je ne suis pas très familier avec le pesage données ted, je ne sais pas comment résoudre les problèmes.

Merci d'avance pour l'aide!

Answer 1

ce code fonctionne

library(survey) 
data(api) 
dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc) 
x <- subset(dclus1 , sch.wide == 'Yes') 
svyby(~api00, ~stype, design=x, svytotal) 

s'il vous plaît modifier votre question en ajoutant un exemple minimal reproductible How to make a great R reproducible example?