search for: otter

Hello, I am performing a BACI analysis with ANOVA using the following glm: fit1<-glm(log(Cucs_m+1)~(BA*Otter)+BA+Otter+ID+Primary, data=b1) The summary(aov(fit1)) shows significance in the interaction; however, now I would like to determine what combinations of BA and Otter are significantly different (each factor has two levels). ID and PRIMARY substrates are categorical and included in the model to he...

...It is possible to produce a working model with the present day version of R so that R can now be directly compared with AD Model Builder for this type of model. I have put the results of the test together with the original Schnute and Richards paper and the working R and AD Model Builder codes on Otter's web site http://otter-rsch.ca/tresults.htm The results are that AD Model builder is roughly 1000 times faster than R for this problem. ADMB takes about 2 seconds to converge while R takes over 90 minutes. This is a simple toy example. Real fisheries models are often hundred of times m...

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 have been using this package for a split plot design along streams (my M.Sc thesis). My response data was best described by a negative binomial dis...

Dear R-users, Half a year ago we put out the R package "glmmADMB" for fitting overdispersed count data. http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html Several people who used this package have requested additional features. We now have a new version ready. The major new feature is that glmmADMB allows Bernoulli responses with logistic and probit links. In addition there is a "ranef.glmm.admb()...

I can't get the versions of sampleenc and sampledec in the manual to work back-to-back. I've tried using speex 1.0.4 and 1.1.6. What am I doing wrong? [gerv@otter]$ file in.wav in.wav: RIFF (little-endian) data, WAVE audio, Microsoft PCM, 16 bit, mono 8000 Hz Strip headers: [gerv@otter]$ sox in.wav in.raw Command as given on manual page: [gerv@otter]$ ./sampleenc in.raw | ./sampledec out.raw Re-add headers: [gerv@otter]$ sox -r 8000 -w -s out.raw out.wav...

...l software for Linux. I believe those people must feel quite confident that they are not covered (or "infected"?) by GPL. Just my $0.02... Andy _____ From: r-help-bounces at stat.math.ethz.ch on behalf of Spencer Graves Sent: Fri 5/19/2006 6:43 PM To: mschwartz at mn.rr.com Cc: otter at otter-rsch.com; r-help at stat.math.ethz.ch; Damien Joly Subject: Re: [R] How can you buy R? [Broadcast] I'd like to know what people think is the meaning of section 2.b of the GPL (http://www.gnu.org/copyleft/gpl.html#SEC1): "You must cause any work that you d...

I'm attempting to calculate a regression in R that I normally use Prism for, because the formula isn't pretty by any means. Prism presents the formula (which is in the Prism equation library as Heterologous competition with depletion, if anyone is curious) in these segments: KdCPM = KdnM*SpAct*Vol*1000 R=NS+1 S=(1+10^(X-LogKi))*KdCPM+Hot a=-1*R b=R*S+NS*Hot+BMax c = -1*Hot*(S*MS+BMax) Y

Hi all, Sorry if the subject is a bit irrelevant or hard to understand, I have no idea on what subject should I put down. My company has a samba domain controller. All users must logon before they can gain access to the fileserver and printers. My home also has a samba domain controller. I use it mainly for file sharing on my home network. Now I am planning to purchase a laptop with wireless

...venience functions, such as predict(), fitted.values(), ... are still missing from this package, but will hopefully be added in later versions (contributions/suggestions are most welcome). Other response distributions than negative binomial or Poisson could easily be added. Download 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...

...he freely available R package glmmADMB can do Adaptive Gaussian Quadrature for this type of model, since it is built using AD Model Builder's random effects module which incorporates this feature. There is now a beta version of the software for people using R on the Mac intel platform. http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html Cheers, Dave -- David A. Fournier P.O. Box 2040, Sidney, B.C. V8l 3S3 Canada Phone/FAX 250-655-3364 http://otter-rsch.com

...to use the glmmADMB software. So far this is the only response we have received. I guess it is up to the R community to decided whether Prof. Ripley speaks for all of you. Cheers, Dave -- David A. Fournier P.O. Box 2040, Sidney, B.C. V8l 3S3 Canada Phone/FAX 250-655-3364 http://otter-rsch.com

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

Hello, I'm running a simulation study of a multilevel model with binary response using the binomial probit link. It is a random intercept and random slope model. GLMMPQL and lmer fail to converge on a *significant* portion of the *generated* datasets, while MlWin gives reasonable estimates on those datasets. This is unacceptable. Does anyone has similar experiences? Regards, Roel de

...in 20 seconds or less (although the more complicated turns out to be overparameterized) Using the random effects module of AD Model Builder we have developed a shared library (Windows dll) that can be called from R via the driver function glmm.admb(). The function can be downloaded from http://otter-rsch.com/admbre/examples/nbmm/nbmm.html The two 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...

...P values for the log-likelihood with the different packages and methods. We thought it would be interesting to revisit this example using the AD Model Builder random effects module which we feel is the most stable software available for this problem at this time. http://otter-rsch.com/admodel.htm You can get Table 2 from Lesaffre et al at http://otter-rsch.com/downloads/other/lesaffre.pdf The data and more information are available from the publisher. http://www.blackwellpublishers.co.uk/rss/Volumes/Cv50p3.htm We considered three q...

...unctions for nonlinear robust regression, > > analogous to e.g. LTS? > > > > Searching google I have found > > 1, an abstract to generalize LTS for nonlinear regression > > models, see: http://smealsearch.psu.edu/1509.html > > 2, an AD-model builder, see: http://otter-rsch.com/admodel/cc1.html > > but no mention of R/S > > > > Thank you in advance > > Best regards > > Christian > > _._._._._._._._._._._._._._._._ > > C.h.i.s.t.i.a.n S.t.r.a.t.o.w.a > > V.i.e.n.n.a A.u.s.t.r.i.a > > _._._._._._._._._._....

Hello, I am using odesolve to simulate a group of people moving through time and transmitting infections to one another. In Matlab, there is a NonNegative option which tells the Matlab solver to keep the vector elements of the ODE solution non-negative at all times. What is the right way to do this in R? Thanks, Jeremy P.S., Below is a simplified version of the code I use to try to do this,

I am using the 'glmmPQL function in the 'MASS' library to fit a mixed effects logistic regression model to simulated data. I am conducting a series of simulations, and with certain simulated datasets, estimation of the random effects logistic regression model unexpectedly terminates. I receive the following error message from R: Error in lme.formula(fixed=zz + arm.long,random=~1 |

I have finally gotten one version of windows 2003 r2 x64 installed as a guest. The bad news is no network. The network tab says it''s a realtek 8139 when I know it is a intel. I tried downloading the intel drivers and they wont install says it''s the wrong hard ware. Is this common that it uses a 8139 as a generic network card? Also the network card shows The

Does someone have code, or point to a source to get it, to model repeated measures zero-inflated poisson data. The data come from a replicated field trial comparing two treatments - a control and a test treatment. Thanks in advance ------------------------------------ Subhash Chandra, DSc Senior Biometrician Primary Industries Research Victoria Department of Primary Industries 1 Ferguson Road