similar to: Labelling certain points on the x-axis

Hi, I've got the following edf: *** x = c(1.6,1.8,2.4,2.7,2.9,3.3,3.4,3.4,4,5.2) F2.5 <- ecdf(x) plot(F2.5, verticals= TRUE, do.p = TRUE, lwd=3, ylab = "", xlab = "", xlim = c(1,5.5)) abline(h= (0:5)*0.2) #mean abline(v=mean(x), lwd=2) mtext(text=expression(bar(x) == 3.07), side=1, adj=0.462, padj=3, cex=1) *** Now I would like to

Hello, I have the following problem: I created an ecdf and a barplot. Unfortunatly, the bars are not where I would like them to be (please see picture below). http://www.nabble.com/file/p12877530/problem.gif That's my code: #------------------------ par(mfrow=c(2,1), mar=c(2,3,3,2)) #ECDF x = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2) F2.5 <- ecdf(x) plot(F2.5,

Hi list members I have a vector of strings x=c("w","ex","ee") And I want to get a logical vector showing the positions where my search string "e" matches the elements partially, i.e. is at least the left-hand part of the target strings, i.e. I want to get a vector FALSE TRUE TRUE. Any ideas? Thanks Steve Powell proMENTE social research research |

Hi, I have just started using R. Now I have the following problem: I want to create an Empirical Cumulative Distribution Function and I only came so far: F10 <- ecdf(x) plot(F10, verticals= TRUE, do.p = TRUE, lwd=3) x=c(1.6,1.8,2.4,2.7,2.9,3.3,3.4,3.4,4,5.2) Now I'd like to use arguments such as xlabs and main but I don't know how to integrate them. I hope someone can help me, I am

Hello, Does anyone know if it's possible to change the text direction of mtext? I'd like to use 'mtext(side=2)' but the text should be 'left-to-right'. I couldn't find an argument for that in the help files. Thanks for any suggestions Tobias -- View this message in context: http://www.nabble.com/changing-text-direction-of-mtext-tf4321931.html#a12307342 Sent from

I have the following code: mtext(text="X[1]") How can I change "[1]" to be subscript? (see picture below) http://www.nabble.com/file/p12241351/subscript.gif subscript.gif Thanks for any help Tobias -- View this message in context: http://www.nabble.com/Subscript-tf4300665.html#a12241351 Sent from the R help mailing list archive at Nabble.com.

Hello, I'm trying to create the following graph: http://www.nabble.com/file/p12974686/bar.gif Is there a command (like pie(...) or barplot(...)) or is this just a special case of eg barplot? I'd be glad if someone could give me a hint on where to start. Regards, Tobias -- View this message in context: http://www.nabble.com/how-to-create-this-graph--tf4546742.html#a12974686 Sent from

Hi, I have an error happened when I use exp() in my Fortran code. My Fortran code is SUBROUTINE f(n,c) REAL, INTENT(IN) :: n(10) REAL, INTENT(OUT) :: c(10) INTEGER :: k DO k=1, 10 c(k) = exp( n(k) ) END DO END I compiled the fortran code as a dll Rcmd SHLIB -o f.dll f.f -L -lacml Then I loaded the dll file to R, and try some number >

Hi. I have the following two ecdf plots in one graph: plot( ecdf(....), do.points=FALSE, verticals=TRUE, main=paste("Ecdf of distances ",DIM,sep=""), col="red" ); lines( ecdf(....), do.points=FALSE, verticals=TRUE ); How do I change the color of the resulting graph? Adding col="red" to either plot or lines results in an error

Hello, I wanted to create a histogram, but somehow I got stuck... The interval limits are: x = 1, 2, 3, 3.5, 4.5, 5, 5.5 The interval widths are therefore: 1, 1, 0.5, 1, 0.5, 0.5 Nothing I tried worked... Can anyone help me please? Thanks Tobias -- View this message in context: http://www.nabble.com/Histogram-tf4315900.html#a12288850 Sent from the R help mailing list archive at Nabble.com.

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

http://r.789695.n4.nabble.com/file/n3567636/ecdfs.jpg ecdfs.jpg http://r.789695.n4.nabble.com/file/n3567636/ecdf_curve.gif ecdf_curve.gif Hello, I have generated a plot of two empirical CDFs (attachment 1). As a result, they are stepwise when plotted. The following code was used: > plot(ecdf(mut), do.points=FALSE, verticals=TRUE, xlim=range(mut, non), > col="red") >

Hi all, given this example #start a<-c(0,70,50,100,70,650,1300,6900,1780,4930,1120,700,190,940, 760,100,300,36270,5610,249680,1760,4040,164890,17230,75140,1870,22380,5890,2430) length(a) b<-c(0,0,10,30,50,440,1000,140,70,90,60,60,20,90,180,30,90, 3220,490,20790,290,740,5350,940,3910,0,640,850,260) length(b) out<-ks.test(log10(a+1),log10(b+1)) # max distance D

Hello, if I plot a horizontal line, e.g., plot(c(1,2),c(1,1),xlim=c(0,3),lwd=2,type="l") or plot(c(1,2),c(1,1),xlim=c(0,3),lwd=4,type="l") then the left end (1st example) or both ends (2nd example) of the lines are not rectangular but slanted on the graphical display (screen). That behavour first occurred when I was trying to plot a stepfun, e.g., y <-

Dear experts, when trying to estimate an kernel density function with density(x) I get the following error message with imported data from either EXCEL or text files: Error in density(spr) : argument must be numeric. Other procedues such as truehist work. If I generate data within R density works fine. Does anybody have an idea? Yours -- Christoph Hanck Wissenschaftliche Hilfskraft

Hi I am conducting 2-sample Kolmogorov Smirnov tests for my Masters project to determine if two independant tree populations have the same size-class distribution or not. The trees have been placed into size-class categories based on their basal diameters. Once I started running the stats on my data, I got confused with the results. Just to show an example of what I was testing I ran stats

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 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 everyone, I have what may appear to be a newbie question, but I have looked everywhere I can think to look and I cannot find an answer. On page 35 of "An Introduction to R" the following command appears: plot(ecdf(eruptions), do.points=FALSE, verticals=TRUE). What is the do.points argument? I know what it does (suppresses printing of the points) but where can I find help on it?

Dear R users/experts, I've heard it's possible to estimate the above kinds of models in R. However, after (an admittedly brief) survey of the packages, I haven't found an obvious candidate. Can you offer any help? (Yes I'm a newby.) Yours sincerely Christoph Hanck Studentische Hilkskraft Lehrstuhl f?r Empirische Wirtschaftsforschung, Prof. Dr. Wilfling