1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
\name{predict.glmrob}
\alias{predict.glmrob}
\title{Predict Method for Robust GLM ("glmrob") Fits}
\description{
Obtains predictions and optionally estimates standard errors of those
predictions from a fitted \emph{robust} generalized linear model (GLM)
object.
}
\usage{
\method{predict}{glmrob}(object, newdata = NULL,
type = c("link", "response", "terms"), se.fit = FALSE,
dispersion = NULL, terms = NULL, na.action = na.pass, \dots)
}
\arguments{
%% the following is +- copy-pasted from predict.glm.Rd:
\item{object}{a fitted object of class inheriting from \code{"glmrob"}.}
\item{newdata}{optionally, a data frame in which to look for variables with
which to predict. If omitted, the fitted linear predictors are used.}
\item{type}{the type of prediction required. The default is on the
scale of the linear predictors; the alternative \code{"response"}
is on the scale of the response variable. Thus for a default
binomial model the default predictions are of log-odds (probabilities
on logit scale) and \code{type = "response"} gives the predicted
probabilities. The \code{"terms"} option returns a matrix giving the
fitted values of each term in the model formula on the linear predictor
scale.
The value of this argument can be abbreviated.
}
\item{se.fit}{logical switch indicating if standard errors are required.}
\item{dispersion}{the dispersion of the GLM fit to be assumed in
computing the standard errors. If omitted, that returned by
\code{summary} applied to the object is used.}
\item{terms}{with \code{type="terms"} by default all terms are returned.
A character vector specifies which terms are to be returned}
\item{na.action}{function determining what should be done with missing
values in \code{newdata}. The default is to predict \code{NA}.}
\item{\dots}{optional further arguments, currently simply passed to
\code{\link{predict.lmrob}()}.}
}
% \details{
% If necessary, more details than the description above ~~
% }
\value{
%% the following is +- copy-pasted from predict.glm.Rd:
%% also correct,here ?
If \code{se = FALSE}, a vector or matrix of predictions.
If \code{se = TRUE}, a list with components
\item{fit}{Predictions}
\item{se.fit}{Estimated standard errors}
\item{residual.scale}{A scalar giving the square root of the
dispersion used in computing the standard errors.}
}
\author{Andreas Ruckstuhl}
\seealso{
\code{\link{glmrob}()} to fit these robust GLM models,
\code{\link{residuals.glmrob}()} and other methods;
\code{\link{predict.lm}()}, the method used for a non-robust fit.
}
\examples{
data(carrots)
## simplistic testing & training:
i.tr <- sample(24, 20)
fm1 <- glmrob(cbind(success, total-success) ~ logdose + block,
family = binomial, data = carrots, subset = i.tr)
fm1
predict(fm1, carrots[-i.tr, ]) # --> numeric vector
predict(fm1, carrots[-i.tr, ],
type="response", se = TRUE)# -> a list
% FIXME: gives a "bad" error -- should rather say "not yet implemented"
% or implement it !
% predict(fm1, carrots[-i.tr, ], interval = "confidence")
% predict(fm1, carrots[-i.tr, ], interval = "prediction")
data(vaso)
Vfit <- glmrob(Y ~ log(Volume) + log(Rate), family=binomial, data=vaso)
newd <- expand.grid(Volume = (V. <- seq(.5, 4, by = 0.5)),
Rate = (R. <- seq(.25,4, by = 0.25)))
p <- predict(Vfit, newd)
filled.contour(V., R., matrix(p, length(V.), length(R.)),
main = "predict(glmrob(., data=vaso))", xlab="Volume", ylab="Rate")
}
\keyword{models}
\keyword{regression}
|
