[BioC] Subsetting data from RGList (Again)

elliott harrison e.harrison at epistem.co.uk
Tue Oct 30 16:17:08 CET 2007


Hi All, 
Sorry to trouble everyone  but I didn't get my subsetting of an RGList
sorted last week.
I've tried loads of different combinations of things that were suggested
but not quite there yet.

MY Data is an RGList object here is a bit
> R[1,1]
An object of class "RGList"
$G
     C:\\Program Files\\Agilent\\GeneSpring
GX\\data\\CRX\\980\\test01_251485020980_S01_GE2-v5_91_0806_2_1_1
[1,]
15662.69
 
$Gb
     C:\\Program Files\\Agilent\\GeneSpring
GX\\data\\CRX\\980\\test01_251485020980_S01_GE2-v5_91_0806_2_1_1
[1,]
73.9323
 
$R
     C:\\Program Files\\Agilent\\GeneSpring
GX\\data\\CRX\\980\\test01_251485020980_S01_GE2-v5_91_0806_2_1_1
[1,]
15662.69
 
$Rb
     C:\\Program Files\\Agilent\\GeneSpring
GX\\data\\CRX\\980\\test01_251485020980_S01_GE2-v5_91_0806_2_1_1
[1,]
73.9323
 
$weights
     C:\\Program Files\\Agilent\\GeneSpring
GX\\data\\CRX\\980\\test01_251485020980_S01_GE2-v5_91_0806_2_1_1
[1,]
0
 
$targets
  SlideNumber   Name
FileName Treatment MasterTreatment
1        9801 Pre1_1 C:\\Program Files\\Agilent\\GeneSpring
GX\\data\\CRX\\980\\test01_251485020980_S01_GE2-v5_91_0806_2_1_1.txt
Pre1             Pre
       Cy3 Cy5       Date
1 WildType Ref 11/09/2007
 
$genes
  Row Col FeatureNum ProbeUID ControlType       ProbeName
GeneName  SystematicName   Status
1   1   1          1        0           1 GE_BrightCorner
GE_BrightCorner GE_BrightCorner Control6
 
$source
[1] "generic"
 
 
My list of probes to select data with is in this object.
> a
           Gene.Probe
Probe1   A_24_P329635
Probe2   A_23_P104098
Probe3   A_23_P419038
Probe4   A_24_P124550
Probe5    A_23_P99044
 
 
 
 
I've tried
A=rownames(a)
RA = R[A,]
 
But don't get any genes
>$genes
[1] Row            Col            FeatureNum     ProbeUID
ControlType    ProbeName      GeneName       SystematicName Status

<0 rows> (or 0-length row.names)
 
I've tried creating a vector using match a load of different ways but it
is always empty
>vec = match(A, R)
>vec = match(A, R$genes$ProbeName)
>vec = match(A, R$genes)
>vec = match(probelist600, R)
>vec = match(probelist600, R$genes)
>vec = match(probelist600, R$genes$ProbeName)
 
 
This is part what ends up in object A
> A
  [1] "Probe1"   "Probe2"   "Probe3"   "Probe4"   "Probe5"   "Probe6"
"Probe7"   "Probe8"   "Probe9"   "Probe10"  "Probe11"  "Probe12"
"Probe13" 

I expected to see my probe names in there.
I'm clearly not getting something. 
Can anyone see where I'm going wrong?
Many Thanks
Elliott


This message has been scanned for viruses by BlackSpider...{{dropped:3}}



More information about the Bioconductor mailing list