-1

Tudo funciona até eu colocar o argumento paired = T, então ele dá a seguinte mensagem: Error in t.test.formula(polinizadores ~ tratamento, paired = T, alternative = "l", : cannot use 'paired' in formula method

Podem me ajudar?

1 Resposta 1

2

A mensagem de erro é bem clara. Não é possível utilizar a sintaxe (ou método) de fórmula para fazer um Teste T Pareado no R. A solução é utilizar a sintaxe tradicional, dada por

t.test(grupo_1, grupo_2, data = dados, paired = TRUE)

Como não foi fornecido um CMR para o problema, resolverei ele utilizando o conjunto de dados sleep, disponível com qualquer instalação do R.

O primeiro passo a ser tomado é colocar os dados no formado wide (largo), como especificado neste link:

library(tidyverse)
library(janitor)
#> 
#> Attaching package: 'janitor'
#> The following objects are masked from 'package:stats':
#> 
#>     chisq.test, fisher.test
sleep_wide <- 
  sleep |> 
  pivot_wider(names_from = group, values_from = extra) |> 
  clean_names()


sleep_wide |> 
  select(-id) |> 
  summarise_all(mean)
#> # A tibble: 1 × 2
#>      x1    x2
#>   <dbl> <dbl>
#> 1  0.75  2.33

Note as médias de 0.75 e 2.33 para ambos os grupos. Desta forma, a diferença x1 - x2 é igual a -1.58, como mostrado no output abaixo.

t.test(sleep_wide$x1, sleep_wide$x2, paired = TRUE)
#> 
#>  Paired t-test
#> 
#> data:  sleep_wide$x1 and sleep_wide$x2
#> t = -4.0621, df = 9, p-value = 0.002833
#> alternative hypothesis: true mean difference is not equal to 0
#> 95 percent confidence interval:
#>  -2.4598858 -0.7001142
#> sample estimates:
#> mean difference 
#>           -1.58

Além disso, como p-valor = 0.002833 < 0.05, a hipótese nula é rejeitada e há diferença significativa entre os grupos.

Mas se fizer questão de usar os dados no formato longo, é possível ajustar um modelo linear misto a estes dados, considerando group como efeito fixo e id como efeito aleatório:

library(lme4)
library(emmeans)

modelo <- lmer(extra ~ group + (1|ID), data = sleep)

contrast(emmeans(modelo, ~ group), "pairwise")
#>  contrast        estimate    SE df t.ratio p.value
#>  group1 - group2    -1.58 0.389  9  -4.062  0.0028
#> 
#> Degrees-of-freedom method: kenward-roger

Note como no output acima temos os mesmos resultados do Teste T Pareado em termos de estimativa, erro padão, graus de liberdade, estatística do teste e p-valor, mostrando que ambos os métodos são equivalentes.

Created on 2024-05-24 with reprex v2.1.0

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .