Complex sampling survey utilisation du package survey

[sylvie ahoussou]

Bonjour tout le monde

Je ne comprends pas quel argument je suis censée utiliser pour la fonction svydesign pour que R comprenne que mon plan de sondage est le suivant:
* sondage à 2 degrés:
- premier degré: ferme, le numéro d'identification des fermes est "num" , et les fermes sont tirées au sort dans la population fpc1
- deuxième degré : animaux, leur numéro d'identification est Id_An et ils sont tirés au sort dans la population Totanim

ma base de données tab1 contient 11792 entrées, la clé primaire (excusez moi pour ce terme "Access") est Id_An

J'ai essayé :
clustot<-svydesign(id=~num+Id_An, fpc=~fpc1+~Totanim, data=tab1)

et j'obtiens le message d'erreur suivant:

> clustot<-svydesign(id=~num+Id_An, fpc=~fpc1+~Totanim, data=tab1)
Erreur dans as.fpc(fpc, strata, ids) :
FPC implies >100% sampling in some strata.

Mais je ne suis pas dans le cas d'un sondage en grappe mais d'un sondage à 2 degrés. Savez vous quels arguments je dois utiliser ?

Merci d'avance

[Renaud Lancelot]

ça ou fé ?

Bonjour tout le monde

Je ne comprends pas quel argument je suis censée utiliser pour la fonction svydesign pour que R comprenne que mon plan de sondage est le suivant:
* sondage à 2 degrés:
- premier degré: ferme, le numéro d'identification des fermes est "num" , et les fermes sont tirées au sort dans la population fpc1
- deuxième degré : animaux, leur numéro d'identification est Id_An et ils sont tirés au sort dans la population Totanim

ma base de données tab1 contient 11792 entrées, la clé primaire (excusez moi pour ce terme "Access") est Id_An

J'ai essayé :

Code : Tout sélectionner

clustot<-svydesign(id=~num+Id_An, fpc=~fpc1+~Totanim, data=tab1)


et j'obtiens le message d'erreur suivant:

> clustot<-svydesign(id=~num+Id_An, fpc=~fpc1+~Totanim, data=tab1)
Erreur dans as.fpc(fpc, strata, ids) :
FPC implies >100% sampling in some strata.

Mais je ne suis pas dans le cas d'un sondage en grappe mais d'un sondage à 2 degrés. Savez vous quels arguments je dois utiliser ?

Merci d'avance

Apparemment, il y a une coquille dans ton code:

Code : Tout sélectionner

fpc=~fpc1+~Totanim

qui devrait peut-être être

Code : Tout sélectionner

fpc = ~ fpc1 + Totanim


C'est tjs une bonne idée d'aérer le code...

Aussi, regarde l'exemple dans l'aide de svydesign:

Code : Tout sélectionner

# two-stage cluster sample
dclus2 <- svydesign(id = ~ dnum + snum, fpc = ~ fpc1 + fpc2, data = apiclus2)

l'argument fpc décrivant la fraction de sondage dans les différentes unités, je ne vois pas pourquoi tu indiquerais la population totale (Totanim) dans cet argument. Si tu veux plus d'aide, il faudrait que tu nous fournisses un petit jeu de données exemple, avec le code générant l'erreur.

Tjen bé raid !

Renaud