pval2FDR: function to do ...
In Senbee/mwt: Moderated Welch Test

Description Usage Arguments Details Value Note Author(s) References See Also Examples

~~ A concise (1-5 lines) description of what the function does. ~~

1	pval2FDR(pval, lim = 0.7)

`pval`	~~Describe `pval` here~~
`lim`	~~Describe `lim` here~~

~~ If necessary, more details than the description above ~~

~Describe the value returned If it is a LIST, use

`comp1`	Description of 'comp1'
`comp2`	Description of 'comp2'

...

~~further notes~~

~Make other sections like Warning with Warning.... ~

~~who you are~~

~put references to the literature/web site here ~

~~objects to See Also as help, ~~~

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function(pval, lim = 0.7)
{
  n= length(pval)
  Fp = rank(pval)/length(pval)
  p0 = sum(pval>lim)/((1-lim)*n)
  p0 = min(p0, 1)
  
  FDRp = p0 * pmin(pval/Fp, 1)
  ord = order(pval)
  FDR.o = FDRp[ord]
  b = rev(cummin(rev(FDR.o)))
  FDR = rep(0, n)
  FDR[ord] = b
  return(FDR)
  }