R CLM Post-Hoc контрастный анализ, если вы используете 4 уровня, а не 2 уровня IV

Для своего исследования я использовал в качестве IV 4 человека (2 м, 2 ж) и 4 одиночные позы (2 закрытые и 2 открытые).

Сначала я запустил свой clm с IV Sex (имеющий m,f) и Pos (закрытый, открытый).

Я использовал контрастное кодирование для Sex и Pos

contrasts(Pos) <- -1*contr.sum(2)
contrasts(Sex) <- -1*contr.sum(2)


Clm1 <- clm(y ~Sex*Pos, data = x)

Во-вторых, я хочу пробежать четырех человек и одиночные позы.

contrasts(SPos) <- -1*contr.sum
contrasts(Per) <- -1*contr.sum

Clm1 <- clm(y ~Per*SPos, data = x)

``` 
But now SPos2, SPos3 and SPos1 may are compared to Pos4 and I do lack one row

```
attr(,"contrasts")
  [,1] [,2] [,3]
1    1    0    0
2    0    1    0
3    0    0    1
4   -1   -1   -1
Levels: 1 2 3 4 
```
I read the solution is a post hoc contrast (-1/3 -1/3 -1/3 1). But I have no idea, how to do it in r.

Im thankful for any coding examples

I found the ScheffeTest() from the  DescTools package, but still I am not skilled enough to use it.


0 ответов

Другие вопросы по тегам