Я пытаюсь выяснить, значительно ли увеличивается количество случаев, когда направление ветра (U) <(- 3,54) в течение> = 5 дней подряд значительно увеличивается со временем на уровне 10%. Ниже приведен код, который я использовал для получения количества экземпляров в год в Austral Summer. Я уже проверил на значительное увеличение с использованием Excel на уровне 5%, но я хотел бы узнать, как сделать это с помощью R и на уровне 10%. Любая помощь будет принята с благодарностью.
####WIND PERIOD CALCULATIONS###########################################################################
#### U <(-3.54)
###Austral summer#######################################################################################
library(dplyr)
AusSum_Wind <- data.frame(year=character(), instances=integer(), stringsAsFactors=FALSE)
RowNum <- 1
for (i in 1993:2015) {
AusSum_Wind[RowNum,1] <- paste(as.character(i), as.character(i+1), sep = "-")
wind %>%
filter((Month >= 10 & Year == i) | (Month <= 3 & Year == (i+1))) %>%
mutate(threshold = U < (-3.54),
group = cumsum(threshold != lag(threshold, default = FALSE))) %>%
group_by(group) %>%
mutate(n_days = n()) %>%
summarise_all(first) %>%
filter(threshold, n_days >= 5) %>%
select(-group, -threshold) -> instances
AusSum_Wind[RowNum,2] <- nrow(instances)
RowNum <- RowNum + 1
}
AusSum_Wind
> AusSum_Wind
year instances
1 1993-1994 1
2 1994-1995 3
3 1995-1996 3
4 1996-1997 1
5 1997-1998 5
6 1998-1999 3
7 1999-2000 4
8 2000-2001 2
9 2001-2002 1
10 2002-2003 0
11 2003-2004 3
12 2004-2005 3
13 2005-2006 1
14 2006-2007 1
15 2007-2008 1
16 2008-2009 0
17 2009-2010 3
18 2010-2011 5
19 2011-2012 1
20 2012-2013 5
21 2013-2014 4
22 2014-2015 3
23 2015-2016 2
>