compareSteadyNoNascent: Identify post-transcriptionally regulated genes from an...

Description Usage Arguments Examples

Description

This function compare exons and introns expression level matrices, from two up to an arbitrary number of samples, in order to identify genes which are oddly regluated, compared to an expected standard behaviour, from the post transcriptional point of view.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
compareSteadyNoNascent(
  inspectIds,
  expressionThreshold = 0.25,
  log2FCThreshold = 2,
  trivialAngle = NaN,
  returnNormScores = FALSE,
  referenceCondition = "median"
)

## S4 method for signature 'INSPEcT_steadyNoNascent'
compareSteadyNoNascent(
  inspectIds,
  expressionThreshold = 0.25,
  log2FCThreshold = 2,
  trivialAngle = NaN,
  returnNormScores = FALSE,
  referenceCondition = "median"
)

Arguments

inspectIds

An object of class INSPEcT_steadyNoNascent

expressionThreshold

A parameter which sets how many log2 fold changes of distance from the median behaviour are imputable to noise.

log2FCThreshold

A parameter which sets the log2 fold change distance from the median behaviour that is imputable to noise.

trivialAngle

A numeric between 0 and 90 to define the standard behavior, if NaN (default) it is computed internally from the data.

returnNormScores

A logical, if TRUE returned the deviations from the standard behavior normalized by the sd.

referenceCondition

The label of the condition to use as reference, if NaN (default) the medians are used.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
data('allcounts', package='INSPEcT')
data('featureWidths', package='INSPEcT')
data('libsizes', package='INSPEcT')

nascentCounts<-allcounts$nascent
matureCounts<-allcounts$mature
conditions<-letters[1:11]
expDes<-rep(conditions,3)

matExp_DESeq2<-quantifyExpressionsFromTrCounts(
      allcounts=matureCounts
      ,libsize=totalLS
      ,exonsWidths=exWdths
      ,intronsWidths=intWdths
      ,experimentalDesign=expDes)

matureInspObj <- newINSPEcT(tpts=conditions,matureExpressions=matExp_DESeq2)

matureInspObj<-compareSteadyNoNascent(inspectIds=matureInspObj
								 ,expressionThreshold=0.25
							 ,log2FCThreshold=.5)
regGenes <- PTreg(matureInspObj)
head(regGenes)
table(regGenes)

INSPEcT documentation built on Nov. 8, 2020, 6:49 p.m.