Thanks a lot, John Fox.. Its working Dileepkumar R Postdoctoral Researcher, Interdisciplinary Centre for Water Research (ICWaR), Indian Institute of Sciences, Bangalore, 560012 Mob: 9953035600, skype id: dileepkunjaai
On Mon, Jul 27, 2020 at 10:27 PM John Fox <j...@mcmaster.ca> wrote: > Dear Dileepkumar R, > > As is obvious from the tick marks, the vertical axis is not log-scaled: > > > log10(99.999) - log10(99.99) > [1] 3.908865e-05 > > log10(99) - log10(90) > [1] 0.04139269 > > > That is, these (approximately?) equally spaced ticks aren't equally > spaced on the log scale. > > The axis is instead apparently (at least approximately) on the logit > (log-odds) scale: > > > library(car) > Loading required package: carData > > logit(99.999) - logit(99.99) > [1] 2.302675 > > logit(99) - logit(90) > [1] 2.397895 > > > You can get a graph close to the one you shared via the following: > > library(car) # repeated so you don't omit it > > logits <- logit(y_values) > > plot(x_value, logits, log="x", axes=FALSE, > + xlim=c(1, 200), ylim=logit(c(10, 99.999)), > + xlab="Precipitation Intensity (mm/d)", > + ylab="Cumulative Probability", > + main="Daily U.S. Precipitation", > + col="magenta") > > axis(1, at=c(1, 2, 5, 10, 20, 50, 100, 200)) > > probabilityAxis(side=2, at=c(10, 30, 50, 90, 99, 99.9, 99.99, > 99.999)/100) > > box() > > Comments: > > This produces probabilities, not percents, on the vertical axis, which > conforms to what the axis label says. Also, the ticks in the R version > point out rather than into the plotting region -- the former is > generally considered better practice. Finally, the graph is not a > histogram as the original title states. > > I hope this helps, > John > > > -------------------------------------------- > John Fox > Professor Emeritus > McMaster University > Hamilton, Ontario, Canada > web: https://socialsciences.mcmaster.ca/jfox/ > > On 7/27/2020 11:56 AM, Dileepkumar R wrote: > > I think the attached sample figure is not visible > > Here is the sample figure: > > > https://drive.google.com/file/d/16Uy3JD0wsEucUv_KOhXCxLZ4U-3wiBTs/view?usp=sharing > > > > sincerely, > > > > > > Dileepkumar R > > > > > > > > > > On Mon, Jul 27, 2020 at 7:13 PM Dileepkumar R <dileepkunj...@gmail.com> > > wrote: > > > >> Dear All, > >> > >> I want to plot a simple cumulative probability distribution graph with > >> like the attached screenshot. > >> But I couldn't fix the y-axis scale as in that screenshot. > >> > >> My data details are follows: > >> > >> y_values > >> > =c(66.78149,76.10846,81.65518,85.06448,87.61703,89.61314,91.20297,92.36884, > >> 93.64070,94.57693,95.23052,95.75163,96.15792,96.58188,96.97933,97.29730, > >> 97.59760,97.91556,98.14520,98.37485,98.57799,98.74580,98.87829,99.06377, > >> 99.16093,99.25808,99.37290,99.45239,99.54072,99.59371,99.62904,99.66437, > >> 99.69970,99.70853,99.72620,99.73503,99.77036,99.79686,99.80569,99.82335, > >> 99.83219,99.84985,99.86751,99.87635,99.87635,99.90284,99.90284,99.90284, > >> 99.91168,99.92051,99.92051,99.93817,99.93817,99.93817,99.95584,99.95584, > >> 99.97350,99.97350,99.97350,99.97350,99.97350,99.97350,99.97350) > >> > >> x_value=seq(63) > >> > >> Thank you all in advance > >> > >> Dileepkumar R > >> > > > > [[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. > > > [[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.