기본 콘텐츠로 건너뛰기

R에서 그래프 저장하기



 R에서 데이터 분석을 하기 위해선 여러 가지 그래프가 필수적입니다. 데이터 시각화는 데이터를 더 쉽게 이해할 수 있도록 도와줍니다. 더 나아가 선형 회귀 분석에서 모델 평가를 위해서는 반드시 잔차도를 비롯한 그래프를 비교 검토해야 합니다. 이 그래프를 저장하는 것은 R 스튜디오에서 기본으로 기능을 지원하기 때문에 R 스튜디오에서 시행해도 되지만, 당연히 R 콘솔에서 저장도 가능합니다. 오늘은 앞서 예제에서 나온 그래프를 저장하는 코드에 대해서 설명하겠습니다. 


 기본적으로 추천하는 그래픽 파일 포맷은 png나 pdf 입니다. jpeg도 물론 지원하지만, 사실 고해상도의 파일을 저장하는 용도로는 그렇게 추천하지 않으며 대개의 R 설명서에서 쓰지 말라고 이야기하고 있습니다. 그래도 R은 기본으로 아래 파일을 지원합니다. 명령어도 그냥 png나 jpeg를 사용할 수 있습니다. ?png 코드로 상세한 내용을 확인할 수 있습니다. 


bmp(filename = "Rplot%03d.bmp",
    width = 480, height = 480, units = "px", pointsize = 12,
    bg = "white", res = NA, family = "", restoreConsole = TRUE,
    type = c("windows", "cairo"), antialias)

jpeg(filename = "Rplot%03d.jpg",
     width = 480, height = 480, units = "px", pointsize = 12,
     quality = 75,
     bg = "white", res = NA, family = "", restoreConsole = TRUE,
     type = c("windows", "cairo"), antialias)

png(filename = "Rplot%03d.png",
    width = 480, height = 480, units = "px", pointsize = 12,
    bg = "white", res = NA, family = "", restoreConsole = TRUE,
    type = c("windows", "cairo", "cairo-png"), antialias)

tiff(filename = "Rplot%03d.tif",
     width = 480, height = 480, units = "px", pointsize = 12,
     compression = c("none", "rle", "lzw", "jpeg", "zip", "lzw+p", "zip+p"),
     bg = "white", res = NA, family = "", restoreConsole = TRUE,
     type = c("windows", "cairo"), antialias)



 여기서는 png 파일로 저장하는 것을 해보겠습니다. 예제는 잔차도입니다. 


library("ggplot2")

set.seed(3311)
diamonds1<-sample 2500="" diamonds="" nrow="" span="">
D1<-diamonds diamonds1="" span="">

model=lm(price~carat, data=D1)
summary(model)

par(mfrow=c(2,2))
plot(model)


 이 코드에서 나온 그래프를 어떻게 저장할 수 있을까요? H 디스크에 R post 폴더에  R graph 폴더를 넣고 residualPlot01이라는 이름의 파일을 만들어 보겠습니다. 


png("H:R post\\R graph\\residualPlot01.png",width = 2000, height = 1500)
par(mfrow=c(2,2))
plot(model)
dev.off()


 우선 png()을 열고 파일의 경로와 이름을 지정해 주어야 합니다. 이때 실제 콘솔에서는 \\ 대신 사선으로 나타난다는 점을 확인할 수 있습니다. 



 그 다음 몇 개의 그림을 하나에 담을 것인지를 결정해주어야 합니다. 여기서는 표준인 2x2 구성으로 했습니다. 다음 plot을 통해 그림을 그려주면 완성이지만, 여기서 디바이스를 닫아야 다음 그림에서 오류가 생기지 않습니다. dev.off()라고 해주면 됩니다. 사용자가 한 가지 더 결정해야 하는 중요한 변수는 해상도로 여기서는 2000x1500 픽셀로 했습니다. jpeg는 잘 쓰지는 않지만 같은 방법으로 저장할 수 있습니다. 




 또 다른 중요한 포맷은 pdf입니다. pdf의 장점은 여러 장의 사진을 하나의 파일에 고품질로 담을 수 있다는 점입니다. 비교를 위해서 두 개의 코드를 사용해 보겠습니다. 


png("H:\\R post\\R graph\\residualPlot02.png",width = 2000, height = 1500)
par(mfrow=c(2,2))
plot(model)
par(mfrow=c(1,1))
ols_cooksd_barplot(model)
ols_dfbetas_panel(model)
ols_dffits_plot(model)
dev.off()

pdf("H:\\R post\\R graph\\residualPlot02.pdf",width = 16, height = 16)
par(mfrow=c(2,2))
plot(model)
par(mfrow=c(1,1))
ols_cooksd_barplot(model)
ols_dfbetas_panel(model)
ols_dffits_plot(model)
dev.off()


 png는 크기에서 픽셀 기준이지만, pdf는 인치 기준이라는 점을 주의해야 합니다. 아무튼 png는 해보면 알겠지만, 마지막 파일만 남고 덮어쓰기가 이뤄집니다. 반면 pdf는 파일을 각 페이지에 저장합니다. 따라서 여러 장의 그래프를 하나의 파일에 저장하고 싶다면 pdf가 가장 유리한 방법입니다. 고품질 출력물의 기준이기도 하기 때문에 pdf는 여러 모로 유용한 포맷입니다. 








 R에서 다양한 pdf를 포함해서 다양한 문서를 만들 수 있는 방법에 대해서는 'R과 Knitr를 활용한 데이터 연동형 문서 만들기'같은 전문 서적을 추전합니다. 반대로 pdf나 문서 파일, 여러 이미지 파일에서 데이터를 추출해서 R에서 처리할 수 있는데, 나중에 다뤄보겠습니다. 

