Como comparar duas amostras de dados com o teste t de R

Se você deseja usar R de t.test () função para comparar os seus dados, você primeiro tem que verificar, entre outras coisas, se ambas as amostras são normalmente distribuídos.

Se você quer saber se a temperatura média difere entre os períodos, o castor é ativo e inativo, você pode fazê-lo com um simples comando:

> T.test (temperatura ~ activ, data = beaver2) Welch Two-Sample t-testdata: temp por activt = -18,5479, df = 80,852, p-valor lt; 2.2e-16alternative hipótese: a verdadeira diferença entre as médias não é igual ao intervalo de confiança de 095 por cento: -0.8927106 estimativas -0.7197342sample: média no grupo 0 significa em grupo 137.0968437.90306

Normalmente, você só pode realizar um teste t de amostras para as quais as variações são aproximadamente iguais. R usa variação de Welch sobre o t-teste, que corrige para variâncias desiguais.

Você ganha um monte de informação aqui:

• A segunda linha dá-lhe a estatística de teste (t para este teste), os graus de liberdade (df), E o valor de p de acordo com. O muito pequeno valor de p indica que as médias das duas amostras diferem significativamente.

• A hipótese alternativa diz o que você pode concluir se o p-valor é inferior ao limite de significância. Geralmente, os cientistas consideram a hipótese alternativa para ser verdade se o p-valor é menor do que 0,05.

• O intervalo de confiança de 95 por cento é o intervalo que contém a diferença entre os meios com 95 por cento de probabilidade, de modo que, neste caso, a diferença entre os meios reside provavelmente entre 0,72 e 0,89.

• A última linha dá-lhe os meios de ambas as amostras.

Você leu a fórmula Temp ~ activ Como # 147-avaliar Temp dentro dos grupos determinados pela activ.# 148- Como alternativa, você pode usar dois vectores separados para as amostras que você deseja comparar e passar tanto para a função, como no exemplo a seguir:

> activetemp lt; - beaver2 $ temp [beaver2 $ activ == 1]> inactivetemp lt; - beaver2 $ temp [beaver2 $ activ == 0]> t.test (activetemp, inactivetemp)