R/st.m0.estimate.r

In mcp.project: Applies different False Discovery Rate controlling procedures

Defines functions st.m0.estimate

Documented in st.m0.estimate

st.m0.estimate=function(sorted,m){
  k=2:m
  m0.k=(m+1-k)/(1-sorted[-1])
  diffs=diff(m0.k,lag=1)
  indicators= diffs>0
  optimal.k=min( k[indicators] )
  m0=min(ceiling(m0.k[optimal.k]),m)
  stopifnot(m0 <=m ,m0 > 0)
  return(m0)
}

Try the mcp.project package in your browser