set.seed(20)
x = seq(from=0, to=20, by=0.1)
err = rnorm(length(x), mean=10, sd=50)
y.obs = 100 + 120*x - 12*x^2 + 0.5*x^3 + err
plot(y.obs ~ x, pch=16)

set.seed(200)
x = seq(from=0, to=200, by=0.1)
err = rnorm(length(x), mean=10, sd=50)
y.obs = 100 + 120*x - 12*x^2 + 0.5*x^3 + err
plot(y.obs ~ x, pch=16)