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 ?