课件：利用多发病例家族之单倍体资料进行遗传相关性研究-.ppt

另一種繪製機率密度函數圖的方法 x=c(-250:250)/50 matplot(x, outer(x, c(0.5,1,2),function(x,sd)dnorm(x,0,sd)), lty=1:3, type=c(l,l,l), ylab=, col=c(black,red,blue)) matlines(x, outer(x, c(-1,2),function(x,mu)dnorm(x,mu,1)), lty=4:5, type=c(l,l,l), ylab=, col=c(green,purple)) legend(1.5,0.8,c(N(0,0.5),N(0,1),N(0,2),N(-1,1),N(-2,1)), lty=1:5,col=c(black,red,blue,green,purple) ) 另一種繪製機率密度函數圖的方法 matplot: 將矩陣(x, y)依據行的資料，畫成圖形。 指令：matplot(x, y, 其餘圖形的設定指令)；在此y是利用outer指令產生。 outer: 利用函數function產生一個矩陣，亦即生成(x, y) 指令：outer(x, y, function) function: 自訂函數。 指令：function(自訂變數名稱) 表示式或公式。 matlines: 增加不同圖形於matplot上面。 matpoints: 同上。 Example of matplot table(iris$Species) # is data.frame with Species factor iS - iris$Species == “setosa” #邏輯變數 iV - iris$Species == versicolor #邏輯變數 op - par(bg = “bisque”) #背景是燒陶的顏色 matplot(c(1, 8), c(0, 4.5), type= n, xlab = Length, ylab = Width, main = “Petal and Sepal Dimensions in Iris Blossoms”) #type=“n”表示座標軸 matpoints(iris[iS,c(1,3)], iris[iS,c(2,4)], pch = c(s, S), col = c(2,4)) matpoints(iris[iV,c(1,3)], iris[iV,c(2,4)], pch = (v,V), col = c(2,4)) legend(1, 4, c(Setosa Petals, Setosa Sepals, Versicolor Petals, Versicolor Sepals), pch = sSvV, col = rep(c(2,4), 2)) 繪製機率密度函數圖 繪製Beta分佈機率密度函數圖形 將Beta(1,1), Beta(1,2), Beta(1, 0.5), Beta(3,2), Beta(4,4), Beta(0.5,0.25) 繪製在同一張圖上 統計圖形 圖形控制 說明 例子 高階繪圖函數 自動生成新的座標軸 plot() hist() 低階繪圖函數 增加點或線至原有的圖上 points() abline() 繪圖參數 控制一頁圖的所有細節 par(mfrow=c(2,1)) par(cex=0.5) 參數類型 說明 例子 高階參數 必須在高階繪圖函數下設定 xlab, type 版面編排參數 可以在高階繪圖函數或是par()函數下設定 cex, col 一般參數 必須透過par()函數才能設定 fig, mfrow 低階繪圖函數 指令 說明 例子 text 可增加文字於設定位置 plot(1:10,type=n) text(2,2,Hi,CWB!) text(locator(1),Hi) box 產生座標軸框框 box() rect 產生新的矩形 rect(-3,-3,3,3) new 增加高階繪圖函數的圖到現有的圖上 plot(1:10,ylab=) par(new=T) plot(10:1,ylab=) axis 可增加座標軸及修改刻度表示顏色等等細節 par(mfrow=c(2,2)) plot(1:10) plot(1:10,axes=F) axis(1) axis(2,lty=2,col=red) axis(3,at=c(5,9),labels=c(A,B)) box() 低階繪圖函數 指令 說明 例子 srt,crt 設定文字或符號旋轉的角度 plot(1:10,type=n) text(2,2,srt=90,srt=90) text(4,4,srt=0,srt=0) cex 設