2010/10/29 Chitra <cbban...@gmail.com> > > £ukasz Rêc³awowicz, > Thanks. > I got this p-value. >
It's exact upper tail (>=), for lower you have to add (<(-test.statistic)) to obtain two-sided p-vaule. > > test() > [1] P-value: 0.00139 > > I still could not figure out how to plot 5000 permuted Qtot vs Itot. > Result of permuted.Qtot<-permn(data$Qtot) is the list class object, so we can plot something very strange like a plot(rep(data$Itot,5040)~unlist(permuted.Qtot)) But I don't think it's what you want to get... -- Mi³ego dnia [[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.