Divide by zero error, even with test data
I've tested installing RiboDiff locally, and installing via conda. I consistently get this error when running on my data and on the provided test data.
/nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/ribodiff/adjlik.py:14: RuntimeWarning: divide by zero encountered in divide n = 1 / disp /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/ribodiff/adjlik.py:15: RuntimeWarning: invalid value encountered in divide p = n / (n + mu) /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/scipy/stats/_discrete_distns.py:179: RuntimeWarning: invalid value encountered in greater_equal return (n > 0) & (p >= 0) & (p <= 1) /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/scipy/stats/_discrete_distns.py:179: RuntimeWarning: invalid value encountered in less_equal return (n > 0) & (p >= 0) & (p <= 1) /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/scipy/optimize/_minimize.py:761: RuntimeWarning: Method 'bounded' does not support relative tole rance in x; defaulting to absolute tolerance. "defaulting to absolute tolerance.", RuntimeWarning) /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/ribodiff/fitdisp.py:35: RuntimeWarning: invalid value encountered in greater idx = np.logical_and(dispRaw > lowerBound, dispRaw < upperBound).nonzero()[0] /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/ribodiff/fitdisp.py:35: RuntimeWarning: invalid value encountered in less idx = np.logical_and(dispRaw > lowerBound, dispRaw < upperBound).nonzero()[0] /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/statsmodels/genmod/generalized_linear_model.py:273: DomainWarning: The identity link function do es not respect the domain of the Gamma family. DomainWarning) /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/ribodiff/fitdisp.py:49: RuntimeWarning: invalid value encountered in less if np.nonzero(dispFitted < 0)[0].size > 0: /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/statsmodels/genmod/generalized_linear_model.py:273: DomainWarning: The identity link function do es not respect the domain of the Gamma family. DomainWarning) /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/ribodiff/adjdisp.py:18: RuntimeWarning: divide by zero encountered in log logprior = (np.log(disp) - np.log(dispFitted)) ** 2 / (2 * varPrior ** 2) /nobackup/cm13ijb/miniconda/envs/RiboDiff/lib/python2.7/site-packages/ribodiff/adjdisp.py:18: RuntimeWarning: invalid value encountered in log logprior = (np.log(disp) - np.log(dispFitted)) ** 2 / (2 * varPrior ** 2)
The resulting outputs appear to be expected. Thanks!