Hola,

Los ejemplos en el libro están basados en la versión 3.10-1 del paquete,
sobre una versión de R 2.7.2:

*The examples in this book used version 3.10-1 of the survey package and
were*
*run in R version 2.7.2. The home page for the survey package ( h t t p : /
/ f a c u l t y .*
*washingt on. edu/t lumley/survey) will have information about any changes
for*
*newer versions of the package as they are released.*

Y la versión que utilizas, del detalle que incluyes, es la más actual 3.24.
Me inclino a pensar que el problema que tienes viene motivado por esta
diferencia en las versiones de los paquetes.

Mi sugerencia es que busques en el histórico del paquete, para bien
instalarte la versión con la que se ilustró el libro o bien que busques en
la página del libro los cambios que habría que introducir en la función
"svycontrast" para que que funcione el ejemplo. La alternativa definitiva es
que le escalaras el problema al propio autor del libro o a los mantenedores
del paquete para que te dirigieran más rápidamente.

Saludos,
Carlos Ortega
www.qualityexcellence.es



On Mon, Jul 4, 2011 at 9:30 PM, Emilio Torres Manzanera <torres@uniovi.es>wrote:

>
> 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<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<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<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@uniovi.es
>
> ______________________________**_________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/**listinfo/r-help-es<https://stat.ethz.ch/mailman/listinfo/r-help-es>
>

	[[alternative HTML version deleted]]

