[BioC] Draft support for HTA 2.0 with oligo
Mark Cowley
m.cowley at garvan.org.au
Mon Feb 10 04:48:37 CET 2014
thanks Benilton!
On 09/02/2014, at 8:47 AM, Benilton Carvalho <beniltoncarvalho at gmail.com>
wrote:
> Hello,
>
> given the several requests I got re: supporting HTA 2.0 chips when using
> oligo, I'm posting this message to bioc (rather than bioc-devel).
>
> I'll add support for HTA on our next release. A draft is ready through the
> devel branch. I may need to make some adjustments prior to the April
> release.
>
> For those interested in giving it a try, you'll need to have a developer
> installation of Bioconductor (please, make sure to use a different path to
> set your library and avoid mixing & matching). Then install oligo v.
> 1.27.3. With the devel bioc installed, you can use the usual:
>
> biocLite('oligo')
>
> You will also need an annotation package specific for HTA, which you can
> download from here: (source version only)
>
> https://dl.dropboxusercontent.com/u/83643/pd.hta.2.0_3.8.0.tar.gz
>
> After downloading the annotation package, install it using:
>
> ## assuming that you loaded R from the directory where the
> ## package was downloaded to.
> install.packages('pd.hta.2.0_3.8.0.tar.gz', repos=NULL)
>
>
> Once both are installed, you should be able to run the standard:
>
> library(oligo)
> cels = list.celfiles()
> raw = read.celfiles(cels)
> pp0 = rma(raw, target='core')
> pp1 = rma(raw, target='probeset')
>
> to obtain your preprocessed data through oligo.
>
> The annotation package will be made available through Bioconductor in the
> next few weeks (devel branch again).
>
> If you're in the mood of building your own annotation package, after
> ensuring you have the aformentioned packages (and versions) installed, also
> install pdInfoBuilder (version 1.27.4). Then, you can use the template
> below to create your package:
>
> ## code to build HTA annot pkg
> library(pdInfoBuilder)
> seed <- new("AffyHTAPDInfoPkgSeed",
> version="3.8.0",
> license="Artistic-2.0",
> pgfFile=dir(patt=".pgf$"),
> clfFile=dir(patt=".clf$"),
> probeFile=dir(patt="probeset.csv$"),
> transFile=dir(patt="transcript.csv$"),
> coreMps=dir(patt=".mps$"),
> geneArray=TRUE,
> author="Benilton Carvalho",
> email="beniltoncarvalho at gmail.com",
> biocViews="AnnotationData",
> genomebuild="hg19",
> organism="Homo sapiens",
> species="Homo sapien",
> url="http://about.me/benilton")
> makePdInfoPackage(seed, destDir=".")
> ## end
>
> If you have any feedback/comment, please contact me.
>
> benilton
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list