Nothing
## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.width=5, fig.height=5 ,fig.align="center"
)
fpath <- ""
## -----------------------------------------------------------------------------
ozone <- na.omit(airquality)
## -----------------------------------------------------------------------------
fit1 <- loess(Ozone~Wind, data=ozone)
plot(Ozone~Wind, data=ozone, xlim=c(1.5, 21), ylim=c(-5,175))
wind <- seq(min(ozone$Wind), max(ozone$Wind), length.out=30)
lines(wind, predict(fit1, data.frame(Wind=wind)))
## ----fig.width=7.3, fig.height=2.7--------------------------------------------
fit2 <- loess(Ozone~Wind+Solar.R, data=ozone)
par(mfrow=c(1,3))
par(mar=c(5,5,3,1))
for (s in quantile(ozone$Solar.R, c(.25,.5,.75))){
plot(Ozone~Wind, data=subset(ozone, Solar.R <= s+20 & Solar.R >= s-20),
xlim=c(1.5, 21),ylim=c(-5,175), main=paste0("Solar.R=",s))
lines(wind, predict(fit2, data.frame(Wind=wind, Solar.R= s)))
}
## ----eval=F-------------------------------------------------------------------
# suppressMessages(library(condvis2))
# condvis(ozone, fit2,sectionvars="Wind", conditionvars="Solar.R")
## ----echo=FALSE, out.width='100%'---------------------------------------------
knitr::include_graphics(paste0(fpath, "ozone.png"))
## -----------------------------------------------------------------------------
fit3 <- loess(Ozone~Wind+Solar.R+Temp, data=ozone)
## ----eval=F-------------------------------------------------------------------
# condvis(ozone, fit3, sectionvars="Wind", conditionvars=c("Solar.R", "Temp"))
## ----echo=FALSE, out.width='100%'---------------------------------------------
knitr::include_graphics(paste0(fpath, "ozone1.png"))
## ----echo=FALSE, out.width='100%'---------------------------------------------
knitr::include_graphics(paste0(fpath, "ozone4.png"))
## ----echo=FALSE, out.width='100%'---------------------------------------------
knitr::include_graphics(paste0(fpath, "ozone3.png"))
## ----eval=F-------------------------------------------------------------------
# library(e1071)
# fit4 <- svm(Ozone~Wind+Solar.R+Temp, data=ozone)
# condvis(ozone, list(loess=fit3,svm=fit4), sectionvars="Wind", conditionvars=c("Solar.R", "Temp"))
## ----echo=FALSE, out.width='100%'---------------------------------------------
knitr::include_graphics(paste0(fpath, "ozone5.png"))
## ----eval=F-------------------------------------------------------------------
# fit5 <- lm(Ozone~Wind+Solar.R+Temp, data=ozone)
# condvis(ozone, fit5,
# sectionvars="Wind", conditionvars=c("Solar.R", "Temp"), predictArgs=list(list(pinterval="confidence")))
#
## ----eval=F-------------------------------------------------------------------
# fitu <- fit3
# class(fitu)<- c("upper", class(fitu))
#
# CVpredict.upper <- function(f, newdata, ...){
# p <- predict(f, newdata, se=T)
# p$fit+ 2*p$se.fit
# }
#
# fitl <- fit3
# class(fitl)<- c("lower", class(fitu))
#
# CVpredict.lower <- function(f, newdata, ...){
# p <- predict(f, newdata, se=T)
# p$fit- 2*p$se.fit
# }
#
# condvis(ozone, list(loess=fit3,lower=fitl,upper=fitu),
# sectionvars="Wind", conditionvars=c("Solar.R", "Temp"),
# linecols=c("red", "blue","blue"))
## ----echo=FALSE, out.width='100%'---------------------------------------------
knitr::include_graphics(paste0(fpath, "ozone6.png"))
## ----eval=F-------------------------------------------------------------------
# library(ks)
# data(iris)
# irisf <- kde(x=iris[,1:3])
#
# condvis(data = iris, model = list(kde=irisf),
# sectionvars= c("Sepal.Length", "Sepal.Width"),
# conditionvars= "Petal.Length", density=T)
#
## ----echo=FALSE, out.width='100%'---------------------------------------------
knitr::include_graphics(paste0(fpath, "iris.png"))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.