R help es - Jul 2011 - Contrastes con el paquete survey (svycontrast)

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

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]]