我是R的新手并且坚持完成任务。如果有人可以提供帮助，我将非常感激！
这是任务： “当使用线性回归模式将每个L基因型标记与表型相关联时，使用线性回归模型计算三个β参数的MLE（ß^）（即，您的R代码必须包括MLE的公式）。 绘制针对L基因型计算的每个参数估计值的直方图（即三个直方图！）“
我有200个人，我的Xd载体和Xa载体是：
> Xd
[1] 0 0 0 0 1 0 0 NA 0 0 0 0 0 0 0 NA 0 NA 0 1 0 0 0 1 0
[26] 0 0 0 0 0 0 0 0 -1 0 0 1 -1 0 NA -1 0 0 0 0 0 0 NA 0 1
[51] 0 0 0 -1 0 0 0 NA 0 0 -1 -1 0 0 0 0 -1 0 0 -1 0 0 0 0 0
[76] 0 1 0 0 0 0 1 0 0 0 NA NA 0 0 0 NA 0 0 NA 0 -1 0 0 0 -1
[101] 0 0 0 NA NA 0 NA 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0
[126] 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 0 0 NA 1 0 0
[151] -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA 0 1 1 0 -1 0 -1
[176] 0 0 0 0 NA 0 0 0 1 0 0 0 0 0 0 0 -1 0 0 0 0 NA 0 0
> Xa
[1] 1 1 1 1 -1 1 1 NA 1 1 1 1 1 1 1 NA 1 NA 1 -1 1 1 1 -1 1
[26] 1 1 1 1 1 1 1 1 -1 1 1 -1 -1 1 NA -1 1 1 1 1 1 1 NA 1 -1
[51] 1 1 1 -1 1 1 1 NA 1 1 -1 -1 1 1 1 1 -1 1 1 -1 1 1 1 1 1
[76] 1 -1 1 1 1 1 -1 1 1 1 NA NA 1 1 1 NA 1 1 NA 1 -1 1 1 1 -1
[101] 1 1 1 NA NA 1 NA 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 1 1 1
[126] 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 NA -1 1 1
[151] -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 NA 1 -1 -1 1 -1 1 -1
[176] 1 1 1 1 NA 1 1 1 -1 1 1 1 1 1 1 1 -1 1 1 1 1 NA 1 1
>
我做的是：
> Xmu=rep(1,200)
> X=cbind(Xmu, Xa, Xd)
然后我收到以下错误 在cbind（Xmu，Xa，Xd）中： 结果行数不是向量长度的倍数（arg 2）
这是什么意思？如何计算我的ß参数的MLE？我会这样做：
Y <- 1 + Xa*1 + Xd*0 + rnorm(200,0,sqrt(1))
betas <- solve(t(X)%*% X) %*% t(X)%*% Y
beta_mu <- betas[1]
beta_a <- betas[2]
beta_d <- betas[3]
此外，“您的代码必须包含MLE公式” - 部分让我感到困惑！？感谢