[R-sig-ME] ordinal regression with MCMCglmm
Kari Ruohonen
kari.ruohonen at utu.fi
Tue Apr 13 13:11:47 CEST 2010
Hi,
I am trying to figure out how to fit an ordinal regression model with
MCMCglmm. The "MCMCglmm Course notes" has a section on multinomial
models but no example of ordinal models. Suppose I have the following
data
> data
resp treat group
1 4 A 1
2 4 A 1
3 3 A 2
4 4 A 2
5 2 A 3
6 4 A 3
7 2 A 4
8 2 A 4
9 3 A 5
10 2 A 5
11 1 B 6
12 1 B 6
13 1 B 7
14 2 B 7
15 2 B 8
16 3 B 8
17 2 B 9
18 1 B 9
19 2 B 10
20 2 B 10
and the "resp" is an ordinal response, "treat" is a treatment and
"group" is membership to a group. Assume I would like to fit an ordinal
model between "resp" and "treat" by having "group" effects as random
effects. How would I specify such a model in MCMCglmm? And how would I
specify the prior distributions?
All help is greatly appreciated.
regards, Kari
More information about the R-sig-mixed-models
mailing list