pdat <- with(data,
             expand.grid(year_begin = seq(min(year_begin), max(year_begin), length = 100),
                         Tree = levels(Tree)))

pdat <- cbind(pdat, fitted = predict(model, newdata = pdat, type = "response"))