R version 3.1.1 (2014-07-10) -- "Sock it to Me"
> bl <- c(140, 138, 150, 148, 135)
> fu <- c(138, 136, 148, 146, 133)
> t.test(fu, bl, alternative = "two.sided", paired = TRUE)
Error in t.test.default(fu, bl, alternative = "two.sided", paired = TRUE) :
data are essentially constant
Depois, altero apenas um caractere no meu conjunto de dados fu:
> fu <- c(138, 136, 148, 146, 132)
e corre ...
> t.test(fu, bl, alternative = "two.sided", paired = TRUE)
Paired t-test
O que estou perdendo aqui?
bl-fu
. Agorasd(bl-fu)
. Se não é óbvio, ainda, fazer estes:dif=bl-fu
em seguida,n=length(dif)
em seguida,mean(dif)/(sd(dif)/sqrt(n))
... você vê agora?Respostas:
Conforme abordado nos comentários, a questão era que as diferenças eram todas 2 (ou -2, dependendo de como você escrever os pares).
Respondendo à pergunta nos comentários:
Bem, isso depende.
Se a distribuição das diferenças realmente fosse normal, essa seria a conclusão, mas pode ser que a suposição de normalidade esteja errada e a distribuição das diferenças nas medições seja realmente discreta (talvez na população que você deseja inferir sobre isso geralmente seja -2) mas ocasionalmente diferente de -2).
De fato, vendo que todos os números são inteiros, parece que talvez seja o caso de discrição.
... nesse caso, não há tanta certeza de que todas as diferenças serão -2 na população - é mais que falta uma evidência na amostra de uma diferença na população significa algo diferente de -2.
(Por exemplo, se 87% das diferenças da população forem -2, há apenas 50-50 chances de que qualquer uma das 5 diferenças da amostra seja diferente de -2. Portanto, a amostra é bastante consistente com a variação de -2 na população)
Mas você também seria questionado sobre a adequação das suposições para o teste t - especialmente em uma amostra tão pequena.
fonte