[R] Can I build an array of regrssion model?

[Thomas Lumley]

On Thu, 19 Dec 2002, Zhongming Yang wrote:

> Thanks,
>
> But why I can't draw regression line with the following code:
>
> rawData = scan("c:/zyang/mass/data/A01/1.PRN",
> what=list(numeric(),numeric()));
> len = length(rawData[[1]]);
> mod = lm(rawData[[2]]~bs(rawData[[1]],10,degree=1));
> plot(rawData[[1]],rawData[[2]],type='l', col="green", xlab="Da/z",
> ylab="m/z");
> abline(mod,col="red");
>

Because that's not what abline does.  It draws a straight line, not a
piecewise linear regression curve.

You might find termplot() useful.

	-thomas