Hi!
I have a row with data looking like this:
GeneID srr1 srr2 srr3 srr4 srr5 srr6 srr7 srr8 srr9 srr10 srr11 srr12 srr13 srr14 srr15 srr16 srr17
ENSG00000223972 0 1 0 1 0 0 0 0 1 0 0 0 1 1 0 0 0
I want to check it for the zero-inflated negative binomial distribution. So I want to use the function zeroinfl from the pscl package. The function requires a formula as its first parametr.
zeroinfl(formula = , data = my_vector, dist = "negbin", EM = TRUE)
What should I put in the formula ?
Ok, thank you. About my data - it present the gene (ENSG00000223972) expression over a number of people - they are named srr1.... The initial idea was to check whether the expression may be discribed by the negative binomial or Poisson distribution.