Question: Compute confidence intervals for predicted values in glm R same as proc genmod SAS
gravatar for besson.andrea
2.9 years ago by
besson.andrea0 wrote:

Hi everyone,

My current goal is to translate a SAS program to R. My issue here is that I can't find the equivalence/how to compute predicted values for a glm model as the ones I get in SAS.

I have found the following explanations in SAS documentation : But I have troubles to translate the formula with a concrete code in R, specially concerning the covariance ...

For now I have computed the CI the following way (I'm working with lists)

  models <- lapply(combia, function (x) {glm(cbind(n,s-n) ~ scm, family = binomial, data = x)})

  # Get confidence intervals
  p <- lapply(models, function (x) {predict(x, type = "link", = TRUE)})

  lowerlogit <- lapply(p, function(x) {x$fit - 1.96*x$})
  upperlogit <- lapply(p, function(x) {x$fit + 1.96*x$})

  borneinf <- lapply(lowerlogit, function(x) {exp(x)/(1+exp(x))})
  bornesup <- lapply(upperlogit, function(x) {exp(x)/(1+exp(x))})

But I don't get the same CI as SAS, and it's a problem for me in the futures steps.

I hope some of you would be able to help me on that, that would be so helpful !!

Thanks in advance

R • 1.2k views
ADD COMMENTlink modified 2.9 years ago by Jean-Karim Heriche24k • written 2.9 years ago by besson.andrea0
gravatar for Jean-Karim Heriche
2.9 years ago by
EMBL Heidelberg, Germany
Jean-Karim Heriche24k wrote:

What you're doing seems OK at first glance. The difference could be in the way the covariance matrix is estimated. This blog post and the referenced SO post might help. Another way could be to use the glm.predict package.

ADD COMMENTlink written 2.9 years ago by Jean-Karim Heriche24k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 925 users visited in the last hour