r语言非整合参数报错怎么解决呀！

Question

这个东西卡在这好久了，矩阵的行列也没错呀，为啥老给我报错，用Matlab可以运行，但在r studio里就不给过，求解答

Answer 1

引自免费微信小程序：皆我百晓生

你的问题出在矩阵乘法上。在R语言中，%*%是矩阵乘法的操作符。出现"非整合参数"的错误通常意味着你试图进行的矩阵乘法操作不符合矩阵乘法的规则。具体来说，第一个矩阵的列数必须等于第二个矩阵的行数。在你的代码中，t(B) %*% B的操作是不正确的，因为转置后的矩阵t(B)的列数（也就是原矩阵B的行数）为5，而矩阵B的列数为2，它们并不匹配。

解决此问题的方法是重新检查你的矩阵乘法操作，确保它们符合矩阵乘法的规则。根据你的代码和目的，我认为你可能想要计算的是线性模型的系数。在R语言中，通常我们会使用 lm() 函数来拟合线性模型，然后获取系数。如果你一定要手动计算，你可以按照下面的步骤进行：

首先确保你的矩阵维度正确，然后计算系数。这里是一个可能的解决方案：

# 创建Y矩阵
Y <- matrix(c(2.86, 3.4, 4, 4.81, 6), ncol = 1)

# 创建B矩阵（包括X和截距）
# 这里假设你有5个观测值，所以我们需要一个4x5的矩阵（加上截距项）
# 我们首先创建一个包含观测值的矩阵（除了截距项），然后添加一列全为1的截距项
X <- matrix(c(-3.68, -6.77, -10.41, -14.75, -20.05), nrow = 5, ncol = 1)
B <- cbind(X, rep(1, 5))  # 添加截距项

# 计算系数
coefficients <- solve(t(B) %*% B) %*% t(B) %*% Y
print(coefficients)

请注意，这里我假设你的线性模型有一个自变量和一个截距项。如果你的模型有多个自变量，你需要相应地调整你的 X 和 B 矩阵。另外，你也可以直接使用R语言的内置函数 lm() 来拟合线性模型并获取系数，例如 lm(Y ~ X)。