Hi everyone,

After following up Alessandro's suggestion of using glmm in the repeated
package it seems that we had our wires crossed.

I'm trying to model some continuous data with a left skew. Originally I
thought that perhaps there was a package that could fit a loglogistic type
model. Alessandro's suggestion using glmm need a binomial which doesnt suit
my data.

How can I model left skewed data in a mixed model?

Thanks for any suggestions in advance,

Alan Haynes




--------------------------------------------------
Email: aghaynes@gmail.com
Mobile: +41794385586
Skype: aghaynes


On Tue, Dec 6, 2011 at 11:36 AM, Alessandro Moscatelli
<moskante@gmail.com>wrote:

> The non-CRAN package "repeated" (function: glmm) allows running a mixed
> model with log log link function. However, the function just fit models
> having a single random effect (the random intercept).
>
> You can download it here:
> http://www.commanster.eu/rcode.html
>
> Alessandro Moscatelli
>
> 2011/12/6 <r-sig-mixed-models-request@r-project.org>
>
> > Send R-sig-mixed-models mailing list submissions to
> >        r-sig-mixed-models@r-project.org
> >
> > To subscribe or unsubscribe via the World Wide Web, visit
> >        https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
> > or, via email, send a message with subject or body 'help' to
> >        r-sig-mixed-models-request@r-project.org
> >
> > You can reach the person managing the list at
> >        r-sig-mixed-models-owner@r-project.org
> >
> > When replying, please edit your Subject line so it is more specific
> > than "Re: Contents of R-sig-mixed-models digest..."
> >
> >
> > Today's Topics:
> >
> >   1. Mixed models with loglog distibution (Alan Haynes)
> >
> >
> > ----------------------------------------------------------------------
> >
> > Message: 1
> > Date: Tue, 6 Dec 2011 11:09:52 +0100
> > From: Alan Haynes <aghaynes@gmail.com>
> > To: r-sig-mixed-models@r-project.org
> > Subject: [R-sig-ME] Mixed models with loglog distibution
> > Message-ID:
> >        <
> CAPdSD+5cgoE4PyOcxmNvR3STsHh4+nzXzymzKm3htNTo8+EpPA@mail.gmail.com
> > >
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > Dear list members,
> >
> > I am currently trying to analyse some data from an decomposition
> experiment
> > where we exclude herbivores and consider the effect of that on
> > decomposition rates.
> >
> > Unfortunately the data produces residuals which are not normally
> > distributed and appear to correspond with a loglog distribution.
> > I have attached some graphical diagnostic plots of the residuals for your
> > perusal (distribution plots of the original data are very similar).
> >
> > This leads me to my question....Is there some way of running a mixed
> model
> > with a loglog distribution?
> > I have found that the 'reliaR' package has the usual d*foo*, r*foo* etc
> > functions but specifying family = Loglog doesnt seem to work
> >
> > Any suggestions?
> >
> >
> > Many thanks in advance,
> >
> > Alan Haynes
> >
> >
> >
> >
> >
> >
> >
> > --------------------------------------------------
> > Email: aghaynes@gmail.com
> > Mobile: +41794385586
> > Skype: aghaynes
> > -------------- next part --------------
> > A non-text attachment was scrubbed...
> > Name: model_residuals.pdf
> > Type: application/pdf
> > Size: 82132 bytes
> > Desc: not available
> > URL: <
> >
> https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20111206/3552a687/attachment.pdf
> > >
> >
> > ------------------------------
> >
> > _______________________________________________
> > R-sig-mixed-models mailing list
> > R-sig-mixed-models@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
> >
> >
> > End of R-sig-mixed-models Digest, Vol 60, Issue 5
> > *************************************************
> >
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

	[[alternative HTML version deleted]]


