Тестирование на значительное увеличение частоты периодов времени, когда ветер дует восточно в течение 5 дней или более подряд - PullRequest
1 голос
/ 04 октября 2019

Я пытаюсь выяснить, значительно ли увеличивается количество случаев, когда направление ветра (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
> 
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...