EBSeq/R/beta.mom.R

   1 beta.mom <-
   2 function(qs.in){
   3         xbar<-mean(qs.in)
   4         s2<-var(qs.in)
   5         term<-(xbar*(1-xbar))/s2
   6         alpha.hat<-xbar*(term-1)
   7         beta.hat<-(1-xbar)*(term-1)
   8         return(c(alpha.hat,beta.hat))
   9 }
  10