[R] Significance of difference in means using Monte Carlo Simulation in R

Lyndz o|@guer@|yndonm@rk429 @end|ng |rom gm@||@com
Wed May 22 04:44:23 CEST 2019


I want to implement a statistical test by the following paper in R:
https://journals.ametsoc.org/doi/full/10.1175/JCLI4217.1

**Details**

The above paper calculates the significance of the difference in means
between two periods :1961-1983 and 1984-2000 of tropical cyclone passage
frequency (not-normally distributed) using Monte Carlo simulation. The
following steps are provided:

>1). First, 9999 randomly sorted 40-yr time series of the typhoon passage
frequency are prepared.

>2). Averages of the former 23-yr values (1961-1983) minus those of the
latter 17-yr values are calculated.

>3). From the rank of the original difference value among 10000 samples,
the significance level is estimated.

Suppose I have the following 40-yr time series.

    > dat<-floor(runif(40, min=0, max=20))
    > dput(dat)
    c(15, 14, 1, 16, 0, 18, 5, 8, 19, 7, 11, 15, 2, 17, 12, 16, 1,
    9, 9, 19, 12, 17, 15, 10, 1, 11, 10, 12, 17, 10, 4, 2, 9, 10,
    9, 5, 13, 6, 0, 17)

**PROBLEMS**

1. I am new in R and I am not sure how to implement this test.
Any suggestion on how to do this in R?
I'll appreciate any help.

*Lyndz*

	[[alternative HTML version deleted]]



More information about the R-help mailing list