petpal 0.5.4__py3-none-any.whl → 0.5.5__py3-none-any.whl

petpal/kinetic_modeling/graphical_analysis.py

@@ -91,6 +91,50 @@ def fit_line_to_data_using_lls_with_rsquared(xdata: np.ndarray,
     return fit_ans[0][0], fit_ans[0][1], r_squared
 
 
+@numba.njit()
+def linear_least_squares_fit_with_stats(xdata: np.ndarray,
+                                        ydata: np.ndarray) -> tuple[float, float, float, float, float]:
+    """Fits a line to the data using least squares and explicitly computes:
+        - Fit R^2
+        - Standard error of the intercept
+        - Standard error of the slope
+
+    Args:
+        xdata (np.ndarray): X-coordinates of the data.
+        ydata (np.ndarray): Y-coordinates of the data.
+
+    Returns:
+        tuple: A tuple containing five float values: the intercept of the fitted line, the slope
+        of the fitted line, the r-squared value, the intercept standard error, and the slope
+        standard error.
+
+    See:
+        - https://mathworld.wolfram.com/LeastSquaresFitting.html
+    """
+    make_2d_matrix = _line_fitting_make_rhs_matrix_from_xdata
+    matrix = make_2d_matrix(xdata)
+    fit_ans = np.linalg.lstsq(matrix, ydata)
+
+    x_mean = np.mean(xdata)
+    y_mean = np.mean(ydata)
+
+    ss_res = fit_ans[1][0]
+    ss_tot = np.sum((y_mean - ydata) ** 2.)
+
+    n = len(xdata)
+
+    sum_square_xdiff = np.sum(xdata**2)-n*x_mean**2
+    sum_square_ydiff = np.sum(ydata**2)-n*y_mean**2
+    sum_square_xydiff = np.sum(xdata*ydata)-n*y_mean*x_mean
+
+    s = np.sqrt((sum_square_ydiff-sum_square_xydiff**2/sum_square_xdiff)/(n-2))
+
+    r_squared = 1.0 - ss_res / ss_tot
+
+    se_intercept = s*np.sqrt(1/n+np.mean(xdata)**2/sum_square_xdiff)
+    se_slope = s*sum_square_xdiff**(-0.5)
+    return fit_ans[0][0], fit_ans[0][1], r_squared, se_intercept, se_slope
+
 @numba.njit()
 def cumulative_trapezoidal_integral(xdata: np.ndarray,
                                     ydata: np.ndarray,

{petpal-0.5.4.dist-info → petpal-0.5.5.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: petpal
-Version: 0.5.4
+Version: 0.5.5
 Summary: PET-PAL (Positron Emission Tomography Processing and Analysis Library)
 Project-URL: Repository, https://github.com/PETPAL-WUSM/PETPAL.git
 Author: Bradley Judge

{petpal-0.5.4.dist-info → petpal-0.5.5.dist-info}/RECORD

@@ -19,7 +19,7 @@ petpal/input_function/idif_necktangle.py,sha256=o5kyAqyT4C6o7zELY4EjyHrkJyX1BWcx
 petpal/input_function/pca_guided_idif.py,sha256=MPB59K5Z5oyIunIWFqFQts61z647xawLNkv8wICrKYM,44821
 petpal/kinetic_modeling/__init__.py,sha256=tW4yRH3TwaXPwKPqdkrbQmSk9hjrF1yRkV_C59PPboQ,382
 petpal/kinetic_modeling/fit_tac_with_rtms.py,sha256=HpK7VWVCCNoSQABY9i28vYpZsMRmvgs4vdcM_ZbdaYE,20971
-petpal/kinetic_modeling/graphical_analysis.py,sha256=uq3cnBwPEcSccQb_vy6cSH1j3kI3GePYMvYehD5HZKs,49456
+petpal/kinetic_modeling/graphical_analysis.py,sha256=a7IOwYnG3Wao2XTjFgsPK563txm7s4lpMbaqUMQ5wUQ,51003
 petpal/kinetic_modeling/parametric_images.py,sha256=sXYracBFUtFyttO-6oiDAldnU8hPN6Y4vKOD1V-DnlE,47301
 petpal/kinetic_modeling/reference_tissue_models.py,sha256=FkLziIgtpA8tOL2gZJFg_nB8VPEBs40T7RsDAJ3nJ-A,39510
 petpal/kinetic_modeling/rtm_analysis.py,sha256=e3EuaHXml4PDALEczwyOPpnThINAGh41UKNlOQHAPqc,25945
@@ -63,8 +63,8 @@ petpal/visualizations/graphical_plots.py,sha256=ZCKUeLX2TAQscuHjA4bzlFm1bACHIyCw
 petpal/visualizations/image_visualization.py,sha256=Ob6TD4Q0pIrxi0m9SznK1TRWbX1Ea9Pt4wNMdRrTfTs,9124
 petpal/visualizations/qc_plots.py,sha256=iaCPe-LWWyM3OZzDPZodHZhP-z5fRdpUgaH7QS9VxPM,1243
 petpal/visualizations/tac_plots.py,sha256=zSGdptL-EnqhfDViAX8LFunln5a1b-NJ5ft7ZDcxQ38,15116
-petpal-0.5.4.dist-info/METADATA,sha256=hUe3oDjdNZ9hW9e8E_nnoo-9-2geYFMHjlWMAVcsLQ4,2478
-petpal-0.5.4.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
-petpal-0.5.4.dist-info/entry_points.txt,sha256=-AAg5GCaTRuwcLR7QWKPuKyBo8mTEMzspRmcBwaPdLo,692
-petpal-0.5.4.dist-info/licenses/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
-petpal-0.5.4.dist-info/RECORD,,
+petpal-0.5.5.dist-info/METADATA,sha256=RPzJJzXXau935KVTE1i0V7VsfqJJNmGyWWcEzROuivc,2478
+petpal-0.5.5.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
+petpal-0.5.5.dist-info/entry_points.txt,sha256=-AAg5GCaTRuwcLR7QWKPuKyBo8mTEMzspRmcBwaPdLo,692
+petpal-0.5.5.dist-info/licenses/LICENSE,sha256=OXLcl0T2SZ8Pmy2_dmlvKuetivmyPd5m1q-Gyd-zaYY,35149
+petpal-0.5.5.dist-info/RECORD,,