Hi Tobias,
Thanks for your answer but it doesn't change the error message
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
Maybe the probleme is that in some farms all the animals were examined (100%
sampling) which is the case in my study sometimes
--------------------------------------------------
From: "Tobias Verbeke" <[EMAIL PROTECTED]>
Sent: Friday, September 12, 2008 10:25 AM
To: "Ahoussou Sylvie" <[EMAIL PROTECTED]>
Subject: Re: [R] Fw: Complex sampling survey _ Use of survey package
Hi Sylvie,
You still use too much tilde signs (~) as Thomas Lumley already pointed
out.
Replace your line
clustot <- svydesign(id=~num+ ~ Id_An, fpc=~fpc1+~Totanim, data=tab1)
with
clustot <- svydesign(id=~num+Id_An, fpc=~fpc1+Totanim, data=tab1)
Hope this helps,
Tobias
P.S.
"~ num + Id_An" is called a formula in R and is omnipresent as
an interface in R functions
Ahoussou Sylvie wrote:
--------------------------------------------------
From: "Ahoussou Sylvie" <[EMAIL PROTECTED]>
Sent: Friday, September 12, 2008 9:48 AM
To: "Thomas Lumley" <[EMAIL PROTECTED]>
Subject: Re: [R] Complex sampling survey _ Use of survey package
Thanks for your answer
I think I made a mistake when I recopied the 5 first rows of my database
here is the table with the comlums of interest
num esp fpc1 Totanim Id_An
2045 G 551 12 10
2046 C 551 68 11
2070 G 551 9 50
2070 S 551 9 51
2070 S 551 9 52
yes Totanim is the total number of animals in the farm and num is the
total number of herds
I keep on obtaining this error message
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.
--------------------------------------------------
From: "Thomas Lumley" <[EMAIL PROTECTED]>
Sent: Friday, September 12, 2008 12:02 AM
To: "Ahoussou Sylvie" <[EMAIL PROTECTED]>
Cc: <r-help@r-project.org>
Subject: Re: [R] Complex sampling survey _ Use of survey package
On Thu, 11 Sep 2008, Ahoussou Sylvie wrote:
Data base = tab1 here are the five first rows of the database (nrow =
11792)
num
esp
Quarters
Totcat
Totshp
Totgt
Tbtpos
fpc1
Totanim
Id_An
10
2045
G
01-Q1
0
0
12
1
551
10
10
11
2046
G
01-Q1
8
0
60
1
551
11
11
50
2070
G
01-Q1
0
3
6
1
551
50
50
51
2070
S
01-Q1
0
3
6
1
551
51
51
52
2070
S
01-Q1
0
3
6
1
551
52
52
This has 11 data values in each row but only 10 variable names. If we
assume the first number doesn't belong it, the fpc1 column lines up
correctly, so I will do that.
It is surprising that Id_An is the same as Totanim for each of these
five records
- First strata = herds
herd's id column = num
sampled in 551 herds (column fpc1)
- Second strata = animals
animal id column= Id_An
sampled in Totanim
I tried :
clustot <- svydesign (id=~num+Id_An, fpc=~fpc1+~Totanim, data=tab1)
This should be
clustot <- svydesign (id=~num+Id_An, fpc=~fpc1+Totanim, data=tab1)
but that doesn't explain the problem
Your svydesign() call looks correct if fpc1 is the total number of
herds in the population and Totanim is the total number of animals in
the particular herd. I am not sure from your explanation whether this
is how the variables are defined.
-thomas
and I get this error message
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.
I don't have any strata with 100% sampling so do you know which
arguments of the fonction svydesign I'm supposed to use ?
Thanks in advance
AHOUSSOU Sylvie
Vétérinaire Epidémiologiste
CIRAD Domaine Duclos
97 170 Petit-Bourg
tel : 05 90 25 59 47
[[alternative HTML version deleted]]
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
