Hola Ivan gracias por tu ayuda, justo ayer a ultima hora hice una
modificacion en el programa q hizo q el programita ande, aunq no entendi
porque lo hizo.
La modificacion fue agregar "*w" al final de cada parametro en la
definicion
de la funcion "grini" q habia creado, cosa q deduje de mirar la
funcion q le
habia robado al ejemplo del manual de boot
ratio <- function(d, w) sum(d$EDAD * w)/sum(d$UNO * w)
En tonces mi funcion paso de ser:
function(d,w){gini(d$MT_YOP_NEW, w=d$FEXP)
a ser:
grini<-function(d,w){gini(d$MT_YOP_NEW**w*, w=d$FEXP**w*)}
Como les decia este cambio hizo q el programa ande, (deje de darme variancia
cero) pero no se cual es su significado. Ahora el tema es q los resultados
son algo distintos del gini original, pero bueno sobre eso seguire mirando.
Creo q lo importante para mi ahora es entender la explicacion, q seguro esta
relacionada con el vector de indices q me dices q boot necesita para
ejecutarse, sobre lo cual me gustaria tener mas informacion.
Saludos
El 29 de septiembre de 2010 08:00,
<r-help-es-request@r-project.org>escribió:
> Envíe los mensajes para la lista R-help-es a
> r-help-es@r-project.org
>
> Para subscribirse o anular su subscripción a través de la WEB
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> O por correo electrónico, enviando un mensaje con el texto "help"
en
> el asunto (subject) o en el cuerpo a:
> r-help-es-request@r-project.org
>
> Puede contactar con el responsable de la lista escribiendo a:
> r-help-es-owner@r-project.org
>
> Si responde a algún contenido de este mensaje, por favor, edite la
> linea del asunto (subject) para que el texto sea mas especifico que:
> "Re: Contents of R-help-es digest...". Además, por favor, incluya
en
> la respuesta sólo aquellas partes del mensaje a las que está
> respondiendo.
>
>
> Asuntos del día:
>
> 1. Re: calcular la variancia de gini por bootstrap (Jorge Ivan Velez)
> 2. Re: Problema con read.spss (guivivi@alumni.uv.es)
> 3. The R Inferno (miguel.angel.rodriguez.muinos@sergas.es)
> 4. Re: The R Inferno (Usuario R)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 28 Sep 2010 20:45:37 -0400
> From: Jorge Ivan Velez <jorgeivanvelez@gmail.com>
> To: jose ramirez costa <joseramirezcosta@gmail.com>
> Cc: R-help-es <r-help-es@r-project.org>
> Subject: Re: [R-es] calcular la variancia de gini por bootstrap
> Message-ID:
>
<AANLkTiks_WC3g4FnLSr6VK=-Ue4Ejsb7FEmf8MvfZY+Q@mail.gmail.com<Ue4Ejsb7FEmf8MvfZY%2BQ@mail.gmail.com>
> >
> Content-Type: text/plain
>
> Buenas noches Jose,
>
> Considera el siguiente ejemplo (modificado de [1]):
>
> # paquete reldist que contiene la funcion gini()
>
> require(reldist)
>
>
> # datos y pesos -- obtenidos como dput(d)
>
> d <- structure(list(x = c(541, 1463, 2445, 3438, 4437, 5401, 6392,
>
> 8304, 11904, 22261), w = c(0.287577520124614, 0.788305135443807,
>
> 0.4089769218117, 0.883017404004931, 0.940467284293845, 0.0455564993899316,
>
> 0.528105488047004, 0.892419044394046, 0.551435014465824, 0.456614735303447
>
> )), .Names = c("x", "w"), row.names = c(NA, -10L),
class = "data.frame")
>
> d
>
>
> # coeficiente gini ponderado
>
> with(d, gini(x, w)) # 0.433399
>
>
> # bootstrap
>
> require(boot)
>
>
> # funcion para calcular el coeficiente gini
>
> # a partir de un data.frame d
>
> # cuyas columnas son x y w
>
> # ---> aqui index es un vector de indices que
>
> # boot necesita para ejecutarse
>
> giniboot <- function(d, index){
>
> d <- d[sample(index),] # no tienes esto
>
> with(d, gini(x, w))
>
> }
>
> res <- boot(d, giniboot, R = 1000)
>
> res
>
>
> # graficos
>
> plot(res)
>
>
> # algunas medidas para el estadistico
>
> summary(res$t)
>
> var(res$t)
>
> quantile(res$t)
>
> Quizas no es propiamente lo que estas buscando, pero te da una mejor idea
> de
> lo que esta fallando en tu programa.
>
> Ahora, si definitivamente no es lo que buscas, ayudanos a ayudarte
> enviandonos una parte de los datos via dput() (ver ejemplo en este
> mensaje),
> asi como la funcion gini que estas usando. Sospecho que el gini que estas
> usando es el mismo de [1], pero es dificil saberlo con la poca informacion
> que nos das.
>
> Un saludo,
> Jorge Ivan Velez
>
> [1] http://rss.acs.unt.edu/Rdoc/library/reldist/html/gini.html
>
>
> 2010/9/28 jose ramirez costa <>
>
> > Hola, paso el mini programita q estoy viendo, lo q me llama la
atencion
> es
> > una parte donde se definen las funciones.
> >
> > Probe primero meter adentro del boots la estadistica a estimar usando
> > directamente gini(varible, pesos) pero no me dejo.
> >
> > Vi q en el ej del manual de boots, siempre define antes la funcion,
> > entonces
> > probe definir antes una funcion haciendo
> > grini<-function(x) {gini(variable, peso)} pero tampoco me dejo.
> >
> > Entonces preste atencion q la funcion q definia en el ejmeplo tenia 2
> > parametros d y w, este es el ejemplo del manual modificado para mi
base:
> > ratio <- function(d, w) sum(d$EDAD * w)/sum(d$UNO * w)
> >
> > Como ves tiene 2 parametros: d y w, entonces le puse eso mismo a mi
> funcion
> > grini: grini<-function(d,w){gini(d$MT_YOP_NEW, w=d$FEXP)}, donde el
w q
> > esta
> > adentro no es el w del parametro de la funcion sino el el parametro
peso
> de
> > la funcion gini.
> >
> > Me hizo el bootstrap pero me da siempre la misma estimacion para todas
> las
> > replicaciones de las muestras, por lo cual la variancia del estimador
por
> > bootstrap da 0.
> >
> > Entonces me di cuenta q algo de ese parametro tenia q ver con eso,
pero
> no
> > se q es.
> >
> > Por lo tanto la pregunta es: q hace ese parametro w en la funcion???
> >
> > Creo q la solucion debe ser algo simple pero falta entender eso, para
> poder
> > solucionarlo.
> >
> > Bueno si tenes tiempo para mirarlo te agradeceria mucho.
> >
> > Saludos
> > *
> > *
> > *El programita:*
> >
> > # usual bootstrap of the ratio of means using the city data
> > parag$UNO<-1
> > grini(parag$MT_YOP_NEW)
> >
> > ratio <- function(d, w) sum(d$EDAD * w)/sum(d$UNO * w)
> > r<-boot(parag, ratio, R=999, sim="ordinary",
stype="w",
> > strata=parag$COMUNA,
> > weights=parag$FEXP)
> > r
> > str(r)
> > View(r$t)
> >
> > grini<-function(d,w){gini(d$MT_YOP_NEW, w=d$FEXP)}
> > g<-boot(parag, grini, R=999, sim="ordinary",
stype="w",
> > strata=parag$COMUNA,
> > weights=parag$FEXP)
> > g
> > str(g)
> > View(g$t)
> >
> > --
> > "El futuro es ilusorio es una trampa q se inventa cualquier
sistema para
> q
> > la gente se acobarde y agache la cabeza, trabaje, produzca y se haga
> > esclava
> > por el PUTO futuro. QUE FUTURO TE ASEGURASTE VOS? En el laburo te
> aseguran
> > q se acabaron los accidentes,la sirrosis, el cancer, el tiro q nos
puede
> > pegar el tipo q entre a robar aca dentro de un minuto?"
> > Fernando Robles DNI 4.131.417
> > Adhiere: José
> >
> >
> >
> >
> > --
> > "El futuro es ilusorio es una trampa q se inventa cualquier
sistema para
> q
> > la gente se acobarde y agache la cabeza, trabaje, produzca y se haga
> > esclava
> > por el PUTO futuro. QUE FUTURO TE ASEGURASTE VOS? En el laburo te
> aseguran
> > q se acabaron los accidentes,la sirrosis, el cancer, el tiro q nos
puede
> > pegar el tipo q entre a robar aca dentro de un minuto?"
> > Fernando Robles DNI 4.131.417
> > Adhiere: José
> >
> > [[alternative HTML version deleted]]
> >
> >
> > _______________________________________________
> > R-help-es mailing list
> > R-help-es@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
> >
>
> [[alternative HTML version deleted]]
>
>
>
> ------------------------------
>
> Message: 2
> Date: Wed, 29 Sep 2010 09:08:07 +0200 (CEST)
> From: <guivivi@alumni.uv.es>
> To: Jorge Ivan Velez <jorgeivanvelez@gmail.com>, R-help-es
> <r-help-es@r-project.org>
> Subject: Re: [R-es] Problema con read.spss
> Message-ID: <5936575528guivivi@alumni.uv.es>
> Content-Type: text/plain; charset="ISO-8859-1"
>
> Buenos días, Jorge Iván, muchas gracias por tu ayuda.
>
> Efectivamente de las tres opciones que me sugerías para solucionar el
> problema, la tercera ha sido la más inmediata.
> En mi base de datos había variables en las había un signo ? donde
> debiera haber una letra acentuada , así que R no las podía leer.
> He quitado estos ? y ya puedo leer la base con read.spss.
>
> Muchas gracias de nuevo por tu ayuda.
>
> Un saludo,
>
> Guillermo
>
> > Buenas tardes Guillermo,
> >
> > El primer mensaje que obtienes simplemente dice que el archivo .sav no
> esta
> > en el directorio del cual lo intentas leer. En el segundo caso, aunque
> > puedes localizar el archivo y comenzar a leerlo, es reportado un error
> > durante ese proceso de lectura. Este error, sin embargo, nada tiene
> que ver
> > con el hecho de que no tengas SPSS disponible en tu equipo. Lo que
esta
> > ocurriendo es que uno de los nombres de las variables en el archivo
> > "guille.sav" comienza con uno de los caracteres no
permitidos (o
> reservados
> > --- &, #, entre otros) en R.
> >
> > En cuanto a que hacer, existen varias posibilidades: 1. corregir la
> > funcion do_read_spss que se que se encuentra en la linea 22 de
read.spss
> > (para ver read.spss, simplemente escribe "read.spss" sin las
comillas
> en la
> > consola del R y luego oprime enter); 2. abrir la base de datos en un
> > programa equivalente, por ejemplo STATISTICA (bajo Windows -- es
> inmediato)
> > y 3., la mas sencilla, a mi modo de ver, solicitarle a quienes
> generaron el
> > archivo original que por favor revisen y cambien los nombres de
aquellas
> > variables cuyos nombres comienzan con uno de los caracteres
reservados.
> >
> > Puedo ayudarte con 2 (off-list) si te parece bien.
> >
> > Un saludo,
> > Jorge Ivan Velez
> >
> >
> >
> > 2010/9/28 <>
> >
> > > Hola,
> > >
> > > Tengo problemas a la hora de cargar un archivo que originalmente
está
> en
> > > SPSS.
> > >
> > > Trabajo con Fedora y uso R mediante el plugin Rgedit en gedit.
> > > Además en Fedora no tengo SPSS, sino PSPP
> > >
> > > Para cargar el archivo he probado estas dos cosas (escribo
también el
> > > mensaje de error que me genera R):
> > >
> > > >library_(foreign)
> > > >read.spss("guille.sav",to.data.frame = T)
> > > Error en read.spss("guille.sav", to.data.frame = T) :
> > > unable to open file: ''No existe el fichero o el
directorio''
> > >
> > > >
read.spss("/home/guillermo/Escritorio/guille.sav",to.data.frame = T)
> > > Error en
read.spss("/home/guillermo/Escritorio/guille.sav", :
> > > error reading system-file header
> > > Además: Mensajes de aviso perdidos
> > > In read.spss("/home/guillermo/Escritorio/guille.sav",
:
> > > /home/guillermo/Escritorio/guille.sav: position 18: Variable
name
> > > begins with invalid character
> > >
> > > He buscado información sobre este problema pero no he encontrado
> > > solución y estoy atascado, ¿alguien me podría ayudar?, parece que
este
> > > problema es bastante común y no sé si depende de la versión del
SPSS o
> > > del PSPP.
> > >
> > > Muchas gracias de antemano por la ayuda y espero que si alguien
tiene
> el
> > > mismo problema, este correo le pueda ser útil para que sigamos
todos
> > > aprendiendo R.
> > >
> > > Un saludo,
> > >
> > > Guillermo
> > >
> > > _______________________________________________
> > > R-help-es mailing list
> > > R-help-es@r-project.org
> > > https://stat.ethz.ch/mailman/listinfo/r-help-es
> > >
> >
>
>
>
> ------------------------------
>
> Message: 3
> Date: Wed, 29 Sep 2010 09:35:39 +0200
> From: <miguel.angel.rodriguez.muinos@sergas.es>
> To: <r-help-es@r-project.org>
> Subject: [R-es] The R Inferno
> Message-ID:
>
<0A5B869BDC0A7445B15F9A4A5379437719280D0568@ASSCCEVS02.sergas.local
> >
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hola a tod@s.
>
> Creo que este paper es interesante. Un poco "sobreactuado" en su
redacción,
> pero puede resultar de utilidad.
> Como dice el propio Patrick Burns: "If you are using R and you think
you''re
> in hell, this is a map for you".
>
> ... que ya sé que no es nuestro caso!
> }:-)
>
> http://www.burns-stat.com/pages/Tutor/R_inferno.pdf
>
>
> No he encontrado ninguna traducción por ahí, así que os lo paso en inglés.
>
>
> Un Saludo,
> ===============================================> Miguel Ángel Rodríguez
Muíños
> Asesoramento en Informática
> Servicio de Epidemioloxía
> Dirección Xeral de Saúde Pública e Planificación
> Consellería de Sanidade
> Xunta de Galicia
> ===============================================>
>
>
>
>
>
>
>
>
> Nota: A información contida nesta mensaxe e os seus posibles documentos
> adxuntos é privada e confidencial e está dirixida únicamente ó seu
> destinatario/a. Se vostede non é o/a destinatario/a orixinal desta mensaxe,
> por favor elimínea. A distribución ou copia desta mensaxe non está
> autorizada.
>
> Nota: La información contenida en este mensaje y sus posibles documentos
> adjuntos es privada y confidencial y está dirigida únicamente a su
> destinatario/a. Si usted no es el/la destinatario/a original de este
> mensaje, por favor elimínelo. La distribución o copia de este mensaje no
> está autorizada.
>
> See more languages: http://www.sergas.es/aviso_confidencialidad.htm
>
> XACOBEO 2010: http://www.xacobeo.org
>
>
>
> ------------------------------
>
> Message: 4
> Date: Wed, 29 Sep 2010 09:49:05 +0200
> From: Usuario R <r.user.spain@gmail.com>
> To: miguel.angel.rodriguez.muinos@sergas.es
> Cc: r-help-es@r-project.org
> Subject: Re: [R-es] The R Inferno
> Message-ID:
>
<AANLkTikekgOACTgCyEpJ0FqyHYf9CmdnCQt6a+QXsv0T@mail.gmail.com<AANLkTikekgOACTgCyEpJ0FqyHYf9CmdnCQt6a%2BQXsv0T@mail.gmail.com>
> >
> Content-Type: text/plain
>
> Hola Miguel Angel,
>
> Esto lo leí hace tiempo, está muy bien, muchas gracias!
> saludos
>
> 2010/9/29 <miguel.angel.rodriguez.muinos@sergas.es>
>
> > Hola a tod@s.
> >
> > Creo que este paper es interesante. Un poco "sobreactuado"
en su
> redacción,
> > pero puede resultar de utilidad.
> > Como dice el propio Patrick Burns: "If you are using R and you
think
> you''re
> > in hell, this is a map for you".
> >
> > ... que ya sé que no es nuestro caso!
> > }:-)
> >
> > http://www.burns-stat.com/pages/Tutor/R_inferno.pdf
> >
> >
> > No he encontrado ninguna traducción por ahí, así que os lo paso en
> inglés.
> >
> >
> > Un Saludo,
> > ===============================================> > Miguel Ángel
Rodríguez Muíños
> > Asesoramento en Informática
> > Servicio de Epidemioloxía
> > Dirección Xeral de Saúde Pública e Planificación
> > Consellería de Sanidade
> > Xunta de Galicia
> > ===============================================> >
> >
> >
> >
> >
> >
> >
> >
> >
> > Nota: A información contida nesta mensaxe e os seus posibles
documentos
> > adxuntos é privada e confidencial e está dirixida únicamente ó seu
> > destinatario/a. Se vostede non é o/a destinatario/a orixinal desta
> mensaxe,
> > por favor elimínea. A distribución ou copia desta mensaxe non está
> > autorizada.
> >
> > Nota: La información contenida en este mensaje y sus posibles
documentos
> > adjuntos es privada y confidencial y está dirigida únicamente a su
> > destinatario/a. Si usted no es el/la destinatario/a original de este
> > mensaje, por favor elimínelo. La distribución o copia de este mensaje
no
> > está autorizada.
> >
> > See more languages: http://www.sergas.es/aviso_confidencialidad.htm
> >
> > XACOBEO 2010: http://www.xacobeo.org
> >
> > _______________________________________________
> > R-help-es mailing list
> > R-help-es@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
>
>
>
> --
> Patricia García González
> r-es.org/
>
> [[alternative HTML version deleted]]
>
>
>
> ------------------------------
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
> Fin de Resumen de R-help-es, Vol 19, Envío 28
> *********************************************
>
--
"El futuro es ilusorio es una trampa q se inventa cualquier sistema para q
la gente se acobarde y agache la cabeza, trabaje, produzca y se haga esclava
por el PUTO futuro. QUE FUTURO TE ASEGURASTE VOS? En el laburo te aseguran
q se acabaron los accidentes,la sirrosis, el cancer, el tiro q nos puede
pegar el tipo q entre a robar aca dentro de un minuto?"
Fernando Robles DNI 4.131.417
Adhiere: José
[[alternative HTML version deleted]]
