similar to: Mixed model with multinomial distribution

Hi to all the people, I'm working with abundance data of some species, but containing too zero values, and the factors are the ones typical in a BACI experiment (Before-and-After-Control-Impact). Thus, these are two fixed factors. As the data does not holds the normality and homogeneity of variances assumptions of clasiccal ANOVA, I'm trying to fit a zero-altered model using the MCMC glmm

Thanks for your answers Stephen and Ben, I hope I am posting on the correct list now. I managed so far to run the multinomial model with random effect with the following command: MCMCglmm(fixed=cbind(Apsy,Mygl,Crle,Crru,Miag,empty) ~ habitat:trait,random=~idh(trait):mesh,family="multinomial12", data=dataA,rcov=~trait:units) (where multiple responses are different species, Habitat

> On Mar 22, 2018, at 1:31 PM, Michelle Kline <michelle.ann.kline at gmail.com> wrote: > > Hi, > > Thanks in advance for any help on this question. I'm running multinomial > models using the MCMCglmm package. The models have 5 outcome variables > (each with count data), and an additional two random effects built into the > models. The issue is that when I use

Hi David, Thanks for your comment. I haven't posted the data because they are unpublished and include human subjects so there are issues with sharing on a list serv, but I thought perhaps someone had encountered a similar problem and would already know the answer. I will reconsider whether my University's ethics approval would allow me to post the data and update the question if I think

Hi, Thanks in advance for any help on this question. I'm running multinomial models using the MCMCglmm package. The models have 5 outcome variables (each with count data), and an additional two random effects built into the models. The issue is that when I use the following code, the summary only gives me results for four of the outcome variables. Here is the code for my model: m3.random

Dear all, I would like to add mixed effects in a multinomial model and I am trying to use MCMCglmm for that. The main problem I face: my data set consits of a trapping data set, where the observation at eah trap (1 or 0 for each species) have been aggregated per traplines. Therefore we have a proportion of presence/absence for each species per trapline. ex: ID_line mesh habitat Apsy Mygl

Hello Amelie, I don't have an answer to your question, but I just wanted to point out this page I noticed recently ( http://hlplab.wordpress.com/2009/05/07/multinomial-random-effects-models-in-r/), which might be helpful. I'm also interested in figuring out how to do a multinomial glmm, so if you find out anything I'd be happy to hear more about it! Based on what I've found so

With binomial/binary responses (0|1) running MCMCglmm with family="multinomial" terminates with Error in if (nJ < 1) { : missing value where TRUE/FALSE needed with family="categorical" there are no errors I have not looked in the code, do I need format the responses TRUE/FALSE , or is this just a bug? -- H?vard Wahl Kongsg?rd

I posted a question earlier concerning a mixed model with binomial link for a split-plot experiment with repeated measures. The answer is the GLMM function in the lme4 package. Now I have to fit a mixed model in which the response variable takes the ordered values 0, 1 or 2 - representing 0, 1 or 2 out of two marked leaves per plant that have died at each census period in the same split-plot

Hi All, The data set that I have is a cluster data, and I want to run a HLM mixed model with multi-level response. Here is my data set: response: - Level (num: 1, 2, 3, 4, 5 - 5 levels) Covariates: - Type (Factor: A, B, C - 3 levels) - yr (num: 2006, 2007, ...) - Male (num: 0=not Male, 1=Male - 2 levels) - Ethnicity (Factor: A, B, H, ..., - 7 levels) - ELL (num: 0, 1, - 2

Hi all, I'm hoping I might be able to get some help with some issues specifying priors for MCMCglmm. I'm trying to fit a gaussian glmm using MCMCglmm to a data set with two (correlated) response variables. The response variables are both logit-transformed proportions (there are a few reasons why I've chosen these with gaussian error over binomal glmm, which I won't go into).

Dear R-users, I would like to fit a glmm with cross-classified random effects with the function MCMCglmm. Something along the lines: model1<-MCMCglmm(response~pred1, random=~re1+re2, data=data) where re1 and re2 should be crossed random effects. I was wondering whether you could tell me specifying cross-classified random effects in MCMCglmm requires a particular syntax? Are there any

1. (Mainly) Statistical issues are generally off topic on this list. You might want to try the r-sig-mixed-models list instead. 2. However, I think a better answer is to seek local statistical expertise in order to have an extended discussion about your research intent in order to avoid producing yet more irreproducible psychological research. Cheers, Bert Bert Gunter "The trouble with

Hi all, I previously emailed about a multinomial model, and after seeking some additional help, realized that since my response/outcome variables are not mutually exclusive, I need to use a multi-response model that is *not* multinomial. I'm now trying to figure out how to specify the priors on the multi-response model. Any help would be much appreciated. My data look like this: X

When I try to use the following code, I get the error message shown. This is quite confusing to me, insofar as family is a recognized argument for MCMCglmm. Can anyone spot an obvious glitch? model1 <-MCMCglmm(fixed = NoRRVpos ~ Year, random = ~County, family="zipoisson", data=Rabies_Project_Init)Error in MCMCglmm(fixed = NoRRVpos ~ Year, random = ~County, family =

Hi Bert, That was distinctly unhelpful, and your outward hostility to a field you obviously don't understand reveals a regrettable level of ignorance. By the way, my research is Anthropology despite my job title. Michelle On Tue, May 1, 2018 at 2:48 PM, Bert Gunter <bgunter.4567 at gmail.com> wrote: > 1. (Mainly) Statistical issues are generally off topic on this list. > You

I am attempted to examine the temporal independence of my data set and think I need an unordered multinomial logistic regression (or logit model) with repeated measures to do so. The data in question is location of chickens. Chickens could be in any one of 5 locations when a snapshot sample was taken. The locations of chickens (bird) in 8 pens (pen) were scored twice a day (AMPM) for 20 days

Dear all, I studied in tank prey fish behavior. Using the design described below (and R code), I want to test the effects of both habitat and predator (and interaction) on prey fish's vertical distribution, which was recorded (with repeated measures) as a categorical variable. I found that package mlogit might fit to my need but I don't know how to specify my complex design in the

Good afternoon, gentlemen! After several days studying and researching on categorical data (various forums with answers from the owner of the library - all incipient) how to interpret the output the function MCMCglmm, come to enlist the help of you, if someone has already worked with MCMCglmm function in the case of variables ordinal dependent. I've read and reread all the pdf's of the

Hi, I ran R (version 2.9.0) CMD build under root in Fedora (9). When it tried to remove "junk files" it removed EVERYTHING in my local account! (See below). Can anyone tell me what happened, and even more importantly if I can I restore what was lost. Panickingly, Jarrod [jarrod at localhost AManal]$ R CMD build MCMCglmm_2.05 * checking for file