similar to: predict function for GLMM

Can GLMM take formula derived from another object? foo <- glm (OVEN ~ h + h2, poisson, dataset) # ok bar <- GLMM (OVEN ~ h + h2, poisson, dataset, random = list (yr = ~1)) #error bar <- GLMM (foo$formula, poisson, dataset, random = list (yr = ~1)) #Error in foo$("formula" + yr + 1) : invalid subscript type I am using R2.1.0, lme4 0.8-2, windows xp. Below is a dataset if you

Members of this mailing list may be interested in the following new book: Beginner's Guide to GLM and GLMM with R. - A frequentist and Bayesian perspective for ecologists - Zuur AF, Hilbe JM and Ieno EN This book is only available from: http://www.highstat.com/BGGLM.htm This book presents Generalized Linear Models (GLM) and Generalized Linear Mixed Models (GLMM) based on both

Is there an easy way to get confidence intervals from "glmm" in Jim Lindsey's library(repeated)? Consider the following slight modification of an example from the help page: > df <- data.frame(r=rbinom(10,10,0.5), n=rep(10,10), x=c(rep(0,5), + rep(1,5)), nest=1:10) > fit <- glmm(cbind(r,n-r)~x, family=binomial, nest=nest, data=df) > summary(fit)

Dear R-users, I wanted to inform you that a new package called R2STATS is available, as a graphical front-end for the glm() and glmer() functions. The GUI is based on the RGTk2 and gWidgets packages by Michael Lawrence and John Verzani, and so requires that the GTK+ library be installed first on your system. This is done "automagically" when installing the RGtk2 package (or the

Dear R-list, I thought that I would let some of you know of a free R package, glmm.ADMB, that can handle mixed models for overdispersed and zero-inflated count data (negativebinomial and poisson). It was built using AD Model Builder software (Otter Research) for random effects modeling and is available (for free and runs in R) at: http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html I

Hello, I am using R to train a logistic regression model and save the resulting model to disk. I am then subsequently reloading these saved objects, and using predict.glm on them in order to make predictions about single-row data frames that are generated in real-time from requests arriving at an HTTP server. The following code demonstrates the sort of R calls that I have in mind: > cases

Hi, I have three questions concerning GLMMs. First, I ' m looking for a measure for the significance of the random variable in a glmm. I'm fitting a glmm (lmer) to telemetry-locations of 12 wildcat-individuals against random locations (binomial response). The individual is the random variable. Now I want to know, if the individual ("TIER") has a significant effect on the model

Greetings R Users! I have a data set of count responses for which I have made repeated observations on the experimental units (stream reaches) over two air photo dates, hence the mixed effect. I have been using Dr. Jim Lindsey's GLMM function found in his "repeated" measures package with the "poisson" family. My problem though is that I don't think the poisson

Could you just use lines(newX, myPred, col=2) -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of Paul Johnson Sent: Thursday, September 23, 2004 10:3 AM To: r help Subject: followup: Re: [R] Issue with predict() for glm models I have a follow up question that fits with this thread. Can you force an overlaid plot

Dear R-list, I try to adjust GLMM on incidence cancer data. Without random effect, in GLM the command is, for example with sex effect, glm(Observed~sex+offset(log(Expected)),family=poisson) because the observed are Poisson distribued with parameter Expected*incidence rate. But know I want to introduce random effect (for example spatial effect) and it seems to me that the "offset" does

Hello, I am using R to train a logistic regression model and save the resulting model to disk. I am then subsequently reloading these saved objects, and using predict.glm on them in order to make predictions about single-row data frames that are generated in real-time from requests arriving at an HTTP server. The following code demonstrates the sort of R calls that I have in mind: > cases

Hi, I know this is probably a basic question... But I don't seem to find the answer. I'm fitting a GLM with a Poisson family, and then tried to get a look at the predictions, however the offset does seem to be taken into consideration: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)), data=data, subset=28:36, family=poisson()) predict

Summary: We need to find a reviewer for our implementation of Ball-Laurus path profiling. It is well known that path profiling generates more precise information about a program's behaviour than edge profiling. We are conducting a research project with the goal of developing a methodology to make feedback-directed optimization (FDO) more sound. We are developing combined profiles that enable

On 27/04/2018 9:25 AM, Hadley Wickham wrote: > Hi all, > > Very surprising (to me!) and mystifying result from predict.glm(): the > predictions vary depending on whether or not I use ns() or > splines::ns(). Reprex follows: > > library(splines) > > set.seed(12345) > dat <- data.frame(claim = rbinom(1000, 1, 0.5)) > mns <- c(3.4, 3.6) > sds <- c(0.24,

We would like to announce the following statistics course; Course: Introduction to Linear mixed effects models, GLMM and MCMC with R When: 10-14 February, 2014 Where: Pousada de juventude parque das nacoes. Rua de Moscavide, Lt 47 ? 101, 1998- 011. Lisbon, Portugal Info: http://www.highstat.com/statscourse.htm Flyer: http://www.highstat.com/Courses/Flyer2014_02SIM_LisbonV2.pdf Kind regards,

We would like to announce the following statistics course: Course: Introduction to Linear Mixed Effects Models and GLMM with R Where:? Lisbon, Portugal When:?? 19-23 February 2018 Course website: http://highstat.com/index.php/courses Course flyer: http://highstat.com/Courses/Flyers/2018/Flyer2018_02LisbonV2.pdf Kind regards, Alain Zuur -- Dr. Alain F. Zuur Highland Statistics Ltd. 9 St

Dear All, I have a query relating to the use of the ?predict? and ?gamm? functions. I am dealing with large (approx. 5000) sets of presence/absence data, which I am trying to model as a function of different of environmental covariates. Ideally my models should include individual and colony as random factors. I have been trying to fit binomial models using the gamm function to achieve this. For

Hi Folks, I'm seeking confirmation of something which is probably true but which I have not managed to find in the documentation. I have a binary response y={0.1} and a variable x and have fitted a probit response to the data with f <- glm( y~x, family=binomial(link=probit) ) and then, with a specified set of x-value X I have used the predict.glm function as p <- predict( f, X,

Hello all R-user I am relative new to the R-environment and also to GLMM, so please don't be irritated if some questions don't make sense. I am using R 2.0.0 on Windows 2000. I investigated the occurrence of insects (count) in different parts of different plants (plantid) and recorded as well some characteristics of the plant parts (e.g. thickness). It is an unbalanced design with 21

Dear All, I know this may be a trivial question. In the past I have used glm to make logistic regressions on data. The output creates an object with the results of the logistic regression. This object can then be used to make predictions. Great. I have a different problem. I need to make predictions from a logistic regression taken from a paper. Thus I need to (by hand) enter the reported odds