logratio {nethet} | R Documentation |
Log-likelihood-ratio statistics used in Differential Network
logratio(x1, x2, x, sig1, sig2, sig, mu1, mu2, mu)
x1 |
data-matrix sample 1 |
x2 |
data-matrix sample 2 |
x |
pooled data-matrix |
sig1 |
covariance sample 1 |
sig2 |
covariance sample 2 |
sig |
pooled covariance |
mu1 |
mean sample 1 |
mu2 |
mean sample 2 |
mu |
pooled mean |
Returns a list with named elements 'twiceLR', 'sig1', 'sig2', 'sig'. 'twiceLR' is twice the log-likelihood-ratio statistic.
n.stadler
x1=matrix(rnorm(100),50,2) x2=matrix(rnorm(100),50,2) logratio(x1,x2,rbind(x1,x2),diag(1,2),diag(1,2),diag(1,2),c(0,0),c(0,0),c(0,0))$twiceLR