# Least absolute deviations regression

### Description

This function is used to fit linear models considering Laplace errors.

### Usage

1 2 |

### Arguments

`formula` |
an object of class |

`data` |
an optional data frame containing the variables in the model. If
not found in |

`method` |
character string specifying the algorithm to use. The default
algorithm is the Barrodale and Roberts algorithm |

`subset` |
an optional expression indicating the subset of the rows of data that should be used in the fit. |

`na.action` |
a function that indicates what should happen when the data contain NAs. |

`control` |
a list of control values for the estimation algorithm to replace
the default values returned by the function |

`model, x, y` |
logicals. If |

`contrasts` |
an optional list. See the |

### Value

an object of class `lad`

representing the linear model fit. Generic
function `print`

, show the results of the fit.

### Author(s)

The design was inspired by the R function `lm`

.

### References

Barrodale, I., and Roberts, F.D.K. (1974).
Solution of an overdetermined system of equations in the L1 norm.
*Communications of the ACM* **17**, 319-320.

Phillips, R.F. (2002).
Least absolute deviations estimation via the EM algorithm.
*Statistics and Computing* **12**, 281-285.

### Examples

1 | ```
lad(stack.loss ~ ., data = stackloss, method = "EM")
``` |