[BioC] HTqPCR - how to read data file with multiple samples
Heidi Dvinge
heidi at ebi.ac.uk
Fri Jan 8 09:40:20 CET 2010
Hello Adam,
originally HTqPCR only read in 1 sample per file. I've recently added the
parameter n.data to readCtData() to indicate the number of sample in each
file being processed (along with fixing a few bugs). I still need to
submit this to BioC though, so you can either wait a few days and get this
from the the devel repository, or if it's urgent I can send you the
corrected readCtData() function off-list.
Cheers
\Heidi
> I have a data file (created by SDS), with 6 samples, like this:
> # Status Sample Detector Task Avg Ct
> 1 Analyzed 1 MammU6-4395470 Endogenous Control 22.733
> 2 Analyzed 2 MammU6-4395470 Endogenous Control 22.582
> 3 Analyzed 3 MammU6-4395470 Endogenous Control 22.656
> 4 Analyzed 4 MammU6-4395470 Endogenous Control 22.785
> 5 Analyzed 5 MammU6-4395470 Endogenous Control 22.684
> 6 Analyzed 7 MammU6-4395470 Endogenous Control 22.411
>
> I tried running:
>
> readCtData("Data.txt", SDS=TRUE, n.features=381, type=5, Ct=6,
> header=TRUE, feature=4)
>
> the HTqPCR package reads the data file as if there were only 1 sample:
>
> An object of class "qPCRset"
> Size: 381 features, 1 samples
> Feature types: Endogenous Control, Target
> Feature names: MammU6-4395470 MammU6-4395470 MammU6-4395470 ...
> Feature classes:
> Feature categories: OK
> Sample names: Data NA NA ...
>
> How do I read the data into HTqPCR so that it recognizes that I have 6
> samples?
>
> ./adam
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> 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