Having trouble with datadist in the Design library. When I set the
datadist option, all of the fitting functions give error messages
similar to "Error in exists(x, envir, mode, inherits) : invalid first
argument". Code below was run in R version 2.5.0 under windows, but get
similar results with R version 6.2, except with additional problem that
it will not fit model with "options(datadist=NULL)" since use of
storage.mode(x) <- "single" seems to now be defunct instead of
deprecated.     

For example
> library(Design)
> set.seed(243)
> temp.df<-data.frame(x1=rnorm(100),x2=rnorm(100))
> temp.df$y<- temp.df$x1 - temp.df$x2 + rnorm(100)/2
> 
> dd<-datadist(temp.df) 
> options(datadist=dd)
> 
> options()$datadist
                        x1          x2          y
Low:effect      -0.6299961 -0.82850974 -1.0273244
Adjust to       -0.1309631 -0.02919389  0.1342749
High:effect      0.7723673  0.70284942  1.2207567
Low:prediction  -1.4721508 -1.76692677 -2.3062141
High:prediction  1.8141205  1.60072989  2.8880009
Low             -2.6437032 -2.60640789 -4.2372689
High             2.8397685  1.95926958  4.8788867
> 
> temp.ols<-ols(y~x1+x2, data=temp.df)
Error in exists(x, envir, mode, inherits) : 
        invalid first argument
> 
> options(datadist=NULL)
> 
> temp.ols<-ols(y~x1+x2, data=temp.df)
Warning messages:
1: use of storage.mode(x) <- "single" is deprecated: use mode<- instead 
2: use of storage.mode(x) <- "single" is deprecated: use mode<- instead 
> temp.ols

Linear Regression Model

ols(formula = y ~ x1 + x2, data = temp.df)

         n Model L.R.       d.f.         R2      Sigma 
       100      256.6          2     0.9231     0.4495 

Residuals:
     Min       1Q   Median       3Q      Max 
-1.27308 -0.26226 -0.01148  0.26467  1.14479 

Coefficients:
              Value Std. Error        t Pr(>|t|)
Intercept -0.001592    0.04509  -0.0353    0.972
x1         1.068831    0.04428  24.1401    0.000
x2        -1.071857    0.04356 -24.6066    0.000

Residual standard error: 0.4495 on 97 degrees of freedom
Adjusted R-Squared: 0.9216 

> version
               _                           
platform       i386-pc-mingw32             
arch           i386                        
os             mingw32                     
system         i386, mingw32               
status                                     
major          2                           
minor          5.0                         
year           2007                        
month          04                          
day            23                          
svn rev        41293                       
language       R                           
version.string R version 2.5.0 (2007-04-23)



Would be grateful for any help. Thanks. 

 

Daniel Holder
Biometrics Research WP53B-120
Merck Research Labs
West Point, PA 19486



------------------------------------------------------------------------------
Notice:  This e-mail message, together with any attachme...{{dropped:17}}

