I have 6 samples of RNA-seq files without any technical replicates (Saprophyte_30,Saprophyte_37,Healed_30,Healed_37,surgery_30,surgery_37). These are A.flavus (fungus)samples which were grown at 30 and 37 deg C of saprophytic condition and its clinical isolates. I am following "Count-based differential expression analysis of RNA sequencing data using R and Bioconductor" protocol paper for my data analysis.
I got the error when I executed following command in DESeq
cds=estimateDispersions(cds,method="blind") res=nbinomTest(cds,"Saphrophyte _30","Saphrophyte _37")
Error in if (dispTable(cds)[condA] == "blind" || dispTable(cds)[condB] == : missing value where TRUE/FALSE needed
Can anyone tell me how to solve this?