[ADMB Users] Transformation of mcmc values

Håkon Holand hakon.holand at bio.ntnu.no
Tue Mar 20 09:02:33 PDT 2012


I was wondering if anyone knows a way of transforming values obtained from a HPD of mcmc runs of a model?

I am running models with a binary response variable (0/1) (family=binomial) and using individual identity as a random variable. I am currently using R 2.14.1 (64-bit) and glmmadmb v.

The intervals looks like this:

> head(HPDinterval(m1))
              lower      upper
beta.1 -487.3376278 -389.89072
beta.2   88.3658266  167.10142
beta.3 -106.7078281  -46.53788
beta.4    1.2919346   34.92751
beta.5   -0.8259944   50.58856
beta.6  -49.8802241  -12.91753

First of all: I know that these values are from the parameters fitted internally, using an orthogonalized version of the original design matrix, not the original coefficients.

The question is:  (if possible) How can I transform these values to, for example, Logit? I would like to have my estimated value and its limits on the same scale at least (and hopefully a "reader friendly" scale).

Any help will be deeply appreciated!

Håkon Holand
Håkon Holand, PhD.Student
Centre for Conservation Biology
Department of Biology
Norwegian University of Science and Technology
NO-7491 Trondheim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.admb-project.org/pipermail/users/attachments/20120320/93e6d8f6/attachment.html>

More information about the Users mailing list