Any p-value significant example

any-pvalue-significant.R

## Load library/set seed
set.seed(1325)
library(qvalue)
library(RSkittleBrewer)
trop = RSkittleBrewer("tropical")

## Set a target p-value
target_pval = 0.8

## Not significant if we 
## test with a set of uniform p-values
otherpvals = runif(1000)
pvals = c(target_pval,otherpvals)
hist(pvals,col=trop[2],xlab="P-value",main="",breaks=20)
qvalue(pvals)$qvalue[1]


## But if we make all the
## other p-values really significant
## then whoa! 

otherpvals = rbeta(1000,1,30)
pvals = c(target_pval,otherpvals)
hist(pvals,col=trop[2],xlab="P-value",main="")
qvalue(pvals)$qvalue[1]