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 <drjimle...@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.panco...@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 <drjimle...@gmail.com> ha >> scritto: >> >> Hi Francesca, >> Do you want something like this? >> >> Jim >> >> On Thu, Dec 5, 2019 at 6:58 PM Francesca <francesca.panco...@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]] ______________________________________________ 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.