Вам не нужно указывать значение $ CONSTANT <- 0, это всего лишь 1 боксплот. Тильда (~) используется для указания разделения. Например, если ваши значения RAIN равны 0/1, вы можете сделать: </p>
dat = data.frame(RAIN=sample(0:1,200,replace=TRUE),
hour = sample(1:12,200,replace=TRUE),
PM2.5=runif(200))
boxplot(PM2.5~ RAIN > 0,data=subset(dat,hour == 12),names=c("RAIN==0","RAIN>0"),xlab="")
Если есть значения, отличные от 0/1 для ДОЖДЬ, делай:
dat = data.frame(RAIN=sample(0:5,200,replace=TRUE),
hour = sample(1:12,200,replace=TRUE),
PM2.5=runif(200))
sub_dat = list(
"RAIN==0"=subset(dat, RAIN == 0 & hour == 12)$PM2.5,
"RAIN>0"=subset(dat, RAIN > 0 & hour == 12)$PM2.5
)
boxplot(sub_dat)