equations for predicting cloglog models

cloglog.R

d <- data.frame(x = rnorm(100) + c(1, 2), y = c(0, 1))

m <- glm(y ~ x, data = d, family = binomial(link = cloglog))


newx <- seq(-5, 5, length.out = 100)
predict.line <- predict(m, newdata = data.frame(x = newx), type = 'response')

plot(y ~ x, d)
lines(predict.line ~ newx)


manual.predict <- 1-exp( -exp(m$coefficients['(Intercept)'] + m$coefficients['x']*newx))

lines(manual.predict ~ newx, col = 'red')