est.hyper: 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	est.hyper(z, D, d12)

`z`	~~Describe `z` here~~
`D`	~~Describe `D` here~~
`d12`	~~Describe `d12` 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(z,D,d12){
  f2 = function(d0,D) var(z) - trigamma(D/2) - trigamma(d0/2)
  lim = f2(100,D)
  if (lim<0) D= d12
  d0.est = uniroot(f2,c(1,100), D=D)$root
  s2.est = exp(mean(z) -digamma(D/2) +digamma(d0.est/2) - log(d0.est/D) )
  return(list(d0=d0.est, s2=s2.est))
  }