Skip to content

Is there an lm() equivalent of panel.lmline()?

2 messages · Alex van der Spek, Greg Snow

Alex van der Spek
# 
panel.lmline returns intercept and slope of y ~ x subsetted to the
combination of conditioning factors given to xyplot in lattice.

for instance:

xyplot(Xvalues ~ log(Qvalues)|Tfac, data = df7, panel = panel.lmline)


I am looking to find the equivalent formulation for lm() proper. If I do
this:

lmcal <- lm(Xvalues ~ log(Qvalues):Tfac, data = df7)

Only one value of the intercept is returned. It appears as all the
'individual' intercepts are summed to form one overall intercept.

Any help much appreciated!

Regards,
Alex van der Spek
Greg Snow
# 
You can use the lmList function in the nlme package to do several
seperate regressions, or use a model that allows for multiple
intercepts.  Possibly Xvalues ~ 0 + log(Qvalues)*Tfac or Xvalues ~ 0 +
Tfac + log(Qvalues):Tfac (assuming Tfac is a factor).
On Tue, Oct 30, 2012 at 5:53 AM, Alex van der Spek <doorz at xs4all.nl> wrote: