R ggplot 无法设置对数轴限制

Question

我试图表现出一些非常扭曲数据的使用直方图分布scale_x_log10从ggplot2 。 我最终将在同一帧中有许多类似的图，因此我需要将轴设置为特定限制（ylim=0,150 和 xlim=0,25000）。 使用这里的示例，我正在尝试：

#histograms by age
PfMSP119_1_5hist <- ggplot(s1_to_5_adj, aes(x = PfMSP119_adj))+
 geom_histogram(bins = 1500) +
   ylim(c(0,150)) +
  scale_x_log10(limit = c(0,25000)) +
  xlab("MFI value") +
  ylab("Frequency") +
  labs(title = "Age 1-5") +
  theme(plot.title = element_text(hjust = 0.5)) +
  theme(panel.grid.minor=element_blank(),
        panel.grid.major=element_blank())

然而，生成的图形是空白的，没有 x 轴截断。

尝试另一种方法，我分别调用scale_x_log10和xlim ：

PfMSP119_1_5hist <- ggplot(s1_to_5_adj, aes(x = PfMSP119_adj))+
 geom_histogram(bins = 1500) +
   ylim(c(0,150)) +
  xlim(c(0,25000)) +
  scale_x_log10() +
  xlab("MFI value") +
  ylab("Frequency") +
  labs(title = "Age 1-5") +
  theme(plot.title = element_text(hjust = 0.5)) +
  theme(panel.grid.minor=element_blank(),
        panel.grid.major=element_blank())

我收到警告：

Scale for 'x' is already present. Adding another scale for 'x', which will replace the existing scale.

生成的图形上的 x 轴未设置为我指定的限制。

数据：

