R/rGenWish.R

rGenWish <-
function(Sigma,df){
    n=nrow(Sigma)
    Y=t(rmvnorm(df,sigma=Sigma))
    S=Y%*%t(Y)
    onevec=matrix(1,n,1)
    D.gen=-2*S+diag(S)%*%t(onevec)+onevec%*%t(diag(S))
    D.gen
    }

Try the rwc package in your browser

Any scripts or data that you put into this service are public.

rwc documentation built on May 2, 2019, 3:34 p.m.