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]]
_______________________________________________
R-sig-ecology mailing list
R-sig-ecology@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