댓글

이 블로그의 인기 게시물

100 테슬라급 자기장 도달

 미국의 로스 알라모스 국립 연구소 (Los Alamos National Laboratory) 에서 과학자들이 지금까지 인간이 개발한 가장 강력한 자기장을 발생시키는 장치 개발에 도전하고 있습니다. 자기장의 세기를 나타내는 방법으로 자기력선의 밀도를 나타내기 위해 단위 면적당 자기력선의 수를 표시하는 단위인 테슬라 (T) 가 있습니다. (1T = 1Wb/㎡  웨버 (Wb) 는 자속의 단위)   의료용으로 사용되는 초전도체를 이용한 MRI 의 경우 1.5 - 3 테슬라급의 강력한 자기장으로 인체 내부를 볼 수 있게 만들지만 과학 연구용으로 이보다 더 강력한 자기장이 필요할 수 있습니다. 최근에 등장한 90 테슬라급 자기장에 이어 이번에 로스 알라모스 국립 연구소에서는 100 테슬라급인 100.75 T 를 실현 했다고 합니다.   이를 구현한 것은 18000 파운드 (8.16 톤 정도) 의 코일과 여기에 에너지를 공급할 1200 메가줄 (Megajoule) 급 모터 제네레이터등의 설비입니다.  (   The 1,200-megajoule motor generator that powers the magnetic pulse.  )  이와 같은 연구를 통해 알아내고자 하는 것은  Quantum Phase transitions and new ultra high field magnetic states Electronic Structure determination Topologically protected states of matter  로 요약할 수 있다고 합니다.   아무튼 수 T 급 MRI 만 해도 자기장의 힘이 엄청난데 100 T 라니 엄청난 자기장이네요. 이는 지구 자기장 세기보다 200만배 강력한 (물론 좁은 범위에서 작용하는 자기장이라 지구 전체...

고대 양서류 이야기 (2) - 악어를 닮은 거대 양서류들

  페름기에는 다양한 양막류가 진화해서 앞서 소개한 육상형 템노스폰딜리는 점차 설 자리를 잃게 됩니다. 하지만 양서류는 본래 자신의 서식지인 물과 습지에서 여전히 번성을 누렸습니다. 당시에는 악어류 같은 대형 양서형 파충류도 없던 시절이었기 때문에 이와 비슷한 생태학적 지위는 여전히 양서류의 몫이었습니다. 이들에 대한 이야기는 제 책인 포식자에서 비교적 간단히 다뤘는데, 오늘은 여기에 대한 보충 설명입니다.  책 정보:  http://book.naver.com/bookdb/book_detail.nhn?bid=13347200 Yes 24:  http://www.yes24.com/24/goods/58772859 11번가:  http://books.11st.co.kr/product/SellerProductDetail.tmall?method=getSellerProductDetail&prdNo=1977867160 알라딘:  http://www.aladin.co.kr/shop/wproduct.aspx?ItemId=134877825 교보문고:  http://www.kyobobook.co.kr/product/detailViewKor.laf?ejkGb=KOR&mallGb=KOR&barcode=9788970447988&orderClick=LAG&Kc= 인터파크 :  http://book.interpark.com/product/BookDisplay.do?_method=detail&sc.prdNo=279593764&sc.saNo=003002003&bid1=search_auto&bid2=detail&bid3=prd_img&bid4=001 영풍문고:  http://www.ypbooks.co.kr/book.yp?bookcd=100843205&gubun=NV ...

이빨이 다시 진화한 개구리

  ( CT scans of Gastrotheca guentheri skulls revealed what appeared to be identical rows of teeth on both the upper and lower jaws, which researchers later confirmed through dissection. Credit: Florida Museum/Daniel Paluh )  개구리는 2억년 전 진화 과정에서 이빨을 잃어버리고 큰 턱과 혀를 이용해 곤충 같은 작은 먹이를 잡아 먹는 방향으로 진화했습니다. 파충류나 포유류 같은 다른 사지류와의 경쟁에서 밀려 양서류가 쇠퇴하고 멸종하던 시기에도 개구리는 여전히 생존할 수 있었던 비결입니다.   이빨이 없는 덕분에 큰 혀를 발사하기 편해졌을지는 모르지만, 이빨이 없으면 종종 불편할 때가 있습니다. 씹는 대신 삼키기 때문에 음식을 씹을 수 없다는 점은 문제되지 않지만, 필사적으로 달아나려는 먹이를 잡기 힘들기 때문입니다. 이런 이유 때문에 일부 개구리는 이빨 같이 보이는 엄니 (fang)을 지니고 있으나 이는 사라진 이빨이 다시 난 것이 아니라 다른 부분이 진화한 것입니다.   이는 돌로의 법칙 ( Dollo's Law )으로 알려져 있습니다. 진화 과정에서 퇴화한 부분이 다시 생겨나지 않으며 대신 필요하면 다른 부분이 진화해 그 역할을 대신한다는 것입니다. 예를 들어 아가미가 사라진 사지 동물은 다시 물에 들어온다고 해도 아가미가 다시 생기진 않습니다. 대신 고래처럼 폐가 커져서 그 기능을 대신하게 됩니다.   하지만 모든 법칙엔 예외가 있기 마련입니다. 남미에서 발견된 멸종 위기 개구리 중 하나인 구엔터 유대류 개구리 ( Gastrotheca guentheri, Guenther's marsupial frog, dentate marsupial frog)는 완전한 형태의 이빨을 지니고 있습니다. 참고로 유대류 개구리라는 명칭은...