bayesHierVar.stHM <-
function(x){
K=ncol(x) #### number of observations
v=K-1
I=nrow(x) ### number of genes
s2est=apply(x,1,var)
s2bar=mean(s2est)
S=var(s2est)*(I-1)
Best=(2/v)/(1+2/v)*(I-1)/I + 1/(1+2/v)*(2/v)*s2bar^2*(I-1)/S
if(Best<1){
var_hier_est=(1-Best)*s2est+Best*s2bar
result=var_hier_est
} else{
result=rep(s2bar,I)
names(result) = names(s2est)
}
return(result)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.