I looked for an answer but could not find it.. My supervisor replaced values of outliers in our expression data set to "NA", using the count data that is corrected for library size (means median-of-ratios). As we want to use the same processing steps, we want to assign these outliers also in my work (network analysis).
Is it possible to insert the identified outliers also in the dds object at the same step in DESeq2 (so after dds = estimateSizeFactors(dds))? In other words, could I assign these outlying values as NA in the count matrix and save this to the dds, so I could work further with it when I transform my data for network analysis? Or should I assign the NAs to the count matrix before creating the dds object? (I plan to use the k-nearest neighbour method after to replace the NAs.)
Thanks a lot,