search for: skaug

...load site: http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html The package is based on the software ADMB-RE, but the full unrestricted R-package is made freely available by Otter Research Ltd and does not require ADMB-RE to run with user supplied data. I you will find this useful, Hans Skaug -- Hans Skaug Department of Mathematics University of Bergen, Norway email: skaug at mi.uib.no ph. (+47) 55 58 48 61

...-RE to run. Versions for Linux and Windows are available. We are still happy to get feedback for users, and to get suggestions for improvement. We have set up a forum at http://www.otter-rsch.ca/phpbb/ for discussions about the software. Regards, Hans _____________________________ Hans Julius Skaug Department of Mathematics University of Bergen Johannes Brunsgate 12 5008 Bergen Norway ph. (+47) 55 58 48 61

...ran code (e.g. >http://www.autodiff.org). One of the difficulties to use them with R is the >need to translate the R code into C of Fortran code, but probably there are >many other problems that I'm not able to see. > >Alberto >-- _____________________________ Hans Julius Skaug Department of Mathematics University of Bergen Johannes Brunsgate 12 5008 Bergen Norway

Dear Dr Skaug and R users, just discovered glmm.admb in R, and it seems a very useful tool. However, I ran into a problem when I compare two models: m1<-glmm.admb(survival~light*species*damage, random=~1, group="table", data=bm, family="binomial", link="logit") m1.1<-glm...

...o models of Booth et al are fit by the commands: glmm.admb(y~Base*trt+Age+Visit,random=~1,group="subject",data=epil2) glmm.admb(y~Base*trt+Age+Visit,random=~Visit,group="subject",data=epil2) I will be happy to receive feedback on the function glmm.admb(). Best regards, Hans Skaug Reference: Booth J.G.; Casella G.; Friedl H.; Hobert J.P, Negative binomial loglinear mixed models. Statistical Modelling, October 2003, vol. 3, no. 3, pp. 179-191

...be demanding since he claims that he wants the program to run on his computer which it apparently does not do at present. Prof. Ripley seems to feel that he is a qualified spokesman for the open source community. I have no idea what the community at large feels about this. What follows is Hans Skaug's post with Prof. Ripley's reply. > On Mon, 22 May 2006, H. Skaug wrote: > >> > About glmmADMB and GPL: >> > >> > We were not very cautious when we put in the GPL statement. >> > What we wanted to say was that the use of glmmADMB is free, and...

7 Jan 09 Hello, I am using R version 2.7.0 in a Windows XP context. I am also using the glmm.admb package (created by Dave Fournier, Hans Skaug, and Anders Nielson) to run mixed-effects negative binomial models. To the best of my knowledge and ability, I have searched and studied the R-help, R-sig-mixed models, and ADMB NBMM for R (through Otter Research Ltd) list servs; R help documentation; and the web in general for information...

...w to integrate AD Model Builder with R, and are going to the JSM in Seattle, the following course may be of interest: Statistical and ecological modeling in AD Model Builder Date: August 5, 2006 (the day before the Joint Statistical Meetings open) For more information see: http://www.mi.uib.no/~skaug/admb/seattleworkshop.html Regards, hans

Dear list, Does there exist R routines for reading output files from Genepop? (GENEPOP is a population genetics software package by Raymond & Rousset; http://genepop.curtin.edu.au/ ) I find several R packages that contain function for writing Genepop input files, but non that does the reverse. Regards, hans [[alternative HTML version deleted]]

the code is list below , i can run this code for laplace approximation but not importance sampling. the code work properly if use impSamp=0 , but error for any other number >0 . could you guys help me out. thanks so much ! source("http://www.mi.uib.no/~skaug/cash/lesaffre_dat.s") library(glmmADMB) fit_glmmADMB_ip <- glmm.admb(y~treat,random=~1,group="subject",impSamp =100,data=lesaffre,family="binomial",link="logit",save.dir= NULL) Initial statistics: 2 variables; iteration 0; function evaluation 0 Function va...

...as _glmm_ADMB_temp_dir_ when programing is running , but this fold mysteriously disappeared after the programing done . searched the whole computer desn't find any clue. I posted my output below. I really appreciate if you can help me with the questions. source("http://www.mi.uib.no/~skaug/cash/lesaffre_dat.s <http://www.mi.uib.no/%7Eskaug/cash/lesaffre_dat.s>") fit_glmmADMB_ip <- glmm.admb(y~treat,random=~1,group="subject",impSamp =100,data=lesaffre,family="binomial",link="logit",save.dir= NULL) Initial statistics: 2 variables; iteratio...

Hi I am running an lmer which works fine with family=poisson mixed.model<-lmer(nobees~spray+dist+flwabund+flwdiv+round+(1|field),family="poisson",method="ML",na.action=na.omit) But it is overdispersed. I tried using family=quasipoisson but get no P values. This didnt worry me too much as i think my data is closer to negative binomial but i cant find any examples of

Do Users of Nonlinear Mixed Effects Models Know Whether Their Software Really Works? Lesaffre et. al. (Appl. Statist. (2001) 50, Part3, pp 325-335) analyzed some simple clinical trials data using a logistic random effects model. Several packages and methods MIXOR, SAS NLMIXED were employed. They reported obtaining very different parameter estimates and P