similar to: color graph in multiple plots

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")

Hi, I would like to overlap the cdf curve for observed and generated data  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))

Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. "z" is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE, xlim=c(1.6,3),ylim=c(1-sum(z>1.6)/length(z), 1)) x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4373,-0.2398), col="red") y

Hi r-users, I would like to add a plot of vertical line segment with arrow from (77,.6) to (77,0) and also a horizontal  line segment with arrow from (0,0.6) to (77,.6) .  So far this is what I have: plot(sq, cdf, type="l", lwd=4,col="blue",xaxs="i",yaxs="i", xlab= "Rainfall (mm)", ylab= "Random no.", main="Random number  and

Hi all, I would like to plot part of the emperical CDF. Suppose the variable is x, I just need the part when x>1,therefore, I am using the following codes. tail <- x>1 plot(ecdf(x[tail]), do.points=FALSE, verticals=TRUE) The "x" value starts from 1, but the yaxs still begins from 0, not the corresponding value when "x" is 1. How can I make it match? Could anyone

Hi r-users,   I have 2 sets of data and I would like to superimpose this cumulative density in one graph.   I know how to put the 2 graphs in one same graph but my problem is the data are different.   > z[1:20]  [1]  2.02347388  3.19514379  0.05188875  1.41333812  3.50249892  4.34272676  6.65639581  5.01533819  4.18207505 [10]  2.87615551  2.28202152  0.49431541  0.06570579  5.68352438

Hi r-users,   I would like to compare the cdf between historical and predicted. My x.obs and x.pre are the frequency data in classes of 0-300. I tried: plot(ecdf(x.obs),ecdf(x.pre),type="l",col="red")   and it gives me: Error in plot.stepfun(x, ..., ylab = ylab, verticals = verticals, pch = pch) :   argument 4 matches multiple formal arguments   Thank you so much for any

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=

hi, i am using par(mrow=c(6,6)) function to get 6x6 plots on one screen. the problem that i am having is that the axis tick labels are far away from the ticks and going into previous plots (see attached figure). i need to know how can i reduce the distance between the ticks and their values (y axis values).the part of the code that i am using (after reading in the data) to create the top row is:

Dear expeRts, I would like to color a certain region in a levelplot. The region for z <= 0.02 should have a dark gray color. Below is a minimal example. It almost does what I want, but The region between z=0.02 and z=1 is also colored in dark gray (instead of just the region for z <= 0.02). How can I solve this? Cheers, Marius ## z values for given x and y f <- function(x)

Hi everyone I want a square scatterplot with abline(0,1) going exactly through the SW and NE corners. By "square" I mean that the plotting region is exactly square, and that the axis limits are identical. x <- 1:20 y <- x+rep(c(-1,1),10) lims <- range(c(x,y)) None of the following do this: plot(x,y) ; abline(0,1) #not square plot(x,y,asp=1);abline(0,1) #diagonal

Full_Name: Fran?ois Birgand Version: 2.9.0 OS: windows xp Submission from: (NULL) (152.1.16.161) When I type this sequence: x11() Ylim<-c(0,14) Xlim<-c(0,14) plot(0,0,xlim=Xlim,ylim=Ylim,col="white",main="",xlab="",ylab="",bty="n",xaxt="n",yaxt="n",xaxs="i",yaxs="i",asp=1)

Dear UseRs, maybe is a silly question: how can I get Empirical CDF values from an object created with ecdf()?? Using print I obtain: Empirical CDF Call: ecdf(t) x[1:57] = 4.1, 4.4, 4.5, ..., 491.3, 671.27 Thanks in advance. Regards, Vito Diventare costruttori di soluzioni Became solutions' constructors "The business of the statistician is to catalyze the scientific

Dear R users, The parameters xaxs and yaxs (function par, package graphics) seem not to work with the function hist (package graphics), even when the parameters xlim and ylim are defined. Is there any way to make yaxs="i" and xaxs="i" work properly with the function hist, mainly to produce histograms that "touch" the horizontal axis? The R documentation and the R

Dear R community, I need to plot the results of some simulations I did using QTL Cartographer. I am plotting LOD scores over three chromosomes. The three plot have to be one next to the other. The procedure I am using is: par(mfrow=c(1,3)) plot(x$x, x$y, ylim=c(0,35), type="l", col="blue", las=1, xaxs="i", yaxs="i", xlab="X Chromosome",

Good Afternoon, I'm trying to create a cdf plot, with the following code. It works well, but I have little doubt, if you can help solve. When I create the plot, like the graph line would still not appear with point #cdf x<-table(Dataset$Apcode) View(s) hist(s) *plot(ecdf(x))* x<-1 37607 2 26625 3 5856 4 25992 5 30585 6 16064 7 9850 .. ... .. 186 52 -- View this message in

Hi I used R for the first time yesterday. I wanted to plot the aliasing effect of sampling a 5.5KHz sinusoid at only 8KHz (below the Nyquist limit). So I wrote a small R script that a) plots 1msec worth of a 5.5KHz sin wave b) plots 1msec of the resulting 2.5KHz alias and c) plots the 8 sampling points on the 5.5KHz source wave. I think I have found a bug. The script is as follows:

Dear all r-users,   I have this code below to draw two squares, small and big square.  I would like to colour the small square with red and the big square with blue for example.  I tried using polygon but fail.  Thank you so much for your help.     par(mar=c(4,4,2,1.2),oma=c(0,0,0,0),xaxs="i", yaxs="i") plot.new() plot.window(xlim= c(0,8),

Hi, I have a vectors x and z, for example, x <- 0:20 z <- round(runif(20,1,7)) y <- 0.5 and I want to display z as an image. However if I then call image() with a vector image(x,y,t(z),zlim=c(1,7),col=heat.colors(7),xlab="Year",ylab="Action",yaxt="n",xaxs="r",yaxs="r") then I get the error Error in image.default(x, y, t(z), zlim =

Full_Name: Steve Ellison Version: 2.4.1 OS: Windows, Linux Submission from: (NULL) (194.73.101.157) bxp() allows specifcation of box locations with at=, but neither adjusts xlim= to fit at nor does it respect xlim provided explicitly. This is because bxp() now includes explicit xlim as c(0.5, n+0.5), without checking for explicitly supplied xlim (or ylim if horizontal). This also prevents