search for: ciplot

...o get my LD50 SE (I would need the SE to be expressed in terms of x). Could anybody 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...

...the Design package. Putting conf.int=.95 as an argument in plot does not work. The default appears to be .99 not .95 as stated in the package Design manual (p. 164). 2. My sweave chuck is below and my output is attached as well as linked here: http://www.sonoma.edu/users/s/stanny/330A/project/ciplot.pdf 3. Alternatively, how do I make the confidence interval text smaller for the second plot (the interval text is cut-off) or color the intervals? 4. I have 14 plots and would like to have 2 per page. Thanks in advance. <<ghg07,fig=T>>= par(mfrow=c(2,1)) s.fghg.new.06 <- s...

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

I am looking to create what I would call a "simple variation" on the boxplot. What I would like to do is to be able to plot the upper and lower confidence limits as the "box" and the 10th and 90th percentile as the whiskers. What I have done is write the code to create a dataframe, the columns of which are the mean, sd, 10th percentile, 90th percentile, lower confidence limit