dput(s1_to_5_adj$PfMSP119_adj)
c(63.83391104, 57.82601353, 61.69849758, 33.1065109, 34.98756265, 
16.55325545, 35.36377301, 20.31535896, 57.18397337, 40.63071792, 
33.47425393, 37.29988296, 24.86658864, 32.51784668, 32.51784668, 
32.51784668, 31.56143942, 30.60503217, 83.20743121, 33.27774601, 
147.2540261, 91.51380154, 95.67351979, 79.0946642, 45.63153704, 
91.26307407, 54.75784444, 35.92885894, 26.47389606, 30.25588121, 
78.94894005, 27.57408649, 16.68957866, 29.75098805, 57.32507454, 
26.84845263, 20.31774794, 54.64094134, 20.58934021, 22.96503332, 
22.17313562, 24.54882872, 26.13262412, 21.38123791, 19.00554481, 
38.01108963, 60.18422524, 49.09765743, 61.76802064, 210.7428022, 
180.8148303, 159.3041005, 101.6304046, 96.64240931, 174.5798362, 
343.1479027, 1070.992956, 292.3111764, 587.5552408, 261.5158518, 
267.7668987, 358.0142608, 1363.627725, 280.659379, 158.5420101, 
247.5349309, 283.4312351, 425.5207724, 200.4210316, 456.2915647, 
48.95313475, 120.267578, 223.008725, 290.6970101, 80.98419823, 
41.23667352, 38.94574722, 54.98223137, 32.0729683, 46.96398929, 
28.63657884, 576.3861393, 392.4025337, 201.5914116, 162.4230268, 
440.554493, 191.8891511, 72.69486729, 54.33379256, 89.93179459, 
53.8284401, 309.2272091, 554.3184044, 24.0321675, 30.89850107, 
35.81701072, 81.37068455, 19.47332622, 19.47332622, 18.77785028, 
15.3004706, 25.25721271, 49.02870703, 14.11432475, 14.11432475, 
40.0154283, 15.56155545, 18.52566125, 22.97181995, 33.34619025, 
16.82050109, 19.11420579, 37.46384334, 45.10952566, 23.70161517, 
198.0231719, 16.05593286, 34.1505622, 26.76665686, 23.99769235, 
22.15171602, 25.84366869, 33.22757403, 28.61263319, 23.99769235, 
23.99769235, 91.37582858, 21.22872785, 23.07470419, 68.94089313, 
288.2721413, 96.47749148, 29.68538199, 30.74557421, 18.47224286, 
22.53613629, 23.64447086, 14.77779429, 19.21113257, 39.94087069, 
58.3751187, 41.47705803, 89.09886539, 32.98155951, 15.22225823, 
49.04949875, 21.14202533, 20.29634431, 34.41740114, 39.00638796, 
6471.028532, 43.59537478, 56.21508853, 65.39306217, 30.97566103, 
25.95917146, 30.43489068, 19.69316455, 27.74945915, 30.43489068, 
20.5883084, 68.0309321, 39.38632911, 42.07176064, 30.43489068, 
26.88069224, 43.68112489, 35.84092299, 35.84092299, 24.64063455, 
22.40057687, 45.99563655, 51.25228073, 47.3097976, 67.02221326, 
42.05315342, 40.73899237, 27.59738193, 53.88060282, 44.68147551, 
77.53550162, 31.53986506, 60.88174922, 26.63576528, 22.83065596, 
80.85857319, 21.87937863, 33.29470661, 24.73321062, 35.67289993, 
62.78430389, 19.97682396, 22.38329865, 29.03779284, 29.03779284, 
45.97650534, 22.61151537, 33.91727306, 28.26439422, 52.00648536, 
47.48418228, 38.43957613, 49.74533382, 25.47063486, 25.47063486, 
21.22552905, 61.02339601, 29.71574066, 27.59318776, 20.16425259, 
20.16425259, 39.26722874, 36.08339938, 23.24303103, 44.46492892, 
53.56002802, 40.42266266, 32.33813013, 25.26416416, 25.19439477, 
35.99199253, 59.98665422, 59.98665422, 41.99065796, 32.39279328, 
31.79292674, 19.8299787, 25.04839415, 19.8299787, 49.05310522, 
205.6055687, 45.92205595, 57.40256993, 19.8299787, 48.78881242, 
26.270899, 40.03184609, 45.03582685, 125.099519, 40.03184609, 
42.73749092, 55.20259244, 37.39530456, 43.48186637, 39.4054414, 
81.52849945, 63.86399123, 76.09326615, 74.86804452, 60.4198254, 
95.88363596, 45.97160628, 89.31626363, 229.2012942, 53.85245307, 
43.34465735, 34.1503361, 48.59855521, 61.96737844, 219.5800584, 
44.45485844, 25.59522153, 110.4635877, 40.20515013, 27, 1333.076644, 
28.38010597, 528.1853056, 160.8206005, 58.3368845, 84.5, 72.52693748, 
320.2722567, 28.38010597, 35.17141806, 65.65331372, 28.13713445, 
65.65331372, 59.79141071, 44.55046288, 24.61999264, 42.30690966, 
43.93409849, 73.22349749, 120.4119736, 58.57879799, 73.22349749, 
49.43107012, 35.83752583, 38.30907934, 34.60174908, 265.6920019, 
22.24398155, 30.27653045, 38.42451986, 51.23269315, 47.57321506, 
34.76504178, 109.7843425, 51.23269315, 153.6980794, 44.8286065, 
23.19293513, 59.14198459, 88.71297689, 80.59544959, 36.52887284, 
40.58763648, 54.50339756, 37.10869621, 47.54551702, 48.75660692, 
48.75660692, 130.6677066, 62.40845686, 48.75660692, 46.80634265, 
45.83121051, 56.55766403, 953.5460808, 52.63079018, 54.17875459, 
32.50725276, 38.69911042, 30.95928834, 35, 45, 48, 47, 79, 11707.95725, 
31, 83, 31, 18, 168, 47, 56, 31, 36, 28.68109761, 35.50993038, 
34.14416383, 35.50993038, 92.59837032, 33.55013417, 28.18211271, 
42.94417174, 33.55013417, 71.12628445, 64.41625761, 48.98319589, 
47, 24, 30, 31, 27, 21, 34, 77.5, 20, 30, 30.5, 34, 33, 196, 
47, 41, 114, 51.64431376, 45.68535448, 78.45963051, 49.657994, 
71.50751135, 51.64431376, 59.5895928, 83.42542991)

Answer 1

我认为您的问题是 log(0) 未定义（或 R 为 -Inf），因此您无法在对数转换轴上将 x 限制设置为 0 而不会出错。

我通常的解决方法是将轴限制设置为 1（因为 log(1) = 0），如下所示。

ggplot(s1_to_5_adj, aes(x = PfMSP119_adj))+
  geom_histogram(bins = 1500) +
  ylim(c(0,150)) +
  scale_x_log10(limit = c(1,25000)) +
  xlab("MFI value") +
  ylab("Frequency") +
  labs(title = "Age 1-5") +
  theme(plot.title = element_text(hjust = 0.5)) +
  theme(panel.grid.minor=element_blank(),
        panel.grid.major=element_blank())