# Plotting profile likelihood curves in R

I am trying to figure out how to plot the profile likelihood curve of a GLM

parameter with 95% pLCI's on the same plot. The example I have been trying

with is below. The plots I am getting are not the likelihood curves that I

was expecting. The y-axis of the plots is tau and I would like that axis

to be the likelihood so that I have a curve that maxes at the parameter

estimate. I am not sure where I find those likelihood values? I may just

be misinterpreting the theory behind this. Thanks for any help you can give.

Max