search for: dmodel

...sting this with a more appropriate subject. Do I need to define limits as the error message seems to suggest? If so, how? 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(dm...

Hi all, Do any packages implement density estimation in a modelling framework? I want to be able to do something like: dmodel <- density(~ a + b, data = mydata) predict(dmodel, newdata) This isn't how sm or KernSmooth or base density estimation works. Are there other packages that do density estimation? Or is there some reason that this is a bad idea. Hadley -- http://had.co.nz/