Search Archives

> [In] all our territories, all nests/"no-nests" of a territory are either occupied or unoccupied. Given your controls are unoccupied trees from within each territory, did a conditional logistic regression work (with territory as stratum indicator)? You need to...

On Mon, 21 Jan 2013, Andreas Karpf wrote: > I searched everywhere but I didn't find what I want, that is why I as the > question here. Does anybody know of a function in R which allows to > estimate ordered...

On Wed, 13 Jun 2012, Yi (Charlie) Chen wrote: > I'm trying to do the Bayesian Inference on variance components of mixed > effect models, is there any suggestion about the package I can start > with? Why not look through the...

On Tue, 15 Dec 2015, Simone Bianchi wrote: > I was suggested to use the package glmmADMB for zero-inflated or hurdle > mixed models. The packages is not available for R 3.2.2. Is there any other > resource I can...

On Sun, 27 Jan 2013, Iasonas Lamprianou wrote: > Dear all, I am running an lmer model using repeated measures (scores on > various school subjects) per student as a dependent variable. The > students are a random effect. My problem is that...

Justin Bem wrote: > I got this message, while using the polr function in MASS > > Error in polr(EQ, don, subset = (cote != 0), method = "probit", na.action = na.omit) : > attempt for find suitable starting values failed > > how can I initialise starting...

Another approach to digest *simple* tables is via lynx -dump <table url> Dumb plotting is available using macanova, which is also related to early versions of S (eg plot(x,y,dumb:T)). -- | David Duffy (MBBS PhD) ,-_|\ | email: davidD at...

I couldn't resist adding a more literal answer unback <- function(x) { chars <- unlist(strsplit(deparse(x),"")) chars <- chars[-c(1,length(chars))] paste(gsub("\\\\","/",chars),collapse="") } unback("\n") | David Duffy (MBBS PhD) ,-_|\ | email: davidD at qimr.edu.au ph...

> Thing is: there is nothing like this for GLMMs. The lowest level of "residual variance" is basically the > distribution variance. > You could think there would be such a thing with a threshold model, but it turns out that the total...

On Mon, 8 Jun 2015, Hans Ekbrand wrote: > If you can solve the problem with nlme() or anything else for that > matter, that's perfectly fine, gamm4 should do this. | David Duffy (MBBS PhD) | email: David.Duffy at qimrberghofer.edu...

> (a) should I use lmer routine? Yes > If so, is it sufficient to call it as lmer(y ~ x1+x2+x3 + (1|Id)? As a start > (b) How can I apply Friedman's ?Multivariate Adaptive Regression Splines? > to such problems...

This is nothing to with mixed modelling per se, but how your contrasts for the fixed effects are set up. If this was fixed effects only, you could fit -1 + factor (no intercept) to give estimates for each level of...

On Wed, 31 Jul 2013, Kristofor Voss wrote: > So when I run the model I first ran (i.e. the one you said is correct), I don't use lme very often, but I believe that if you use weights...

> The modeling question regards my depiction of time. The amount of elapsed time > between participants and testing points is not equal equal. Students in the same classroom > will be measured at equal time points, but while one class might have...

If you want to double check your results, consider making the full haplodiploid kinship matrix and running GenABEL::polygenic_hglm (or hglm directly). | David Duffy (MBBS PhD) ,-_|\ | email: davidD at qimr.edu.au ph: INT+61+7+3362-0217 fax...

Paul Meagher <paul at datavore.com> wrote: > > I am looking at developing a user modelling type app (new data points coming > in and wanting to dynamically update regression co-efficients for each user) > which could be viewed as a type...

On Wed, 28 Nov 2012, Paul Johnson wrote: > This might interest some of the more genetic mixed modellers... > > http://www.ncbi.nlm.nih.gov/pubmed/22398122 http://arxiv.org/abs/1112.5006 | David Duffy (MBBS PhD) ,-_|\ | email: davidD at qimr...

Of course the really annoying thing is that this particular problem runs fine in lme() ;) David Duffy. (PS. though not using other packages eg Karin Meyer's WOMBAT using PX-EM or AIREML) -- | David Duffy (MBBS PhD) ,-_|\ | email: davidD at...

Maybe this is old hat, but folks might have seen https://amstat.tandfonline.com/doi/abs/10.1198/jcgs.2009.06127?journalCode=ucgs20 where they fit poisson and binomial GLMMs for 33M records. Just have to develop and code a...

Only traits 2, 7 and 8 seem to have much relationship with tim, and there are a few outliers in there. Did the multitrait model run OK when you just use the mean trait value across occasions? | David Duffy (MBBS...