This is a consolation, because I cannot get it in ggplot either!
Thanks for the code!
F.--------------------------------
Il giorno 5 dic 2019, alle ore 11:17, Jim Lemon <drjimlemon at gmail.com> ha scritto:
Sorry it's not ggplot, I couldn't work that one out.
# using the data frame structure that Rui kindly added
# and perhaps Rui can work out how to do this in ggplot
DBcomplete<-DB[complete.cases(DB),]
library(plotrix)
png("fp.png")
par(mfrow=c(1,2))
density11_0<-density(DBcomplete$delta11_L[DBcomplete$participation1==0])
density2_0<-density(DBcomplete$delta2_L[DBcomplete$participation1==0])
plot(0,xlim=c(-30,50),ylim=c(0,max(density11_0$y)),type="n",
xlab="delta",ylab="density",main="participation == 0")
plot_bg("lightgray")
grid(col="white")
polygon(density11_0,col="#ff773344")
polygon(density2_0,col="#3377ff44")
density11_1<-density(DBcomplete$delta11_L[DBcomplete$participation1==1])
density2_1<-density(DBcomplete$delta2_L[DBcomplete$participation1==1])
plot(0,xlim=c(-30,50),ylim=c(0,max(density11_1$y)),type="n",
xlab="delta",ylab="density",main="participation == 1")
plot_bg("lightgray")
grid(col="white")
polygon(density11_1,col="#ff773344")
polygon(density2_1,col="#3377ff44")
par(cex=0.9)
legend(5,0.11,c("delta11_L","delta2_L"),fill=c("#ff773344","#3377ff44"))
dev.off()
Jim
On Thu, Dec 5, 2019 at 9:14 PM Francesca <francesca.pancotto at gmail.com> wrote:
Exactly. I was trying to remelt data in the right way, but I could not get there yet. Can you suggest me this code?
Thanks a lot
F.
----------------------------------
Il giorno 5 dic 2019, alle ore 11:11, Jim Lemon <drjimlemon at gmail.com> ha scritto:
Hi Francesca,
Do you want something like this?
Jim
On Thu, Dec 5, 2019 at 6:58 PM Francesca <francesca.pancotto at gmail.com> wrote:
Hi, sorry for bothering again.
I was wondering how I can reshape the data, if in your code,
I would like to have only two panels, where in the panel with Participation =0, I represent delta11_L of participation1==0
and delta2_L of participation2==0, and in the right panel, I want Participation=1, but representing together
delta11_L of participation1==1, and delta2_L of participation2==1.
I get messed up with the joint melting of participation, which determines the facet, but then I cannot assign the proper fill to the density plots which depend on it, and on the other hand I would like to have in the same plot with mixed participation.
I hope it is clear.
Nonetheless, the previous plot is useful to understand something I had not thought about.
Thanks again for your time.
F.
----------------------------------
<fp.png>
[[alternative HTML version deleted]]