adaptivetesting 1.1.5__tar.gz → 1.1.6__tar.gz

{adaptivetesting-1.1.5 → adaptivetesting-1.1.6}/.github/workflows/python-test.yml

@@ -16,7 +16,7 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-        python-version: ["3.10", "3.11", "3.12"]
+        python-version: ["3.12", "3.13", "3.14"]
 
     steps:
     - uses: actions/checkout@v3
@@ -42,7 +42,7 @@ jobs:
     strategy:
       fail-fast: false
       matrix:
-          python-version: ["3.10", "3.11", "3.12"]
+          python-version: ["3.12", "3.13", "3.14"]
       
     steps:
       - uses: actions/checkout@v3

{adaptivetesting-1.1.5 → adaptivetesting-1.1.6}/CHANGELOG.md

@@ -1,3 +1,6 @@
+# Version 1.1.6
+- Bug Fig #49 (optimization interval issue in BayesModal)
+
 # Version 1.1.5
 - Bug Fix #46 (numerical stability issue concerning Bayesian ability estimation)
 
@@ -10,4 +13,4 @@
 # Version 1.1.2
 
 - The `ItemPool` class now correctly honors the `id` field of items. (#30, #31)
-- Examples and the readme file have been updated accordingly.
+- Examples and the readme file have been updated accordingly.

{adaptivetesting-1.1.5 → adaptivetesting-1.1.6}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: adaptivetesting
-Version: 1.1.5
+Version: 1.1.6
 Summary: adaptivetesting is a Python package that can be used to simulate and evaluate custom CAT scenarios as well as implement them in real-world testing scenarios from a single codebase
 Project-URL: Homepage, https://github.com/condecon/adaptivetesting
 Project-URL: Documentation, https://github.com/condecon/adaptivetesting/tree/main/docs
@@ -12,7 +12,7 @@ Keywords: computerized-adaptive-testing,item-response-theory,psychology,statisti
 Requires-Python: >=3.12
 Requires-Dist: matplotlib>=3.9.0
 Requires-Dist: numpy>=2.0.0
-Requires-Dist: pandas>=2.2.0
+Requires-Dist: pandas>=2.3.0
 Requires-Dist: scipy>=1.12.0
 Requires-Dist: tqdm>=4.66.0
 Description-Content-Type: text/markdown

{adaptivetesting-1.1.5 → adaptivetesting-1.1.6}/adaptivetesting/math/estimators/__bayes_modal_estimation.py

@@ -62,7 +62,8 @@ class BayesModal(IEstimator):
                 self.c,
                 self.d,
                 self.response_pattern,
-                self.prior
+                self.prior,
+                self.optimization_interval
             )
         # else, we have to calculate the full posterior distribution
         # because the optimizers do not correctly identify the maximum of the function
@@ -91,7 +92,7 @@ class BayesModal(IEstimator):
                 ])
 
                 # add prior
-                unmarginalized_posterior = lik_values * self.prior.pdf(mu)
+                unmarginalized_posterior = lik_values * np.log(self.prior.pdf(mu))
                 # find argmin and return mu
                 estimate_index = np.argmin(unmarginalized_posterior)
                 return float(mu[estimate_index].astype(float))

{adaptivetesting-1.1.5 → adaptivetesting-1.1.6}/adaptivetesting/math/estimators/__functions/__bayes.py

@@ -1,6 +1,5 @@
 import numpy as np
 from scipy.optimize import minimize_scalar, OptimizeResult # type: ignore
-from .__estimators import likelihood
 from ..__prior import Prior
 from ....models.__algorithm_exception import AlgorithmException
 from .__estimators import probability_y0, probability_y1
@@ -19,15 +18,15 @@ def maximize_posterior(
 
     Args:
         a (np.ndarray): item parameter a
-    
+
         b (np.ndarray): item parameter b
-        
+
         c (np.ndarray): item parameter c
-        
+
         d (np.ndarray): item parameter d
-        
+
         response_pattern (np.ndarray): response pattern (simulated or user generated)
-        
+
         prior (Prior): prior distribution
 
         optimization_interval (Tuple[float, float]): interval used for the optimization function
@@ -39,18 +38,39 @@ def maximize_posterior(
         p1 = probability_y1(mu, a, b, c, d)
         p0 = probability_y0(mu, a, b, c, d)
 
-        log_likelihood = np.sum((response_pattern * np.log(p1 + 1e-300)) + \
-                    ((1 - response_pattern) * np.log(p0 + 1e-300)))
+        log_likelihood = np.sum(
+            (response_pattern * np.log(p1 + 1e-300)) + (
+                (1 - response_pattern) * np.log(p0 + 1e-300)
+            ))
         log_prior = np.log(prior.pdf(mu) + 1e-300)
-    
+
         return log_likelihood + log_prior
-    
+
     result: OptimizeResult = minimize_scalar(lambda mu: -log_posterior(mu),
                                              bounds=optimization_interval,
                                              method="bounded") # type: ignore
-    
+
     if not result.success:
         raise AlgorithmException(f"Optimization failed: {result.message}")
-    
+
     else:
         return float(result.x)
+
+
+def likelihood(mu: np.ndarray,
+               a: np.ndarray,
+               b: np.ndarray,
+               c: np.ndarray,
+               d: np.ndarray,
+               response_pattern: np.ndarray
+               ):
+    
+    p1 = probability_y1(mu, a, b, c, d)
+    p0 = probability_y0(mu, a, b, c, d)
+
+    log_likelihood = np.sum(
+        (response_pattern * np.log(p1 + 1e-300)) + (
+            (1 - response_pattern) * np.log(p0 + 1e-300)
+        ))
+
+    return -np.exp(log_likelihood)

{adaptivetesting-1.1.5 → adaptivetesting-1.1.6}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name="adaptivetesting"
-version="1.1.5"
+version="1.1.6"
 description="adaptivetesting is a Python package that can be used to simulate and evaluate custom CAT scenarios as well as implement them in real-world testing scenarios from a single codebase"
 readme = "README.md"
 authors = [
@@ -9,7 +9,7 @@ authors = [
 requires-python = ">=3.12"
 dependencies = [
     "numpy>=2.0.0",
-    "pandas>=2.2.0",
+    "pandas>=2.3.0",
     "scipy>=1.12.0",
     "tqdm>=4.66.0",
     "matplotlib>=3.9.0"