0
$\begingroup$

I want to compute the estimate of $\beta$ for a linear model $Y = X\beta + \varepsilon $ with $$\varepsilon \sim N_d(0, \sigma^2V),$$ where $V$ is a $d\times d$ definitive posive, symmetric matrix.

It is straightforward to generalize the theory of ordinary linear models (where $\varepsilon \sim N_d(0, \sigma^2I_d))$, to this general case. In fact, it is sufficient to consider the Cholesky decomposition of the matrix $V$ and then transform the variables (details, for instance, here).

Now that the theory is clear, how can I specify and solve this problem in R?

Improve this question
$\endgroup$
7
  • 1
    $\begingroup$ You can use the gls function in the nlme library cran.r-project.org/web/packages/nlme . Technically this is called generalized least squares (GLS), the errors are unequal or correlated but still normal. In generalized linear models (GLM), the errors follow some distribution other than the normal distribution. $\endgroup$
    – StatGrrl
    Commented Feb 10, 2019 at 12:57
  • $\begingroup$ Thank you for the package suggestion and the clarification of the terminology $\endgroup$
    – Nisba
    Commented Feb 10, 2019 at 13:16
  • $\begingroup$ If V is a diagonal matrix then the errors have unequal variances but are independent - this is weighted least squares (WLS).The gls function allows for adding a correlation structure and/or weights. But if all you need is WLS, the lm function allows for this. $\endgroup$
    – StatGrrl
    Commented Feb 10, 2019 at 13:18
  • $\begingroup$ I need to consider a correlation structure $\endgroup$
    – Nisba
    Commented Feb 10, 2019 at 13:21
  • $\begingroup$ Is V known or must it be estimated? $\endgroup$
    – StatGrrl
    Commented Feb 10, 2019 at 13:31

0

Reset to default

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.