[R] Q about package Icens: change the color of the shading in plot function

[Mendolia, Franco]

Hello!

I want to use the Icens package for analyzing interval-censored data. This code from the manual gives me what I want.

library(Icens)
data(cosmesis)
csub1 <- subset(cosmesis, subset=Trt==0, select=c(L,R))
e1 <- VEM(csub1)
plot(e1)

However, I would like to change the color of the shading from green to something less green, say gray. Any ideas how I could do that? I looked at par, but I wasn't able to find what I need.

Thanks,
Franco