I have RNAseq data in the form of normalized counts. The counts themselves follow an over-dispersed negative binomial, and I would like to generate random data which replicate these distributions. For example, I've subsampled a set of the counts ~4200 n, and the descriptives are as follows,
Statistic Std. Error Mean 46.28 32.6 Lower -17.64 Upper 110.19 Median 0 Variance 4493154.99 STD 2119.71 Minimum 0 Maximum 136974 Skewness 63.86 0.04 Kurtosis 4123.38 0.08
I've tried the following:
data <- rnbinom(4200,46.28,0.9)
It seems to give a distribution similar to my data, but cuts off the upper bound range. So I'm not exactly sure if this is a robust simulation.
I do apologize for what is likely a rudimentary question. I'm not exactly a skilled statistician, so it's somewhat difficult to follow some of the documentation on these processes. Would anyone be kind enough to explain this to me in the way a novice would understand?