Re: [R] question regarding an error message from survey package

2007-11-24 Thread James Reilly
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

[R] question regarding an error message from survey package

2007-11-24 Thread eugen pircalabelu
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)