ggplot2에서 facet_wrap 및 scales = "free"로 개별 축 제한 설정


116

예측 된 값과 잔차의 플롯과 함께 예측 된 값과 실제 값을 나란히보기 위해 패싯 플롯을 만들고 있습니다. shiny다른 훈련 매개 변수를 사용하여 모델링 노력의 결과를 탐색하는 데 도움 이 될 것 입니다. 85 %의 데이터로 모델을 훈련시키고 나머지 15 %에 대해 테스트 한 다음이를 5 번 반복하여 매번 실제 / 예측 값을 수집합니다. 잔차를 계산 한 후 data.frame다음과 같이 보입니다.

head(results)
       act     pred       resid
2 52.81000 52.86750 -0.05750133
3 44.46000 42.76825  1.69175252
4 54.58667 49.00482  5.58184181
5 36.23333 35.52386  0.70947731
6 53.22667 48.79429  4.43237981
7 41.72333 41.57504  0.14829173

내가 원하는 것 :

  • 측면 음모에 의해 사이드 predactpredresid
  • x / y 범위 / 한계 predact동일, 이상적으로는 min(min(results$act), min(results$pred))~까지max(max(results$act), max(results$pred))
  • x / y 범위 / 한계 predvs. 내가 실제 vs. 예측 플롯에 수행하는 작업에 영향을받지 resid 않습니다 . x예측 된 값과 y잔차 범위에 대한 플로팅 은 괜찮습니다.

두 플롯을 나란히보기 위해 데이터를 녹입니다.

library(reshape2)
plot <- melt(results, id.vars = "pred")

이제 플롯 :

library(ggplot2)
p <- ggplot(plot, aes(x = pred, y = value)) + geom_point(size = 2.5) + theme_bw()
p <- p + facet_wrap(~variable, scales = "free")

print(p)

그것은 내가 원하는 것에 매우 가깝습니다.

여기에 이미지 설명 입력

내가 원하는 것은 실제 대 예측에 대한 x 및 y 범위가 동일 할 것으로 예상되는 것이지만이를 지정하는 방법을 잘 모르겠 으며 예측 대 잔차 플롯에 대해 수행 할 필요 가 없습니다 . 범위는 완전히 다릅니다.

나는 모두이 같은 추가 시도 scale_x_continousscale_y_continuous:

min_xy <- min(min(plot$pred), min(plot$value))
max_xy <- max(max(plot$pred), max(plot$value))

p <- ggplot(plot, aes(x = pred, y = value)) + geom_point(size = 2.5) + theme_bw()
p <- p + facet_wrap(~variable, scales = "free")
p <- p + scale_x_continuous(limits = c(min_xy, max_xy))
p <- p + scale_y_continuous(limits = c(min_xy, max_xy))

print(p)

그러나 그것은 min()잔여 가치 를 선택 합니다.

여기에 이미지 설명 입력

내가 가진 마지막 아이디어 중 하나는 녹이기 전에 최소값 actpred변수 의 값을 저장 한 다음 녹은 데이터 프레임에 추가하여 표시되는 패싯을 지정하는 것입니다.

head(results)
       act     pred       resid
2 52.81000 52.86750 -0.05750133
3 44.46000 42.76825  1.69175252
4 54.58667 49.00482  5.58184181
5 36.23333 35.52386  0.70947731

min_xy <- min(min(results$act), min(results$pred))
max_xy <- max(max(results$act), max(results$pred))

plot <- melt(results, id.vars = "pred")

plot <- rbind(plot, data.frame(pred = c(min_xy, max_xy),
  variable = c("act", "act"), value = c(max_xy, min_xy)))

p <- ggplot(plot, aes(x = pred, y = value)) + geom_point(size = 2.5) + theme_bw()
p <- p + facet_wrap(~variable, scales = "free")

print(p)

포인트도 표시되는 것을 제외하고는 내가 원하는 것을 수행합니다.

여기에 이미지 설명 입력

이와 같은 작업에 대한 제안 사항이 있습니까?


이 아이디어 를 추가하는 것을 보았지만 비트 geom_blank()를 지정하고 aes()제대로 작동하는지 또는 geom_point()히스토그램 사용에 해당하는 것이 무엇인지 잘 모르겠습니다 aes(y = max(..count..)).


