gammaHat: calculate gamma hat

In mnbem: Matrix Normal Block EM

Usage

gammaHat(vmat, w)

Arguments

vmat
w

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function( vmat, w )
{
    m <- nrow( w ); n <- nrow( vmat ) / m; mm <- ( 0:( n - 1 ) ) * m; ghat <- matrix( 0, n, n )
    for( j in 1:m ) {
        for( l in 1:m ){ 
            ghat <- ghat + w[ j, l ] * vmat[ mm + j, mm + l ]
        }
    }
return(ghat)
  }

mnbem documentation built on May 2, 2019, 6:16 p.m.