朱莉娅中没有截距的单变量回归

Question

我可以使用以下方法进行多元回归：

using GLM
using Random
obs = 1000000
X = rand(obs)
y = X .+ rand(Normal(),obs) .+ 7
fit(LinearModel, rand(obs,2), y)

但是我似乎无法像

那样进行单变量回归

fit(LinearModel, X, y)

由于X不是矩阵，导致方法错误。

我可以使用以下命令对截距进行回归：

fit(LinearModel, hcat(ones(obs),X), y)

如何在Julia中不进行截距的情况下进行单变量回归？

Answer 1

如果您要问的是，最简单的方法是使用以下方法将X转换为Matrix：

fit(LinearModel, hcat(X), y)

或

fit(LinearModel, reshape(X, obs, 1), y)

在这种情况下也可以简单地写：

X\y

将在不使用GLM的情况下为您提供所需的答案（但在这种情况下，您只会得到估计值）。