Hello, I'm very new in this field. We used single-cell RT-PCR on neurons to quantify the expression of 60 genes (Fluidigm) and then try to map them to their electric phenotypes .

As expected , we got a lot of "zeros" values for the genes that are expressed less in average. I 'd like to fit the distribution of zeros count (y axis) vs the mean value of expression with a negative binomial as advised in the article (Droplet scRNA-seq is not zero-inflated, svenson ,nature biotech, 2020 link I tried to use their code ( https://github.com/vals/Blog/tree/master/171116-zero-inflation) but I did not well understand how to transfer my values into the Anndata format. I tried many thing but I 'm stuck Is there a simpler way to do it ? I'm rather familiar with R Should I use the raw count values or the log normalized ones? Thanks a lot