[R-sig-ME] lme4: plotting profile density (not Zeta) manually not by lattice
Simon Harmel
@|m@h@rme| @end|ng |rom gm@||@com
Mon Oct 26 01:54:58 CET 2020
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]]
More information about the R-sig-mixed-models
mailing list