[R-sig-ME] lme4: plotting profile density (not Zeta) manually not by lattice

[Simon Harmel]

Dear All,

I'm trying to plot the sampling distributions of my model parameters using `
densityplot()` from the `lattice` package but lattice often throws an error
even if one of the estimate's density distribution is highly skewed or
funny-looking.

Is there a better package or a better way (even manually) to plot the
densities (not Zeta) from a `profile()` call?

library(lattice)
library(lme4)

hsb <- read.csv('
https://raw.githubusercontent.com/rnorouzian/e/master/hsb.csv')
m31 <- lmer(math ~ ses*meanses + (ses | sch.id), data = hsb)

p <-  profile(m31)
densityplot(p)
Error in UseMethod("predict") :
  no applicable method for 'predict' applied to an object of class "NULL"

	[[alternative HTML version deleted]]