Lyndz
2019-May-22 02:44 UTC
[R] Significance of difference in means using Monte Carlo Simulation in R
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 passagefrequency are prepared.>2). Averages of the former 23-yr values (1961-1983) minus those of thelatter 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]]
Bert Gunter
2019-May-22 14:39 UTC
[R] Significance of difference in means using Monte Carlo Simulation in R
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 at 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 at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]