[BioC] error with pd.ragene1.0.st.v1 and oligo

[Ian B Stanaway]

Hi Benilton,

I am also trying to use this probe id database to assign gene symbols
to probe ids and it seems that most of the probe ids in the 'correct'
db are not shared with this chip.  Do you think there is a revision
problem with this array db also?

ragene10stprobeset.db

probeset.list <-featureNames(exp_data_rma)
xtmp <- ragene10stprobesetSYMBOL
mappedprobes <- mappedkeys(xtmp)
genes.symbols <- as.list(xtmp[mappedprobes])
xtmp <- paste(probeset.list, mappedprobes, genes.symbols,
genes.refseqNMids, collapse="\n")
write.table(t(xtmp), file="nm_id_lookup", row.names=F,
col.names=FALSE, sep="\t", quote=FALSE)

then when I look at the list the 2 for columns should be matching
probe ids but they don't match:

10700965 10703255 Zswim5 c("NM_001107963", "NP_001101433")
10700966 10703266 RGD1561339 c("XM_577732", "XP_577732")
10700967 10703267 RGD1564858 c("XM_233450", "XP_233450")
10700968 10703268 RGD1564858 c("XM_233450", "XP_233450")
10700969 10703303 RGD1560972 c("XM_577704", "XP_577704")
10700970 10703304 RGD1560808 c("XM_577734", "XP_577734")
10700971 10703311 Smoc2 c("NM_001106215", "NP_001099685")
10700972 10703312 Smoc2 c("NM_001106215", "NP_001099685")
10700973 10703313 Smoc2 c("NM_001106215", "NP_001099685")
10700974 10703314 Smoc2 c("NM_001106215", "NP_001099685")


#if I leave this in it gobbles up memory
gene.info.tmp <- data.frame(probeset.list=probeset.list,
genes.symbols=genes.symbols, stringsAsFactors=F)

thanks,

Ian


On Wed, Feb 12, 2014 at 11:46 AM, Ian B Stanaway <bard at uw.edu> wrote:
> Many Thanks!
>
> On Wed, Feb 12, 2014 at 11:39 AM, Benilton Carvalho
> <beniltoncarvalho at gmail.com> wrote:
>> Ian, I'm working on a fix, which should be available in the next couple of
>> days. Benilton
>>
>> On Feb 12, 2014 4:25 PM, "Ian B Stanaway" <bard at uw.edu> wrote:
>>>
>>> Hi,
>>>
>>> I should also note that if I remove the hist and boxplot sections of
>>> the script if proceeds and calculates all the normalizations, design,
>>> contrast.matrix and outputs topTables just fine.  It seems the naming
>>> in the graphical functions underlying code is off.
>>>
>>> best,
>>>
>>> Ian
>>>
>>> On Tue, Feb 11, 2014 at 1:59 PM, Ian [guest] <guest at bioconductor.org>
>>> wrote:
>>> >
>>> > Hi,
>>> >
>>> > I am analyzing some rat arrays. I get this error:
>>> >
>>> > Error in validObject(.Object) :
>>> >   invalid class "GeneFeatureSet" object: featureNames differ between
>>> > assayData and featureData
>>> > Calls: hist ... assayData<- -> assayData<- -> .init_NChannelSet ->
>>> > validObject
>>> >
>>> > I am pretty sure my code is close to correct as I put together this
>>> > script to analyze the pd.rat230.2 arrays for another project and it worked
>>> > great.
>>> >
>>> > Now I have switched to a project that used the pd.ragene.1.0.st.v1
>>> > instead and now I am getting this error and am in need of figuring it out
>>> > again.
>>> >
>>> > Any pointers would be great.
>>> >
>>> > Thanks,
>>> >
>>> > Ian
>>> >
>>> >  -- output of sessionInfo():
>>> >
>>> >> datadir <- "in_vitro_rat"
>>> >> setwd(datadir)
>>> >>
>>> >> celFiles <- list.celfiles()
>>> >> exp_data <- read.celfiles(celFiles, pkgname="pd.ragene.1.0.st.v1")
>>> > Platform design info loaded.
>>> > Reading in : 10.CEL
>>> > Reading in : 11.CEL
>>> > Reading in : 12.CEL
>>> > Reading in : 13.CEL
>>> > Reading in : 14.CEL
>>> > Reading in : 15.CEL
>>> > Reading in : 16.CEL
>>> > Reading in : 17.CEL
>>> > Reading in : 18.CEL
>>> > Reading in : 19.CEL
>>> > Reading in : 1.CEL
>>> > Reading in : 20.CEL
>>> > Reading in : 21.CEL
>>> > Reading in : 22.CEL
>>> > Reading in : 23.CEL
>>> > Reading in : 24.CEL
>>> > Reading in : 2.CEL
>>> > Reading in : 3.CEL
>>> > Reading in : 4.CEL
>>> > Reading in : 5.CEL
>>> > Reading in : 6.CEL
>>> > Reading in : 7.CEL
>>> > Reading in : 8.CEL
>>> > Reading in : 9.CEL
>>> >> png(file="exp.hist.png",width = 500, height = 600, type = "cairo")
>>> >> hist(exp_data)
>>> > Error in validObject(.Object) :
>>> >   invalid class "GeneFeatureSet" object: featureNames differ between
>>> > assayData and featureData
>>> > Calls: hist ... assayData<- -> assayData<- -> .init_NChannelSet ->
>>> > validObject
>>> > Execution halted
>>> >
>>> >
>>> > --
>>> > Sent via the guest posting facility at bioconductor.org.