This function takes the output of get_mle called with h2.method = "grid" and return.profile = TRUE and draws a line plot of the profile log-likelihood against \(h^2\).

plot_profile_h2(
  fit,
  xlab = "h2",
  ylab = "Profile log-likelihood",
  main = "",
  mark.max = TRUE,
  ...
)

Arguments

fit

A list returned by get_mle that contains a $profile element (i.e., fitted with h2.method = "grid" and return.profile = TRUE).

xlab

A character string specifying the x-axis label. Default "h2".

ylab

A character string specifying the y-axis label. Default "Profile log-likelihood".

main

A character string specifying the plot title. Default "".

mark.max

A logical scalar. If TRUE (default), draw a vertical dashed line at the \(h^2\) value that maximizes the profile log-likelihood.

...

Further graphical parameters passed to plot.

Value

The profile data frame (invisibly).