nethet (version 1.4.2)

logratio: Log-likelihood-ratio statistics used in DiffNet

Description

Log-likelihood-ratio statistics used in Differential Network

Usage

logratio(x1, x2, x, sig1, sig2, sig, mu1, mu2, mu)

Arguments

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

Value

• Returns a list with named elements 'twiceLR', 'sig1', 'sig2', 'sig'. 'twiceLR' is twice the log-likelihood-ratio statistic.

Examples

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