Кажется, что для решения задачи используется функция решить две матрицы. Может кто-нибудь подробно объяснить, что он делает с двумя матрицами? Почему возникла необходимость в c (rep (1, n)?
demo.data <-data.frame(y=c(1.55,0.42,1.29,0.73,0.76,-1.09,1.41,-0.32),x1=c(1.13,-0.73,0.12,0.52,-0.54,-1.15,0.20,-1.09),x2=c(1,0,1,1,0,1,0,1))
demo.data
Y <- matrix(demo.data$y)
Y
n <- nrow(demo.data)
n
X <- matrix(c(rep(1,n),demo.data$x1,demo.data$x2),nrow=n,ncol=3)
X
BETA.HAT <- solve(t(X)%*%X)%*%t(X)%*%Y
BETA.HAT