다음은 재생할 데이터입니다 (용해 전 실제 값, 예측 값 및 잔여 값).

> dput(results)
structure(list(act = c(52.81, 44.46, 54.5866666666667, 36.2333333333333, 
53.2266666666667, 41.7233333333333, 35.2966666666667, 30.6833333333333, 
39.25, 35.8866666666667, 25.1, 29.0466666666667, 23.2766666666667, 
56.3866666666667, 42.92, 41.57, 27.92, 23.16, 38.0166666666667, 
61.8966666666667, 37.41, 41.6333333333333, 35.9466666666667, 
48.9933333333333, 30.5666666666667, 32.08, 40.3633333333333, 
53.2266666666667, 64.6066666666667, 38.5366666666667, 41.7233333333333, 
25.78, 33.4066666666667, 27.8033333333333, 39.3266666666667, 
48.9933333333333, 25.2433333333333, 32.67, 55.17, 42.92, 54.5866666666667, 
23.16, 64.6066666666667, 40.7966666666667, 39.0166666666667, 
41.6333333333333, 35.8866666666667, 25.1, 23.2766666666667, 44.46, 
34.2166666666667, 40.8033333333333, 24.5766666666667, 35.73, 
61.8966666666667, 62.1833333333333, 74.6466666666667, 39.4366666666667, 
36.6, 27.1333333333333), pred = c(52.8675013282404, 42.7682474758679, 
49.0048248585123, 35.5238560262515, 48.7942868566949, 41.5750416040131, 
33.9548164913007, 29.9787449128663, 37.6443975781139, 36.7196211666685, 
27.6043278172077, 27.0615724310721, 31.2073056885252, 55.0886903524179, 
43.0895814712768, 43.0895814712768, 32.3549865881578, 26.2428426737583, 
36.6926037128343, 56.7987490221996, 45.0370788180147, 41.8231642271826, 
38.3297859332601, 49.5343916620086, 30.8535641206809, 29.0117492750411, 
36.9767968381391, 49.0826677983065, 54.4678549541069, 35.5059204731218, 
41.5333417555995, 27.6069075391361, 31.2404889715121, 27.8920960978598, 
37.8505531149324, 49.2616631533957, 30.366837650159, 31.1623492639066, 
55.0456078770405, 42.772538591063, 49.2419293590535, 26.1963523976241, 
54.4080781796616, 44.9796700541254, 34.6996927469131, 41.6227713664027, 
36.8449646519306, 27.5318686661673, 31.6641793552795, 42.8198894266632, 
40.5769177148146, 40.5769177148146, 29.3807781312816, 36.8579132935989, 
55.5617033901752, 55.8097119335638, 55.1041728261666, 43.6094641699075, 
37.0674887276681, 27.3876960746536), resid = c(-0.0575013282403773, 
1.69175252413213, 5.58184180815435, 0.709477307081826, 4.43237980997177, 
0.148291729320228, 1.34185017536599, 0.704588420467079, 1.60560242188613, 
-0.832954500001826, -2.50432781720766, 1.98509423559461, -7.93063902185855, 
1.29797631424874, -0.169581471276786, -1.51958147127679, -4.43498658815778, 
-3.08284267375831, 1.32406295383237, 5.09791764446704, -7.62707881801468, 
-0.189830893849219, -2.38311926659339, -0.541058328675241, -0.286897454014273, 
3.06825072495888, 3.38653649519422, 4.14399886836018, 10.1388117125598, 
3.03074619354486, 0.189991577733821, -1.82690753913609, 2.16617769515461, 
-0.088762764526507, 1.47611355173427, -0.268329820062384, -5.12350431682565, 
1.5076507360934, 0.124392122959534, 0.147461408936991, 5.34473730761318, 
-3.03635239762411, 10.1985884870051, -4.18300338745873, 4.31697391975358, 
0.0105619669306023, -0.958297985263961, -2.43186866616734, -8.38751268861282, 
1.64011057333683, -6.36025104814794, 0.226415618518729, -4.80411146461488, 
-1.1279132935989, 6.33496327649151, 6.37362139976954, 19.5424938405001, 
-4.17279750324084, -0.467488727668119, -0.254362741320246)), .Names = c("act", 
"pred", "resid"), row.names = c(2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 
10L, 11L, 12L, 13L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 
24L, 25L, 26L, 28L, 29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 
38L, 39L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 
52L, 54L, 55L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L
), class = "data.frame")

