# NOT RUN {
n=20
s=.2
set.seed(10) #just to get reproducible results..
#covariates
x<-seq(.1,10,l=n)
z<-rnorm(n)
#response
mu<- 10+.5*x- z #linear regression function
y<-rlnorm(n, log(mu)-s^2/2, s) #data..
o<- lognlm(y~x+z) #the model
summary(o, sandw=TRUE)
# }
Run the code above in your browser using DataLab