similar to: Highlight overlapping area between two curves

Dear R helpers, I would like to have abline, for a lm model for example, lying within data range. Do you know how to get it? Thank in advance Nguyen D Nguyen #CODE x<- rnorm(200, 35,5) y<- rnorm(200, 0.87,0.12) plot(y~x, xlim=c(0,50), pch=17, bty="l") abline(lm(y~x)) # I would like abline is between min(x) and max(x) [[alternative HTML version deleted]]

Dear all, I found something strange when calculating sin of pi value sin(pi) [1] 1.224606e-16 pi [1] 3.141593 sin(3.141593) [1] -3.464102e-07 Any help and comment should be appreciated. Regards Nguyen ____________________________ Nguyen Dinh Nguyen Garvan Institute of Medical Research Sydney, Australia

Dear R helpers I have a data set sth like this: set.seed(123);dat <- data.frame(ID= c(rep(1,2),rep(2,3), rep(3,3), rep(4,4), rep(5,5)), var1 =rnorm(17, 35,2), var2=runif(17,0,1)) dat ID var1 var2 1 1 33.87905 0.02461368 2 1 34.53965 0.47779597 3 2 38.11742 0.75845954 4 2 35.14102 0.21640794 5 2 35.25858 0.31818101 6 3 38.43013

Dear all, I have a stat question that may not be related to R, but I would like to have your advice. I have just read a medical paper in which the authors report the 1-p (where p is the cumulative survival probability from the Kaplan Meier curve) as incidence of disease. Specifically, the study followed ~12000 women on drug A and ~20000 women on drug B for 12 months. During that period

Dear all, A friend of mine, who just installed Rcmdr package. When calling the package, the error as following comes up. And actually, I've checked in CRAN packages, the tcltk package doesn't exist anymore, but tcltk2. Any help is appreciated Regards Nguyen D Nguyen Garvan Institute of Medical Research Sydney, Australia library(Rcmdr) Loading required package: tcltk Loading Tcl/Tk

Dear Christian, Every single time check update package, ?cwhmisc? always requests updating. I?m aware that the package was latest updated in CRAN on 20Nov2008. Is there anything wrong with my R library or somethingelse? I use R 2.8.1 on Window XP service pack 2 Regards Nguyen Garvan Institute of Medical Research Sydney, Australia

Hi everyone. My problem: I've 4 distributions: A, B, C and D: #A mA=16.6 sA=3.0 #B mB=18.9 sB=3.2 #C mC=20.3 sC=2.1 #D mD=24 sD=0.8 ###Graphiques ensembles plot(function(x) dnorm(x,mA,sA),0,40,col="orange",ylim=c(0,0.5)) plot(function(x) dnorm(x,mB,sB),0,40,add=T,col="green") plot(function(x) dnorm(x,mC,sC),0,40,add=T,col="blue")

Dear Prof Gray and everyone, As our package developers discussed about incompatibility between Design and survival packages, I faced another problem with cmprsk- a survival dependent packacge. The problem is exactly similar to what happened to the Design package that when I just started running cuminc function, R was suddenly closed. These incidents suggest that maybe many other survival

Hello, I have finished setting up the load balancer with IPROUTE ... also patch the kernel to support DGD and now it''s working fine with the valuable guide at LARTC website, Julian Anastasov, and the kind people in this mailing list. Now I would like to launch a web server and a ftp server to the public but I''m stuck into a problem and really need your help. Currently internal

Dear Prof Harrell and everyone, My PC: Window XP service pack 3 and service pack 2 R version 2.8.1 and 2.9 alpha For the last 3 days, after updating R, my two computers have been facing problems when running existing and runable R commands that involves with Design package I attempt to use 'survest', but I failed all the times with R (both 2.8.1 and 2.9 alpha) being shut down

Hi List, I would like to plot multiple curves (parametric density curves) in one plot. For example: # parameters for three normal density curves parms = data.frame(ID=c(1,2,3),mu=c(50,55,60),sigma=c(10,12,15)) # I can easily draw three normal density curves using curve(): curve(dnorm(x,mean=parms$mu[1],sd=parms$sigma[1]),from=0, to=150, ylab="density", col="red")

Hi, I am trying to plot two distribution graph on one plot. But I dont know how. I set them to the same x, y limit, even same x, y labels. Code: x1=rnorm(25, mean=0, sd=1) y1=dnorm(x1, mean=0, sd=1) x2=rnorm(25, mean=0, sd=1) y2=dnorm(x2, mean=0, sd=1) plot(x1, y1, type='p', xlim=range(x1,x2), ylim=range(y1, y2), xlab='x', ylab='y') plot(x2, y2, type='p',

Sorry about the subject --- On Thu, 24/9/09, KABELI MEFANE <kabelimefane@yahoo.co.uk> wrote: From: KABELI MEFANE <kabelimefane@yahoo.co.uk> Subject: Re: [R] Multiply Normal Curves To: R-help@r-project.org Date: Thursday, 24 September, 2009, 11:48 AM R -helpers   i have been trying to do this problem without must success,i managed to do a graph for x, but it is not what i want to

I might be missing something but I thought that AUC was a measure for comparing ROC curves, so there is nothing else needed to "compare" them. The larger AUC is the higher correlation of 2 variables compared. No other measures or calculations are needed. Jarek Tuszynski -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On

Dear all, I try to compare the performances of several parameters to diagnose lameness in dogs. I have several ROC curves from the same dataset. I plotted the ROC curves and calculated AUC with the ROCR package. I would like to compare the AUC. I used the following program I found on R-help archives : From: Bernardo Rangel Tura Date: Thu 16 Dec 2004 - 07:30:37 EST

Hi, is there a way to plot distribution curves (say normal or chi sq etc) from within R? For example I looked up the *chisq family of functions but I'm not sure as to how I would use them to generate a plot of the chi sq distribution (for arbitrary d.o.f). Thanks, ------------------------------------------------------------------- Rajarshi Guha <rajarshi at presidency.com>

Dear R Users, Do you have a sample code for developing a nomogram with competing-risks? Any help is appreciated. Kind regards, ND Nguyen

I did some serach around. It seems that ROC curve computation is not supported on R. Anyone has some leads? Thanks. Jason ===== Jason G. Liao, Ph.D. Division of Biometrics UMDNJ School of Public Health 335 George Street, Suite 2200 New Brunswick, NJ 08903-2688 phone (732) 235-9748, fax (732) 235-9777 http://www.geocities.com/jg_liao __________________________________________________ Great

Hi all, I've got this graphics question which really should be easy. I want to shade an area between bounds under a curve. A suitable beginning seems to be the following: > plot(dnorm,-4,4) > segments(-4,0,4,0) > segments(-2,0,-2,dnorm(-2)) > segments(2,0,2,dnorm(2)) It is the area between -2 and 2 which I want to shade (or something similar). Hints anyone? Robert

Hello, list why not add the smart proposal by Greg Snow as a built-in function in {stats}, just changing the "x234" and "newc" lines to allow for more distributions to be generated ? Or do I miss an already existing function to do that ? Regards. Olivier # slight modification of the original code by Greg Snow [mailto:Greg.Snow at intermountainmail.org] # on April 04, 2007