能不能再上面写的代码上加上决定系数,均方误差,总偏差平方和,残差平方和的计算 import pandas as pd import numpy as np from scipy.linalg import inv # 读取数据集 data = pd.read_csv('your_dataset.csv') # 提取自变量和因变量 X = data.iloc[:, :-1].values Y = data.iloc[:, -1].values # 增加常数列 X = np.concatenate([np.ones((X.shape[0], 1)), X], axis=1) # 计算回归系数 beta = np.dot(np.dot(inv(np.dot(X.T, X)), X.T), Y) # 输出结果 print('Coefficients:', beta)
时间: 2024-02-13 08:05:56 浏览: 81
快速计算数组均方根误差RMSE 和RPD的程序
5星 · 资源好评率100%
当然可以,下面是加上这些计算的代码示例,包括计算决定系数,均方误差,总偏差平方和和残差平方和:
```python
import pandas as pd
import numpy as np
from scipy.linalg import inv
from sklearn.metrics import r2_score, mean_squared_error
# 读取数据集
data = pd.read_csv('your_dataset.csv')
# 提取自变量和因变量
X = data.iloc[:, :-1].values
Y = data.iloc[:, -1].values
# 增加常数列
X = np.concatenate([np.ones((X.shape[0], 1)), X], axis=1)
# 计算回归系数
beta = np.dot(np.dot(inv(np.dot(X.T, X)), X.T), Y)
# 输出结果
print('Coefficients:', beta)
# 计算决定系数
r2 = r2_score(Y, np.dot(X, beta))
print("R-squared:", r2)
# 计算均方误差
mse = mean_squared_error(Y, np.dot(X, beta))
print("Mean squared error:", mse)
# 计算总偏差平方和
tss = sum((Y - np.mean(Y)) ** 2)
print("Total sum of squares:", tss)
# 计算残差平方和
rss = sum((Y - np.dot(X, beta)) ** 2)
print("Residual sum of squares:", rss)
```
其中,`r2_score` 函数用于计算决定系数,`mean_squared_error` 函数用于计算均方误差。`tss` 和 `rss` 分别表示总偏差平方和和残差平方和,它们的计算方式同上一段回答中所述。
阅读全文