similar to: Glimmix and glm

Hello. We are trying to install R on our network, and I wanted to ask if there is a user license agreement. I will be grateful if somebody can send me a link to it; if one exists. Thank you very much ******************************************* Antonio Paredes USDA- Center for Veterinary Biologics Biometrics Unit 510 South 17th Street, Suite 104 Ames, IA 50010 (515) 232-5785 [[alternative

Hello everyone. I would like to know if there is any R library to perform adaptive rejection sampling. Thanks in advance. ******************************************* Antonio Paredes USDA- Center for Veterinary Biologics Biometrics Unit 510 South 17th Street, Suite 104 Ames, IA 50010 (515) 232-5785 [[alternative HTML version deleted]]

Hello everyone. I want to write some likelihood functions to use with optim. For example AR-structures, I have written a couple of functions that appear to get the job done. But, since I am new to R, I wanted to ask if there are any references in this topic. Thanks in advance, ******************************************* Antonio Paredes USDA- Center for Veterinary Biologics Biometrics Unit

Hi, I am trying to fit a model with lmer in R and proc glimmix in SAS. I have simplified my code but I am surprised to see I get different results from the two softwares. My R code is : lmer(y~age_cat + (1|cat),data=fic,family=binomial(link = "logit"), NaGQ=1) My SAS code is : ods output Glimmix.Glimmix.ParameterEstimates=t_estimates; proc glimmix data=tab_psi method=laplace;

I'm fitting generalized linear mixed models to using several fixed effects (main effects and a couple of interactions) and a grouping factor (site) to explain the variation in a dichotomous response variable (family=binomial). I wanted to compare the output I obtained using PROC GLIMMIX in SAS with that obtained using lmer in R (version 2.6.1 in Windows). When using lmer I'm specifying

CENTER FOR VETERINARY BIOLOGICS UNITED STATES DEPARTMENT OF AGRICULTURE Join our team of statisticians and data scientists in 2018 We bring a broad perspective to our consulting and analytical work with scientists at the USDA National Centers for Animal Health. Our Statistics Section is looking for a good communicator who is comfortable dealing with the analytical aspects of practical

Hello everyone, I wanted to ask what will be a good editor to write R scripts in Fedora 7. Tony [[alternative HTML version deleted]]

I recently came across a flyer from REvolution Computing, and I wanted to ask if this is R going private? Tony. [[alternative HTML version deleted]]

Hello everyone, I'll to request some input on what is available for use as an R/Excel interface; any help will be appreciated. Tony. [[alternative HTML version deleted]]

I'm trying to fit a generalized mixed effects model to a data set where each subject has paired categorical responses y (so I'm trying to use a binomial logit link). There are about 183 observations and one explanatory factor x. I'm trying to fit something like: (lmer(y~x+(1|subject))) I also tried fitting the same type of model using glmmPQL from MASS. In both cases, I get a

Dear R users, I've been trying to find an R package which does the PL estimation of spatial GLMMs especially with the negative binomial model. so it would be something similar to the "proc GLIMMIX" with the PL method in SAS. I've looked up some possible packages related to GLMMs, but it doesn't seem to be anyone using the PL estimation. Thanks for your help! Fei He UCR

Hello everyone. Currently I am running R 2.2.1 (windows), and I will like to update to 2.3. I wanted to ask if it is possible to update without having to removed 2.2.1; or do I first need to delete 2.2.1? Thank you very much. Tony [[alternative HTML version deleted]]

Hello everyone, I will like to know if there is an R package to compute exact confidence intervals for the ratio of two binomial proportions. Tony. [[alternative HTML version deleted]]

Dear R Users, May be this message should be directy send to Douglas Bates ... I just want to know if I can use the AIC value given in the output of an lmer model to classify my logistic models. I heard that the AIC value given in GLIMMIX output (SAS) is false because it come from a calculation based on pseudo-likelyhood. Is it the same for lmer ??? thanks, Arnaud Arnaud MOSNIER Biologiste

Dear I am analysing my data wit a mixed model. I used SAS but I want to redo the same analysis in R. Here the SAS code and what I wrote in R. It seems to work but the results are not the same. I don't know how to specify the class variable in R or specify the variance matrix. Can you please help me? Thanks Jurgen ## SAS: proc glimmix data=trend method=RSPL; class pid; model mdrfinal

Dear R-users, I'm new to R, so my apologies if this question doesn't make sense. I've tried the following model in lmer, and it works perfectly: model<-lmer(aphids~densroot+zone+(1|zone/site), family=quasipoisson) But if I try the exact same model with a different variable, totmas, the model looks as follows: model<-lmer(aphids~totmas+zone+(1|zone/site), family=quasipoisson)

Dear R-users, I am currently trying to switch from SAS to R, and am not very familiar with R yet, so forgive me if this question is irrelevant. If I try to find the significance of the fixed factor "spikes" in a generalized linear mixed model, with "site" nested within "zone" as a random factor, I compare following two models with the anova function:

Statistical Programmer - USDA Center for Veterinary Biologics The USDA Center for Veterinary Biologics is seeking a statistical programmer with good R skills. The position is a two-year term position in the CVB Statistics Section with the possibility of becoming permanent. Applicants must be United States citizens. The position calls for a thorough knowledge of R and aptitude in both programming

Hi All, I am making the switch to R and uncertain which of the several packages for mixed models is appropriate for my analysis. I am waiting for Pinheiro and Bates' book to arrive via inter-library loan, but it will be a week or more before it arrives. I am trying to fit a generalized linear mixed model of survival data (successes/trials) as a function of several categorical fixed and

Statistical Programmer - USDA Center for Veterinary Biologics The USDA Center for Veterinary Biologics is seeking a statistical programmer with good R skills. The position is a two-year term position in the CVB Statistics Section with the possibility of becoming permanent. Applicants must be United States citizens. The position calls for a thorough knowledge of R and aptitude in both programming