blob: efdea3d502d01cfdca997e87260a12b0feb5a6fd (
plain) (
blame)
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
|
\name{calplot}
\alias{calplot}
\alias{calplot.default}
\alias{calplot.lm}
\title{Plot calibration graphs from univariate linear models}
\description{
Produce graphics of calibration data, the fitted model as well
as prediction and confidence bands.
}
\usage{
calplot(object, xlim = "auto", ylim = "auto",
xlab = "Concentration", ylab = "Response", alpha=0.05)
}
\arguments{
\item{object}{
A univariate model object of class \code{\link{lm}} or
\code{\link[MASS:rlm]{rlm}}
with model formula \code{y ~ x} or \code{y ~ x - 1}.
}
\item{xlim}{
The limits of the plot on the x axis.
}
\item{ylim}{
The limits of the plot on the y axis.
}
\item{xlab}{
The label of the x axis.
}
\item{ylab}{
The label of the y axis.
}
\item{alpha}{
The confidence level for the confidence and prediction bands.
}
}
\value{
A plot of the calibration data, of your fitted model as well as lines showing
the confidence limits as well as the prediction limits.
}
\examples{
# Example of a Calibration plot for a weighted regression
data(massart97ex3)
attach(massart97ex3)
yx <- split(y,factor(x))
s <- round(sapply(yx,sd),digits=2)
w <- round(1/(s^2),digits=3)
weights <- w[factor(x)]
m <- lm(y ~ x,w=weights)
calplot(m)
}
\author{
Johannes Ranke
\email{jranke@uni-bremen.de}
\url{http://www.uft.uni-bremen.de/chemie/ranke}
}
\keyword{regression}
|
