# NOT RUN {
data(paris)
plot(pm10~pm25, data=paris)
o<-lm(pm10~pm25, data=paris) #negative intercept! it's meaningless..
o1<-lognlm(pm10 ~ pm25, data=paris, lik=TRUE)
abline(coef=coef(o), col=2)
abline(coef=coef(o1), col=3)
# }
Run the code above in your browser using DataLab