On Tue, Aug 7, 2012 at 8:05 AM, Dominic Roye <dominic.roye at gmail.com>
wrote:> 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.