ISoLDE.Rcheck/tests_i386/testsIsolde.Rout
R version 3.5.1 Patched (2018-07-24 r75005) -- "Feather Spray"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: i386-w64-mingw32/i386 (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(ISoLDE)
>
> targetfile <- system.file("extdata", "target_file.txt", package = "ISoLDE")
> rawfile <- system.file("extdata", "rawASRcounts_file.txt", package = "ISoLDE")
> normfile <- system.file("extdata", "normASRcounts_file.txt", package = "ISoLDE")
>
> rawASRcounts <- readRawInput(raw_file = rawfile, del = "tab", rownames = TRUE, colnames = TRUE)
Reading raw input file...
Done
> normASRcounts <- readNormInput(norm_file = normfile, del = "tab", rownames = TRUE, colnames = TRUE, dec = ".")
Reading normalized input file...
Done
> target <- readTarget(target_file = targetfile, asr_counts = rawASRcounts, del = "tab")
Reading target file...
Done
>
> res_filterT <- filterT(rawASRcounts = rawASRcounts, normASRcounts = normASRcounts, target = target, bias="parental")
Filtering your data...
Filtering threshold is 3
Done
> filteredASRcounts <- res_filterT$filteredASRcounts
> resiso<-isolde_test(bias = "parental", asr_counts = filteredASRcounts, target = target, graph = TRUE, text = FALSE, ext = "png", nboot = 3000,pcore=1)
Warning: pcore is a percentage between 0 and 100, be sure you want to use 1% of your available cores.
ISoLDE is searching for allele-specific gene expression due to parental biases...
Done
According to the satisfying number of replicates the full method can be applied. See the vignette for more details
Testing...
This step can last for a few minutes. Please be patient :)
Note: the full method uses a bootstrap step which means some results might change from one test to another.
Open MP enabled
Have Cores : 20
UseCores : 1
nu0= 0.200000
Open MP enabled
Have Cores : 20
UseCores : 1
lambda0 : 1.600000
Open MP enabled
Have Cores : 20
UseCores : 1
nu0= 0.100000
Open MP enabled
Have Cores : 20
UseCores : 1
lambda0 : 1.600000
Open MP enabled
Have Cores : 20
UseCores : 1
nu0= 0.100000
Open MP enabled
Have Cores : 20
UseCores : 1
lambda0 : 1.600000
Done
49 allele specific expressed (ASE) genes found and written in file ISoLDE_result_ASE_10-17-2018_02-54-29.tsv.
6092 biallelic (BA) genes found and written in file ISoLDE_result_BA_10-17-2018_02-54-29.tsv.
36 Undetermined (UN) genes found and written in file ISoLDE_result_UN_10-17-2018_02-54-29.tsv.
Among these undetermined (UN) genes, 20 genes are coherent.
These genes correspond to the ones with the flag column set on "FLAG_consistency" of "FLAG_significance" in the undetermined genes file. See the vignette for more information on flags.
0 genes have been filtered due to few reads.
> res_filterT <- filterT(rawASRcounts = rawASRcounts, normASRcounts = normASRcounts, target = target, bias="strain")
Filtering your data...
Filtering threshold is 3
Done
> filteredASRcounts <- res_filterT$filteredASRcounts
> resiso<-isolde_test(bias = "strain", asr_counts = filteredASRcounts, target = target, graph = TRUE, text = FALSE, ext = "png", nboot = 3000, method = "threshold")
ISoLDE is searching for allele-specific gene expression due to strain biases...
The predefined threshold method is about to proceed, but we strongly recommand the bootstrap method because ISoLDE detected more than 2 replicates in each cross.
Testing...
Done
308 allele specific expressed (ASE) genes found and written in file ISoLDE_result_ASE_10-17-2018_02-54-31.tsv.
4622 biallelic (BA) genes found and written in file ISoLDE_result_BA_10-17-2018_02-54-31.tsv.
1330 Undetermined (UN) genes found and written in file ISoLDE_result_UN_10-17-2018_02-54-31.tsv.
Among these undetermined (UN) genes, 913 genes are coherent.
These genes correspond to the ones with the flag column set on "FLAG_consistency" of "FLAG_significance" in the undetermined genes file. See the vignette for more information on flags.
0 genes have been filtered due to few reads.
>
> proc.time()
user system elapsed
158.25 0.60 158.85
|
ISoLDE.Rcheck/tests_x64/testsIsolde.Rout
R version 3.5.1 Patched (2018-07-24 r75005) -- "Feather Spray"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64 (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(ISoLDE)
>
> targetfile <- system.file("extdata", "target_file.txt", package = "ISoLDE")
> rawfile <- system.file("extdata", "rawASRcounts_file.txt", package = "ISoLDE")
> normfile <- system.file("extdata", "normASRcounts_file.txt", package = "ISoLDE")
>
> rawASRcounts <- readRawInput(raw_file = rawfile, del = "tab", rownames = TRUE, colnames = TRUE)
Reading raw input file...
Done
> normASRcounts <- readNormInput(norm_file = normfile, del = "tab", rownames = TRUE, colnames = TRUE, dec = ".")
Reading normalized input file...
Done
> target <- readTarget(target_file = targetfile, asr_counts = rawASRcounts, del = "tab")
Reading target file...
Done
>
> res_filterT <- filterT(rawASRcounts = rawASRcounts, normASRcounts = normASRcounts, target = target, bias="parental")
Filtering your data...
Filtering threshold is 3
Done
> filteredASRcounts <- res_filterT$filteredASRcounts
> resiso<-isolde_test(bias = "parental", asr_counts = filteredASRcounts, target = target, graph = TRUE, text = FALSE, ext = "png", nboot = 3000,pcore=1)
Warning: pcore is a percentage between 0 and 100, be sure you want to use 1% of your available cores.
ISoLDE is searching for allele-specific gene expression due to parental biases...
Done
According to the satisfying number of replicates the full method can be applied. See the vignette for more details
Testing...
This step can last for a few minutes. Please be patient :)
Note: the full method uses a bootstrap step which means some results might change from one test to another.
Open MP enabled
Have Cores : 20
UseCores : 1
nu0= 0.100000
Open MP enabled
Have Cores : 20
UseCores : 1
lambda0 : 1.600000
Open MP enabled
Have Cores : 20
UseCores : 1
nu0= 0.100000
Open MP enabled
Have Cores : 20
UseCores : 1
lambda0 : 1.600000
Open MP enabled
Have Cores : 20
UseCores : 1
nu0= 0.100000
Open MP enabled
Have Cores : 20
UseCores : 1
lambda0 : 1.600000
Done
47 allele specific expressed (ASE) genes found and written in file ISoLDE_result_ASE_10-17-2018_02-56-00.tsv.
6095 biallelic (BA) genes found and written in file ISoLDE_result_BA_10-17-2018_02-56-00.tsv.
35 Undetermined (UN) genes found and written in file ISoLDE_result_UN_10-17-2018_02-56-00.tsv.
Among these undetermined (UN) genes, 15 genes are coherent.
These genes correspond to the ones with the flag column set on "FLAG_consistency" of "FLAG_significance" in the undetermined genes file. See the vignette for more information on flags.
0 genes have been filtered due to few reads.
> res_filterT <- filterT(rawASRcounts = rawASRcounts, normASRcounts = normASRcounts, target = target, bias="strain")
Filtering your data...
Filtering threshold is 3
Done
> filteredASRcounts <- res_filterT$filteredASRcounts
> resiso<-isolde_test(bias = "strain", asr_counts = filteredASRcounts, target = target, graph = TRUE, text = FALSE, ext = "png", nboot = 3000, method = "threshold")
ISoLDE is searching for allele-specific gene expression due to strain biases...
The predefined threshold method is about to proceed, but we strongly recommand the bootstrap method because ISoLDE detected more than 2 replicates in each cross.
Testing...
Done
308 allele specific expressed (ASE) genes found and written in file ISoLDE_result_ASE_10-17-2018_02-56-03.tsv.
4622 biallelic (BA) genes found and written in file ISoLDE_result_BA_10-17-2018_02-56-03.tsv.
1330 Undetermined (UN) genes found and written in file ISoLDE_result_UN_10-17-2018_02-56-03.tsv.
Among these undetermined (UN) genes, 913 genes are coherent.
These genes correspond to the ones with the flag column set on "FLAG_consistency" of "FLAG_significance" in the undetermined genes file. See the vignette for more information on flags.
0 genes have been filtered due to few reads.
>
> proc.time()
user system elapsed
91.00 0.68 91.71
|