This is one way. I read your two columns of data into two variables, x and y. Then:
test <- replicate(5000, cor(x, sample(y))) > hist(test) Kjetil On Thu, Oct 28, 2010 at 9:49 AM, Chitra <cbban...@gmail.com> wrote: > >> Dear R experts, >>I am sorry for my inability. >>I have the following dataset: > > > Qtot Itot >>1 73 684 >>2 64 451 >>3 71 378 >>4 65 284 >>5 47 179 >>6 31 117 >>7 19 69 >> >>Now I need to perform Monte Carlo Pertutation test underlaying the > following condition. >> >> >>Condition >> >>In order to choose randomly (5000 times) for the Qtot in the interval 0 to > the maximum possible value of >> >>Qtot for each observed value of the independent variable (Itot), Monte > Carlo tests needed to be perform to test the null hypothesis. The >>maximum observed value is either the calculated size of the regional pool > (Interpolated value), or the measured size of the actual >> >>pool. Each time correlation coefficient (r) between independent and > dependent variable has to be calculated in order to achieve >> >>the empirical distribution of r for the null hypothesis conditions. >> >>Now please help me how to do? >>I have also intalled combinat library in R. But could not perform by > myself. >>I will be highly grateful. >>Thanks. >> > --> > Chitra Bahadur Baniya (PhD) > Lecturer > Central Department of Botany > Tribhuvan University > Kirtipur > Kathmandu, Nepal > > -- > View this message in context: > http://r.789695.n4.nabble.com/Please-help-me-about-Monte-Carlo-Permutation-tp3017131p3017131.html > Sent from the R help mailing list archive at Nabble.com. > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > 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. > ______________________________________________ R-help@r-project.org mailing list 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.
