> Call:
> lm(formula = curr_data[[tmin_col]] ~ curr_data[[year_col]] -
> 1 | curr_data[[month_col]])
First, this is not a sensible formula for lm; lm() does not use '|' to denote
grouping. It would be a valid formula for xyplot, in which | specifies grouping
variables. In lm(), '|' is simply being treated as 'or', which is why one of
your coefficients is for ' curr_data[[year_col]] | curr_data[[month_col]]TRUE'
Second, you should not normally need things like curr_data[[month_col]], either
in lm or xyplot. If curr_data is a data frame, things like
lm(tmin ~ year, data=curr_data)
should work.
Third, 'nested' models in lm use the nesting operator '/', not '|'. So if you
want 12 lines with separate intercept and gradient from an lm, you need (with
month as a factor and the default intercept suppressed)
lm(tmin~month+year/month -1, data=curr_data) #-1 suppresses the intercept and
provides a zero-based intercept for each month
This gives you 12 'month' intercepts and one gradient per month. If you wanted
a common intercept you'd do
lm(tmin~ year/month, data=curr_data)
But beware; the coefficients in both cases cannot be interpreted as a simple
gradient and intercept for each month: if I recall correctly, the gradients for
month2 and on are modelled as an additive increment on the first month
gradient. Use predict() if you want an easy way to predict a value for a given
(possibly fractional) time of year and 'month'.
[Incidentally I don’t immediately see why that is a sensible thing to do - this
fits a monthly summary against a numeric year. But I'm going to assume you know
what you want there.]
Finally, though, this model will not help you much with lattice as there's no
_simple_ way of putting those lines on different panels in a lattice plot. If
you just want a line on each of 12 panels, that's much easier. You can use the
panel() function with panel.lmline to put it there. For example, if you want to
plot a line over the data, use
xyplot(tmin~year|month, curr_data,
panel=function(x, y, ...) {
panel.xyplot(x, y, ...)
panel.lmline(x, y, ...)
}
)
S Ellison
*******************************************************************
This email and any attachments are confidential. Any use, copying or
disclosure other than by the intended recipient is unauthorised. If
you have received this message in error, please notify the sender
immediately via +44(0)20 8943 7000 or notify postmas...@lgcgroup.com
and delete this message and any copies from your computer and network.
LGC Limited. Registered in England 2991879.
Registered office: Queens Road, Teddington, Middlesex, TW11 0LY, UK
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
