From 465c0f3a1721081b0ca26eea72e0b113ecda03be Mon Sep 17 00:00:00 2001 From: wbalmer Date: Wed, 3 Jan 2024 14:18:33 -0500 Subject: [PATCH] correct chi2 calculation --- backtracks/backtracks.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/backtracks/backtracks.py b/backtracks/backtracks.py index dff77a1..eac9c29 100644 --- a/backtracks/backtracks.py +++ b/backtracks/backtracks.py @@ -198,7 +198,7 @@ def __init__(self, target_name: str, candidate_file: str, nearby_window: float = self.stationary_params = [self.ra0, self.dec0, 0, 0, 0, self.rao, self.deco, self.pmrao, self.pmdeco, self.paro, self.radvelo] # Compute useful chi2 value self.stationary_loglike = self.loglike(self.stationary_params) - self.stationary_chi2_red = np.log(-2.*self.stationary_loglike)/(len(self.epochs)-self.ndim) + self.stationary_chi2_red = -2.*self.stationary_loglike/((2*(len(self.epochs)-1))-self.ndim) jd_start, jd_end, number = ephem_open() print('[BACKTRACK INFO]: Opened ephemeris file') @@ -582,7 +582,7 @@ def fit(self, dlogz=0.5, npool=4, dynamic=False, nlive=200, mpi_pool=False, resu # Compute useful chi2 value self.median_loglike = self.loglike(self.run_median) - self.median_chi2_red = np.log(-2.*self.median_loglike)/(len(self.epochs)-self.ndim) + self.median_chi2_red = -2.*self.median_loglike/((2*len(self.epochs))-self.ndim) return self.results @@ -604,7 +604,7 @@ def load_results(self, fileprefix: str = './'): # recompute useful chi2 value self.median_loglike = self.loglike(self.run_median) - self.median_chi2_red = np.log(-2.*self.median_loglike)/(len(self.epochs)-self.ndim) + self.median_chi2_red = -2.*self.median_loglike/((2*len(self.epochs))-self.ndim) def generate_plots( self,