Hi R-mailing list I would have the following set-up below with a simplified data-frame. Through a loop which includes certain criteria for the densities I would like to plot the different density-distributions in the same plot. Of course I hope I don't do any mistakes with all the indexes of the dataframe. All I would like to have is the different densities in the same plot with a general background density d.
In a second step there might be also the need to combine all the different densities within the loop to one single distribution. I would be grateful for any help on this. Cheers > a [1] 1 4 7 11 16 21 23 > b [1] 15 17 13 23 14 11 12 > c [1] 0 15 0 0 0 0 0 > d [1] 0 0 7 0 16 0 0 > e [1] 0 0 0 11 0 21 0 > f [1] 0 0 123 124 122 121 0 q<- data.frame(a,b,c,d,e,f) a b c d e f 1 1 15 0 0 0 0 2 4 17 15 0 0 0 3 7 13 0 7 0 123 4 11 23 0 0 11 124 5 16 14 0 16 0 122 6 21 11 0 0 21 121 7 23 12 0 0 0 0 indx<-which(q[,"a"]==q[,"d"]) indy<-which(q[,"a"]==q[,"e"]) s<-[c(indx,indy),] std<-sd(S[,"f"]) t<-round(max(q[,"a"]),0) u<-round(min(q[,"a"]),0) colnames(v)<-c("a") D<-q[,"a"] d<-density(D) pdf('L') for (h in length(s[["a"]])) { for (i in length(S[["a"]])) { inda<-which(q[h,"a"]==q[h,"d"]) indb<-which(q[i,"a"]==q[i,"e"]) R<-(s[c(inda),"f"]+s[c(indb),"f"])/2 indh<-which(abs(q[,"b"]-R)<=(2*std)) indi<-which(abs(q[,"c"]-R)<=(2*std)) A<-q[indh,c("a","b")] colnames(A)<- c("a","g") B<-q[indy,c("a","c")] colnames(B)<- c("a","g") C<-rbind(A,B) c<-density(C[,"a"]) plot(d, col="blue",xlim=c(u,t),ylim=c(0,0.5)) polygon(d, col="blue",border="blue") lines(c, col="red") } } dev.off() -- Niklaus Hürlimann Université de Lausanne Institut de Minéralogie et Géochimie L'Anthropole CH-1015 Lausanne Suisse E-mail: niklaus.hurlim...@unil.ch Tel:+41(0)21 692 4452 [[alternative HTML version deleted]]
______________________________________________ R-help@r-project.org mailing list 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.