x <- rexp(100)
y <- 1/x + rnorm(100)
summary(data.frame(y=y,x=x))