Skip to content

problem with read.table

6 messages · x, Talita Perciano, Nordlund, Dan (DSHS/RDA) +1 more

x
# 
I forgot to show that my file indeed contains x, y values:

x	y
1	8.2
2	18
3	-17.6
4	-3.4
5	6.8
6	11.3
7	1.5
8	11
9	-3.3
10	8.1
--- On Wed, 4/15/09, x <to_rent_2000 at yahoo.com> wrote:

            

      
      
      
    

  
  


  


    

    

      
      

        
      

    

    

      
      

        


  

        

      x
    

    

      
      #
    

  


  

      
attach() fixed the problem. Thanks!!

        
--- On Wed, 4/15/09, Talita Perciano <talitaperciano at gmail.com> wrote:

        

            

      
      
      
    

        

      
      
    

            

      
      
      
    

  
  


  


      

    

    

      
      

        


  

        

      Nordlund, Dan (DSHS/RDA)
    

    

      
      #
    

  


  

          

      
      
      
    

        
After 

    data_model.df = read.table("./verify.txt", header=TRUE, nrows=10)

you can attach the dataframe and refer to the individual variables in the dataframe by name

    attach(data_model.df)

or you can refer to the variables as

    data_model.df$x
    data_model.df$y 

or

    data_model.df[1]
    data_model.df[2]

You might want to read the "Introduction to R" manual, in particular, section 6.

Hope this is helpful,

Dan

Daniel J. Nordlund
Washington State Department of Social and Health Services
Planning, Performance, and Accountability
Research and Data Analysis Division
Olympia, WA  98504-5204

  
  


  


      

    

    
12 days later

    

      
      

        


  

        

      x
    

    

      
      #
    

  


  

      
Hi all,

Do I need to define limits as the error message seems to suggest? The error message, my code, the output and the first few lines of my data are all below.

Thank you!

"Error in Getlim(at, allow.null = TRUE, need.all = TRUE) : 
  variable dmodel.df does not have limits defined in fit or with datadist"

My code:
==================
library(Hmisc); library(Design); library(lattice);  

dmodel.df = read.table("./data_cub3.txt", header=TRUE, nrows=100)
f <- ols(dmodel.df$y1 ~ rcs(dmodel.df$x1,3) )
print(f)

dd <- datadist(dmodel.df$x1)
options(datadist="dd")
describe(dmodel.df)

print( Function(f) )
plot(dmodel.df$x1, dmodel.df$y1)
plot(f, add=TRUE, col="blue", pch=2)

Output:
====================
Linear Regression Model

ols(formula = dmodel.df$y1 ~ rcs(dmodel.df$x1, 3))

         n Model L.R.       d.f.         R2      Sigma 
       100      501.6          2     0.9934      45128 

Residuals:
     Min       1Q   Median       3Q      Max 
-68644.8 -31355.6   -849.9  31823.3 154196.6 

Coefficients:
           Value Std. Error      t  Pr(>|t|)
Intercept  35925    12789.6  2.809 0.0060121
dmodel.df  -1620      422.9 -3.832 0.0002260
dmodel.df' 25202      523.9 48.102 0.0000000

Residual standard error: 45130 on 97 degrees of freedom
Adjusted R-Squared: 0.9932 

Error in Getlim(at, allow.null = TRUE, need.all = TRUE) : 
  variable dmodel.df does not have limits defined in fit or with datadist

Sample data:
=========================
config	benchmark	x1	noise	y1
1	verify2	1	0.72	282
2	verify2	2	1.6	256

  
  


  


      

    

    

      
      

        


  

        

      David Winsemius
    

    

      
      #
    

  


  

      
On Apr 28, 2009, at 6:19 PM, x wrote:

        

            

      
      
      
    

        
Try instead:

dd <- datadist(dmodel.df)
f <- ols( y1 ~ rcs(x1,3), data=dmodel.df )

            

      
      
      
    

        

      
      
    

        
David Winsemius, MD
Heritage Laboratories
West Hartford, CT