exportPattern("^[[:alpha:]]+")
importFrom("Matrix",t,solve,diag,colSums,crossprod, tcrossprod)

