[R] strangeness in Predict() {rms}

[Bill.Venables at csiro.au]

This works without a glitch on my linux system (info below).  You might try upgrading your R to 2.10.1, perhaps.

> sessionInfo()
R version 2.10.1 (2009-12-14) 
x86_64-unknown-linux-gnu 

locale:
 [1] LC_CTYPE=en_AU.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_AU.UTF-8        LC_COLLATE=en_AU.UTF-8    
 [5] LC_MONETARY=C              LC_MESSAGES=en_AU.UTF-8   
 [7] LC_PAPER=en_AU.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_AU.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] splines   grid      stats     graphics  grDevices utils     datasets 
[8] methods   base     

other attached packages:
[1] rms_2.1-0       plyr_0.1.9      Design_2.3-0    Hmisc_3.7-0    
[5] survival_2.35-9

loaded via a namespace (and not attached):
[1] ASOR_0.1       cluster_1.12.1 lattice_0.18-3 tcltk_2.10.1   tools_2.10.1  
>  


Bill Venables
CSIRO/CMIS Cleveland Laboratories


-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Dylan Beaudette
Sent: Wednesday, 17 February 2010 10:05 AM
To: r-help at r-project.org
Subject: [R] strangeness in Predict() {rms}

Hi,

Running the following example from ?Predict() throws an error I have never 
seen before:

set.seed(1)
x1 <- runif(300)
x2 <- runif(300)
ddist <- datadist(x1,x2); options(datadist='ddist')
y  <- exp(x1+ x2 - 1 + rnorm(300))
f  <- ols(log(y) ~ pol(x1,2) + x2)
p1 <- Predict(f, x1=., conf.type='mean')

Error in paste(nmc[i], "=", if (is.numeric(x)) format(x) else x, sep = "") : 
  cannot coerce type 'closure' to vector of type 'character'
In addition: Warning message:
In is.na(v) : is.na() applied to non-(list or vector) of type 'closure'

Here is the output from sessionInfo()

R version 2.9.0 (2009-04-17) 
i686-pc-linux-gnu 

locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] grid      splines   stats     graphics  grDevices utils     datasets 
[8] methods   base     

other attached packages:
[1] plyr_0.1.9         mgcv_1.5-5         RColorBrewer_1.0-2 nlme_3.1-94       
[5] rms_2.1-0          Hmisc_3.7-0        survival_2.35-6    lattice_0.17-25   

loaded via a namespace (and not attached):
[1] cluster_1.12.0


Any ideas?
Thanks!

Dylan


-- 
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341

______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.