[BioC] URGENT Help required: Getting this error with GAGE analysis
Jung [guest]
guest at bioconductor.org
Wed Jan 25 19:56:56 CET 2012
Hi there,
I am getting this error on R console. The format of the input files, Micro_array_dataset.txt and Gene_set.gmt is shown below:
Micro_array_dataset.txt sample shown below, it's a tab delimited file:
Gene TCGA_AB_2803 TCGA_AB_2805 TCGA_AB_2806 TCGA_AB_2807 TCGA_AB_2808 TCGA_AB_2810
TSPAN6 0.87040619 0 0 1.49516441 0.145067698 0
TNMD 0.14858841 0.037147103 0 0 0 0
DPM1 45.07516892 63.28462838 64.9535473 41.44256757 46.61233108 42.66807432
SCYL3 12.62339374 22.15584843 25.60428336 13.98088962 17.0583196 4.036902801
Gene_set.gmt is a Gene Matrix Transposed file format required for GSEA analysis, sample shown below:
chr5q23 na PTMAP2 FTMT PRR16 MGC32805
chr16q24 na LOC642452 SNORD68 RPL13 TRAPPC2L
chr8q24 na LYPD2 KIFC2 EIF2C2 ADCY8
chr13q11 na LOC645626 ZMYM2 FGF9 XPO4
chr7p21 na DNAH11 LOC221814 MYL7 CRS
Any kind of help to get rid of the error shown below will be really appreciated.
Thanks,
Jung
-- output of sessionInfo():
> Micro_array_data <- readExpData(file = "Micro_array_dataset.txt")
> Gene_set <- readList("Gene_set.gmt")
> Reference_condition <- c(1,3,5)
> Target_condition <- c(2,4,6)
> A1_compare_un <- gage(Micro_array_data, Gene_set, ref = Reference_condition, samp = Target_condition)
Error in saaPrep(exprs, ref = ref, samp = samp, same.dir = same.dir, compare = compare, :
exprs needs to be a numeric matrix or vector
> # Essential_member_genes <- essGene(Gene_set, Micro_array_data,ref = NULL)
> # Non_redundant_significant_gene_set_list <- esset.grp()
>
--
Sent via the guest posting facility at bioconductor.org.
More information about the Bioconductor
mailing list