CallT2Dpeak_qvalue: CallT2Dpeak_qvalue

In chenweng1991/RePACT: What the Package Does (One Line, Title Case)

CallT2Dpeak_qvalue

Description

This function is to run regression (linear/logistic) based on the number of PCs, and characteristics of samples.

Usage

CallT2Dpeak_qvalue(
  cellvsPeak.m.aggr = ATAConCCA.betaT2D.tjct.3nd.10bin.ob$cellvsPeak.m.aggr,
  depths = ATAConCCA.betaT2D.tjct.3nd.10bin.ob$depths,
  index = ATAConCCA.betaT2D.tjct.3nd.10bin.ob$index,
  qcut = 0.1,
  slopecut1 = 0.5,
  slopecut2 = -0.5,
  doscale = T
)

Arguments

cellvsPeak.m.aggr,	gene~traj_bin counts matrix
depths	total fragments of cells within each traj_bin
index	avg pseudo-index of each bin
qcut	qvalue cutoff to measure the significance of glm(expr~index)
slopecut1	UP slope cutoff, slope from glm(expr~index), UP genes
slopecut2	DN slope cutoff, slope from glm(expr~index), DN genes
doscale	if build glm model based on scaled_expression and scaled pseudo-inddex

Value

The function return a list: "PCvariance", "PCanfpheno", "object.raw.withinfo", "model", "reg.plot.2d", "model.para"

Examples

betaT2D.diffGene.20bin.PCA <- CallT2Dpeak_qvalue(beta.RNA.PCA.20bin.ob$cellvsPeak.m.aggr, depths=beta.RNA.PCA.20bin.ob$depths, index=beta.RNA.PCA.20bin.ob$index, qcut=0.2,slopecut1=0.3,slopecut2=-0.3,doscale=T)

chenweng1991/RePACT documentation built on Aug. 28, 2023, 6:28 p.m.