[BioC] concentrations

Mohammad Esad-Djou shahrgol at web.de
Wed Jun 1 17:14:20 CEST 2005


Hallo all,
thank you very much for your suggestions.

> Next time, please learn to ask the question more clearly (with example
> of the desired output) to avoid the rest of guessing (e.g. what do you
> mean by "correctly finding numeric values").

Yes.  You are right. I try to reformulate my question. I would like to compute concentrations for different ProbeSet. 
I used SpikeIn example in 

Description of affy
Laurent Gautier, Rafael Irizarry, Leslie Cope, and Ben Bolstad
April 21, 2005

page 22:

> data(SpikeIn)
> pms <- pm(SpikeIn)
> mms <- mm(SpikeIn)
> par(mfrow = c(1, 2))
> concentrations <- matrix(as.numeric(sampleNames(SpikeIn)), 20,
+ 12, byrow = TRUE)
> matplot(concentrations, pms, log = "xy", main = "PM", ylim = c(30,
+ 20000))
> lines(concentrations[1, ], apply(pms, 2, mean), lwd = 3)
> matplot(concentrations, mms, log = "xy", main = "MM", ylim = c(30,
+ 20000))
> lines(concentrations[1, ], apply(mms, 2, mean), lwd = 3)


I wrote:
>library(affy)
>data.raw <- ReadAffy(filenames="./R/ME_cel/Expt1_R1.CEL", ....)
>pset1 <- probeset(data.raw,geneNames(data.raw)[1])

Test:
>pset1
$"1007_s_at"
ProbeSet object:
  id=1007_s_at
  pm= 16 probes x  4  chips

and:
> SpikeIn
ProbeSet object:
  id=AFFX-BioB-5_at
  pm= 20 probes x  12  chips

R program calls both objects ProbeSet, but I cannot go same way as above example: (e.g.)

> pms <- pm(pset1)
Error in pm(pset1) : No direct or inherited method for function "pm" for this call

> sampleNames(pset1)
Error in sampleNames(pset1) : No direct or inherited method for function "sampleNames" for this call

> sampleNames(SpikeIn)
 [1] "0.50"   "0.75"   "1.00"   "1.50"   "2.00"   "3.00"   "5.00"   "12.50" 
 [9] "25.00"  "50.00"  "75.00"  "150.00"


My main question:
1. How can I compute concentrations for different ProbeSet?

2. What didn't I consider for ProbeSet (SpikeIn and pset1)?

Thanks,
Mohammad Esad-Djou



More information about the Bioconductor mailing list