用R中的delta方法估计线性回归的截距标准误差

# 生成模拟数据
set.seed(123)
n <- 100
x <- runif(n, min = 0, max = 10)
y <- 2 + 3 * x + rnorm(n, mean = 0, sd = 1)

# 拟合线性回归模型
model <- lm(y ~ x)

# 计算方差-协方差矩阵
vcov_matrix <- vcov(model)

# 提取截距的方差
intercept_variance <- vcov_matrix[1, 1]

# 计算截距的标准误差
intercept_se <- sqrt(intercept_variance)

# 输出截距的标准误差
print(intercept_se)

summary(model)$coefficients[1, "Std. Error"]