R/logproby.R

Defines functions logproby

logproby <-
function(y,XX,ZZ,aa,bb,s2,nu){
  df=nu
  n=length(y)
  U=aa*ZZ+bb*XX+diag(n)
  S=s2*U
  log.prob=dmt(y,mean=rep(0,n),df=df,S=S,log=TRUE)
  log.prob
}

Try the PAGWAS package in your browser

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

PAGWAS documentation built on May 30, 2017, 6:40 a.m.