RiboDiff icon indicating copy to clipboard operation
RiboDiff copied to clipboard

Divide by zero error, even with test data

Open IsabelBirds opened this issue 6 years ago • 0 comments

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!

IsabelBirds avatar Nov 11 '19 13:11 IsabelBirds