R语言建模线性回归

一、

a.

R 复制代码
# 给定的 (x, y) 数据
x <- c(2, 9, 10, 7)
y <- c(3, 13, 12, 11)

# 线性回归模型 y = a + bx
model1 <- lm(y ~ x)
summary(model1)  # 查看回归结果

# 提取系数 a 和 b
a <- coef(model1)[1]
b <- coef(model1)[2]

# 预测值
y_pred <- predict(model1)

# 计算 RMSE
rmse <- sqrt(mean((y - y_pred)^2))
cat("RMSE for model 1:", rmse, "\n")

# 计算 R²
r_squared <- summary(model1)$r.squared
cat("R² for model 1:", r_squared, "\n")

b.

R 复制代码
# b
# 给定的 (x, y) 数据
x <- c(2, 9, 10, 7)
y <- c(3, 13, 12, 11)
# 定义 z = log(x)
z <- log(x)

# 拟合 z 的线性回归模型
model_log <- lm(y ~ z)

# 打印模型参数 a 和 b
summary(model_log)

# 计算预测值
y_pred_log <- predict(model_log)

# 计算 RMSE
rmse_log <- sqrt(mean((y - y_pred_log)^2))
cat("RMSE (log):", rmse_log, "\n")

# 计算 R^2
r_squared_log <- summary(model_log)$r.squared
cat("R-squared (log):", r_squared_log, "\n")

c.

二、

a.

题目给出了三组数据点 (x_{11}, x_{12}, y_1), (x_{21}, x_{22}, y_2), (x_{31}, x_{32}, y_3):

• (x_{11}, x_{12}, y_1) = (0, 1, 1)

• (x_{21}, x_{22}, y_2) = (3, -2, 9)

• (x_{31}, x_{32}, y_3) = (-1, 4, -2)

问题 a 要求估计多元线性回归模型的 R^2 ,但是不需要进行计算。

R^2 是决定系数,表示模型解释的方差比例。 R^2 = 1 表示模型完美拟合数据,而 R^2 = 0 表示模型无法解释任何方差。

因为只有 3 个数据点,而我们有两个自变量 x_1 和 x_2 ,意味着回归模型将完美拟合这些数据点。因此,在这种情况下,模型的 R^2 值将是 1,因为我们可以通过回归完全解释数据的变化。

b.

R 复制代码
# 构造 X 矩阵和 y 向量
X <- matrix(c(1, 1, 1, 1,   # 截距项
              0, 3, -1, 2,  # x1
              1, -2, 4, 2), # x2
            nrow=4, byrow=TRUE)

y <- c(1, 9, -2, 2)

# 计算回归系数
beta <- solve(t(X) %*% X) %*% t(X) %*% y
beta
相关推荐
侃侃_天下11 小时前
最终的信号类
开发语言·c++·算法
echoarts11 小时前
Rayon Rust中的数据并行库入门教程
开发语言·其他·算法·rust
Aomnitrix12 小时前
知识管理新范式——cpolar+Wiki.js打造企业级分布式知识库
开发语言·javascript·分布式
每天回答3个问题12 小时前
UE5C++编译遇到MSB3073
开发语言·c++·ue5
伍哥的传说12 小时前
Vite Plugin PWA – 零配置构建现代渐进式Web应用
开发语言·前端·javascript·web app·pwa·service worker·workbox
小莞尔13 小时前
【51单片机】【protues仿真】 基于51单片机八路抢答器系统
c语言·开发语言·单片机·嵌入式硬件·51单片机
我是菜鸟0713号13 小时前
Qt 中 OPC UA 通讯实战
开发语言·qt
JCBP_13 小时前
QT(4)
开发语言·汇编·c++·qt·算法
Brookty13 小时前
【JavaEE】线程安全-内存可见性、指令全排序
java·开发语言·后端·java-ee·线程安全·内存可见性·指令重排序
百锦再14 小时前
[特殊字符] Python在CentOS系统执行深度指南
开发语言·python·plotly·django·centos·virtualenv·pygame