Hi,
I am trying to plot the conditional effects of a gamm fit with brms, the categorical predictor (ECO_GROUP2) has 7 levels so the resulting plot is very busy. Is there a way to only plot a subset of the levels (i.e. Lodges & Marshes)?
fit<-brm(Use~s(Month, by=ECO_GROUP2)+(1+ECO_GROUP2|Animal_ID,
data=data,
family=“bernoulli”,
warmup=1000,
iter=10000,
chains=3,
cores=3)
plot(conditional_effects(fit, method=“posterior_inpred”), ask=FALSE)
ECO_GROUP2 is categorical with 7 levels and Month is numeric
Thanks!