Included EBSeq for downstream differential expression analysis

[rsem.git] / EBSeq / R / beta.mom.R

diff --git a/EBSeq/R/beta.mom.R b/EBSeq/R/beta.mom.R
new file mode 100644 (file)
index 0000000..269996d
--- /dev/null
+++ b/EBSeq/R/beta.mom.R
@@ -0,0 +1,10 @@
+beta.mom <-
+function(qs.in){
+       xbar<-mean(qs.in)
+       s2<-var(qs.in)
+       term<-(xbar*(1-xbar))/s2
+       alpha.hat<-xbar*(term-1)
+       beta.hat<-(1-xbar)*(term-1)
+       return(c(alpha.hat,beta.hat))
+}
+