plotastrodata 1.9.13__tar.gz → 1.9.14__tar.gz

{plotastrodata-1.9.13/plotastrodata.egg-info → plotastrodata-1.9.14}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: plotastrodata
-Version: 1.9.13
+Version: 1.9.14
 Summary: plotastrodata is a tool for astronomers to create figures from FITS files and perform fundamental data analyses with ease.
 Home-page: https://github.com/yusukeaso-astron/plotastrodata
 Download-URL: https://github.com/yusukeaso-astron/plotastrodata

{plotastrodata-1.9.13 → plotastrodata-1.9.14}/plotastrodata/__init__.py

@@ -1,4 +1,4 @@
 import warnings
 
 warnings.simplefilter('ignore', FutureWarning)
-__version__ = '1.9.13'
+__version__ = '1.9.14'

{plotastrodata-1.9.13 → plotastrodata-1.9.14}/plotastrodata/noise_utils.py

@@ -152,11 +152,17 @@ class Noise:
         if 'pbcor' in self.sigma:
             bounds.append([0.1, 2])
         # curve_fit does not work for this fitting.
-        # 0.01 in sigma is set only to search the best-fit parameters.
-        # Thus, this sigma does not justify the parameter errors.
-        fitter = EmceeCorner(bounds=bounds, model=model,
-                             xdata=self.hbin, ydata=self.hist,
-                             sigma=np.max(self.hist) * 0.01)
+        # For binned data, sigma^2 is the expected number of data in each bin.
+
+        def logl(p):
+            dh = self.hbin[1] - self.hbin[0]
+            ndata = len(np.ravel(self.data)) * dh
+            numobs = self.hist * ndata
+            numexp = model(self.hbin, *p) * ndata
+            chi2 = np.sum((numobs - numexp)**2 / numexp.clip(1, None))
+            return -0.5 * chi2
+
+        fitter = EmceeCorner(bounds=bounds, logl=logl)
         fitter.fit(**_kw)
         self.popt = fitter.popt
         self.mean = float(self.popt[1] * self.s0 + self.m0)

{plotastrodata-1.9.13 → plotastrodata-1.9.14/plotastrodata.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: plotastrodata
-Version: 1.9.13
+Version: 1.9.14
 Summary: plotastrodata is a tool for astronomers to create figures from FITS files and perform fundamental data analyses with ease.
 Home-page: https://github.com/yusukeaso-astron/plotastrodata
 Download-URL: https://github.com/yusukeaso-astron/plotastrodata