manuscript/emmeans.iv.pkgcompare.R

In mac-theobio/effects: isolated confidence intervals and bias-corrected predictor effects

library(shellpipes)
library(dplyr)
library(emmeans)

loadEnvironments()

emmeans_iv <- (emmeans(fit, specs=~age
		, type="response"
		, at=list(age=unique(focal))
		, weights="proportional"
		, nesting=NULL
	)
	|> data.frame()
	|> rename(
		fit=emmean
		, se=SE
		, lwr=lower.CL
		, upr=upper.CL
	)
	|> mutate(model="emmeans")
	|> select(-df)
)
head(emmeans_iv)

saveVars(emmeans_iv)