title: “Coursera Regresion model” author: “Makasani Christopher” date: “2023/09/11” output: html_document
Executive summary
library(ggplot2)
data(mtcars)
summary(mtcars)
## mpg cyl disp hp
## Min. :10.40 Min. :4.000 Min. : 71.1 Min. : 52.0
## 1st Qu.:15.43 1st Qu.:4.000 1st Qu.:120.8 1st Qu.: 96.5
## Median :19.20 Median :6.000 Median :196.3 Median :123.0
## Mean :20.09 Mean :6.188 Mean :230.7 Mean :146.7
## 3rd Qu.:22.80 3rd Qu.:8.000 3rd Qu.:326.0 3rd Qu.:180.0
## Max. :33.90 Max. :8.000 Max. :472.0 Max. :335.0
## drat wt qsec vs
## Min. :2.760 Min. :1.513 Min. :14.50 Min. :0.0000
## 1st Qu.:3.080 1st Qu.:2.581 1st Qu.:16.89 1st Qu.:0.0000
## Median :3.695 Median :3.325 Median :17.71 Median :0.0000
## Mean :3.597 Mean :3.217 Mean :17.85 Mean :0.4375
## 3rd Qu.:3.920 3rd Qu.:3.610 3rd Qu.:18.90 3rd Qu.:1.0000
## Max. :4.930 Max. :5.424 Max. :22.90 Max. :1.0000
## am gear carb
## Min. :0.0000 Min. :3.000 Min. :1.000
## 1st Qu.:0.0000 1st Qu.:3.000 1st Qu.:2.000
## Median :0.0000 Median :4.000 Median :2.000
## Mean :0.4062 Mean :3.688 Mean :2.812
## 3rd Qu.:1.0000 3rd Qu.:4.000 3rd Qu.:4.000
## Max. :1.0000 Max. :5.000 Max. :8.000
dim(mtcars)
## [1] 32 11
Regression Analysis Starting regression on mpg and am.
am_model <- lm(mpg ~ am,mtcars)
summary(am_model)
##
## Call:
## lm(formula = mpg ~ am, data = mtcars)
##
## Residuals:
## Min 1Q Median 3Q Max
## -9.3923 -3.0923 -0.2974 3.2439 9.5077
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 17.147 1.125 15.247 1.13e-15 ***
## am 7.245 1.764 4.106 0.000285 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 4.902 on 30 degrees of freedom
## Multiple R-squared: 0.3598, Adjusted R-squared: 0.3385
## F-statistic: 16.86 on 1 and 30 DF, p-value: 0.000285
Explaining more of variables in the data and variance
fullmodel <- lm(mpg ~ .,data =mtcars)
summary(fullmodel)
##
## Call:
## lm(formula = mpg ~ ., data = mtcars)
##
## Residuals:
## Min 1Q Median 3Q Max
## -3.4506 -1.6044 -0.1196 1.2193 4.6271
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 12.30337 18.71788 0.657 0.5181
## cyl -0.11144 1.04502 -0.107 0.9161
## disp 0.01334 0.01786 0.747 0.4635
## hp -0.02148 0.02177 -0.987 0.3350
## drat 0.78711 1.63537 0.481 0.6353
## wt -3.71530 1.89441 -1.961 0.0633 .
## qsec 0.82104 0.73084 1.123 0.2739
## vs 0.31776 2.10451 0.151 0.8814
## am 2.52023 2.05665 1.225 0.2340
## gear 0.65541 1.49326 0.439 0.6652
## carb -0.19942 0.82875 -0.241 0.8122
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 2.65 on 21 degrees of freedom
## Multiple R-squared: 0.869, Adjusted R-squared: 0.8066
## F-statistic: 13.93 on 10 and 21 DF, p-value: 3.793e-07
stepmodel<- step(fullmodel)
## Start: AIC=70.9
## mpg ~ cyl + disp + hp + drat + wt + qsec + vs + am + gear + carb
##
## Df Sum of Sq RSS AIC
## - cyl 1 0.0799 147.57 68.915
## - vs 1 0.1601 147.66 68.932
## - carb 1 0.4067 147.90 68.986
## - gear 1 1.3531 148.85 69.190
## - drat 1 1.6270 149.12 69.249
## - disp 1 3.9167 151.41 69.736
## - hp 1 6.8399 154.33 70.348
## - qsec 1 8.8641 156.36 70.765
## <none> 147.49 70.898
## - am 1 10.5467 158.04 71.108
## - wt 1 27.0144 174.51 74.280
##
## Step: AIC=68.92
## mpg ~ disp + hp + drat + wt + qsec + vs + am + gear + carb
##
## Df Sum of Sq RSS AIC
## - vs 1 0.2685 147.84 66.973
## - carb 1 0.5201 148.09 67.028
## - gear 1 1.8211 149.40 67.308
## - drat 1 1.9826 149.56 67.342
## - disp 1 3.9009 151.47 67.750
## - hp 1 7.3632 154.94 68.473
## <none> 147.57 68.915
## - qsec 1 10.0933 157.67 69.032
## - am 1 11.8359 159.41 69.384
## - wt 1 27.0280 174.60 72.297
##
## Step: AIC=66.97
## mpg ~ disp + hp + drat + wt + qsec + am + gear + carb
##
## Df Sum of Sq RSS AIC
## - carb 1 0.6855 148.53 65.121
## - gear 1 2.1437 149.99 65.434
## - drat 1 2.2139 150.06 65.449
## - disp 1 3.6467 151.49 65.753
## - hp 1 7.1060 154.95 66.475
## <none> 147.84 66.973
## - am 1 11.5694 159.41 67.384
## - qsec 1 15.6830 163.53 68.200
## - wt 1 27.3799 175.22 70.410
##
## Step: AIC=65.12
## mpg ~ disp + hp + drat + wt + qsec + am + gear
##
## Df Sum of Sq RSS AIC
## - gear 1 1.565 150.09 63.457
## - drat 1 1.932 150.46 63.535
## <none> 148.53 65.121
## - disp 1 10.110 158.64 65.229
## - am 1 12.323 160.85 65.672
## - hp 1 14.826 163.35 66.166
## - qsec 1 26.408 174.94 68.358
## - wt 1 69.127 217.66 75.350
##
## Step: AIC=63.46
## mpg ~ disp + hp + drat + wt + qsec + am
##
## Df Sum of Sq RSS AIC
## - drat 1 3.345 153.44 62.162
## - disp 1 8.545 158.64 63.229
## <none> 150.09 63.457
## - hp 1 13.285 163.38 64.171
## - am 1 20.036 170.13 65.466
## - qsec 1 25.574 175.67 66.491
## - wt 1 67.572 217.66 73.351
##
## Step: AIC=62.16
## mpg ~ disp + hp + wt + qsec + am
##
## Df Sum of Sq RSS AIC
## - disp 1 6.629 160.07 61.515
## <none> 153.44 62.162
## - hp 1 12.572 166.01 62.682
## - qsec 1 26.470 179.91 65.255
## - am 1 32.198 185.63 66.258
## - wt 1 69.043 222.48 72.051
##
## Step: AIC=61.52
## mpg ~ hp + wt + qsec + am
##
## Df Sum of Sq RSS AIC
## - hp 1 9.219 169.29 61.307
## <none> 160.07 61.515
## - qsec 1 20.225 180.29 63.323
## - am 1 25.993 186.06 64.331
## - wt 1 78.494 238.56 72.284
##
## Step: AIC=61.31
## mpg ~ wt + qsec + am
##
## Df Sum of Sq RSS AIC
## <none> 169.29 61.307
## - am 1 26.178 195.46 63.908
## - qsec 1 109.034 278.32 75.217
## - wt 1 183.347 352.63 82.790
summary(stepmodel)
##
## Call:
## lm(formula = mpg ~ wt + qsec + am, data = mtcars)
##
## Residuals:
## Min 1Q Median 3Q Max
## -3.4811 -1.5555 -0.7257 1.4110 4.6610
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 9.6178 6.9596 1.382 0.177915
## wt -3.9165 0.7112 -5.507 6.95e-06 ***
## qsec 1.2259 0.2887 4.247 0.000216 ***
## am 2.9358 1.4109 2.081 0.046716 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 2.459 on 28 degrees of freedom
## Multiple R-squared: 0.8497, Adjusted R-squared: 0.8336
## F-statistic: 52.75 on 3 and 28 DF, p-value: 1.21e-11
Inference
t.test(mpg ~am,mtcars)
##
## Welch Two Sample t-test
##
## data: mpg by am
## t = -3.7671, df = 18.332, p-value = 0.001374
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -11.280194 -3.209684
## sample estimates:
## mean in group 0 mean in group 1
## 17.14737 24.39231
Appendix and Plots
boxplot(mpg~am,data = mtcars,col =(c("Yellow","Red")),ylab = "Miles per Gallon",xlab = "Transmission Type")
Potting using base graphics
par(mfrow=c(2,2))
plot(stepmodel)
data(mtcars)
pairs(mtcars,mapping = aes(colour=am))
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy.coords(x, y), type = type, ...): "mapping" is not a
## graphical parameter
## Warning in plot.window(...): "mapping" is not a graphical parameter
## Warning in plot.xy(xy, type, ...): "mapping" is not a graphical parameter
## Warning in title(...): "mapping" is not a graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter
## Warning in axis(side = side, at = at, labels = labels, ...): "mapping" is not a
## graphical parameter