Displaying 1 result from an estimated 1 matches for "2001syr".
Did you mean:
2001se
2005 Jun 14
2
ordinary polynomial coefficients from orthogonal polynomials?
...from
lm(billions ~ a+b*decade+c*decade^2+d*decade^3)
but that gives: "Error in eval(expr, envir, enclos) :
Object "a" not found"
> decade <- c(1950, 1960, 1970, 1980, 1990)
> billions <- c(3.5, 5, 7.5, 13, 40)
> # source: http://www.ipcc.ch/present/graphics/2001syr/large/08.17.jpg
>
> pm <- lm(billions ~ poly(decade, 3))
>
> plot(decade, billions, xlim=c(1950,2050), ylim=c(0,1000),
main="average yearly inflation-adjusted dollar cost of extreme weather
events worldwide")
> curve(predict(pm, data.frame(decade=x)), add=TRUE)
&...