R中的线性回归模型矩阵

Question

考虑模型 ，在哪里

$ X $的矩阵是

如何将此矩阵输入R？

Answer 1

假设x = 1,2，...，N，N = 10并且y = 2x +ε，ε〜N（0,1），那么您将编写如下内容：

N = 10;
set.seed(123)
x =  1:N
e = rnorm(N)  
y = 2*x + e;
mod <- lm( y ~x);
Xmatrix = matrix( c(rep(1,N), x), ncol=2)

请参阅以下有关R中矩阵和矩阵计算的链接， 以获取有关此问题的更多详细信息。

Answer 2

您可以将该矩阵输入R并执行Y〜X，也可以省略截距并执行Y〜$ {X} _i $ -1，其中$ X_i $只是一个向量。