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-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Dylan Beaudette Sent: Wednesday, 17 February 2010 10:05 AM To: r-help@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@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. ______________________________________________ R-help@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.