Hi Setsuko,
I haven't looked at this in a long time, and I can't seem to find the CDF file
that I used locally on my computer. The information regarding the antigenomic
probesets should be readily available in the files you get from Affymetrix
though, if you want to custom build a CDF file.
To be honest, when I did this >2 years ago, there was no "official" GCRMA
implementation for Gene 1.0 ST arrays. Maybe Bioconductor has this now for
these arrays?
I've cc'd the aroma.affymetrix mailing list in the hope that someone has a
better solution for you.
Cheers,
Mark
On Jun 21, 2011, at 8:28 AM, Setsuko Sahara wrote:
> Dear Mark
>
> I found your previous e-mail regarding an application of gcrma to data for
> Gene ST arrays. Do you happen to have a chance to provide us of antigenomic
> probesets somewhere? Or do you recommend to try your previous scripts?
>
> Sincerely,
>
> setsu
>
>
>
> Mark Robinson
> Sun, 29 Mar 2009 23:12:09 -0700
>
> Hi Mario.
>
> I have made some modifications to the reading of probe_tab files and
> to the computing of affinities so that this procedure can run now,
> either as you have done below by choosing lowly expressed probes, or
> (perhaps preferably) by using the 'antigenomic' probes on the array:
>
> library(aroma.affymetrix)
> verbose <- Arguments$getVerbose(-30); timestampOn(verbose)
>
> cdf<-AffymetrixCdfFile$fromChipType("HuGene-1_0-st-
> v1",verbose=verbose,tags="PD")
> cs<-AffymetrixCelSet$fromName("tissues",cdf=cdf,verbose=verbose)
>
> bcGc <- GcRmaBackgroundCorrection(cs,
> type="affinities",indicesNegativeControl = negativeControlIndices)
> csGBC <- process(bcGc,verbose=verbose)
>
> controlIndices <- which(!isPm(cdf))
> bc <-
> GcRmaBackgroundCorrection
> (cs,type="affinities",indicesNegativeControl=controlIndices)
> csBC <- process(bc,verbose=verbose,force=TRUE)
>
> I needed to make a CDF file that contained these antigenomic probesets
> as they are not present in the binary-converted CDF files I created
> previously. I will make available these CDFs once I can test
> everything.
>
> Unfortunately, I do not have a good way of testing that my
> modifications are doing exactly the right thing, as I am also not
> intimately familiar with the gcrma model/code. To be honest, I don't
> know of anyone that has successfully run gcrma on these chips or Exon
> ST chips, Bioconductor or otherwise. Do you? If so, please let me
> know.
>
> These changes can be made available in the next release or possibly
> earlier with a patch, but I just want to test the changes first.
>
> Cheers,
> Mark
>
> Setsuko Sahara
>
>
>
>
------------------------------
Mark Robinson, PhD (Melb)
Epigenetics Laboratory, Garvan
Bioinformatics Division, WEHI
e: mrobin...@wehi.edu.au
e: m.robin...@garvan.org.au
p: +61 (0)3 9345 2628
f: +61 (0)3 9347 0852
------------------------------
______________________________________________________________________
The information in this email is confidential and intended solely for the
addressee.
You must not disclose, forward, print or use it without the permission of the
sender.
______________________________________________________________________
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest
version of the package, 2) to report the output of sessionInfo() and
traceback(), and 3) to post a complete code example.
You received this message because you are subscribed to the Google Groups
"aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to aroma-affymetrix@googlegroups.com
To unsubscribe and other options, go to http://www.aroma-project.org/forum/
