example.R

# Create some data
x <- seq(0, 10, by=0.01)
y <- exp(1)^-x
df <- data.frame(x=x, y=y)

# Add some noise.
df$y <- jitter(df$y, amount=.1)

print("Plotting original data...")
plot(df)
print("Press [enter] to continue")
line <- readline()

print("Training model based on log(y) ~ x.")
model <- lm(log(y) ~ x)
print(summary(model))
plot(model)
print("Press [enter] to continue")
line <- readline()

print("Predictions...")
df$pred <- predict(model, list(df$x))
plot(df$y)
lines(exp(df$pred),col="red",lwd=2)
print("Press [enter] to continue")
line <- readline()