arrayQualityMetrics with count table of RNA-Seq data - how to make ExpressionSetobject or ExpressionSetIllumina object?
0
0
Entering edit mode
Mabeyer • 0
@mabeyer-21106
Last seen 5.4 years ago

Hello, I have the output of HtSeq (=count tables as .txt files) as input for R and want to use arrayQualityMetrics. Is this somehow possible? I found this in the vignette: "If you are working on one colour arrays other than Affymetrix genechips, you can load your data into Bioconductor as an ExpressionSetobject, or if you work with Illumina data and the beadarray package, as an ExpressionSetIllumina object"

I have no idea if it is possible to receive either an ExpressionSetobject or an ExpressionSetIllumina object.

count_table <- read.table("all_counts_merged2.txt", header = TRUE)"
head(count_table)
#                  id   C     T1    T2   T3
# 1 rna-XM_011330573.3    0     0     0     0
# 2 rna-XM_011330574.2    0     0     0     0
# 3 rna-XM_011330575.3    0     0     0     0
# 4 rna-XM_011330577.3    0     0     0     0
# 5 rna-XM_011330578.3  133   374   382   246
# 6 rna-XM_011330579.3    0     0     0     0

I also generated a count Matrix:

countMatrix <- as.matrix(count_table[2:5]) 
rownames(countMatrix) <- count_table$id
head(countMatrix)
#                                 id   C     T1    T2   T3 
# rna-XM_011330573.3    0     0     0     0
# rna-XM_011330574.2    0     0     0     0
# rna-XM_011330575.3    0     0     0     0
# rna-XM_011330577.3    0     0     0     0
# rna-XM_011330578.3  133   374   382   246
# rna-XM_011330579.3    0     0     0     0

PS: I handle time course data (at the moment only one replicate) C=0 h, T1=24h, T2=48h, T3=120h (h=hours)

sessionInfo()
R version 3.6.0 (2019-04-26)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 10 x64 (build 17134)

Matrix products: default

Random number generation:
 RNG:     Mersenne-Twister 
 Normal:  Inversion 
 Sample:  Rounding 

locale:
[1] LC_COLLATE=German_Germany.1252  LC_CTYPE=German_Germany.1252    LC_MONETARY=German_Germany.1252 LC_NUMERIC=C                    LC_TIME=German_Germany.1252    

attached base packages:
[1] parallel  stats4    stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
 [1] genefilter_1.66.0           gridExtra_2.3               RColorBrewer_1.1-2          pheatmap_1.0.12             limma_3.40.2                gplots_3.0.1.1             
 [7] arrayQualityMetrics_3.40.0  ggplot2_3.2.0               DESeq2_1.24.0               SummarizedExperiment_1.14.0 DelayedArray_0.10.0         BiocParallel_1.17.18       
[13] matrixStats_0.54.0          Biobase_2.44.0              GenomicRanges_1.36.0        GenomeInfoDb_1.20.0         IRanges_2.18.1              S4Vectors_0.22.0           
[19] BiocGenerics_0.30.0        

loaded via a namespace (and not attached):
Error in x[["Version"]] : subscript out of bounds
In addition: Warning messages:
1: In FUN(X[[i]], ...) : DESCRIPTION von Paket 'rlang' fehlt oder zerstört
2: In FUN(X[[i]], ...) : DESCRIPTION von Paket 'dplyr' fehlt oder zerstört
3: In FUN(X[[i]], ...) :
  DESCRIPTION von Paket 'beadarray' fehlt oder zerstört
arrayQualityMetrics RNA-Seq ExpressionSetobject ExpressionSetIllumina • 847 views
ADD COMMENT

Login before adding your answer.

Traffic: 706 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6