R/Gibbs4.R

In BAGS: A Bayesian Approach for Geneset Selection

Gibbs4<-
function(noRow,noCol,iter,GrpSzs,YMu,L0,V0,L0A,V0A,MM,AAPi,ApriDiffExp,result1,result2,result3)
     {
output=.C("Gibbs4",as.integer(noRow),as.double(noCol),as.integer(iter),as.double(GrpSzs),as.double(YMu),as.double(L0),as.double(V0),as.double(L0A),as.double(V0A),as.double(MM),as.double(AAPi),as.double(ApriDiffExp),as.double(result1),as.double(result2),as.double(result3))
      Resultado1=matrix(output[[13]],noRow,iter)
      Resultado2=matrix(output[[14]],noRow,iter)
      Resultado3=matrix(output[[15]],noRow,iter)
      
      Resultado=list(Resultado1,Resultado2,Resultado3)
      names(Resultado)=c("alfa.1","alfa.2","alfa.3")
      return(Resultado)
     }