search for: datobs

Hi r-users,   I would like to overlap a smooth line on the histogram.  I tried using spline but it does not work. sq     <- seq(0,900,by=50) sq.50  <- as.character(sq) datobs <- sum_pos ## first, plot histogram histo <- hist(datobs,breaks=sq,freq=F) ## extract counts from histogram and calculate the probability in ## each specified interval. Also check the length of the interval is ## the same with the length of the fitted model. hct   <- c(0,histo$counts)/l...

Hi r-users, I would like to overlap 2 ecdf plots.  I tried this below and it gives me two plots of ecdf but just both just in black. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs)) lines(ecdf(gam_sum_gen)) Then I try to add colors etc and also the legend but fail. par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot(ecdf(datobs),main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", pch = pch,vert...

...oes anybody encounter the same problem when we overlap histogram and density     that the density line seem to shift to the right a little bit?           If you do have the same problem, what should we do to correct that?           Thank you.           par(mar=c(4,4,2,1.2),oma=c(0,0,0,0))     hist(datobs,prob=TRUE, main ="Volume of a catchment from four     stations",col="yellowgreen", cex.axis=1,     xlab="rainfall",ylab="Relative frequency", ylim= c(0,.003), xlim=c(0,1200))     lines(density(dd), lwd=3,col="red")     #legend("topright"...

...ata  Here is my code: plot(cdf,main ="CDF of the sum for winter season-Hume",cex.axis=1.2,xlab="Rainfall (mm)", xaxs="i",yaxs="i",col=c("black","red"), lty=c(1,1),ylab="Cumulative probability", xlim=c(0,800),lwd=1) lines(ecdf(datobs)) legend("topright", legend = c("observed","fitted"),        col = c("black","red"), pch=c(NA,NA), lty = c(1, 1),        lwd=c(3,3),bty="n", pt.cex=2) It gives me the plots but it is not smooth.  How do I adjust so that I will get two sm...

Hi r-users, I would like my axes to intersect at (0,0).  I tried  xaxs="i",yaxs="i" but it does not change anything.  I hope anybody can help me with this problem.  Here is my code. hist(datobs, prob=TRUE, main ="PDF of the sum of two stations",col="yellowgreen", cex.axis=1.2, xlab="Rainfall (mm)", ylab="Relative frequency", ylim= c(0,.008), xlim=c(0,600),xaxs="i",yaxs="i") lines(density(dd), lwd=3,col="red") legend(&...