Displaying 3 results from an estimated 3 matches for "multinomial12".
Did you mean:
multinomial
2012 Sep 12
7
multinomial MCMCglmm
...0
15 028S6E 28 hedge 7 0 0 1 0 0 0 0 0
0 0
empty
11 28
12 28
13 24
14 21
15 22
When I run the following:
> test1 <-
MCMCglmm(fixed=cbind(Apsy,Mygl,Crle,Crru,Miag,Miar,Mimi,Mumu,Misu,Soar,Somi,empty)~habitat,random=~mesh,family="multinomial12",data=metalSmA[,c(2,9,23:34)],rcov=~us(trait):units)
I got some error when running regarding the variance structure:
> "ill-conditioned G/R structure: use proper priors if you haven't or
rescale data if you have"
I guess that the problem comes from the nature of my observ...
2012 Nov 06
1
Multinomial MCMCglmm
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 the explanatory
variable and Mesh the random effect)
The main question I am facing now is:
Why the multinomial model fit does not provide by default a parameter
estimate for each response category but a...
2012 Sep 12
0
R-help Digest, Vol 115, Issue 12
...0 0 0
> 0 0
> empty
> 11 28
> 12 28
> 13 24
> 14 21
> 15 22
>
> When I run the following:
>
> > test1 <-
>
> MCMCglmm(fixed=cbind(Apsy,Mygl,Crle,Crru,Miag,Miar,Mimi,Mumu,Misu,Soar,Somi,empty)~habitat,random=~mesh,family="multinomial12",data=metalSmA[,c(2,9,23:34)],rcov=~us(trait):units)
>
> I got some error when running regarding the variance structure:
>
> > "ill-conditioned G/R structure: use proper priors if you haven't or
> rescale data if you have"
>
> I guess that the problem com...