[BioC] lumiR returns ExpressionSet object instead of LumiBatch

Liu X [guest] guest at bioconductor.org
Fri Jun 28 04:30:24 CEST 2013


Hello,

It is my first time trying to use lumi package to analyze beadchip data by following the instruction in vignettes (dated 2013/4/4). However I run into problems. Trying to use lumiR to read in a beadstudio export file with following commands:

library(lumi)
library(illumi
sampfile <- 'path/to/sample_probe_profile.txt'
lumi.sample<-lumiR(sampfile, lib.mapping="lumiMouseIDMapping", annotationColumn = c('PROBE_ID', 'SYMBOL', 'ILMN_GENE', 'CHROMOSOME', 'DEFINITION', 'SYNONYMS') )

It runs without error, but class(lumi.sample) tells me that it is an ExpressionSet object not LumiBatch. And of course I cannt apply lumi methods to analyze it.

What wrong with my steps? 

Thanks!

Here is the first 15 lines of the probe_profile file:
PROBE_ID        SYMBOL  NT1.AVG_Signal  NT1.Detection Pval      NT2.AVG_Signal  NT2.Detection Pval      NT3.AVG_Signal  NT3.Detection Pval      R1.AVG_Signal   R1.Detection PvalR2.AVG_Signal   R2.Detection Pval       R3.AVG_Signal   R3.Detection Pval       M1_2.AVG_Signal M1_2.Detection Pval     HG1_2.AVG_Signal        HG1_2.Detection Pval    M1.AVG_Signal    M1.Detection Pval       M2.AVG_Signal   M2.Detection Pval       M3.AVG_Signal   M3.Detection Pval       HG1.AVG_Signal  HG1.Detection Pval      HG2.AVG_Signal  HG2.Detection Pval       HG3.AVG_Signal  HG3.Detection Pval      NT1_2.AVG_Signal        NT1_2.Detection Pval    R1_2.AVG_Signal R1_2.Detection Pval     SEARCH_KEY      ILMN_GENE       CHROMOSOME       DEFINITION      SYNONYMS
ILMN_2896528    0610006I08Rik   2607.672        0.00    2527.176        0.00    2774.688        0.00    1997.495        0.00    2124.155        0.00    2544.899        0.00    3062.838 0.00    3381.947        0.00    2546.513        0.00    3017.973        0.00    2987.024        0.00    3376.653        0.00    2492.156        0.00    3577.373        0.00     2659.114        0.00    2353.836        0.00    ILMN_220369     0610006I08RIK   19      Mus musculus RIKEN cDNA 0610006I08 gene (0610006I08Rik), mRNA.
ILMN_2721178    0610006I08Rik   335.6214        0.00    309.3783        0.00    311.5328        0.00    255.2664        0.00    253.7414        0.00    304.8765        0.00    363.3045 0.00    416.1963        0.00    285.1985        0.00    417.5847        0.00    356.8044        0.00    371.7544        0.00    308.4683        0.00    533.7448        0.00     304.7859        0.00    259.3861        0.00    ILMN_220369     0610006I08RIK   19      Mus musculus RIKEN cDNA 0610006I08 gene (0610006I08Rik), mRNA.
ILMN_3033922    0610007C21Rik   1435.676        0.00    1407.494        0.00    1813.653        0.00    1138.474        0.00    1199.863        0.00    1369.412        0.00    1186.747 0.00    1436.43 0.00    1138.299        0.00    1344.861        0.00    1320.021        0.00    1378.295        0.00    1026.762        0.00    1909.174        0.00    1361.688 0.00    1232.837        0.00    ILMN_244468     0610007C21RIK   5       Mus musculus RIKEN cDNA 0610007C21 gene (0610007C21Rik), transcript variant 1, mRNA.    AI316792; HSPC013; Apr3; p18
ILMN_3092673    0610007C21Rik   11323.47        0.00    9251.417        0.00    11511.11        0.00    9007.652        0.00    8195.785        0.00    9388.674        0.00    9282.038 0.00    9625.286        0.00    8375.805        0.00    9982.474        0.00    9821.961        0.00    9076.46 0.00    7963.4  0.00    9219.595        0.00    10254.2 0.00     9168.383        0.00    ILMN_229637     0610007C21RIK   5       Mus musculus RIKEN cDNA 0610007C21 gene (0610007C21Rik), transcript variant 2, mRNA.    AI316792; HSPC013; Apr3; p18
ILMN_2816356    0610007P08Rik   25.62827        0.09398496      45.53978        0.02380952      0.5504513       0.3784461       23.20286        0.0989975       7.337133        0.2568922        16.80947        0.1566416       18.40724        0.1478697       31.25397        0.0764411       0.1861096       0.3759398       10.07403        0.2255639       16.067   0.1503759       21.37516        0.1027569       6.433805        0.2393484       24.07508        0.08646616      24.31498        0.08646616      26.45578        0.08270676       ILMN_196322     0610007P08RIK   13      Mus musculus RIKEN cDNA 0610007P08 gene (0610007P08Rik), transcript variant 2, mRNA.    RAD26L; 1700019D06Rik
ILMN_2808939    0610007P14Rik   712.1257        0.00    624.3039        0.00    779.4008        0.00    695.1705        0.00    671.4506        0.00    787.93  0.00    756.97860.00     476.1805        0.00    574.3113        0.00    746.5471        0.00    743.4769        0.00    422.7865        0.00    376.804 0.00    620.1172        0.00    607.81810.00     753.4628        0.00    ILMN_239772     0610007P14RIK   12      Mus musculus RIKEN cDNA 0610007P14 gene (0610007P14Rik), mRNA.  AU019315; C77855; 1190004E09Rik; ORF11
ILMN_2634564    0610007P22Rik   595.3325        0.00    451.3885        0.00    580.5104        0.00    371.2965        0.00    492.7807        0.00    479.4933        0.00    549.2484 0.00    494.6838        0.00    401.5907        0.00    437.3422        0.00    539.7564        0.00    471.3762        0.00    331.5264        0.00    516.2159        0.00     535.5032        0.00    477.19  0.00    ILMN_213192     0610007P22RIK   17      Mus musculus RIKEN cDNA 0610007P22 gene (0610007P22Rik), mRNA.  AL022759; MGC113771; 1110014B11Rik
ILMN_2737647    0610008C08Rik   -0.5213413      0.4085213       7.743881        0.254386        4.068563        0.3157895       6.312949        0.2593985       7.170675        0.2606516        10.64135        0.2155388       25.66834        0.0877193       7.385303        0.2593985       -5.414614       0.5413534       0.8812458       0.3822055       14.31114 0.1716792       -3.741079       0.5012531       13.13283        0.1390978       -2.912084       0.4724311       -13.54885       0.8233083       -6.836047       0.6052632ILMN_221600     0610008C08RIK   X       Mus musculus RIKEN cDNA 0610008C08 gene (0610008C08Rik), mRNA.  RP23-272D10.2; MGC130105; MGC130106; 1110019O03Rik
ILMN_2734484    0610009B22Rik   1660.084        0.00    1600.733        0.00    1853.014        0.00    1533.854        0.00    1641.063        0.00    1688.237        0.00    1758.677 0.00    1718.719        0.00    1528.311        0.00    1867.37 0.00    1768.41 0.00    1607.629        0.00    1211.871        0.00    1692.573        0.00    1825.8470.00     1744.845        0.00    ILMN_221361     0610009B22RIK   11      Mus musculus RIKEN cDNA 0610009B22 gene (0610009B22Rik), mRNA.  RP23-79E13.7
ILMN_2952292    0610009D07Rik   161.6432        0.00    150.4215        0.00    216.1236        0.00    169.6933        0.00    141.3559        0.00    208.3837        0.00    183.1769 0.00    113.102 0.00    145.3057        0.00    174.5983        0.00    196.8766        0.00    98.69936        0.00    85.1308 0.00    135.3256        0.00    207.43860.00     199.6864        0.00    ILMN_218698     0610009D07RIK   12      Mus musculus RIKEN cDNA 0610009D07 gene (0610009D07Rik), mRNA.  6030419K15Rik; AV001342; Sf3b14
ILMN_2699078    0610009D07Rik   -13.64666       0.7719298       -17.76805       0.9285714       -13.25819       0.7243108       -7.510522       0.6127819       -21.9306        0.9699248        -17.45373       0.8696742       -19.45185       0.9360902       -9.678173       0.6303258       -5.608812       0.5501253       -10.17932       0.6729323       -17.24747        0.8746867       -16.90044       0.9172933       -8.415766       0.7230576       -11.58076       0.7481203       -20.62218       0.9761904       -9.215716       0.6654136        ILMN_218698     0610009D07RIK   12      Mus musculus RIKEN cDNA 0610009D07 gene (0610009D07Rik), mRNA.  6030419K15Rik; AV001342; Sf3b14
ILMN_1213681    0610009O20Rik   177.286 0.00    146.7101        0.00    148.6533        0.00    136.2325        0.00    194.6553        0.00    193.7644        0.00    165.37730.00     167.7943        0.00    163.8992        0.00    221.8292        0.00    197.5723        0.00    203.4886        0.00    139.8207        0.00    300.4458        0.00    163.7285 0.00    179.5151        0.00    ILMN_220613     0610009O20RIK   18      Mus musculus RIKEN cDNA 0610009O20 gene (0610009O20Rik), mRNA.  2700004E22Rik
ILMN_2735413    0610010D20Rik   7.146307        0.2706767       7.9561  0.2468672       24.70525        0.1102757       5.685019        0.2719298       10.20794        0.2155388-4.846704       0.5213032       -3.845815       0.4874687       -6.052074       0.5238096       -6.612185       0.5827067       1.393941        0.3759398       -11.49365       0.735589 1.478938        0.3558897       5.659925        0.2506266       9.050623        0.226817        -2.104225       0.4398496       -15.24391       0.8558897       ILMN_221425      0610010D20RIK   19      Mus musculus RIKEN cDNA 0610010D20 gene (0610010D20Rik), mRNA.
ILMN_2735415    0610010D20Rik   -11.26077       0.7067669       -0.4047703      0.4035088       6.900579        0.2807018       0.7097692       0.3834586       -7.749537       0.6215539        -16.24752       0.8446115       17.96398        0.1503759       5.702412        0.2819549       5.526793        0.2869674       11.87896        0.2017544       1.067933 0.3734336       -10.01922       0.7205514       -0.05420115     0.4110276       -3.374624       0.4837093       12.72869        0.1691729       -8.716571       0.6453634ILMN_221425     0610010D20RIK   19      Mus musculus RIKEN cDNA 0610010D20 gene (0610010D20Rik), mRNA.




 -- output of sessionInfo(): 

R version 3.0.1 (2013-05-16)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

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

other attached packages:
 [1] TxDb.Hsapiens.UCSC.hg19.knownGene_2.9.2  GenomicFeatures_1.12.2                  
 [3] lumiMouseIDMapping_1.10.0                lumi_2.12.0                             
 [5] BSgenome.Scerevisiae.UCSC.sacCer1_1.3.17 sagenhaft_1.30.0                        
 [7] SparseM_0.99                             edgeR_3.2.3                             
 [9] limma_3.16.5                             hexbin_1.26.2                           
[11] BSgenome_1.28.0                          ShortRead_1.18.0                        
[13] latticeExtra_0.6-24                      RColorBrewer_1.0-5                      
[15] lattice_0.20-15                          Rsamtools_1.12.3                        
[17] Biostrings_2.28.0                        rtracklayer_1.20.2                      
[19] GenomicRanges_1.12.4                     IRanges_1.18.1                          
[21] gaschYHS_1.0.11                          org.Hs.eg.db_2.9.0                      
[23] org.Sc.sgd.db_2.9.1                      BiocInstaller_1.10.2                    
[25] GO.db_2.9.0                              RSQLite_0.11.4                          
[27] DBI_0.2-7                                AnnotationDbi_1.22.6                    
[29] Biobase_2.20.0                           BiocGenerics_0.6.0                      

loaded via a namespace (and not attached):
 [1] affy_1.38.1           affyio_1.28.0         annotate_1.38.0       beanplot_1.1         
 [5] biomaRt_2.16.0        bitops_1.0-5          colorspace_1.2-2      hwriter_1.3          
 [9] illuminaio_0.2.0      KernSmooth_2.23-10    MASS_7.3-26           Matrix_1.0-12        
[13] matrixStats_0.8.1     mclust_4.1            methylumi_2.6.1       mgcv_1.7-24          
[17] minfi_1.6.0           multtest_2.16.0       nleqslv_2.0           nlme_3.1-109         
[21] nor1mix_1.1-4         preprocessCore_1.22.0 R.methodsS3_1.4.2     RCurl_1.95-4.1       
[25] reshape_0.8.4         siggenes_1.34.0       splines_3.0.1         stats4_3.0.1         
[29] survival_2.37-4       tools_3.0.1           XML_3.95-0.2          xtable_1.7-1         
[33] zlibbioc_1.6.0       

--
Sent via the guest posting facility at bioconductor.org.



More information about the Bioconductor mailing list