Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/influence.mlm.R

This collection of functions is designed to compute regression deletion
diagnostics for multivariate linear models following Barrett & Ling (1992)
that are close analogs of
methods for univariate and generalized linear models handled by the
`influence.measures`

in the stats package.

In addition, the functions provide diagnostics for deletion of
subsets of observations of size `m>1`

.

1 2 3 4 5 6 7 8 |

`model` |
An |

`do.coef` |
logical. Should the coefficients be returned in the |

`m` |
Size of the subsets for deletion diagnostics |

`x` |
An |

`FUN` |
For |

`funnames` |
logical. Should the |

`...` |
Other arguments passed to methods |

`digits` |
Number of digits for the print method |

`influence.mlm`

is a simple wrapper for the computational function, `mlm.influence`

designed to provide an S3 method for class `"mlm"`

objects.

There are still infelicities in the methods for the `m>1`

case in the current implementation.
In particular, for `m>1`

, you must call `influence.mlm`

directly, rather than using
the S3 generic `influence()`

.

`influence.mlm`

returns an S3 object of class `inflmlm`

, a list with the following components

`m` |
Deletion subset size |

`H` |
Hat values, |

`Q` |
Residuals, |

`CookD` |
Cook's distance values |

`L` |
Leverage components |

`R` |
Residual components |

`subsets` |
Indices of the observations in the subsets of size |

`labels` |
Observation labels |

`call` |
Model call for the |

`Beta` |
Deletion regression coefficients– included if |

Michael Friendly

Barrett, B. E. and Ling, R. F. (1992).
General Classes of Influence Measures for Multivariate Regression.
*Journal of the American Statistical Association*, **87**(417), 184-191.

`influencePlot.mlm`

, `mlm.influence`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | ```
# Rohwer data
Rohwer2 <- subset(Rohwer, subset=group==2)
rownames(Rohwer2)<- 1:nrow(Rohwer2)
Rohwer.mod <- lm(cbind(SAT, PPVT, Raven) ~ n+s+ns+na+ss, data=Rohwer2)
# m=1 diagnostics
influence(Rohwer.mod)
# try an m=2 case
res2 <- influence.mlm(Rohwer.mod, m=2, do.coef=FALSE)
res2.df <- as.data.frame(res2)
head(res2.df)
scatterplotMatrix(log(res2.df))
influencePlot(Rohwer.mod, id.n=4, type="cookd")
# Sake data
Sake.mod <- lm(cbind(taste,smell) ~ ., data=Sake)
influence(Sake.mod)
influencePlot(Sake.mod, id.n=3, type="cookd")
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.