Skip to contents

listDeByTool: makes a data.frame with genes (lines) and all executed DE tool. Each tool that identify gene as DE, are marked by 1.

Usage

listDeByTool(consexpressionList, geneNames, deList)

Arguments

consexpressionList

A list generated by consexpressionR function

geneNames

name of all genes in experiment (list)

deList

List with names of all genes analysed in experiment

Value

data.frame with genes (lines) and DE tools (column), where 1 shows a DE indication by tool.

Examples

data(gse95077)
treats = c("BM", "JJ")
cons_result <- runExpression(numberReplics = 3, groupName = treats,
                              rDataFrameCount = gse95077,
                              sepCharacter = ",",
                              experimentName = "test_cons",
                              outDirPath = "." )
#> Using classic mode.
#> Getting annotation of the Homo Sapiens...
#> Using reference genome 38.
#> Calculating gene expression values...
#> RQ fit ......
#> SQN 
#> 
#>  
#>  ===== ERROR: KnowSeq execution is failed === 
#>  Error in mclustBIC(data = structure(c(-3.52446659876917, -2.42801002792396, : could not find function "mclustBIC"
#>  
#> Removing intercept from test coefficients
#> 
#>  ------------ limma executed!
#> [1] "Computing (M,D) values..."
#> [1] "Computing probability of differential expression..."
#> 
#>  ------------ NOISeq executed! 
#> Warning: `expect_is()` was deprecated in the 3rd edition.
#>  Use `expect_type()`, `expect_s3_class()`, or `expect_s4_class()` instead
#> Warning: `expect_is()` was deprecated in the 3rd edition.
#>  Use `expect_type()`, `expect_s3_class()`, or `expect_s4_class()` instead
#> Warning: `expect_is()` was deprecated in the 3rd edition.
#>  Use `expect_type()`, `expect_s3_class()`, or `expect_s4_class()` instead
#> Initial number of DE patterns = 2
#> Final number of DE patterns = 2
#> 
#>  ------------ EBSeq executed!
#> 
#>  
#>  ===== ERROR: in DESeq2 execution is failed === 
#>  === Contrast group is: [  ] 
#>  === Control group is: [  ] 
#>  ERROR!! This values are not defined in groupName: [ BM ]
#> 
#>  
#>  ===== ERROR: in DESeq2 execution is failed === 
#>  === Contrast group is: [  ] 
#>  === Control group is: [  ] 
#>  ERROR!! This values are not defined in groupName: [ JJ ]
#> 
#>  ------------ DESeq2 executed!
#> Estimating sequencing depths...
#> Resampling to get new data matrices...
#> perm= 1
#> perm= 2
#> perm= 3
#> perm= 4
#> perm= 5
#> perm= 6
#> perm= 7
#> perm= 8
#> perm= 9
#> perm= 10
#> perm= 11
#> perm= 12
#> perm= 13
#> perm= 14
#> perm= 15
#> perm= 16
#> perm= 17
#> perm= 18
#> perm= 19
#> perm= 20
#> perm= 21
#> perm= 22
#> perm= 23
#> perm= 24
#> perm= 25
#> perm= 26
#> perm= 27
#> perm= 28
#> perm= 29
#> perm= 30
#> perm= 31
#> perm= 32
#> perm= 33
#> perm= 34
#> perm= 35
#> perm= 36
#> perm= 37
#> perm= 38
#> perm= 39
#> perm= 40
#> perm= 41
#> perm= 42
#> perm= 43
#> perm= 44
#> perm= 45
#> perm= 46
#> perm= 47
#> perm= 48
#> perm= 49
#> perm= 50
#> perm= 51
#> perm= 52
#> perm= 53
#> perm= 54
#> perm= 55
#> perm= 56
#> perm= 57
#> perm= 58
#> perm= 59
#> perm= 60
#> perm= 61
#> perm= 62
#> perm= 63
#> perm= 64
#> perm= 65
#> perm= 66
#> perm= 67
#> perm= 68
#> perm= 69
#> perm= 70
#> perm= 71
#> perm= 72
#> perm= 73
#> perm= 74
#> perm= 75
#> perm= 76
#> perm= 77
#> perm= 78
#> perm= 79
#> perm= 80
#> perm= 81
#> perm= 82
#> perm= 83
#> perm= 84
#> perm= 85
#> perm= 86
#> perm= 87
#> perm= 88
#> perm= 89
#> perm= 90
#> perm= 91
#> perm= 92
#> perm= 93
#> perm= 94
#> perm= 95
#> perm= 96
#> perm= 97
#> perm= 98
#> perm= 99
#> perm= 100
#> Number of thresholds chosen (all possible thresholds) = 8
#> Getting all the cutoffs for the thresholds...
#> Getting number of false positives in the permutation...
#> 
#>  ------------ SAMSeq executed!
#> 
#>  ------------ DIFERENTIAL EXPRESSION ANALYSIS WAS COMPLETE!
expDef_result <- expressionDefinition(resultTool = cons_result, groups = treats)
deByTool <- listDeByTool(consexpressionList = cons_result,
                         geneNames = row.names(gse95077),
                         deList = expDef_result)