[BioC] AB Rat Genome Survey Microarray
Iain Gallagher
iaingallagher at btopenworld.com
Thu Feb 16 20:03:34 CET 2012
For those that come afterwards.
Firstly I had the structure of the data file wrong - in the first two columns put the probeID and geneID, then create 3 columns for each array consisting of the Signal, S/N and Flags where the column names are in the form:
Signal GSM492839 S/N GSM492839 Flags GSM492839 Signal GSM517967 S/N GSM517967 Flags GSM517967... etc etc
Here each GSM.... is the name of an array (these were downloaded from GEO).
Secondly the 'Data file' as defined in the ABarray has to be comma delimited (at least on linux - I haven't tried other
platforms) not tab delimited (as it says it can be in the vignette). Trying to use tab delimited data results in an error about there being a mismatch in the number of Signal and S/N columns.
The package can now at least read the data and I'll see how I progress from here.
HTH
best
i
> sessionInfo()
R version 2.14.1 (2011-12-22)
Platform: x86_64-pc-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_GB.utf8 LC_NUMERIC=C
[3] LC_TIME=en_GB.utf8 LC_COLLATE=en_GB.utf8
[5] LC_MONETARY=en_GB.utf8 LC_MESSAGES=en_GB.utf8
[7] LC_PAPER=C LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_GB.utf8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] ABarray_1.22.0
loaded via a namespace (and not attached):
[1] Biobase_2.14.0 MASS_7.3-17 multtest_2.10.0 splines_2.14.1
[5] survival_2.36-12 tcltk_2.14.1 tools_2.14.1
----- Original Message -----
From: Iain Gallagher <iaingallagher at btopenworld.com>
To: "sunya at appliedbiosystems.com" <sunya at appliedbiosystems.com>
Cc: bioconductor <bioconductor at stat.math.ethz.ch>
Sent: Sunday, 12 February 2012, 9:19
Subject: [BioC] AB Rat Genome Survey Microarray
Hello List & Yongming Andrew Sun
I've been tasked with the analysis of some AB Rat Genome Survey Microarray data. The data I have are text files with the following columns:
Assay_Name Row Column Probe_ID Probe_Type Gene_ID X Y Assay_Normalized_Signal Signal CL_Sig CL_Raw SDEV CV S_N CL_Sig_Error CL_Raw_SDEV Flags Sample_Name
I was planning to use the ABarray package for this analysis but reading through the vignette I'm not following the examples. From the vignette:
"Data File. The rows of the file represent probes. The first column is probeID (or Probe Name), the second
column is geneID, next set of columns should contain Signal, S/N and FLAGS for all samples... If array name in the header of data file is present and
arrayName is defined in experiment design file, arrayName will be used to distinguish which column is for which
hybridization sample."
So, I take my array files and create a large text file representing all arrays (20) in long format.
> arrays <- read.table('allArrays.csv', header=T, sep='\t')
> head(arrays)
arrayName Probe_ID Gene_ID Signal S.N Flags
1 GSM517967 20693443 rCG63297 264.31 -0.88 1
2 GSM517967 20693548 rCG43972.1 3258.81 15.17 0
3 GSM517967 20693561 rCG44662 1200.68 -0.22 1
4 GSM517967 20693609 rCG59171 30762.77 40.18 0
5 GSM517967 20693611 rCG58434 838.39 2.13 0
6 GSM517967 20693655 rCG21764 5331.88 20.66 0
...
537135 GSM524823 22425691 rCG48508 535.35 1.27 0
537136 GSM524823 22425699 rCG40754 518.45 -0.26 1
537137 GSM524823 22425720 rCG31597 975.71 3.62 0
537138 GSM524823 22425784 rCG61959.1 331.31 -1.62 1
537139 GSM524823 22425829 rCG38412 27281.49 45.11 0
537140 GSM524823 22426049 rCG34097 2172.00 3.65 0
I create a design file as suggested in the vignette:
> des <- read.table('design.csv', header=T, sep='\t')
> head(des)
sampleName arrayName phenotype
1 LV_LRT1 GSM492839 LRT
2 LV_LRT2 GSM517967 LRT
3 LV_LRT3 GSM517968 LRT
4 LV_LRT4 GSM517969 LRT
5 LV_LRT5 GSM517970 LRT
6 LV_LRS1 GSM518311 LRS
Then, trying to create the expressionSet object:
> eset <- ABarray('allArrays.csv', 'design.csv', 'phenotype')
Using arrayName to match experiment with signal in file: allArrays.csv
Error in `[.data.frame`(pd, , group) : undefined columns selected
Hmm, could someone let me know where I'm going wrong?
Thanks & best
Iain
_______________________________________________
Bioconductor mailing list
Bioconductor at r-project.org
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list