0
library(ggstatsplot)
library(MASS)
library(dplyr)

immer <- immer %>% mutate(production_difference = Y2 - Y1) %>% 
  select(Y1, Y2,production_difference)

head(immer)

t.test(immer$production_difference)

gghistostats(
  data = immer, ## data from which variable is to be taken
  x = production_difference, ## numeric variable
  test.value = 0 
)  

Por qué el confidence interval de gghistostats no coincide con el confidence interval de t.test? El p-value y el t-statistic coinciden.

Lo único que se me ocurre es q el confidence interval sea del tamaño del efecto de Cohen que está entre 0 y 1. Y q el signo negativo indique que hay una disminución. Será eso? Allí si concidiría por que gghistostats indica que el tamaño del efecto es -0.59 y el confidence interval de gghistostats es (-0.97, -0.21)

P.S. El default de gghistostats es parametric que usa la mean. En el ejemplo uso el default: https://indrajeetpatil.github.io/ggstatsplot/reference/gghistostats.html

introducir la descripción de la imagen aquí

5
  • El método t.test utiliza la distribución t de Student para calcular el intervalo de confianza para la media de una muestra cuando se desconoce la varianza poblacional. El intervalo de confianza se calcula a partir del t-statistic y del nivel de confianza deseado (por ejemplo, 95%). Por otro lado, el gráfico de gghistostats muestra un intervalo de confianza en torno a la mediana de los datos. El intervalo de confianza se calcula a partir de la distribución normal y el nivel de confianza deseado. Es posible que los intervalos de confianza no coincidan debido a estas diferencias en los métodos de
    – John
    el 29 dic. 2022 a las 6:17
  • Gracias por tu respuesta. El default de gghistostats es parametric que usa la mean: indrajeetpatil.github.io/ggstatsplot/reference/…
    – Beginner
    el 29 dic. 2022 a las 6:24
  • Además el t-statistic y el p-value del t.test coinciden con los de gghistostats.
    – Beginner
    el 29 dic. 2022 a las 6:36
  • Lo único que se me ocurre es q el confidence interval sea del tamaño del efecto de Cohen que está entre 0 y 1. Y q el signo negativo indique que hay una disminución. Será eso? Allí si concidiría por que gghistostats indica que el tamaño del efecto es -0.59 y el confidence interval de gghistostats es (-0.97, -0.21)
    – Beginner
    el 29 dic. 2022 a las 6:38
  • Estoy casi seguro q gghistostats está usando el mismo t.test. Sin embargo, espero que estés en lo correcto. De nuevo, muchas gracias por tu opinión.
    – Beginner
    el 29 dic. 2022 a las 6:43

1 respuesta 1

1

Simplemente gghistostats no retorna los intervalos de confianza de t.test() (Ver), los únicos coeficientes que aprovecha de t.test() son:

[1] "mu"          "statistic"   "df.error"   
[4] "p.value"     "method"      "alternative"

Para los intervalos de confianza aplica una de estas dos posibles funciones para calcularlos, de acuerdo al parámetro effsize.type, hedges_g() y cohens_d().

Puedes verificar que debería coincidir con alguno de estos intervalos:

effectsize::hedges_g(
  x       = immer$production_difference,
)

effectsize::cohens_d(
  x       = immer$production_difference,
)
3

Tu Respuesta

Al pulsar en “Publica tu respuesta”, muestras tu consentimiento a nuestros términos de servicio, política de privacidad y política de cookies

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.