Question: fRMA for Gene ST arrays
0
gravatar for maria
6.0 years ago by
maria10
maria10 wrote:
> > > Kamila, > > You can preprocess Gene ST arrays using the latest version of frma -- > just load the data using the oligo package. And you can obtain GNUSE > values using the GNUSE function on the preprocessed data object. > > A barcode implementation for the Gene ST arrays is not currently > available, but it is in the works -- hopefully it will be ready before > the fall BioC release. > > Best, > Matt > > On Wed, Apr 3, 2013 at 8:41 AM, Naxerova, Kamila > <naxerova at="" ...=""> wrote: > > Dear list, > > > > to preprocess Affymetrix whole transcript arrays, use of oligo and xps is recommended over the affy > package. I am interested in using the GNUSE algorithm to assign quality scores to some arrays, and would > also like to try out the barcode function in the frma package. I am now wondering... is frma compatible with > whole transcript arrays? > > > > Thanks a lot. > > Kamila > > _______________________________________________ > > Bioconductor mailing list > > Bioconductor at ... > > https://stat.ethz.ch/mailman/listinfo/bioconductor > > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > Hi, First, I apologize if I posted this message at the wrong place. I tried to run barcode for the hugene.1.0.st.v1 and got an output for only 4000 probes. The problem seems that when barcode maps the probe names from "hugene.1.0.st.v1frmavecs" to my data the overlap is only of 4000 although tau and mu contains about 37000 items. I wonder if I did a mistake in my coding (I'm new in microarray analysis) or if this is really what happens. (see script below). I also tried to edit the code of barcode to add a new platform and I assigned the vector hugene.1.0.st.v1frmavecs to the pkg. (pkg <- "hugene.1.0.st.v1frmavecs") but it does not help. Thank you in advance for your help Maria library("oligo") library(pd.hugene.1.0.st.v1) celFiles = list.celfiles() affyRaw= read.celfiles(celFiles) library("hugene.1.0.st.v1frmavecs") data("hugene.1.0.st.v1frmavecs") data("hugene.1.0.st.v1barcodevecs") library("frma") eset=frma(affyRaw,input.vecs=hugene.1.0.st.v1frmavecs, summarize="median_polish" ) b=barcode(exprs(eset),platform=NULL,mu=hugene.1.0.st.v1barcodevecs$mu, tau=hugene.1.0.st.v1barcodevecs$tau, output="p-value") R version 2.15.2 (2012-10-26) Platform: i386-w64-mingw32/i386 (32-bit) attached base packages: [1] stats graphics utils methods base other attached packages: [1] frma_1.10.0 hugene.1.0.st.v1frmavecs_1.0.0 [3] pd.hugene.1.0.st.v1_3.8.0 RSQLite_0.11.4 [5] DBI_0.2-7 oligo_1.22.0 [7] Biobase_2.18.0 oligoClasses_1.20.0 [9] BiocGenerics_0.4.0 loaded via a namespace (and not attached): [1] affxparser_1.30.2 affy_1.36.1 affyio_1.26.0 [4] BiocInstaller_1.8.3 Biostrings_2.26.3 bit_1.1-11 [7] codetools_0.2-8 ff_2.2-12 foreach_1.4.1 [10] GenomicRanges_1.10.7 grDevices_2.15.2 IRanges_1.16.6 [13] iterators_1.0.6 MASS_7.3-23 parallel_2.15.2 [16] preprocessCore_1.20.0 splines_2.15.2 stats4_2.15.2 [19] tools_2.15.2 zlibbioc_1.4.0
microarray probe oligo frma • 1.2k views
ADD COMMENTlink modified 6.0 years ago by Matthew McCall830 • written 6.0 years ago by maria10
Answer: fRMA for Gene ST arrays
0
gravatar for Matthew McCall
6.0 years ago by
United States
Matthew McCall830 wrote:
Maria, It looks like you version of frma is quite old (the current release version is 1.14.0). I'd start by updating the relevant packages. What are the dimensions of the preprocessed expression values: dim(exprs(eset)) And the dimensions of the barcode values: dim(b) Best, Matt On Thu, Jan 16, 2014 at 10:29 AM, maria <maria.pacheco at="" uni.lu=""> wrote: >> >> >> Kamila, >> >> You can preprocess Gene ST arrays using the latest version of frma -- >> just load the data using the oligo package. And you can obtain GNUSE >> values using the GNUSE function on the preprocessed data object. >> >> A barcode implementation for the Gene ST arrays is not currently >> available, but it is in the works -- hopefully it will be ready before >> the fall BioC release. >> >> Best, >> Matt >> >> On Wed, Apr 3, 2013 at 8:41 AM, Naxerova, Kamila >> <naxerova at="" ...=""> wrote: >> > Dear list, >> > >> > to preprocess Affymetrix whole transcript arrays, use of oligo and xps > is recommended over the affy >> package. I am interested in using the GNUSE algorithm to assign quality > scores to some arrays, and would >> also like to try out the barcode function in the frma package. I am now > wondering... is frma compatible with >> whole transcript arrays? >> > >> > Thanks a lot. >> > Kamila >> > _______________________________________________ >> > Bioconductor mailing list >> > Bioconductor at ... >> > https://stat.ethz.ch/mailman/listinfo/bioconductor >> > Search the archives: > http://news.gmane.org/gmane.science.biology.informatics.conductor >> > > > Hi, > First, I apologize if I posted this message at the wrong place. I tried to > run barcode for the hugene.1.0.st.v1 and got an output for only 4000 probes. > The problem seems that when barcode maps the probe names from > "hugene.1.0.st.v1frmavecs" to my data the overlap is only of 4000 although > tau and mu contains about 37000 items. I wonder if I did a mistake in my > coding (I'm new in microarray analysis) or if this is really what happens. > (see script below). I also tried to edit the code of barcode to add a new > platform and I assigned the vector hugene.1.0.st.v1frmavecs to the pkg. > (pkg <- "hugene.1.0.st.v1frmavecs") but it does not help. > Thank you in advance for your help > > Maria > > > > > > > library("oligo") > library(pd.hugene.1.0.st.v1) > > celFiles = list.celfiles() > > affyRaw= read.celfiles(celFiles) > > library("hugene.1.0.st.v1frmavecs") > data("hugene.1.0.st.v1frmavecs") > data("hugene.1.0.st.v1barcodevecs") > > > library("frma") > > eset=frma(affyRaw,input.vecs=hugene.1.0.st.v1frmavecs, > summarize="median_polish" ) > > b=barcode(exprs(eset),platform=NULL,mu=hugene.1.0.st.v1barcodevecs$mu, > tau=hugene.1.0.st.v1barcodevecs$tau, output="p-value") > > > > > R version 2.15.2 (2012-10-26) > Platform: i386-w64-mingw32/i386 (32-bit) > > > attached base packages: > [1] stats graphics utils methods base > > other attached packages: > [1] frma_1.10.0 hugene.1.0.st.v1frmavecs_1.0.0 > [3] pd.hugene.1.0.st.v1_3.8.0 RSQLite_0.11.4 > [5] DBI_0.2-7 oligo_1.22.0 > [7] Biobase_2.18.0 oligoClasses_1.20.0 > [9] BiocGenerics_0.4.0 > > loaded via a namespace (and not attached): > [1] affxparser_1.30.2 affy_1.36.1 affyio_1.26.0 > [4] BiocInstaller_1.8.3 Biostrings_2.26.3 bit_1.1-11 > [7] codetools_0.2-8 ff_2.2-12 foreach_1.4.1 > [10] GenomicRanges_1.10.7 grDevices_2.15.2 IRanges_1.16.6 > [13] iterators_1.0.6 MASS_7.3-23 parallel_2.15.2 > [16] preprocessCore_1.20.0 splines_2.15.2 stats4_2.15.2 > [19] tools_2.15.2 zlibbioc_1.4.0 > > _______________________________________________ > 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 -- Matthew N McCall, PhD 112 Arvine Heights Rochester, NY 14611 Cell: 202-222-5880
ADD COMMENTlink written 6.0 years ago by Matthew McCall830
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 419 users visited in the last hour