metacountregressor 0.1.312__py3-none-any.whl → 0.1.313__py3-none-any.whl

metacountregressor/metaheuristics.py

@@ -615,10 +615,16 @@ class DifferentialEvolution(object):
                     obj_trial = self._obj_fun.get_fitness(trial, self.pf.get_objective_is_multi())
                     if not self.pf.get_objective_is_multi():
                         average_iteration += obj_trial[self._obj_fun._obj_1]
-                        logger(i, obj_trial, None, True, self.get_instance_name(), 1)
-                    else:
-                        logger(i, obj_trial, None, True, self.get_instance_name(), self.pf.get_objective_is_multi())
+                        try:
 
+                            logger(i, obj_trial, None, True, self.get_instance_name(), 1)
+                        except:
+                            pass
+                    else:
+                        try:
+                            logger(i, obj_trial, None, True, self.get_instance_name(), self.pf.get_objective_is_multi())
+                        except:    
+                            pass
                 except Exception as e:
                     print('why is there an exception')
                     print(e)
@@ -637,7 +643,7 @@ class DifferentialEvolution(object):
 
                         iterations_without_improvement = 0
                         self._population[j] = obj_trial
-
+                        
                         logger(self.it_process, obj_trial, self._population, True,
                                self.instance_number + '/population_logger_strict_non_pareto.csv', 1)
                         logger(self.it_process, obj_trial, self._pareto_population, True,

metacountregressor/solution.py

@@ -3557,6 +3557,9 @@ class ObjectiveFunction(object):
                                varnames if x not in correlation]
             return
 
+    
+
+
     def _chol_mat(self, correlationLength, br, Br_w, correlation):
         # if correlation = True correlation pos is randpos, if list get correct pos
         
@@ -3601,13 +3604,15 @@ class ObjectiveFunction(object):
             try:
                 rv_val = chol[chol_count] if is_correlated else br_w[rv_count]
             except:
+                print('exception here start') 
                 print(self.rdm_cor_fit, 'rdm_cor_fit')
                 print(self.rdm_fit, 'rdm_fit')
                 print('varnames', varnames)
-                print(br, 'br')
+                print('br', br, 'br')
                 print(Br_w, 'Br_w')
                 print(chol, 'chol')
                 print(br_w, 'br_w')
+                print('exception here end')
             chol_mat_temp[rv_count_all, rv_count_all] = rv_val
             rv_count_all += 1
             if is_correlated:
@@ -4869,10 +4874,10 @@ class ObjectiveFunction(object):
                 n_coeff = self.get_param_num(dispersion)
                 Kf_a, Kr_a, Kr_c, Kr_b_a, Kchol_a, Kh = self.get_num_params()
                 if Kchol_a != Kchol:
-                    print('hold qhy')
+                    print('this should not hh qhy')
 
-                if Kr_b != Kr_b_a:
-                    print('hold qhy')
+                if Kr_b != Kr_a:
+                    print('hold qhy this should never happen')
 
 
 
@@ -4886,13 +4891,14 @@ class ObjectiveFunction(object):
 
                     )
                 Bf = betas[0:Kf]  # Fixed betas
-
-
-           # Bf_new, Br_new, Br_std_new, Br_rema = self.extract_parameters(betas, Kf, Kr, Kchol_a, Kr_b_a)
+            TEST_ME = False
+            if not TEST_ME:
+                Bf, br, br_std, Br_rema = self.extract_parameters(betas, Kf, Kr, Kchol_a, Kr_b_a)
 
 
             Vdf = dev.np.einsum('njk,k -> nj', Xdf, Bf, dtype=np.float64)  # (N, P)
-            br = betas[Kf:Kf + Kr]
+            if TEST_ME:
+                br = betas[Kf:Kf + Kr]
 
 
 
@@ -4910,9 +4916,9 @@ class ObjectiveFunction(object):
             # Kchol_a + Krb_a
             #its grabbing from the
 
+            if TEST_ME:
 
-
-            brstd = betas[Kf + Kr:Kf + Kr + Kr_b + Kchol]
+                brstd = betas[Kf + Kr:Kf + Kr + Kr_b + Kchol]
 
             # initialises size matrix
             proba = []  # Temp batching storage
@@ -7209,6 +7215,20 @@ class ObjectiveFunction(object):
         indices = [i for i, name in enumerate(self._characteristics_names) if name in names and isinstance(name, str)]
         return indices
 
