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) options(datadist=d) fit =

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)

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 Department

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

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

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

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