x1 = df.iloc[:, 0] x2 = df.iloc[:, 1] x3 = df.iloc[:, 2] y = df.iloc[:, 3] 拟合相关性

在这段代码中，`df` 是一个 Pandas DataFrame（数据框），`iloc` 是用来选择DataFrame中的特定行和列的索引方法。

`x1 = df.iloc[:, 0]`, `x2 = df.iloc[:, 1]`, 和 `x3 = df.iloc[:, 2]` 分别选取了 DataFrame 的前三列，从左到右，第一列为 x1，第二列为 x2，第三列为 x3。

`y = df.iloc[:, 3]` 则选取了 DataFrame 的第四列作为 y 值。

接下来，如果要进行拟合相关性分析，通常指的是计算这些变量之间的线性或非线性关系。这可以通过以下步骤实现：

1. **数据检查**：确保 `x1, x2, x3` 和 `y` 中没有缺失值或异常值。
2. **计算相关系数**：使用 `pandas` 的 `corr()` 或者 `scipy.stats.pearsonr()` 函数来计算 x1, x2, x3 对 y 的皮尔逊相关系数。
- 如果是直接的线性相关性分析，`pd.DataFrame.corr()` 可能足够。
- 如果需要更详细的统计描述（如偏斜度、比例等），可能需要用到 `scipy.stats.describe()` 和 `sklearn.metrics.pairwise_distances()`
3. **绘制散点图**：用 Matplotlib 或 Seaborn 绘制散点图，直观地查看各对变量之间的关系。
4. **回归分析**：如果感兴趣的是预测 y 基于 x1, x2, x3，可以用 `statsmodels.formula.api.ols` 或 `sklearn.linear_model.LinearRegression` 进行简单线性回归或其他模型。

假设我们想计算 x1, x2, x3 对 y 的皮尔逊相关系数并输出结果，代码可能是这样的：

import pandas as pd
import numpy as np
from scipy.stats import pearsonr

# 假设df已经包含了所有数据
corr_matrix = df.iloc[:, :4].corr()

# 计算x1, x2, x3与y的相关系数
x1_y_corr = corr_matrix.loc['x1', 'y']
x2_y_corr = corr_matrix.loc['x2', 'y']
x3_y_corr = corr_matrix.loc['x3', 'y']

print(f"x1 与 y 的相关系数: {x1_y_corr}")
print(f"x2 与 y 的相关系数: {x2_y_corr}")
print(f"x3 与 y 的相关系数: {x3_y_corr}")

# 使用pearsonr进一步获取协方差和显著性水平
(x1_p_value, _) = pearsonr(df['x1'], df['y'])
(x2_p_value, _) = pearsonr(df['x2'], df['y'])
(x3_p_value, _) = pearsonr(df['x3'], df['y'])

print("皮尔逊检验的 p 值:")
print(f"x1: {x1_p_value}, x2: {x2_p_value}, x3: {x3_p_value}")