+    def sanity_check(self, Xr):
+        a =len(self.rdm_cor_fit)+len(self.rdm_fit) != Xr.shape[2]
+        if a:
+            print('why')
+            print('The number of random effects does not match the data')
+            print(Xr.shape)
+            print(self.rdm_cor_fit)
+            print(self.rdm_fit)
+            
+            #raise Exception('The number of random effects does not match the data')
+        else:
+            return True
+    
+    
     """
     This function is named 'makeRegression'. It takes several parameters; 'self', 'model_nature', 'layout', '*args', and '**kwargs'.
     The purpose of this function is to execute a regression based on the training and testing datasets.
@@ -7315,7 +7335,7 @@ class ObjectiveFunction(object):
             Xr_cor = df_tf[:, :, indices3]
             # FIXME not sure if this is the right way orientatied
             Xr = np.concatenate((Xr, Xr_cor), axis=2)
-
+            self.sanity_check(Xr)
         model_nature['Xr'] = Xr
         if self.is_multi:
             Xr_test = df_test[:, :, indices2]

{metacountregressor-0.1.312.dist-info → metacountregressor-0.1.313.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: metacountregressor
-Version: 0.1.312
+Version: 0.1.313
 Summary: Extensive Testing for Estimation of Data Count Models
 Home-page: https://github.com/zahern/CountDataEstimation
 Author: Zeke Ahern

{metacountregressor-0.1.312.dist-info → metacountregressor-0.1.313.dist-info}/RECORD

@@ -6,15 +6,15 @@ metacountregressor/halton.py,sha256=jhovA45UBoZYU9g-hl6Lb2sBIx_ZBTNdPrpgkzR9fng,
 metacountregressor/helperprocess.py,sha256=komn_EjCShFg-_8INh1Q06GdWWN4nxhHOgJ3V1X6aJI,26309
 metacountregressor/main.py,sha256=tGOm8DdbdyDf316qIxDAre6l6GzfJIWYNYIBaSeIemI,23685
 metacountregressor/main_old.py,sha256=eTS4ygq27MnU-dZ_j983Ucb-D5XfbVF8OJQK2hVVLZc,24123
-metacountregressor/metaheuristics.py,sha256=P0Xjlvhp1cEwZFACrqeeets6x8BK7F2iDyu1OfS4bog,107010
+metacountregressor/metaheuristics.py,sha256=Qhu06AdF8HwEpyN-Q-6deG4b4kr7BKslgL8i2T7t1aI,107234
 metacountregressor/pareto_file.py,sha256=whySaoPAUWYjyI8zo0hwAOa3rFk6SIUlHSpqZiLur0k,23096
 metacountregressor/pareto_logger__plot.py,sha256=mEU2QN4wmsM7t39GJ_XhJ_jjsdl09JOmG0U2jICrAkI,30037
 metacountregressor/setup.py,sha256=5UcQCCLR8Fm5odA3MX78WwahavxFq4mVD6oq0IuQvAY,936
 metacountregressor/single_objective_finder.py,sha256=jVG7GJBqzSP4_riYr-kMMKy_LE3SlGmKMunNhHYxgRg,8011
-metacountregressor/solution.py,sha256=k6SkxlVAR6MoPZi5zJ0xh39f-6UGxqZX1bL_XRb2MPk,319289
+metacountregressor/solution.py,sha256=biAgfV5X4EXhzw1vGkNec53C-6ph02WpUWPtCQJfHKc,320016
 metacountregressor/test_generated_paper2.py,sha256=pwOoRzl1jJIIOUAAvbkT6HmmTQ81mwpsshn9SLdKOg8,3927
-metacountregressor-0.1.312.dist-info/licenses/LICENSE.txt,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-metacountregressor-0.1.312.dist-info/METADATA,sha256=2SPbRZ94EwJOVLWVbFdHjbmHXYa53AEmrIpmjcmESMA,23581
-metacountregressor-0.1.312.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
-metacountregressor-0.1.312.dist-info/top_level.txt,sha256=zGG7UC5WIpr76gsFUpwJ4En2aCcoNTONBaS3OewwjR0,19
-metacountregressor-0.1.312.dist-info/RECORD,,
+metacountregressor-0.1.313.dist-info/licenses/LICENSE.txt,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+metacountregressor-0.1.313.dist-info/METADATA,sha256=kUz8SzD7wzpQxrwHbBSqlIqA7ZRnkJClAMvxYlZVrns,23581
+metacountregressor-0.1.313.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+metacountregressor-0.1.313.dist-info/top_level.txt,sha256=zGG7UC5WIpr76gsFUpwJ4En2aCcoNTONBaS3OewwjR0,19
+metacountregressor-0.1.313.dist-info/RECORD,,