Re: [R] Troubles plotting lrm output in Design Library

2008-05-30 Thread Prof Brian Ripley
On Fri, 30 May 2008, Frank E Harrell Jr wrote: Prof Brian Ripley wrote: On Fri, 30 May 2008, Patrizio Frederic wrote: Thanks for the solution. But I wish someone would tell us why the solution makes sense. If you have changed your R environment in any way from the defaults please let us kn

Re: [R] Troubles plotting lrm output in Design Library

2008-05-30 Thread Frank E Harrell Jr
Prof Brian Ripley wrote: On Fri, 30 May 2008, Patrizio Frederic wrote: Thanks for the solution. But I wish someone would tell us why the solution makes sense. If you have changed your R environment in any way from the defaults please let us know. Frank -- Frank E Harrell Jr Professor and

Re: [R] Troubles plotting lrm output in Design Library

2008-05-30 Thread Prof Brian Ripley
On Fri, 30 May 2008, Patrizio Frederic wrote: Thanks for the solution. But I wish someone would tell us why the solution makes sense. If you have changed your R environment in any way from the defaults please let us know. Frank -- Frank E Harrell Jr Professor and Chair School of

Re: [R] Troubles plotting lrm output in Design Library

2008-05-30 Thread Frank E Harrell Jr
Patrizio Frederic wrote: Thanks for the solution. But I wish someone would tell us why the solution makes sense. If you have changed your R environment in any way from the defaults please let us know. Frank -- Frank E Harrell Jr Professor and Chair School of Medicine

Re: [R] Troubles plotting lrm output in Design Library

2008-05-30 Thread Patrizio Frederic
> Thanks for the solution. But I wish someone would tell us why the solution > makes sense. If you have changed your R environment in any way from the > defaults please let us know. > > Frank > > -- > Frank E Harrell Jr Professor and Chair School of Medicine > Depa

Re: [R] Troubles plotting lrm output in Design Library

2008-05-30 Thread Frank E Harrell Jr
Patrizio Frederic wrote: 2008/5/29 Frank E Harrell Jr <[EMAIL PROTECTED]>: Patrizio Frederic wrote: dear Harrell, thank you for quick reply and suggestions. I still have the problem: library(Design) x = rnorm(100) y = runif(100)<(exp(x)/(1+exp(x))) y = 0*y+1*y d = datadist(x,y) options

Re: [R] Troubles plotting lrm output in Design Library

2008-05-30 Thread Patrizio Frederic
2008/5/29 Frank E Harrell Jr <[EMAIL PROTECTED]>: > Patrizio Frederic wrote: >> >> dear Harrell, >> thank you for quick reply and suggestions. I still have the problem: >> >> library(Design) >> x = rnorm(100) >> y = runif(100)<(exp(x)/(1+exp(x))) >> y = 0*y+1*y >> d = datadist(x,y) >> optio

Re: [R] Troubles plotting lrm output in Design Library

2008-05-29 Thread Frank E Harrell Jr
Patrizio Frederic wrote: dear Harrell, thank you for quick reply and suggestions. I still have the problem: library(Design) x = rnorm(100) y = runif(100)<(exp(x)/(1+exp(x))) y = 0*y+1*y d = datadist(x,y) options(datadist="d") fit = lrm(y~x) # works fine, but plot(fit) #produce the error

Re: [R] Troubles plotting lrm output in Design Library

2008-05-29 Thread Patrizio Frederic
dear Harrell, thank you for quick reply and suggestions. I still have the problem: library(Design) x = rnorm(100) y = runif(100)<(exp(x)/(1+exp(x))) y = 0*y+1*y d = datadist(x,y) options(datadist="d") fit = lrm(y~x) # works fine, but plot(fit) #produce the error message Error in value.chk

Re: [R] Troubles plotting lrm output in Design Library

2008-05-29 Thread Frank E Harrell Jr
Patrizio Frederic wrote: Dear R-helpers, I'm having a problem in using plot.design in Design Library. Tho following example code produce the error: n <- 1000# define sample size set.seed(17) # so can reproduce the results age<- rnorm(n, 50, 10) blood.pressure <- rnorm(n, 120, 15

[R] Troubles plotting lrm output in Design Library

2008-05-29 Thread Patrizio Frederic
Dear R-helpers, I'm having a problem in using plot.design in Design Library. Tho following example code produce the error: > n <- 1000# define sample size > set.seed(17) # so can reproduce the results > age<- rnorm(n, 50, 10) > blood.pressure <- rnorm(n, 120, 15) > cholesterol<