Hist function

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120807/aac0c8c8/attachment.pl>
Hello everyone,

First i explain my aim. I would like to calculate the frecuency of
temperature (datos.mx1, columns 3-6)  for each month with the thresholds of
table lf.med and lc.med.

numero1 <- hist(subset(datos.mx1[,3],
datos.mx1$Mes==1),plot=FALSE,breaks=c(min(datos.mx1[,3]),lf.med[1,2],lc.med[1,2],max(datos.mx1[,3])))
# First month and only the first weather station (column 3)

I now have to apply the above command for all 12 months at all the weather
stations. I've been trying to get with a for-loop, but without success.

How can i do it the best way?
I'd suggest the tapply() function. If you could use dput() to make a
small reproducible example of your data, I'll help you set up the
exact syntax, but I think it's something like this:

tapply(datos.mx1, datos.mx1$Mes, function(x) hist(x[,3], plot = FALSE,
breaks = c(min(x[,3], lf.med[1,2], lc.med[1,2], max(x[,3]))))

Best,
Michael
After this command, i'm only interested in the part of $counts. I hope
someone can help me. Any help is greatly appreciated!!!  Thank you,
everybody`!!

datos.mx1
          Fecha Mes Serra.da.Foladoira Santiago Sergude Rio.Do.Sol
1    01.01.2006   1                4.1     14.4     9.3        5.7
2    02.01.2006   1                5.5      9.8    11.0        7.6
3    03.01.2006   1                9.9     13.0    12.7        9.9
4    04.01.2006   1                9.7     13.3    12.7        9.7
5    05.01.2006   1               10.2     12.2    13.1        8.0
6    06.01.2006   1                4.2     14.2    10.2        5.4
7    07.01.2006   1                1.0      6.9     3.2        3.2
.
.
.......................etc.

lf.med
   Mes Serra.da.Foladoira  Santiago   Sergude Rio.Do.Sol
1    1           7.833702  8.811945  8.122162   7.899513
2    2           7.766978  8.592693  8.041483   7.887327
3    3           7.864609  8.861457  8.385405   8.080646
4    4           8.204098  8.986149  8.803744   8.438141
5    5           8.716766  9.317366  9.385637   8.977087
6    6           9.377651  9.874497 10.194340   9.686213
7    7           9.838470 10.055430 10.447219   9.991618
8    8          10.046204 10.099868 10.520872  10.098330
9    9           9.780911 10.038629 10.156048   9.851299
10  10           9.425308  9.779886  9.646881   9.496192
11  11           8.417168  9.346191  8.791033   8.549046
12  12           7.968706  8.795800  8.038208   8.017603

lc.med
   Mes Serra.da.Foladoira Santiago  Sergude Rio.Do.Sol
1    1           21.29353 26.62562 23.71784   22.04786
2    2           21.97787 27.07213 25.68340   23.39020
3    3           22.69749 27.42612 26.97889   24.33741
4    4           24.49683 29.21054 29.76726   26.59280
5    5           25.24233 29.83056 30.75132   27.28277
6    6           26.46742 30.70485 33.02328   28.77892
7    7           27.72731 30.24450 34.28162   29.97767
8    8           27.64356 30.84850 34.52161   30.00666
9    9           27.21998 31.18390 33.19776   29.54618
10  10           25.35830 28.83195 29.63192   26.42831
11  11           22.37895 26.38226 25.69239   23.47575
12  12           21.44064 26.09142 23.99676   22.26208

        [[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.