R/MeanA.SGB.R

In SGB: Simplicial Generalized Beta Regression

Defines functions MeanA.SGB

Documented in MeanA.SGB

MeanA.SGB <-
function(shape1,scale,shape2){
    pdig <- digamma(shape2)/shape1
    np <- length(shape2)
    EAu <- scale%*%diag(exp(pdig),nrow=np,ncol=np)
    EAu <- EAu/rowSums(EAu) 
    colnames(EAu) <- colnames(scale)
    return(EAu)
  }

Try the SGB package in your browser