[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