similar to: How to label certain values on a curve?

Hi All, I'm trying to figure out the cumulative incidence curve in R in some limited time. I found in package "cmprsk", the command "plot.cuminc" can get this curve. But I noticed that there is no mark for the censored time there, comparing with the KM curve by "plot.survfit". Here are my codes (attached is the data): ----------------

For a given random variable rv, for instance, rv = rnorm(1000), I plot its density curve and calculate some quantiles: plot(density(rv)) P10P50P90 = = quantile(rv,probs = c(10,50,90)/100) I would like to color the area between P10 and P90 and under the curve and mark the P50 on the curve. > rv = rnorm(1000) > plot(density(rv)) > P10P50P90 = = quantile(rv,probs = c(10,50,90)/100) Could

I am attempting to wrap the histogram function in my own custom function, so that I can quickly generate some standard plots. A part of what I want to do is to draw a normal curve over the histogram: > x <- rnorm(1000) > hist(x, freq=F) > curve(dnorm(x), lty=3, add=T) (for normal use, x would be a vector of empirical values, but the rnorm() function works for testing) That

Hi, Excuse me for asking this silly question. But I really couldn't understand why cov() and ccov() don't work for my calculation of covariance matrix. a <- matrix(1:8, 2, 4) a [,1] [,2] [,3] [,4] [1,] 1 3 5 7 [2,] 2 4 6 8 > ccov(a) Error in solve.default(cov, ...) : Lapack routine dgesv: system is exactly singular I also tried colume bind, but it

Dear useRs, my question could be very basic for which i apologize in advance. Each column of a matrix with dimensions 365 rows and 37 columns was drawn against another matrix of dimensions 365 rows and 1 column. with that i was able to draw 37 curves on the same axis. now i want to draw an average curve of these 37 curves on the same axis in such a way that all the curves (average and 37

Hello, ? What I really want to do is to add a rejection region in the form of a long rectangle to a density plot I have drawn.? I am getting? 2 plots.? How can I add rectangle to first plot?? see code below. First section works fine.? It just is not quite what I want. # NORMAL DISTRIBUTION PLOT OF RAW DATA WITH UPPER CRITICAL LEVEL - ok xcrit=144.1 # *** single-sample Upper one-tailed

Hi, I wonder if anyone could advise me with this: I've been trying to make a standard curve in R with lm() of some standards from a spectrophotometer, so as I can express the curve as a formula, and so obtain values from my treated samples by plugging in readings into the formula, instead of trying to judge things by eye, with a curve drawn by hand. It is a curve and so I used the

Hello, Can anyone suggest a simple way to generate a Kaplan-Meier plot with 2 survfit objects, just like this one:? https://drive.google.com/file/d/1fEcpdIdE2xYtA6LBQN9ck3JkL6-goabX/view?usp=sharing Suppose I have 2 survfit objects: fit1 is for the curve on the left (survtime has been truncated to the cutoff line: year 5), fit2 is for the curve on the right (minimum survival time is at the

Hi, I have a problem about graphics. I would like to plot two graphs: a barplot and curve. Here is the code: > barplot(dpois(0:45,20),xlim=c(0,45),names=0:45) > curve(dnorm(x,20,sqrt(20)),from=0,to=45,add=T) Both graphs are drawn in the same figure, however the scale in both graphs dooes not match. For some reason the second plot is shifted to left. I think there is a problem

Hello: Using the built-in dataset aml as an example: data(aml) If I use instead dummy variables: aml$x1 = (aml$x=="maintained")aml$x2 = (aml$x=="unmaintained") and I want to plot the survival curve using x1, x2, and I just want the 2 levels, rather than 4 curves from: fit <- survfit(Surv(time, status) ~ x1+x2, data=aml) plot(fit) I guess because there are 2 levels

On 05/08/2010 8:46 PM, David B. Thompson, Ph.D., P.E., D.WRE, CFM wrote: > On Aug 5, 2010, at 5:20 PM, Duncan Murdoch wrote: > >> On 05/08/2010 7:18 PM, David B. Thompson, Ph.D., P.E., D.WRE, CFM wrote: >>> I'm running R 2.11.1 (MacBook Pro and OS X 10.6.4) and am trying to set a line type in the contour() function. What I did was: >>> >>>> >

R-Masters, I need to produce high resolution line plots and place labels on the curves. It seems that cex must be high relative to the other cex values in order to produce sufficiently large & legible tick labels at high resolutions. But high cex values cause the curve labels to become gigantic when using Hmisc. I've struggled and searched the archives, but cannot find a way of

hello, i managed to get CI's for my curve - but now I need the intervall for a certain y point (y_tenth) of the curve.. can anyone help me with this? #####data: por<-data.frame(list(structure(list(run = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L), .Label = c("1", "3", "4"), class = "factor"), press = c(15, 21, 24, 29.5, 15, 21,

Dear all, I want to let R automatically download available files from a website in certain folders. Since the files I need, if exist, are in a certain directory, I used download.file() function with a predesigned directory. However, some of the folders do not have the file I want. When this occurs, my program returns error and terminates. So I wonder if there is any way that I can check if a

Hello Everyone I am new to R I have drawn indifference curves using the program below (Contour Plot) u <- function(x, y) x^0.5 + y^0.5 x <- seq(0, 1000, by=1) y <- seq(0, 1000, by=1) a <- c(10, 20, 30) contour(x, y, outer(x, y, u),levels=a,col="blue") Now can any body please tell me how to draw Heat maps and that too on the same indifference curve plot (contour)

Hi, I'm attempting to calculate the 0.25 and 0.97 quantiles for tree height (0-50 meters) against tree age (0-300 years) and I am running into some difficulty with the plotted grafic. I've run the examples in the quantreg help and can get those to work properly and by plugging in my data I can also get the lines plotted on my dataset. Unfortunately I'm running into a problem with the

Hi, I have a query about the maximum length of vector that can be plotted in one go in a postscript driver. Try the following code (in 1.7.0; version details below): t <- seq(from=0, to=4*pi, length=200000) y <- sin(t) postscript(file="o.ps") plot(t, y, type="l") dev.off() If I view the postscript file o.ps in "gv", it takes many seconds before eventually

Hello, I'm trying to fit a logistic curve to data but I'm having a hard time discovering how. Every tutorial I've come across either assumes the logistic curve has 0<y<1 or assumes I have multiple categories of data I simply have two vectors, a and b, of equal length with no missing data, and I suspect they follow a logistic curve. The vectors are a<-c(39609, 39643,

Dear listers, As an exercise, I am trying to fit a logistic model with nlme. Blue tit pulli (youngs) were weighted occasionnally (for field reasons) along time in 17 nestboxes. Individuals where not idenfied but their age was known. This means that for a given age several measurements were done but individuals could not be identified from a time to the other. This makes repeated values for

Hi all My client want a shared folder so that only two people can access and write to it. I added the following to my smb.conf [private] writeble = yes guest ok = no path = /home/samba/private valid users = wendy, pierre write list = wendy, pierre force user = nobody All i need yo do now is make the folder "private" to nobody. How would i go on by doing that and would this work.