Como realizar o teste de humor mediano em r

O teste de mediana de humor é utilizado para comparar as medianas de dois ou mais grupos independentes.

A função median_test da biblioteca de moedas pode ser usada para realizar este teste em R, que usa a seguinte sintaxe:

median_test(resposta~grupo, dados)

Ouro:

• resposta: um vetor de valores de resposta
• grupo: um vetor de valores de agrupamento
• dados: um quadro de dados contendo os vetores de resposta e de grupo

O exemplo a seguir ilustra como usar esta função para realizar o teste de humor mediano em R.

Exemplo: teste de mediana de humor em R

Suponha que um professor queira saber se dois métodos de estudo diferentes produzem ou não notas diferentes em testes entre os alunos de sua turma. Para testar isso, ela pede aleatoriamente a 10 alunos que usem um método de estudo e a outros 10 alunos que usem outro. Após duas semanas, cada aluno faz o mesmo exame.

Ela decide usar o teste de mediana de Mood para determinar se a pontuação média do exame difere entre os dois grupos.

Etapa 1: Crie o quadro de dados.

 #createdata
method = rep(c('method1', 'method2'), each=10)
score = c(75, 77, 78, 83, 83, 85, 89, 90, 91, 97, 77, 80, 84, 84, 85, 90, 92, 92, 94, 95)
examData = data.frame(method, score)

#viewdata
examData

    method score
1 method1 75
2 method1 77
3 method1 78
4 method1 83
5 method1 83
6 method1 85
7 method1 89
8 method1 90
9 method1 91
10 method1 97
11 method2 77
12 method2 80
13 method2 84
14 method2 84
15 method2 85
16 method2 90
17 method2 92
18 method2 92
19 method2 94
20 method2 95

Etapa 2: execute o teste de humor mediano.

 #load the coin library
library(corner)

#perform Mood's Median Test
median_test(score~method, data = examData)

#output
	Asymptotic Two-Sample Brown-Mood Median Test

data: score by method (method1, method2)
Z = -0.43809, p-value = 0.6613
alternative hypothesis: true mu is not equal to 0

O valor p do teste é 0,6613 . Como este valor não é inferior a 0,05, não rejeitamos a hipótese nula. Não temos evidências suficientes para afirmar que existe uma diferença estatisticamente significativa nas notas medianas dos exames entre os dois grupos.

Por padrão, esta função atribui uma pontuação 0 às observações exatamente iguais à mediana. No entanto, você pode especificar esse valor como 0,5 ou 1 usando o argumento mid.score .

Por exemplo, o código a seguir executa exatamente o mesmo teste de humor da mediana, mas atribui um valor de 0,5 às observações iguais à mediana:

 #perform Mood's Median Test
median_test(score~method, mid.score="0.5" , data = examData)

#output
	Asymptotic Two-Sample Brown-Mood Median Test

data: score by method (method1, method2)
Z = -0.45947, p-value = 0.6459
alternative hypothesis: true mu is not equal to 00

O valor p do teste acabou sendo 0,6459 , que é ligeiramente inferior ao valor p anterior de 0,6613 . No entanto, a conclusão do teste ainda é a mesma: não temos evidências suficientes para afirmar que as notas medianas dos exames entre os dois grupos sejam significativamente diferentes.