iProMix.LRT: Calculate the likelihood ratio test statistics (as well as...
In JiayiJi/iProMix: Cell-Type Specific Association Analysis for Bulk Profiling Data

Description Usage Arguments Value Examples

Calculate the likelihood ratio test statistics (as well as chisq p-value) for iProMix

iProMix.LRT(
  y,
  x,
  cov = NULL,
  pi,
  reduce1 = c(2, 1),
  reduce2 = NULL,
  tuningPar = 1e-06,
  diffNum = 1e-04,
  numitersNum = 200
)

`y`	The quantitative measure (e.g. protein/expression) of a gene
`x`	The quantitative measure of anther gene (e.g. ACE2 protein levels) that we would like to know their cell-type specific dependency with Y
`cov`	The covariates for adjustment. Their impact on the mean value of X and Y are adjusted
`pi`	The proportion of cell type 1
`reduce1`	A index of the row and column of the variance-covariance matrix that should be forced to be zero in cell type 1. Default is NULL.
`reduce2`	A index of the row and column of the variance-covariance matrix that should be forced to be zero in cell type 2. Default is NULL.
`tuningPar`	Default is 1e-8. It is used in the embedded graphic lasso procedure for estimating correlation. A larger tuningPar can be selected if one is interested in penalized estimates.
`diffNum`	Default 0.0001. The convergency criterion for EM algorithm.
`numitersNum`	Default 200. The number of iterations in EM algorithm.

list with 10 elements. It contains

`var1:`	The estimated (y,x) variance of cell type 1
`var2:`	The estimated (y,x) variance of cell type 2
`mu1:`	The estimated (y,x) mean function of cell type 1
`mu2:`	The estimated (y,x) mean function of cell type 2
`cor.score1:`	The estimated X-Y correlation in cell type 1
`cor.score2:`	The estimated X-Y correlation in cell type 2
`full.ll:`	The estimated log likelihood function from the full model
`reduced.ll:`	The estimated log likelihood function from the reduced model
`LRT:`	The likelihood ratio test statistics for iProMix
`LRT.pvalue:`	The chisq p-value for iProMix

library(iProMix)
set.seed(111)
y <- rnorm(100,10,1)
x <- rnorm(100,10,1)
pi <- runif(100)
iProMix.LRT_result <- iProMix.LRT(y = y, x = x, pi = pi, reduce1=c(2,1), reduce2=NULL)