R-codes for Calculating Fluctuation Assay Results and 95% Confidence Intervals Based on Ma-Sandri-Sarkar Maximum Likelihood

Abstract

Abstract The Luria-Delbrück fluctuation assay is an essential experiment in calculating mutation rates, especially in genetic and mutation research. Its reliability and accuracy have made it the go-to method for numerous researchers. In this article, we provide R-codes that statistically analyze the assay results more easily and offer the most challenging codes for calculating 95% confidence intervals based on the gold standard method “Ma-Sandri-Sarkar Maximum Likelihood.” Recently, the maximization of the likelihood function through optimization functions in R can be a challenging task. The recursive format of the likelihood function is known to cause memory stack issues. Our findings indicate that utilizing a non-recursive version of the function can increase the tractability of the maximization process. With these codes, future scientists can unlock valuable statistical insights related to the biological mechanisms that drive genetic variation and can, therefore, contribute to developing novel therapeutic interventions and innovative solutions to various biological and medical challenges.