[BioC] Calculating probe affinities directly from sequence data - GCRMA
Steve Pederson
stephen.pederson at adelaide.edu.au
Sun Sep 27 09:43:44 CEST 2009
Hi,
I'm working on some alternative ideas for fitting exon arrays & was
wondering if there is a way to calculate theoretical probe affinities
directly from sequence data such as:
> myseq<-paste(sample(rep(c("A","C","T","G"),10),25),collapse="")
> myseq
[1] "ACTCACTTCGCTAAAACATTGGCTA"
I can find the spline coefficients using "affinity.spline.coefs", but it
looks like you can only perform the calculations via
> compute.affinities.local(myAffyBatch)
Currently, you can't specify exon arrays as an AffyBatch, but the
sequences for each probe are available in a probe_tab file.
Thanks,
Steve
More information about the Bioconductor
mailing list