常用的R绘图程式6RGraphics.PDF

第 6 章: 常用的 R 繪圖程式 6: R Graphics R 內建許多圖形工具函式, 這些圖形工具可以顯示各種統計繪圖並且自建一些全新的圖. 以下指令, 可以先看看 R 的圖形示範. demo(graphics) demo(image) demo(lattice) R 圖形工具函式之特色是用同一個繪圖函式, 對不同的類型物件, 可以作出不同的圖形. 圖形工具函式 既可互動式使用, 也可以批次處理使用. 在許多情況下, 互動式使用是最有效的. 圖形工具函式所產生之 圖形結果, 無法指定成物件, 必須送到 圖形裝置 (graphic device), 圖形裝置可以是一個視窗或某一特 定格式之圖形檔案. R 有一系列 圖形引數 (graphical parameters), 這些圖形引數可以修改或制定所 需的圖形環境. 繪圖指令可以分成了三個基本的類型: 1. 高階繪圖 (high-levlel plotting functions): 指令在圖形裝置上產生一個新的圖區, 它可以包 括坐標軸, 標籤, 標題等等. 2. 低階繪圖 (low-level plotting functions): 指令會在一個已經存在的繪圖上, 加上其它的圖 形元素, 如額外的點, 線與標籤等等. 3. 互動式繪圖 (interactive graphics functions): 指令允許互動式地用其他設備 (如滑鼠) 在一 個已經存在的繪圖上, 加上圖形資訊或者萃取圖形中的資訊. 在 R 中, grid 套件有一個獨立的強大之圖形系統, 根據 grid 套件, 所建構的另一個較容易使用之 lat- tice 套件, 可以產生多重漂亮起專業的統計繪圖, 類似 S-PLUS 裏的 Trellis 系統. 6.1 繪圖簡介 在此先示範一些常見的 R 繪圖函式與使用方法, 如圖 6.1. 許多時候, 開始用 R 繪圖函式會有一些 困惑, 但這些困惑通常會很快消失. 注意: 統計繪圖是一專業學科, 好的統計繪圖是要多次嘗試才能得到 較好得結果, 可參考些基本工具書, 如 Cleveland (1993, 1994), Tufte (1990, 1997), Wainer (1984, 1994, 1996) c 林建甫 (October 27, 2006) 2 6.1 繪圖簡介 x.norm-rnorm(100,mean=0,sd=1) x.group-sample(c(1:4), 100, replace=TRUE) x.cat-sample(c(1:8), 100, replace=TRUE, prob=c(8:1)) x.cat.tab-table(x.cat) barplot(x.cat.tab) pie(x.cat.tab) hist(x.norm) boxplot(x.norm) qqnorm(x.norm) plot(density(x.norm)) 以上圖形指令, 執行結果如圖 6.1. 例 1: NIDDM 臨床試驗比較合成胰島素與一般胰島素對降低血糖和總膽固醇的效應 一位研究者進行一個臨床試驗, 針對成人第二型糖尿病 (NIDDM), 在同時控制其他的因素, 包括 性別和發生糖尿病之年齡, 之後, 比較 2 種不同的治療: 新的治療藥物為 “合成胰島素 (Syn. Insulin)” 與 “一般胰島素 (Reg. Insulin)”, 對降低血糖和總膽固醇的效應. 總共 32 位病患, 以性別和發生糖尿 病之年齡分成 4 個階層, 隨機地分配到合成胰島素治療組與一般胰島素控制 (Reg. Insulin), 性別分 層為男與女, 發生糖尿病之年齡分層為 20-39 歲 和 40-59 歲. 測量基線 (baseline) 之空腹血糖和總膽 固醇 (FBS1, TC1), 與測量治療 12 個月之後的空腹血糖和總膽固醇 (FBS2, TC2). 研究的第 1 個 目地