Hi Ani,
Not sure about this, but is this what you want?
ann10<-aggregate(mon10[3:12],list(mon10$Group.1),which.max)
Jim
On Wed, Dec 16, 2020 at 4:55 PM ani jaya <gaaauul at gmail.com>
wrote:>
> Dear R-Help,
>
> I have a data frame containing monthly maxima of rainfall in 10
> locations for 30 year and want to look at a month where an annual
> maxima happens. I can get the annual maxima using aggregate. I try to
> extract the month using the code below.
>
> ann10<-aggregate(mon10[3:12],list(mon10$Group.1),max,na.rm=T)
> when<-matrix(NA,30,10)
> dum<-list()
> for(j in 1:10){
> for(i in 1:30){
> dum<-which(mon10[,j+2]==ann10[i,j+1])
> when[i,j]<-mon10$Group.2[dum[1]]
> }
> }
>
> I am aware that the code above is wrong when multiple max value
> detected (so I use dum[1] there). Any elegant solution how to find the
> month?? Thank you
>
> dput(head(mon10,36))
>
> structure(list(Group.1 = structure(c(1L, 1L, 1L, 1L, 1L, 1L,
> 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,
> 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label =
c("1985",
> "1986", "1987", "1988", "1989",
"1990", "1991", "1992", "1993",
> "1994", "1995", "1996", "1997",
"1998", "1999", "2000", "2001",
> "2002", "2003", "2004", "2005",
"2006", "2007", "2008", "2009",
> "2010", "2011", "2012", "2013",
"2014"), class = "factor"), Group.2 > structure(c(1L,
> 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L,
> 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 1L, 2L, 3L, 4L, 5L, 6L, 7L,
> 8L, 9L, 10L, 11L, 12L), .Label = c("01", "02",
"03", "04", "05",
> "06", "07", "08", "09",
"10", "11", "12"), class = "factor"),
> V1 = c(0, 0, 0, 0, 0, 0, 1, 13, 0, 0, 0, 0, 175, 123, 76,
> 112, 35, 52, 14, 28, 63, 68, 30, 67, 107, 46, 44, 22, 33,
> 25, 10, 5, 50, 10, 32, 59), V2 = c(167, 50, 48, 114, 128,
> 37, 54, 5, 23, 36, 10, 43, 175, 123, 110, 112, 58, 12, 18,
> 171, 63, 68, 30, 79, 189, 40, 35, 45, 23, 10, 0, 1, 1, 5,
> 34, 71), V3 = c(63.0054, 68.9842, 1, 115.9904, 5.796, 53.0929,
> 63.0092, 18.3692, 0, 30.1854, 12.214, 55.0721, 91.9728, 64.9951,
> 28.431, 21.6815, 38.4651, 48.1583, 1, 1, 1, 7.654, 50.1237,
> 86.9713, 69.9766, 14.2489, 45.2282, 12.4218, 18.5532, 23.9453,
> 38.4967, 0, 0, 1.8495, 41.3493, 69.98), V4 = c(112, 94, 70,
> 18, 6.1, 55, 84.1, 7.1, 0, 4.1, 24.9, 26.9, 42, 38, 38, 62,
> 26, 50, 8, 0, 6, 40, 18, 77, 30, 78, 31, 28, 13, 27, 50,
> 0, 0, 0, 29, 84), V5 = c(114.2545, 96.2488, 71.9131, 14.299,
> 22.2288, 55.4017, 1.971, 3.69, 0, 13.301, 57.5355, 27.3004,
> 39.6649, 34.3766, 34.3278, 63.1532, 21.4376, 49.5373, 4.556,
> 0, 2.7328, 36.8717, 66.3924, 79.1253, 25.1753, 80.1422, 26.2089,
> 23.232, 9.4119, 22.3628, 49.4622, 0, 0, 0, 24.1889, 86.2128
> ), V6 = c(65, 35, 83, 46, 50, 26, 28, 10, 6, 13, 90, 97,
> 178, 138, 79, 47, 30, 73, 5, 0, 1, 24, 55, 46, 80, 62, 56,
> 57, 21, 18, 0, 0, 0, 9, 36, 175), V7 = c(110, 100, 100, 63,
> 167, 40, 5, 7, 76, 0, 15, 17, 55.1, 116.1, 57.9, 21.1, 0.3,
> 103.1, 14, 0, 5.1, 5.1, 1, 20.1, 31, 26.9, 0.5, 0, 27.9,
> 41, 0, 6.1, 0, 0, 0, 0), V8 = c(78, 87, 37, 98, 0, 0, 0,
> 0, 0, 0, 0, 0, 64, 90, 89, 73, 104, 70, 101, 15, 0, 0, 0,
> 0, 0, 1, 38, 97, 114, 28, 0, 0, 0, 0, 0, 0), V9 = c(29, 17,
> 72, 29, 16, 25, 23, 43, 37, 10, 15, 7, 14, 77, 54, 43, 22,
> 58, 63, 27, 12, 12, 16, 9, 22, 33, 60, 54, 40, 37, 35, 12,
> 1, 4, 11, 49), V10 = c(63.1161, 11.9425, 35.3226, 34.0615,
> 92.9325, 34.0785, 39.009, 21.5103, 17.9209, 50.9037, 32.7218,
> 52.2346, 63.1381, 99.946, 48.0861, 63.1813, 43.7723, 86.9132,
> 39.1252, 49.6261, 40.7285, 112.9464, 49.8403, 66.5253, 21.5789,
> 12.0566, 89.927, 71.7187, 37.816, 5.3817, 1, 0, 0, 4.636,
> 40.7964, 88.9239)), row.names = c(1L, 31L, 61L, 91L, 121L,
> 151L, 181L, 211L, 241L, 271L, 301L, 331L, 2L, 32L, 62L, 92L,
> 122L, 152L, 182L, 212L, 242L, 272L, 302L, 332L, 3L, 33L, 63L,
> 93L, 123L, 153L, 183L, 213L, 243L, 273L, 303L, 333L), class =
"data.frame")
>
> sessionInfo()
> R version 4.0.2 (2020-06-22)
> Platform: x86_64-w64-mingw32/x64 (64-bit)
> Running under: Windows 10 x64 (build 18363)
>
>
> Best,
> Ani
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.