[R-es] Media Ponderada

[Freddy López]

Hola José.

A mí no me queda muy clara la pregunta.

La frase:

> Pero no soy capaz de hallar la media ponderada en funcion del arte

me confunde un poco.

A ver, la media, cuando la calculas ponderada, te retornará un valor
(estoy pensando en voz alta). Y esta se calcula con un vector que está
pesado por otro vector (numérico, naturalmente) de igual tamaño.

Lo que se me ocurre es que estás pidiendo la media estratificada a
través de cada valor del factor arte, que, me parece, es otra cosa.

Si este es el caso, pienso que una solución (ojalá alguien tenga una
más elegante) puede ser crear un conjunto de datos que no esté
resumido y luego calcular la media a través de arte.

Siendo 'datos':

> datos
    arte talla ejemplares
1 pareja    30          6
2   beta    31         18
3 pareja    32         18
4   baca    33         12
5    gnr    57         12
6    tnf    35          6

quizá pueda servir:

datos2<-c()
for( j in 1:nrow(datos))
{
     i<-1
     while(i<=ejemplares[j])
    {
          datos2<-rbind(datos2,datos[j,])
     i<-i+1
     }
}

(comprobar con:
>  mean(datos2$talla)
[1] 36.16667
> weighted.mean(datos$talla,datos$ejemplares)
[1] 36.16667
> )

y luego

> by(datos2$talla,datos2$arte,mean)

que te dará la media estratificada por arte. Seguro se puede pulir la
salida, pero eso es harina de otro costal...

Salud.


2011/7/29 <r-help-es-request en r-project.org>
>
> Envíe los mensajes para la lista R-help-es a
>        r-help-es en 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 en r-project.org
>
> Puede contactar con el responsable de la lista escribiendo a:
>        r-help-es-owner en 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. Media Ponderada (jose cebrian)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 29 Jul 2011 09:53:45 +0100 (BST)
> From: jose cebrian <pepeceb en yahoo.es>
> To: "R-help-es en r-project.org" <R-help-es en r-project.org>
> Subject: [R-es] Media Ponderada
> Message-ID:
>        <1311929625.95075.YahooMailNeo en web29608.mail.ird.yahoo.com>
> Content-Type: text/plain
>
> Hola a todos.
> Ya he escrito sobre esto pero vuelvo a no ser capaz de resolver una cosa
>
> Tengo este archivo
>
>
>
> año   mes      fecha  puerto           barco   origen       arte
>
>  2007 enero 2007-01-04 Ribeira  juan        Norte     PAREJA
>  2007 enero 2007-01-04 Ribeira    luis        -Norte     BETA
>  2007 enero 2007-01-04 Ribeira  pedro     Norte     PAREJA
>  2007 enero 2007-01-04 Ribeira    marcosNorte     BACA
>  2007 enero 2007-01-04 Ribeira       nolo    Norte     GNR
> 2007 enero 2007-01-04 Ribeira      bein     Norte     TNF
>
>  p.desem  talla ejemplares
>    1485        30          6
>    1485         31         18
>     1485         32         18
>     1485         33         12
>    1485         57         12
>     1485         35          6
>
>
>
> Puedo hallar la Media Ponderada de la Talla en funcion de los ejemplares (weighted.mean (talla, ejemplares))
> Pero no soy capaz de hallar la media ponderada en funcion del arte
>
> Tal que me quede algo asi
>
>   BACA     BETA   PAREJA
> 40.12        28.66    42.27
>
> Esto me sale con      tapply(talla+ejemplares,arte, weighted.mean, na.rm=TRUE)
>
>
> Lo que pasa es que esto es la media y yo quiero la media ponderada
>
> Bueno, gracias por leerme
> Saludos
> Jocedo
>
>        [[alternative HTML version deleted]]
>
>
>
> ------------------------------
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
> Fin de Resumen de R-help-es, Vol 29, Envío 30
> *********************************************



--
«But Gwindor answered: 'The doom lies in yourself, not in your name.'»

JRR Tolkien