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
相关推荐
weixin_3077791340 分钟前
PySpark实现ABC_manage_channel逻辑
开发语言·python·spark
??? Meggie2 小时前
【Python】保持Selenium稳定爬取的方法(防检测策略)
开发语言·python·selenium
酷爱码4 小时前
如何通过python连接hive,并对里面的表进行增删改查操作
开发语言·hive·python
画个大饼4 小时前
Go语言实战:快速搭建完整的用户认证系统
开发语言·后端·golang
喵先生!5 小时前
C++中的vector和list的区别与适用场景
开发语言·c++
Thomas_YXQ5 小时前
Unity3D Lua集成技术指南
java·开发语言·驱动开发·junit·全文检索·lua·unity3d
xMathematics6 小时前
计算机图形学实践:结合Qt和OpenGL实现绘制彩色三角形
开发语言·c++·qt·计算机图形学·cmake·opengl
yuanManGan8 小时前
C++入门小馆: 深入了解STLlist
开发语言·c++
北极的企鹅888 小时前
XML内容解析成实体类
xml·java·开发语言
BillKu8 小时前
Vue3后代组件多祖先通讯设计方案
开发语言·javascript·ecmascript