Hola
¿Alguno tiene un código para sumarizar medidas?
my.df <- data.frame( "ID" = c( rep("A", 3 ),
rep("B",2),"C", rep("D",3) ),
"Obs.1" = rnorm( 9,0,1 ), "Obs.2" = rnorm( 9,0,3 ) )
Algo que quede mas bonito que esto:
by( my.df, my.df$ID, mean )
Gracias un saludo!
--
Patricia García González
[[alternative HTML version deleted]]
Hola Patricia. Prueba do.call(rbind, by( my.df[,-1], list(ID=my.df$ID), mean)) observa los cambios en la llamada a by. Un saludo Gregorio 2011/5/24 Usuario R <r.user.spain@gmail.com>> Hola > > ¿Alguno tiene un código para sumarizar medidas? > > my.df <- data.frame( "ID" = c( rep("A", 3 ), rep("B",2),"C", rep("D",3) ), > "Obs.1" = rnorm( 9,0,1 ), "Obs.2" = rnorm( 9,0,3 ) ) > > > > > Algo que quede mas bonito que esto: > > by( my.df, my.df$ID, mean ) > > > Gracias un saludo! > -- > Patricia García González > > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es > >-- Dr. Gregorio R. Serrano Dpto. Economía Cuantitativa (UCM) Voz:+34 91394 2361 Fax:+34 91394 2591 http://www.grserrano.es [[alternative HTML version deleted]]
Muchas gracias! El 24 de mayo de 2011 11:14, Gregorio R. Serrano <grserrano@ccee.ucm.es>escribió:> Hola Patricia. > > Prueba > > do.call(rbind, by( my.df[,-1], list(ID=my.df$ID), mean)) > > observa los cambios en la llamada a by. > > Un saludo > Gregorio > > 2011/5/24 Usuario R <r.user.spain@gmail.com> > >> Hola >> >> ¿Alguno tiene un código para sumarizar medidas? >> >> my.df <- data.frame( "ID" = c( rep("A", 3 ), rep("B",2),"C", rep("D",3) ), >> "Obs.1" = rnorm( 9,0,1 ), "Obs.2" = rnorm( 9,0,3 ) ) >> >> >> >> >> Algo que quede mas bonito que esto: >> >> by( my.df, my.df$ID, mean ) >> >> >> Gracias un saludo! >> -- >> Patricia García González >> >> [[alternative HTML version deleted]] >> >> >> _______________________________________________ >> R-help-es mailing list >> R-help-es@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-help-es >> >> > > > -- > Dr. Gregorio R. Serrano > Dpto. Economía Cuantitativa (UCM) > Voz:+34 91394 2361 > Fax:+34 91394 2591 > http://www.grserrano.es >-- Patricia García González [[alternative HTML version deleted]]
Otra opcion con plyr. library(plyr) ddply(my.df, .(ID), numcolwise(mean)) Luciano El 24 de mayo de 2011 06:21, Usuario R <r.user.spain@gmail.com> escribió:> Muchas gracias! > > El 24 de mayo de 2011 11:14, Gregorio R. Serrano > <grserrano@ccee.ucm.es>escribió: > > > Hola Patricia. > > > > Prueba > > > > do.call(rbind, by( my.df[,-1], list(ID=my.df$ID), mean)) > > > > observa los cambios en la llamada a by. > > > > Un saludo > > Gregorio > > > > 2011/5/24 Usuario R <r.user.spain@gmail.com> > > > >> Hola > >> > >> ¿Alguno tiene un código para sumarizar medidas? > >> > >> my.df <- data.frame( "ID" = c( rep("A", 3 ), rep("B",2),"C", rep("D",3) > ), > >> "Obs.1" = rnorm( 9,0,1 ), "Obs.2" = rnorm( 9,0,3 ) ) > >> > >> > >> > >> > >> Algo que quede mas bonito que esto: > >> > >> by( my.df, my.df$ID, mean ) > >> > >> > >> Gracias un saludo! > >> -- > >> Patricia García González > >> > >> [[alternative HTML version deleted]] > >> > >> > >> _______________________________________________ > >> R-help-es mailing list > >> R-help-es@r-project.org > >> https://stat.ethz.ch/mailman/listinfo/r-help-es > >> > >> > > > > > > -- > > Dr. Gregorio R. Serrano > > Dpto. Economía Cuantitativa (UCM) > > Voz:+34 91394 2361 > > Fax:+34 91394 2591 > > http://www.grserrano.es > > > > > > -- > Patricia García González > > [[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]]
Hola, ¿Que tal utilizar aggregate? aggregate(cbind(Obs.1, Obs.2) ~ ID, data=my.df, FUN=mean) Saludos. Oscar. ------------------- Oscar Perpiñán Lamigueiro Dpto. Ingeniería Eléctrica EUITI-UPM http://procomun.wordpress.com