궁금한 점이 있습니다. 동일한 그래프에 실제와 잔차를 표시하지 않는 이유는 무엇입니까?
Ricardo Saporta 2013 년

2
플롯을 별도로 생성 한 다음 grid.arrange.
joran 2013-08-04

@RicardoSaporta 예를 들어 링크 할 수있는 Google 이미지가 있습니까? 용융 후 데이터를 사용하여 ggplot(plot, aes(x = pred, y = value)) + geom_point()패싯없이 수행 할 것을 제안하고 있습니까? 비 무작위성 / 왜곡을 감지하기 어렵게 만들기 위해 잔차의 척도를 실제로 축소하지 않을까요?
Hendy

1
내 다른 의견은 패싯 코드 적다 는 것입니다 ... 나는 녹인 다음에서 variable생성 한 값으로 플롯하고 패싯해야 했습니다 melt(). 그런 다음 다시 lapply다양한 조합을 구성하기 위해 만든 목록에 저장할 수 있다고 가정 합니다. 입력 해 주셔서 감사합니다. grid솔루션 을 만들고 싶다면 대답을 받아 들일 수 있지만 그것이 우리가 취하는 경로라면 다른 grid기반 솔루션 과 중복 될 수도 있습니다 .
Hendy

1
@joran과 나는 사람들에게 거의 항상 레이아웃을 엉망으로 만드는 사용 하지 말라고 조언 하고 있습니다 grid.arrange. gtable의 오랜 버그가 해결 되었으면합니다.
baptiste

답변:


115

다음은 더미 geom_blank레이어 가있는 코드입니다.

range_act <- range(range(results$act), range(results$pred))

d <- reshape2::melt(results, id.vars = "pred")

dummy <- data.frame(pred = range_act, value = range_act,
                    variable = "act", stringsAsFactors=FALSE)

ggplot(d, aes(x = pred, y = value)) +
  facet_wrap(~variable, scales = "free") +
  geom_point(size = 2.5) + 
  geom_blank(data=dummy) + 
  theme_bw()

여기에 이미지 설명 입력


11
이것에 대한 좋은 변형은 expand_limits(pred=range_act, value=range_act)사용 geom_blank하지만 사용하기 더 간단합니다.
eregon

6
이것은 제한을 확장 할뿐 (축소하지는 않음) 범위를 줄이는 방법이 있습니까? @baptiste
Indranil Gayen에게

32

나는 당신이 원하는 것을 이해하지 못하지만 내가 이해 한 것을 바탕으로

x 스케일이 동일 해 보이지만 동일하지 않은 것은 y 스케일입니다. 이것은 스케일을 = "free"로 지정했기 때문입니다.

scales = "free_x"를 지정하여 x 만 자유로울 수 있습니다 (이 경우 pred가 정의에 따라 동일한 범위를 갖는 것과 동일 함)

p <- ggplot(plot, aes(x = pred, y = value)) + geom_point(size = 2.5) + theme_bw()
p <- p + facet_wrap(~variable, scales = "free_x")

나를 위해 일했다, 그림을 참조하십시오

여기에 이미지 설명 입력

나는 당신이 그것을 너무 어렵게 만들고 있다고 생각합니다. 최소와 최대가있는 공식을 기반으로 한도를 정의 한 것을 기억하는 것 같습니다.면 처리되면 해당 값 만 사용했다고 생각하지만 코드를 찾을 수 없습니다


7

coord_cartesian 명령으로 범위를 지정하여 원하는 y 축 범위를 설정할 수도 있습니다. 이전 게시물에서와 같이 scales = free_x를 사용합니다.

p <- ggplot(plot, aes(x = pred, y = value)) +
     geom_point(size = 2.5) +
     theme_bw()+
     coord_cartesian(ylim = c(-20, 80))
p <- p + facet_wrap(~variable, scales = "free_x")
p

여기에 이미지 설명 입력

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.