metacountregressor 0.1.311__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

@@ -3549,37 +3549,30 @@ class ObjectiveFunction(object):
 
         return inv_hess
 
+    def _correlated_pos(self, correlation, varnames):
+        if (isinstance(correlation, list)):
+            correlationpos = [varnames.tolist().index(x) for x in
+                              varnames if x in correlation]
+            uncorrelatedpos = [varnames.tolist().index(x) for x in
+                               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
-        dont_run = 0
-        if dont_run:
-            correlationpos = []
-            varnames = ['a', 'b', 'c', 'd', 'e', 'f']
-            varnames = np.asarray(varnames)
-            randvars = ['c', 'd', 'e', 'f']
-            correlation = ['e', 'f']
-            if randvars:
-                # Position of correlated variables within randvars
-                correlationpos = [varnames.tolist().index(x)
-                                  for x in varnames if x in randvars]
-            if (isinstance(correlation, list)):
-                self.correlationpos = [varnames.tolist().index(x) for x in
-                                       varnames if x in correlation]
-                self.uncorrelatedpos = [varnames.tolist().index(x) for x in
-                                        varnames if x not in correlation]
-            # if correlation = True correlation pos is randpos, if list get correct pos
-            correlationpos = []
-        else:
-            varnames = self.none_handler(
+        
+        
+        
+        varnames = self.none_handler(
                 self.rdm_fit) + self.none_handler(self.rdm_cor_fit)
 
-        if (isinstance(correlation, list)):
+        
             # Kchol, permutations of specified params in correlation list
-            Kchol = int((len(correlation) *
+        Kchol = int((len(correlation) *
                          (len(correlation) + 1)) / 2)
-        else:
-            # i.e. correlation = True, Kchol permutations of random paramaters
-            Kchol = int((len(br) * (len(br) + 1)) / 2)
+        
 
         # creating cholesky matrix for the variance-covariance matrix
         # all random variables not included in correlation will only
@@ -3608,7 +3601,18 @@ class ObjectiveFunction(object):
                 is_correlated = True
             else:
                 is_correlated = False
-            rv_val = chol[chol_count] if is_correlated else br_w[rv_count]
+            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, '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:
@@ -4870,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')
 
 
 
@@ -4887,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]
 
 
 
@@ -4911,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
@@ -7210,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.
@@ -7316,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.311.dist-info → metacountregressor-0.1.313.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: metacountregressor
-Version: 0.1.311
+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.311.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=h-sHkShA-P5oWJuowCRt0J2wixcxCgoL_tlyM2NCBxQ,319693
+metacountregressor/solution.py,sha256=biAgfV5X4EXhzw1vGkNec53C-6ph02WpUWPtCQJfHKc,320016
 metacountregressor/test_generated_paper2.py,sha256=pwOoRzl1jJIIOUAAvbkT6HmmTQ81mwpsshn9SLdKOg8,3927
-metacountregressor-0.1.311.dist-info/licenses/LICENSE.txt,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-metacountregressor-0.1.311.dist-info/METADATA,sha256=U0e8zrNVJziARj-x5z-ezThGsPGvwBwLY_mr4gVjZ5Q,23581
-metacountregressor-0.1.311.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
-metacountregressor-0.1.311.dist-info/top_level.txt,sha256=zGG7UC5WIpr76gsFUpwJ4En2aCcoNTONBaS3OewwjR0,19
-metacountregressor-0.1.311.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,,