[BioC] re incomplete analysis in Deseq
Julian [guest]
guest at bioconductor.org
Mon Mar 5 11:50:37 CET 2012
I'm using deseq with 454 data and it worked for one set of data but the same script is failing me the second time around with a different set of experimental data.
The input data is a matrix of counts of 454 seqs per sample, I have 36 pre and 36 post samples.
When I run the estmateSizeFactors I get all my samples as NA.
Any ideas why?
-- output of sessionInfo():
> cds <- estimateSizeFactors(cds) # Estimates size factors based on the count data
> sizeFactors( cds )
X01_MA_1 X02_MA_10 X03_MA_100 X04_MA_102 X05_MA_11 X06_MA_13 X07_MA_14 X08_MA_15 X09_MA_17 X10_MA_18 X11_MA_19
NA NA NA NA NA NA NA NA NA NA NA
X12_MA_2 X13_MA_20 X14_MA_22 X15_MA_23 X16_MA_24 X17_MA_25 X18_MA_4 X19_MA_47 X20_MA_5 X21_MA_69 X22_MA_7
NA NA NA NA NA NA NA NA NA NA NA
X23_MA_71 X24_MA_73 X25_MA_75 X26_MA_77 X27_MA_79 X28_MA_8 X29_MA_81 X30_MA_83 X31_MA_86 X32_MA_88 X33_MA_9
NA NA NA NA NA NA NA NA NA NA NA
X34_MA_90 X35_MA_92 X36_MA_94 X37_MA_96 X38_MA_98 X39_MA_101 X40_MA_103 X41_MA_26 X42_MA_27 X43_MA_29 X44_MA_30
NA NA NA NA NA NA NA NA NA NA NA
X45_MA_31 X46_MA_33 X47_MA_34 X48_MA_36 X49_MA_37 X50_MA_40 X51_MA_41 X52_MA_42 X53_MA_43 X54_MA_44 X55_MA_45
NA NA NA NA NA NA NA NA NA NA NA
X56_MA_46 X57_MA_49 X58_MA_50 X59_MA_52 X60_MA_54 X61_MA_55 X62_MA_70 X63_MA_72 X64_MA_74 X65_MA_76 X66_MA_78
NA NA NA NA NA NA NA NA NA NA NA
X67_MA_80 X68_MA_82 X69_MA_84 X70_MA_87 X71_MA_89 X72_MA_91 X73_MA_93 X74_MA_95 X75_MA_97 X76_MA_99
NA NA NA NA NA NA NA NA NA NA
--
Sent via the guest posting facility at bioconductor.org.
More information about the Bioconductor
mailing list