[R] Error in Getlim()
x
to_rent_2000 at yahoo.com
Wed Apr 29 00:19:33 CEST 2009
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
More information about the R-help
mailing list