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

Bert Gunter bgunter@4567 @end|ng |rom gm@||@com
Wed May 22 16:39:02 CEST 2019

```This list is for Help and usually expects you to first show us your own
efforts. You need to do your own homework and spend some time with one or
more of the many good R tutorials on the Web. Some suggestions canbe found
here, though you will find many more by simple web searches:

https://www.rstudio.com/online-learning/#r-programming

Cheers,
Bert

Bert Gunter

"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )

On Wed, May 22, 2019 at 12:56 AM Lyndz <olagueralyndonmark429 using gmail.com>
wrote:

> 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]]
>
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help