## I recommend using lattice for this task.
## First I show the example from my book and package (HH).
## Then I use this on your example.
library(HH) ## Package supporting Heiberger and Holland,
## Statistical Analysis and Data Display (Second
edition, 2015)
HHscriptnames(4) ## Filename on your computer for script for all
Chapter 4 examples
## this is Chunk 23
###################################################
### code chunk number 23: grap.tex:1953-1981
###################################################
bwdata <- data.frame(Y=(rt(80, df=5)*5 + rep(c(20,25,15,22, 22,28,16,14), 10)),
week=ordered(rep(c(1:4, 1:4), each=10)),
treatment= rep(c("A", "B"), each=40))
position(bwdata$week) <- c(1, 2, 4, 8)
levels(bwdata$week) <- c(1, 2, 4, 8)
bwdata$week.treatment <- with(bwdata, interaction(treatment, week))
position(bwdata$week.treatment) <-
as.vector(t(outer(c(1, 2, 4, 8), c(-.18,.18), "+")))
BR <- likertColor(2, colorFunctionOption="default")[2:1]
## uses panel.bwplot.intermediate.hh to control position and colors
## hhpdf("bwplotposition.pdf", width=7, height=5)
bwplot(Y ~ week.treatment, data=bwdata,
panel=panel.bwplot.intermediate.hh, xlim=c(0, 9),
box.width=.25,
pch=c(17, 16), col=BR,
xlab="Week", ylab=list(rot=0),
scales=list(x=list(at=position(bwdata$week), tck=1)),
key=list(
text=list(c("A","B"), col=BR),
points=list(pch=c(17, 16), col=BR),
space="top", columns=2, between=1, border=TRUE,
title="Treatment", cex.title=.9)) +
layer(panel.abline(h=37, col="gray60", lty=3, lwd=2))
## hhdev.off()
## The placement features provided by panel.bwplot.intermediate.hh are
## 1. The two boxes at each time position are clearly distinguished
## from boxes at other time positions.
##
## 2. Times do not need to be evenly spaced.
## Now your sample data and lattice code for your desired graph
## Script to demonstrate what I am trying to do.
##
## Simulate some data:
Year <- factor(rep(4:8,each=50,times=2))
Type <- rep(c("National","Local"),each=250)
M0 <- 1300+50*(0:4)
set.seed(42)
M1 <- M0 + runif(5,-100,-50)
X0 <- rnorm(250,rep(M0,each=50),150)
X1 <- rnorm(250,rep(M1,each=50),100)
DemoDat <- data.frame(Year=Year,Score=c(X0,X1),Type=Type)
if (FALSE) { ## Rolf Turners original code
## Grouped boxplots:
library(ggplot2)
print(ggplot(data=DemoDat) +
geom_boxplot(aes(x=Year, y=Score, color=Type,fill=Type),
position=position_dodge(1),alpha=0) +
theme_minimal() +
scale_fill_discrete(name="National v. Local") +
ylim(700,2100))
}
DemoDat$Year.Type <- with(DemoDat, interaction(Year, Type))
position(DemoDat$Year.Type) <-
as.vector(t(outer(c(4, 5, 6, 7, 8), c(-.18, .18), "+")))
RB <- c("red", "black")
SYT <-
bwplot(Score ~ Year.Type, data=DemoDat,
panel=panel.bwplot.intermediate.hh,
xlim=c(-.1, 9.1),
ylim=c(690, 2110),
box.width=.22,
col=RB,
xlab="Year", ylab=list(rot=0),
scales=list(x=list(at=0:9, tck=1),
y=list(at=seq(700, 2100, 100), tick=1)),
par.settings=list(box.dot=list(pch="|"),
plot.symbol=list(pch="-", col=RB, cex=1.5)),
key=list(
text=list(levels(DemoDat$Type), col=RB, cex=.8),
lines=list(col=RB), size=1.5,
columns=2, between=.5, between.columns=.6,
space="right", border=FALSE,
title="\nNational v. Local", cex.title=.9),
main="Matches specifications"
)
SYT
update(SYT, main="Outliers made invisible, not recommended",
par.settings=list(plot.symbol=list(cex=0)))
## Rich
On Sat, Jul 28, 2018 at 7:04 PM, Rolf Turner <r.tur...@auckland.ac.nz> wrote:
>
> On 29/07/18 02:54, Jeff Newmiller wrote:
>
>> 1) I don't know... it looks to me like you did not run my code.
>
>
> Aaaarrrgghhh. I *thought* I had, but instead left "fill=Type" inside the
> aes() call and neglected to add fill=NULL outside this call. Duhhhh!!! It's
> tough being mentally challenged, let me assure you.
>
>> I have included a complete reprex below... try it out in a fresh session.
>> If you still get the problem, check your sessionInfo package versions
>> against mine.
>
>
> Yep. Works like a charm.
>
>> 2) This still smells like your fill parameter is inside the aes function
>> with Type as value. This causes a legend to be created, and since that
>> legend has a different name ("Type") than the colour scale, they are
>> separated. Confirm that you are using fill outside the aes function (because
>> you don't want fill to depend on the data) and have the constant NULL as
>> value (so it won't generate any fill graphical representation).
>
>
> Yeah. Well. Duhhh. I'm a retread.
>>
>>
>> 3) I missed that... the ylim()/scales_y_continuous(breaks=) limits
>> constrain which data are included as input into the graph. The
>> coord_cartesian function forces the limits as desired.
>
>
> Bewdy, ta.
>
>>
>> 4) While showing outliers is a standard semantic feature of boxplots
>> whether produced by ggplot or lattice or base or non-R solution,
>
>
> Indeed. But the client is always right! :-)
>
>> you can please the client by making the outliers transparent.
>
>
> And your code shows me how! Which I need. Bewdy, ta.
>
>
>> There is a link to the generated image below.
>>
>> ################
>> # Simulate some data:
>> Type <- rep( c( "National", "Local" ), each = 250 )
>> M0 <- 1300+50*(0:4)
>> set.seed( 42 )
>> M1 <- M0 + runif( 5, -100, -50 )
>> X0 <- rnorm( 250, rep( M0, each = 50 ), 150 )
>> X1 <- rnorm( 250, rep( M1, each = 50 ), 100 )
>>
>> library(ggplot2)
>> Year <- factor( rep( 4:8, each = 50, times = 2)
>> , levels = 0:8 )
>> DemoDat <- data.frame( Year = Year
>> , Score = c( X0, X1 )
>> , Type = Type
>> )
>>
>> ggplot( data = DemoDat
>> , aes( x = Year
>> , y = Score
>> , color = Type
>> )
>> , fill = NULL
>> ) +
>> geom_boxplot( position = position_dodge( 1 )
>> , outlier.alpha = 0
>> ) +
>> theme_minimal() +
>> scale_colour_manual( name = "National v. Local"
>> , values = c( "red", "black" ) ) +
>> scale_x_discrete( drop = FALSE ) +
>> scale_y_continuous( breaks=seq( 700, 2100, 100 ) ) +
>> coord_cartesian( ylim = c( 700, 2100 ) )
>>
>> # 
>
>
> Looks perfect. Thanks *HUGELY* for your patience with my stupidity.
>
> <SNIP>
>
> cheers,
>
> Rolf
>
> --
> Technical Editor ANZJS
> Department of Statistics
> University of Auckland
> Phone: +64-9-373-7599 ext. 88276
>
> ______________________________________________
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
