[BioC] Kooky probe names (followup)

Crispin Miller CMiller at picr.man.ac.uk
Wed Jun 11 15:01:45 MEST 2003


Hi,
More information about the probe-names problems I'm having with probe names becoming X12324.at instead of 1234_at


Firstly, this is the version of R I'm using:

>R.version
         _
platform i686-pc-linux-gnu
arch     i686
os       linux-gnu
system   i686, linux-gnu
status   Beta
major    1
minor    7.1
year     2003
month    06
day      08
language R


I've managed to reproduce the wierdness with the following minimal script:
I start off with two dataframes:
class(first);
class(second);
> class(first);
[1] "data.frame"
> class(second);
[1] "data.frame"

containing some numbers calculated for each probeset (these are both generated using identical functions, just comparing different groups of chips):

first[1:2,]
second[1:2,]

> first[1:2,]
          X1        X2        X3       X4
1007_s_at  0 0.1218871 0.2664602 1.202853
1053_at    0 0.2757304 0.4460427 1.362298
> second[1:2,]
          X1        X2          X3        X4
1007_s_at  0 0.6255407 -0.08225843 -1.058674
1053_at    0 0.6037322 -0.19470686 -1.144492

I then find entries where X2<0.05 in each data.frame:
> first[first$X2<0.5,][1:2,]
          X1        X2        X3       X4
1007_s_at  0 0.1218871 0.2664602 1.202853
1053_at    0 0.2757304 0.4460427 1.362298
> second[second$X2<0.5,][1:2,]
           X1         X2          X3        X4
X121.at     0 0.09278053 -0.19887301 -1.147801
X1255.g.at  0 0.41586962  0.09071178  1.064895

Weird!

Now, I play a bit:
> second[1:2,][1:2,]
          X1        X2          X3        X4
1007_s_at  0 0.6255407 -0.08225843 -1.058674
1053_at    0 0.6037322 -0.19470686 -1.144492
> second[1:2,][1:3,]
           X1        X2          X3        X4
X1007.s.at  0 0.6255407 -0.08225843 -1.058674
X1053.at    0 0.6037322 -0.19470686 -1.144492
NA         NA        NA          NA        NA
>

Does this give a clue as to what's going on?

Cheers,
Crispin
 
--------------------------------------------------------

 
This email is confidential and intended solely for the use of th... {{dropped}}



More information about the Bioconductor mailing list