R question - Comment extraire les valeurs d'attributs de l'objet bystat et les placer dans des variables

Question

J'utilise la fonction bystat du paquet Hmisc dans R. Comment puis-je extraire des valeurs d'attribut et les placer dans des variables. Par exemple, je veux calculer la moyenne et le SD pour la variable aaf et les mettre dans une matrice de données.

t <- with(d.aaf,bystats(y=aaf,plot_bid,fun=function(x) { 
     c(Mean = round(mean(x),digits=2),SD = round(sd(x),digits=2)) 
    })) 

> str(t) 
bystats [1:121, 1:3] 5 5 5 5 5 4 5 5 3 4 ... 
- attr(*, "dimnames")=List of 2 
    ..$ : chr [1:121] "P00000000006001288020278" "P00000000006001288085814" 
        "P00000000006001288151350" "P00000000006001288216886" ... 
    ..$ : chr [1:3] "N" "Mean" "SD" 
- attr(*, "heading")= chr "function(x) { c(Mean = round(mean(x),digits=2), 
              SD = round(sd(x),digits=2)) } 
          of aaf by plot_bid" 
- attr(*, "byvarnames")= chr "plot_bid" 

La façon dont je le fais est d'abord convertir « t » en dataframe, que je ne pense pas est très efficace. Merci pour vos suggestions.

Answer 1

Vous pouvez utiliser ddply à partir du package plyr qui est directement transmis à une trame de données.

library(plyr) 
    t<-ddply(d.aaf, "plot_bid", summarise, mean=round(mean(aaf),2), SD=round(sd(aaf),2))  
    SD<-t$SD 
    mean<-t$mean