curve fitting problems in R -
i have quick question here new r language. not sense code right curve fitting: poly fit checks out, log , exponential fit unsure of (e fit def wrong). help appreciated , thanks. -sam
polyfit <- function(xs,ys,degree) { #polynomial fitting info set grade 3 fit3 <- lm(ys~poly(xs,degree,raw=true)) xx <- seq(0,160, length=50) plot(xs,ys,pch='@') lines(xx, predict(fit3, data.frame(xs=xx)), col="green") } logfit <- function(xs,ys) { #graph info set log(x), y logestimate = lm(ys ~ log(xs)) plot(xs,ys,pch='@') lines(xs,predict(logestimate),col='green') } efit <- function(xs,ys) { logestimate = lm(log(ys) ~ xs) plot(xs,ys,pch='@') lines(xs,predict(logestimate),col='green') } for clarification, xs = x points, ys = y points
no, this:
efit <- function(xs,ys) { expestimate = lm(log(ys) ~ xs) plot(xs, ys,pch='@') lines(xs, exp(predict(expestimate)),col='green') } i changed model name create distinct logfit function, though it's internal function. in model, fit log(y), predict log(y), if want see how prediction stack against real y, need untransform them, , exp(log(y)) = y.
r
No comments:
Post a Comment