Thanks Bob / Dylan:

Bob:
I am concerned with the possibility of increase the risk of false 
positives because two tests are actually performing.

What do you think of multiplying p-value by 2 (the number of tests) to 
control risk?

set.seed(10)
rnorm(999)->null.dist
obs<-2.5
alpha<-0.05
pval=mean(obs>null.dist)
min(pval, 1-pval)*2

Dylan:
If tried with this, but it seems that wilcox.test need two variables of 
the same length

wilcox.test(obs, null.dist,alternative=c("two.sided"),paired=T,conf.int=T)
Error en wilcox.test.default(obs, null.dist, alternative = 
c("two.sided"),  :
   'x' and 'y' must have the same length

Best wishes
Tano
-- 
*Cayetano Gutiérrez Cánovas*
Aquatic Ecology Group <http://www.um.es/ecoaqua>
Universidad de Murcia (Spain)

website <http://www.um.es/ecoaqua/personal/Gutierrez-Canovas_en.html>
Scholar Google Citations 
<http://scholar.google.com/citations?user=Yowrt8gAAAAJ>
Lessons from nature <http://lessonsfromnature.org/es/>
skype: tano.gc
twitter: @LdN_Esp <https://twitter.com/LdN_Esp>

	[[alternative HTML version deleted]]

