multiNFSSEMiPALM2 {fssemR}

multiNFSSEMiPALM2

Description

Implementing NFSSEM algorithm for network inference. If Xs is identify for different conditions, multiNFSSEMiPALM will be use, otherwise, please use multiNFSSEMiPALM2 for general cases

Usage

multiNFSSEMiPALM2(
  Xs,
  Ys,
  Bs,
  Fs,
  Sk,
  sigma2,
  lambda,
  rho,
  Wl,
  Wf,
  p,
  maxit = 100,
  inert = inert_opt("linear"),
  threshold = 1e-06,
  verbose = TRUE,
  sparse = TRUE,
  trans = FALSE,
  B2norm = NULL,
  strict = FALSE
)

Arguments

Xs

eQTL matrices

Ys

Gene expression matrices

Bs

initialized GRN-matrices

Fs

initialized eQTL effect matrices

Sk

eQTL index of genes

sigma2

initialized noise variance from ridge regression

lambda

Hyperparameter of lasso term in NFSSEM

rho

Hyperparameter of fused-lasso term in NFSSEM

Wl

weight matrices for adaptive lasso terms

Wf

weight matrix for columnwise l2 norm adaptive group lasso

p

number of genes

maxit

maximum iteration number. Default 100

inert

inertial function for iPALM. Default as k-1/k+2

threshold

convergence threshold. Default 1e-6

verbose

Default TRUE

sparse

Sparse Matrix or not

trans

Fs matrix is transposed to k x p or not. If Fs from ridge regression, trans = TRUE, else, trans = FALSE

B2norm

B2norm matrices generated from ridge regression. Default NULL.

strict

Converge strictly or not. Default False

Value

fit List of NFSSEM model

Bs

coefficient matrices of gene regulatory networks

Fs

coefficient matrices of eQTL-gene effect

mu

Bias vector

sigma2

estimate of covariance in SEM