[R-es] Contrastes con el paquete survey (svycontrast)

[Emilio Torres Manzanera]

Estimados usuarios:

   Estoy intentando reproducir el ejemplo 6.4 de Thomas Lumley. Complex  
Survey. Editorial Wiley. 2010 (ver la página en google:  
http://books.google.es/books?id=9wN0HbBEmb0C&pg=PA114&dq=thomas+lumley+scottish+household+survey&hl=es&ei=HwwSToaXHsSO8gPwsLDGDg&sa=X&oi=book_result&ct=result&resnum=1&ved=0CCoQ6AEwAA#v=onepage&q&f=false),  
pero no consigo acertar cómo se configura el contraste. ¿A alguno de  
ustedes le sucede lo mismo?
Gracias por adelantado
Emilio


## Figure 6.4 pág. 114 Thomas Lumley. Complex Survey.
## Aquí se puede ver la página
##  
http://books.google.es/books?id=9wN0HbBEmb0C&pg=PA114&dq=thomas+lumley+scottish+household+survey&hl=es&ei=HwwSToaXHsSO8gPwsLDGDg&sa=X&oi=book_result&ct=result&resnum=1&ved=0CCoQ6AEwAA#v=onepage&q&f=false
library(survey)
load(url("http://faculty.washington.edu/tlumley/svybook/shs.rda")) # Baja  
un fichero de 1.5M
m <- svyglm( intuse ~ I(age-18)*sex , design = shs, family =  
quasibinomial())
m2 <- svyglm( intuse ~ (pmin(age, 35) + pmax(age, 35))*sex , design = shs,  
family = quasibinomial())
summary(m) # OK. Coincide con el libro
summary(m2) # OK. Coincide con el libro
svycontrast(m2, quote('pmin(age, 35)' +'pmin(age, 35):sexfemale' )) #  
¡Falla!
## Error en deriv.default(expr, names(datalist)) :
##   expresión inválida en 'FindSubexprs'
##
sessionInfo() # Mi configuración
## R version 2.12.0 (2010-10-15)
## Platform: i486-pc-linux-gnu (32-bit)

## locale:
##  [1] LC_CTYPE=es_ES.UTF-8       LC_NUMERIC=C
##  [3] LC_TIME=es_ES.UTF-8        LC_COLLATE=es_ES.UTF-8
##  [5] LC_MONETARY=C              LC_MESSAGES=es_ES.UTF-8
##  [7] LC_PAPER=es_ES.UTF-8       LC_NAME=C
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C
## [11] LC_MEASUREMENT=es_ES.UTF-8 LC_IDENTIFICATION=C

## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base

## other attached packages:
## [1] survey_3.24

## loaded via a namespace (and not attached):
## [1] tools_2.12.0
-- 
-------------------------------------------------
Emilio Torres Manzanera
Fac. de Comercio - Universidad de Oviedo
c/ Luis Moya 261, E-33203 Gijón (Spain)
Tel. 985 182 197 email: torres en uniovi.es