Van dos opciones, una usando una librería y otra a mano (derivada del código en esa librería). El problema es que no me coinciden los estadísticos que calculo usando la distribución con lo que obtienes en tu prueba, ver al final.
sjPlot
En este caso la librería sjPlot::
tiene la función dist_t()
toma como argumentos los grados de libertad, el valor t
de la prueba y (opcionalmente) el nivel alfa de la decisión (criterio de corte del p-value). Grafica solamente pruebas en una sola cola.
En tu caso podrías usar:
dist_t(deg.f = 3430, #Los grados de libertad de la distribución
p = 1) + #Para pintar toda el área de rechazo
geom_vline(xintercept = 1.15) #Para marcar el valor t que salió de tu test

Es muy simple y esa librería permite hacer lo mismo con diferentes distribuciones. Eso sí, tienes unas dependencias muy pesadas (carga muchísimos otros paquetes en la instalación).
A mano
Si buscas una solución más personalizada https://github.com/strengejacke/sjPlot/blob/master/R/sjPlotDist.R está el código de la función, que podrías adaptar a tus necesidades (sobre todo el cálculo de los cuantiles)
En este caso se grafica una prueba t en ambas colas de las distribución.
De hecho revisando ese código hice esta adaptación al problema:
# Crear los datos de una dist t con df grados de libertad
library(tidyverse)
df <- 3430.8 #Grados de libertad en la distribución T
t_max <- stats::qt(p = 0.00001, df = df, lower.tail = F) #Valor de T para un evento muy poco probable. Indica hasta dónde vamos a dibujar la curva.
t_alfa <- 1-qt(p = 0.05, df = df) # El valor T para el p-value que usamos como alfa, asumo que 0.05
p_value <- pt(q = 1.1584, df = df, lower.tail = F)*2 #Solo para verificar que corresponde con el resultado de la prueba
t_observado <- 1.1584
para_t <- tibble(
x = seq(-t_max, t_max, length.out = 1000 * t_max), #Creo una secuencia usando el intervalo de mínimo y máximo
y = dt(x, df= df), #Encuentro la densidad en la distribución T para esa secuencia
alfa = ifelse(x >= t_alfa | x <= -t_alfa, "Es posible rechazar H0", #Esto es para colorear después
"No es posible rechazar H0"))
ggplot(para_t, aes(x = x, y = y )) +
geom_ribbon(aes(ymax = y, ymin = 0, fill = alfa)) +
geom_vline(xintercept = t_observado) + #Dibujo una línea en el valor T observado
geom_text(x = t_observado-0.1, #Lo muevo para que no se superponga con la línea
y = 0.2, #Como a la mitad de y, ajustar a ojo
label = "T observado = 1.15",
angle = 90) +
theme(legend.position = "bottom")
