# Example 1
n <- 300
error <- rglg(n,0,1,1)
y <- 0.5 + error
fit <- glg(y~1,data=as.data.frame(y))
deviance_residuals(fit)
# Example 2
n <- 300
error <- rglg(n,0,1,1)
x <- runif(n,-3,3)
y <- 0.5 + 2*x + error
fit <- glg(y~x,data=as.data.frame(y,x))
deviance_residuals(fit)
Run the code above in your browser using DataLab