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)

请参阅Matrices and matrix computations in R上的以下链接，了解有关此事的更多详情。

Answer 2

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