search for: limfo

...@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), 1, function(x)sum(x^2))) list(hat=pval, se=se, x=X[,xcol]) } limfo <- cilim(obj, xcol) hat <- limfo$hat se <- limfo$se x <- limfo$x upper <- hat+se lower <- hat-se ord <- order(x) plot(x, hat, yaxt="n", type="l", xlab=nam, ylab="") rug(x) lines(x[ord], lower[ord]) lines(x[or...

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