Início / coding / Perguntas / 79307500
Accepted
Daniel Lobo
Asked: 2024-12-25 19:47:13 +0800 CST

Definindo limites para o eixo x e eixo y para o gráfico QQ

  • 772

Eu tenho abaixoQ-Q plot

library(ggplot2)
ggplot(data.frame('Val' = 
c(-0.6434733395, -0.0000935279, -1.2176148541, -0.4817091973, 1.1364799871, -0.8810002309, -0.8383323744, -0.7285213185, 1.6299068100, 1.0358721646, -1.5631088786, -0.6265271885,-1.1331727142, -0.7129360871,
1.6266852445, -0.3328868518, -0.9638501136, 0.1088177257, -0.6034307366, -0.3342069674, 0.0932110885, 0.1088177257, 1.9742990366, 0.6584310789, 0.7201386429, 0.9314809258, 1.1455416757)), aes(sample = Val)) + 
stat_qq() + stat_qq_line() + scale_x_continuous(limits = c(-3, 3)) + scale_y_continuous(limits = c(-3, 3))

Aqui eu coloquei os limites de x axis and y axiscomo c(-3, 3)manualmente. Se eu tivesse colocado c(-2, 2), já que with c(-3, 3) Q-Q plotnão ocupa a janela de plotagem inteira, então a linha diagonal não seria visível.

Então, eu queria saber se há alguma maneira de colocar esse limite (simétrico) programaticamente, com base apenas nos dados, para que meu gráfico ocupe toda a janela do gráfico, independentemente dos dados, e ainda veja a linha diagonal desenhada corretamente?

1 respostas

  1. Best Answer

    sim! Você pode calcular quantis teóricos e reais usando qqnorm()e então encontrar o valor absoluto máximo entre todos os quantis. Arredonde e então use como use este valor simetricamente para ambos os eixos:

    library(ggplot2)

# your data
data <- data.frame('Val' = c(-0.6434733395, -0.0000935279, -1.2176148541, -0.4817091973, 
                             1.1364799871, -0.8810002309, -0.8383323744, -0.7285213185, 
                             1.6299068100, 1.0358721646, -1.5631088786, -0.6265271885,
                             -1.1331727142, -0.7129360871, 1.6266852445, -0.3328868518, 
                             -0.9638501136, 0.1088177257, -0.6034307366, -0.3342069674, 
                             0.0932110885, 0.1088177257, 1.9742990366, 0.6584310789, 
                             0.7201386429, 0.9314809258, 1.1455416757))


theoretical_q <- qqnorm(data$Val, plot.it = FALSE)$x # Calculate theoretical quantiles for the data

actual_q <- qqnorm(data$Val, plot.it = FALSE)$y # Get actual data quantiles


max_abs_value <- max(abs(c(theoretical_q, actual_q))) # Find the maximum absolute value among both theoretical and actual quantiles
# Round up to the next 0.5 to get a nice limit
limit <- ceiling(max_abs_value * 2) / 2

# plot with calculated symmetric limits
ggplot(data, aes(sample = Val)) + 
  stat_qq() + 
  stat_qq_line() + 
  scale_x_continuous(limits = c(-limit, limit)) + 
  scale_y_continuous(limits = c(-limit, limit))

    Fora: fora

    Ou simplesmente deixe a escala de lado completamente:

    ggplot(data, aes(sample = Val)) + 
  stat_qq() + 
  stat_qq_line()

    sem escala

    • 0

relate perguntas