Eva, The printed output can be dumped to a file using command sink(). See its help page for info.
The oecosimu result is the output of the basic statistic function amended with item "oecosimu" with structure: > str(mod$oecosimu) List of 6 $ z : num -66.5 $ pval : num 0.01 $ simulated : num [1, 1:99] 159621 160463 161865 161720 162472 ... $ method : chr "r00" $ statistic : Named num 54904 ..- attr(*, "names")= chr "statistic" $ alternative: chr "two.sided" You can list and save itmes like mod$oecosimu$statistic, mod$oecosimu$z, mod$oecosimu$pval which all are combined into printed ouput. You can use commands like write.table, write.csv to save these results for external software. Please note the 'alternative' algorithm: it seems to me that your hypotheses may one-sided. Cheers, Jari Oksanen Dr Oksanen I'm checking if the simulated number of checkerboard units of a species pair is larger or smaller than the observed number of checkerboard units in the original dataset. If it is larger I list the species pair as exhibiting cooccurrence and if it is smaller i list it as exhibiting competition. The procedure I'm following (which I'm not sure is correct) is that I use the probabilities i get from the output at the 0.1 level. As I understand, this means that I'm checking at the 0,05 level in each tail of the distribution. Am I wrong? Thank's again Eva -- View this message in context: http://r.789695.n4.nabble.com/lost-in-vegan-package-tp2400145p2403096.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.