[BioC] Two class paired analysis using siggenes
Mary Winn
mary.e.winn at gmail.com
Thu Jul 24 23:52:18 CEST 2008
I am using siggenes to do a two class paired analysis of 36
microarrays. I do not appear to have any problems but upon
summarizing the results I receive "NA" for fold change. I believe I
read somewhere in the vingette that fold change is only calculated for
two class unpaired analysis. Is this true? If not, what must I do in
order to calculate the fold change in a two class paired analysis?
Thanks,
Mary
University of California, San Diego
> K<-18
> cl<-c((1:18),-1:-18)
> mydata<-
as
.matrix
(read.table("filter.txt",header=TRUE,sep="\t",row.names=1,as.is=TRUE))
> sam.out<-sam(mydata,cl,rand=123456)
> sum.sam.out<-summary(sam.out,5.8)
>sum.sam.out
SAM Analysis for the Two-Class Paired Case
s0 = 0
Number of permutations: 100
MEAN number of falsely called variables is computed.
Delta: 5.8
cutlow: -7.086
cutup: 6.704
p0: 0.146
Identified Genes: 7988
Falsely Called Genes: 0
FDR: 0
Identified Genes (using Delta = 5.8):
Row d.value stdev rawp q.value R.fold
Name
1 15497 60.02 0.0715 0 0 NA
1439466_s_at
2 7512 52.20 0.0501 0 0 NA
1427036_a_at
3 11424 45.64 0.0501 0 0 NA
1433603_at
4 22076 42.52 0.0526 0 0 NA
1451019_at
5 24816 42.11 0.0720 0 0 NA
1455269_a_at
6 13851 40.63 0.0775 0 0 NA
1436983_at
7 15488 40.49 0.0757 0 0 NA
1439455_x_at
8 25290 40.40 0.0370 0 0 NA
1455975_x_at
9 6774 40.08 0.0814 0 0 NA
1426025_s_at
10 2060 39.92 0.0629 0 0 NA
1418594_a_at
More information about the Bioconductor
mailing list