sklearn解决药学类分析化学第26页第12题

 真是没想到，我居然有使用sklearn完成分析化学作业的一天。

# -*- coding: UTF-8 -*-
# Created on 2023/10/15-17:55
def sub_problem_1_and_3():
    from sklearn.linear_model import LinearRegression

    c = [[0], [10], [20], [30], [40], [50]]
    A = [0.009, 0.135, 0.261, 0.383, 0.519, 0.653]

    reg = LinearRegression().fit(c, A)

    ferrous_concentration = (0.350 - reg.intercept_) / reg.coef_[0]

    print(f"回归方程式为：A = {reg.intercept_:.3f} + {reg.coef_[0]:.3f} * c")
    print(f"亚铁离子浓度为 {ferrous_concentration:.3f}")


def sp2():
    from scipy.stats import pearsonr
    import numpy as np

    c = [0, 10, 20, 30, 40, 50]
    A = [0.009, 0.135, 0.261, 0.383, 0.519, 0.653]

    corr_coef, _ = pearsonr(c, A)
    print(f"皮尔逊相关系数为 {corr_coef:f}")


sub_problem_1_and_3()

sp2()