[R] Influence.measures
Peter B. Mandeville
mandevip at uaslp.mx
Wed Jun 23 11:00:19 CEST 1999
I am using rw0641 with Windows 98. To list just the influential
repetitiones that result from "influence.measures", I am using the input
result <- lm(y~x)
and the code from the example in the help for "influence.measures"
INFLM <- function(result){
inflm <- influence.measures(result)
which(apply(inflm$is.inf,1,any))
}
It works fine up to now with the exception of Huber's data published in
Atkinson 1985:8 (Plots, Transformations and Regression)
Observation x y
1 -4 2.48
2 -3 0.73
3 -2 -0.04
4 -1 -1.44
5 0 -1.32
6 10 0.00
which is not very well excepted. The command
> inflm <- influence.measures(result)
reports
Warning: NaNs produced in function "qf"
and
> which(apply(inflm$is.inf,1,any))
reports
Error: names attribute must be the same length as the vector
What should I do?
Thank you very much.
Peter B.
--
Peter B. Mandeville mandevip at deimos.tc.uaslp.mx
Jefe del Depto. de Informática y Bioestadística rpe1531 at pasteur.fmed.uaslp.mx
Facultad de Medicine Tel: 48 26-23-45 ext. 232
Universidad Autónoma de San Luis Potosí Fax: 48 28-23-52
Av. V. Carranza 2405
Col. Los Filtros
Apartado Postal 145
San Luis Potosí, S.L.P.
78210 México
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list