# NOT RUN {
vehicles_obs <- vehicles[vehicles@data$obs1pred0 == 1,]
vehicles_obs$transheavy <- bxcxt(vehicles_obs$heavy, vehicles_obs$width)[[1]]
qqnorm(vehicles_obs$heavy)
qqnorm(vehicles_obs$transheavy)
lambda <- bxcxt(vehicles_obs$heavy, vehicles_obs$width)[[2]]
# }
Run the code above in your browser using DataLab