Kim Milferstedt
2006-Apr-19 17:14 UTC
[R] Hmisc + summarize + quantile: Why only quantiles for first variable in data frame?
Hi, I'm working on a data set that contains a couple of factors and a number of dependent variables. From all of these dependent variables I would like to calculate mean, standard deviation and quantiles. With the function FUN I get all the means and stdev that I want but quantiles are only calculated for the first of the dependent variables (column 8 in the summarize command). What do I have to do differently in order to get all the quantiles that I want? Thanks, Kim sgldm2 <- read.table("E:/analysistemp/060412_test_data2.txt", header=T) attach(sgldm2) names(sgldm2) FUN <- function(x)c(Mean=mean(x,na.rm=TRUE), STDEV=sd(x,na.rm=TRUE), Quantile=quantile(x, probs= c(0.25,0.50,0.75),na.rm=TRUE)) ordering <- llist(time_h_f, Distance_f) resALL <- summarize(sgldm2[,8:10], ordering, FUN)
Frank E Harrell Jr
2006-Apr-19 18:48 UTC
[R] Hmisc + summarize + quantile: Why only quantiles for first variable in data frame?
Kim Milferstedt wrote:> Hi, > > I'm working on a data set that contains a couple of factors and a > number of dependent variables. From all of these dependent variables > I would like to calculate mean, standard deviation and quantiles. > With the function FUN I get all the means and stdev that I want but > quantiles are only calculated for the first of the dependent > variables (column 8 in the summarize command). What do I have to do > differently in order to get all the quantiles that I want? > > Thanks, > > Kim > > sgldm2 <- read.table("E:/analysistemp/060412_test_data2.txt", header=T) > attach(sgldm2) > names(sgldm2) > > FUN <- function(x)c(Mean=mean(x,na.rm=TRUE), > STDEV=sd(x,na.rm=TRUE), Quantile=quantile(x, probs= > c(0.25,0.50,0.75),na.rm=TRUE)) > ordering <- llist(time_h_f, Distance_f) > > resALL <- summarize(sgldm2[,8:10], ordering, FUN)Please read the documentation and see the examples. The first argument to summarize is a matrix or vector and if a matrix, FUN must use matrix operations if you want column-by-column results. FH> > ______________________________________________ > R-help at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html >-- Frank E Harrell Jr Professor and Chair School of Medicine Department of Biostatistics Vanderbilt University