# best seen with: x11(7.5, 4.7)
# just one variable with marginal densities
PlotMarDens( y=d.pizza$temperature, x=d.pizza$delivery_min, grp=1
, xlab="delivery_min", ylab="temperature", col=SetAlpha("brown", 0.4)
, pch=15, lwd=3
, panel.first= grid(), args.legend=NA
, main="Temp ~ delivery"
)
# use a group variable
PlotMarDens( y=d.pizza$temperature, x=d.pizza$delivery_min, grp=d.pizza$area
, xlab="delivery_min", ylab="temperature", col=c("brown","orange","lightsteelblue")
, panel.first=list( grid() )
, main = "temperature ~ delivery_min | area"
)
# reset layout
par(mfrow=c(1,1))
Run the code above in your browser using DataLab