Your calibration model has one linear predictor, so it only has two
parameters, but you have specified three population totals. A slightly
different setup should work; try again with:
st <- factor(1:3, levels=c(3, 1:2))
pop<-c('(Intercept)'=100, st1=10, st2=20)
James Reilly
Department of Statis
Hi R-users!
Can anyone tell me what this error message refers to exactly
> a<-c('m','f','m')
> b<-c(1:3)
> st<-c(1:3)
> data.frame(a,b,st)
a b st
1 m 1 1
2 f 2 2
3 m 3 3
> data.frame(a,b,st)->l
> p<-svydesign(ids=~1, data=l)
> pop<-c('(intercept)'=100, st1=10,st2=20)
> calibrate(p,~st, pop)
2 matches
Mail list logo