[BioC] Most stable gene pairs in array experiment
Marcelo Laia
marcelolaia at gmail.com
Tue Oct 20 14:30:02 CEST 2009
2009/10/20 Michael Dondrup <Michael.Dondrup at bccs.uib.no>:
> Hi,
>
> you can try something like this or use two for loops:
>
>
>> apply (mygenes, 1 , function(row) { apply (mygenes, 1, function(x) {
>> var(row-x) } ) } )
>
> geneA geneB geneC geneD geneE
> geneA 0.00000 0.04108 0.06397 0.12217 0.08233
> geneB 0.04108 0.00000 0.15543 0.12807 0.09365
> geneC 0.06397 0.15543 0.00000 0.08628 0.08903
> geneD 0.12217 0.12807 0.08628 0.00000 0.01517
> geneE 0.08233 0.09365 0.08903 0.01517 0.00000
>
>
> Cheers
> Michael
Hi,
I am following the discussion and I'm finding very interesting. Congratulations!
After this, I could compare the two genes, two-by-two, and I could
conclude that the pair with minor variance are the two most stable
genes of all?
Is this genes appropriated for qPCR internal control? Or am I totally
wrong here?
Thank you very much!
--
Marcelo Luiz de Laia
Universidade do Estado de Santa Catarina
UDESC - www.cav.udesc.br
Lages - SC - Brazil
Linux user number 487797
More information about the Bioconductor
mailing list