From: <matgopa1_at_umbc.edu>

Date: Mon 20 Feb 2006 - 16:29:51 EST

*>
*

*> summarystat<-function(x)
*

*>
*

*> test<-summarystat(x);test
*

*>
*

R-help@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 Received on Mon Feb 20 16:37:05 2006

Date: Mon 20 Feb 2006 - 16:29:51 EST

Hello all,

I have a simple function which calculates summary statistics of a dataset in terms of a factor (say area).

*> x = data.frame(Area = c(rep("cleanup", 5), rep("ref", 5)), TcCB =
*

c(rnorm(5)+2, rnorm(5)));x

Area TcCB

1 cleanup 2.5829747 2 cleanup 2.6796868 3 cleanup 2.5437094 4 cleanup 2.8453616 5 cleanup 1.1789683 6 ref 1.0140391 7 ref -0.8433729 8 ref 0.6512422 9 ref 0.2341083 10 ref -0.2688026

+ { + no.samples<-by(x,x$Area,function(x) length(x$TcCB)) + mean<-by(x,x$Area,function(x) mean(x$TcCB)) + quantile<-by(x,x$Area,function(x) summary(x$TcCB)) + stdev<-by(x,x$Area,function(x) sd(x$TcCB)) + final<-do.call("cbind",c(quantile,mean,stdev,no.samples)) + return(final) + }

cleanup ref cleanup ref cleanup ref cleanup ref Min. 1.179 -0.8434 2.36614 0.1574428 0.6737748 0.736001 5 5 1st Qu. 2.544 -0.2688 2.36614 0.1574428 0.6737748 0.736001 5 5 Median 2.583 0.2341 2.36614 0.1574428 0.6737748 0.736001 5 5 Mean 2.366 0.1574 2.36614 0.1574428 0.6737748 0.736001 5 5 3rd Qu. 2.680 0.6512 2.36614 0.1574428 0.6737748 0.736001 5 5 Max. 2.845 1.0140 2.36614 0.1574428 0.6737748 0.736001 5 5

Now the results are arranged as per quantile function.

When the results are printed, I would like to have the results for mean, stdev, no.samples, quantiles one after the other with the function names for the two factors namely cleanup and reference. Can somebody help in doing so? I did refer to a previous thread on formatting results from function. But my case is little different.

Moreover, i would like to know, if there is any R command which produces this general summary statistics as above?

Thanks for your time.

Mathangi

R-help@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 Received on Mon Feb 20 16:37:05 2006

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:42:36 EST
*