R help - Mar 2008 - summary of summaries

Hi, 

 

I have a few hundreds files with numerical information of different length
but with the same column structure. I use the following code to get summary
statistics

 

fplist <- list.files(pattern=".*analysis")

for (fp in fplist){

            x2 <- read.delim(fp)

            summary(x2)

}

 

Summary gives something like:

summary (x2)

       V1              V2            V3              V4                  V5
V6                V7               V8               V9        

 Min.   :11.00   4:9    : 11   Min.   : 1.00   Min.   : 19128823   Min.   :
19812573   Min.   :   0.00   Min.   :  1.00   Min.   :0.1538   Min.   :
1.000  

 1st Qu.:13.00   2:9    :  9   1st Qu.: 4.50   1st Qu.: 48952449   1st Qu.:
49085044   1st Qu.:  14.33   1st Qu.:  3.00   1st Qu.:0.7223   1st Qu.:
5.500  

 Median :20.00   3:10   :  9   Median :10.00   Median : 63129783   Median :
63218749   Median :  61.54   Median : 10.00   Median :1.0069   Median :
8.000  

 Mean   :25.68   4:8    :  9   Mean   :12.20   Mean   : 72876334   Mean   :
73051728   Mean   : 175.39   Mean   : 32.64   Mean   :1.0522   Mean   :
9.644  

 3rd Qu.:34.00   3:8    :  7   3rd Qu.:23.00   3rd Qu.: 96700016   3rd Qu.:
96713037   3rd Qu.: 264.88   3rd Qu.: 40.00   3rd Qu.:1.2586   3rd
Qu.:12.000  

 Max.   :72.00   4:13   :  6   Max.   :23.00   Max.   :203343850   Max.
:204024777   Max.   :1440.30   Max.   :289.00   Max.   :2.7690   Max.
:27.000  

                 (Other):300


      V10       

 Min.   : 6.00  

 1st Qu.:10.00  

 Median :14.00  

 Mean   :18.03  

 3rd Qu.:24.00  

 Max.   :52.00  

 

Now I want to get a summary statistics of the summaries. Obviously I would
need to store the summary results inarrays. But I was wondering if that
makes sense and how an expert would solve this? I don't want to create an
array for each column if possible.

 

Thanks for the kind support,

 

Bernd


	[[alternative HTML version deleted]]

Try this:

summary(do.call(rbind, lapply(list(df1, df2), sapply, summary)))

On 17/03/2008, Bernd Jagla <baj2107 at columbia.edu>
wrote:
> Hi,
>
>
>
>  I have a few hundreds files with numerical information of different length
>  but with the same column structure. I use the following code to get
summary
>  statistics
>
>
>
>  fplist <- list.files(pattern=".*analysis")
>
>  for (fp in fplist){
>
>             x2 <- read.delim(fp)
>
>             summary(x2)
>
>  }
>
>
>
>  Summary gives something like:
>
>  summary (x2)
>
>        V1              V2            V3              V4                  V5
>  V6                V7               V8               V9
>
>   Min.   :11.00   4:9    : 11   Min.   : 1.00   Min.   : 19128823   Min.  
:
>  19812573   Min.   :   0.00   Min.   :  1.00   Min.   :0.1538   Min.   :
>  1.000
>
>   1st Qu.:13.00   2:9    :  9   1st Qu.: 4.50   1st Qu.: 48952449   1st
Qu.:
>  49085044   1st Qu.:  14.33   1st Qu.:  3.00   1st Qu.:0.7223   1st Qu.:
>  5.500
>
>   Median :20.00   3:10   :  9   Median :10.00   Median : 63129783   Median
:
>  63218749   Median :  61.54   Median : 10.00   Median :1.0069   Median :
>  8.000
>
>   Mean   :25.68   4:8    :  9   Mean   :12.20   Mean   : 72876334   Mean  
:
>  73051728   Mean   : 175.39   Mean   : 32.64   Mean   :1.0522   Mean   :
>  9.644
>
>   3rd Qu.:34.00   3:8    :  7   3rd Qu.:23.00   3rd Qu.: 96700016   3rd
Qu.:
>  96713037   3rd Qu.: 264.88   3rd Qu.: 40.00   3rd Qu.:1.2586   3rd
>  Qu.:12.000
>
>   Max.   :72.00   4:13   :  6   Max.   :23.00   Max.   :203343850   Max.
>  :204024777   Max.   :1440.30   Max.   :289.00   Max.   :2.7690   Max.
>  :27.000
>
>                  (Other):300
>
>
>       V10
>
>   Min.   : 6.00
>
>   1st Qu.:10.00
>
>   Median :14.00
>
>   Mean   :18.03
>
>   3rd Qu.:24.00
>
>   Max.   :52.00
>
>
>
>  Now I want to get a summary statistics of the summaries. Obviously I would
>  need to store the summary results inarrays. But I was wondering if that
>  makes sense and how an expert would solve this? I don't want to create
an
>  array for each column if possible.
>
>
>
>  Thanks for the kind support,
>
>
>
>  Bernd
>
>
>         [[alternative HTML version deleted]]
>
>  ______________________________________________
>  R-help at r-project.org mailing list
>  https://stat.ethz.ch/mailman/listinfo/r-help
>  PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
>  and provide commented, minimal, self-contained, reproducible code.
>

-- 
Henrique Dallazuanna
Curitiba-Paran?-Brasil
25? 25' 40" S 49? 16' 22" O