##pdf("nit.pdf") nit0 <- read.table("http://www.stat.ufl.edu/~winner/data/fertyld.dat", col.names=c("nit", "nit20", "yield")) x2=nit0$nit20^2 nit1 <- data.frame(y=nit0$yield, x=nit0$nit20, x2=x2) attach(nit1) nit1 nit.mod2 <- lm(y ~ x + x2) summary(nit.mod2) nit.mod3 <- nls(y ~ M - a*(R^x), start=c(M=100, a=75, R=0.5)) summary(nit.mod3) nit.mod4 <- nls(y ~ M*(1-C*exp(k*x)), start=c(M=100, C=.75, k=-1)) summary(nit.mod4) nit2 <- subset(nit1,x>0) nit.mod1 <- nls(y ~ a*(x^b), start=c(a=1, b=1), data=nit2) summary(nit.mod1) ## dev.off()