Displaying 2 results from an estimated 2 matches for "cilim".
Did you mean:
calim
2010 Jan 07
1
LD50 and SE in GLMM (lmer)
...y help me with this?
And one more question: is it correct to use pnorm (where John Maindonald used exp(hat)/(1+exp(hat)))?
Thanks so much in advance!
Linda
Previous post by John Maindonald:
ciplot <- function(obj=model4, data=zerotest,xcol=2,nam="hours"){
cilim <- function(obj, xcol){
b <- fixef(obj)
vcov <- summary(obj)@vcov
X <- unique(model.matrix(obj))
hat <- X%*%b
pval <- exp(hat)/(1+exp(hat)) # NB, designed for logit link
U <- chol(as.matrix(summary(obj)@vcov))
se <- sqrt(apply(X%*%t(U),...
2008 May 07
3
predict lmer
Hi,
I am using lmer to analyze habitat selection in wolverines using the
following model:
(me.fit.of <-
lmer(USED~1+STEP+ALT+ALT2+relM+relM:ALT+(1|ID)+(1|ID:TRKPT2),data=vdata,
control=list(usePQL=TRUE),family=poisson,method="Laplace"))
Here, the habitat selection is calaculated using a so-called discrete
choice model where each used location has a certain number of
alternatives