主页 / coding / 问题 / 79307500
Accepted
Daniel Lobo
Asked: 2024-12-25 19:47:13 +0800 CST

设置 QQ 图的 x 轴和 y 轴的限制

  • 772

我有以下Q-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))

x axis and y axis这里我手动设置了 的界限c(-3, 3)。如果我设置了c(-2, 2)，由于c(-3, 3) Q-Q plot不占据整个绘图窗口，所以对角线将不可见。

所以我想知道是否有任何方法可以以编程方式设置这样的限制（对称），仅基于数据，以便我的绘图可以占据整个绘图窗口而不管数据如何，并且仍然可以看到正确绘制的对角线？

1 个回答

  1. Best Answer

    是的！您可以使用计算理论和实际分位数qqnorm()，然后找出所有分位数中的最大绝对值。将其四舍五入，然后将其用作对称地在两个轴上使用此值：

    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))

    出去： 出去

    或者干脆完全不进行缩放：

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

    无尺度

    • 0

相关问题