# NOT RUN {
lhspec <- bspec(lh)
# compare distributions visually:
par(mfrow=c(2,1))
plot(lhspec)
plot(two.sided(lhspec))
par(mfrow=c(1,1))
# ...and numerically:
print(cbind("frequency"=lhspec$freq,
"median-1sided"=quantile(lhspec,0.5),
"median-2sided"=quantile(two.sided(lhspec),0.5)))
# }
Run the code above in your browser using DataLab