一、R預(yù)設(shè)調(diào)色板
這一系列函數(shù)有5個，即：
rainbow
heat.colors
terrain.colors
topo.colors
cm.colors
在R環(huán)境里面輸入問號（?）和上面任一函數(shù)名就可以獲得這5個函數(shù)的用法說明。這些函數(shù)最少需要一個參數(shù)，n，表示要得到顏色的數(shù)量。n在系統(tǒng)允許范圍內(nèi)沒有限制。下面用彩虹色調(diào)色板函數(shù)rainbow產(chǎn)生的顏色繪一個色盤：

> setwd("D:/"); n=1000
> png("rainbow.disc.png", bg = "transparent")
> par(mar = c(0,0,0,0))
> pie(rep(1,times=n),labels="",col=rainbow(n),border=rainbow(n))
> dev.off()運行后回在D盤根目錄下得到一個rainbow.disc.png文件，圖形如下：

五個調(diào)色板產(chǎn)生的顏色對比如下：

> par(mfrow=c(5,1)); par(mar=c(0.1,0.1,2,0.1));par(xaxs="i", yaxs="i")
> n <- 10000
> barplot(rep(1,times=n),col=rainbow(n),border=rainbow(n),axes=FALSE, main="Rainbow colors");box()
> barplot(rep(1,times=n),col=heat.colors(n),border=heat.colors(n),axes=FALSE, main="heat.colors");box()
> barplot(rep(1,times=n),col=terrain.colors(n),border=terrain.colors(n),axes=FALSE,main="terrain.colors"); box()
> barplot(rep(1,times=n),col=topo.colors(n),border=topo.colors(n),axes=FALSE, main="topo.colors");box()
> barplot(rep(1,times=n),col=cm.colors(n),border=cm.colors(n),axes=FALSE, main="cm.colors"); box()

當(dāng)然你可能不需要這么多顏色，改變n值就可以了：

> n <- 10
> barplot(rep(1,times=n),col=rainbow(n),border=rainbow(n),axes=FALSE, main="Rainbow colors");box()
> barplot(rep(1,times=n),col=heat.colors(n),border=heat.colors(n),axes=FALSE, main="heat.colors");box()
> barplot(rep(1,times=n),col=terrain.colors(n),border=terrain.colors(n),axes=FALSE,main="terrain.c

x <- rnorm(10000,0,1)
hist(x, col="grey", prob=TRUE)
lines(density(x), col="red")

首先可以利用fitdistr函數(shù)求得weibull分布的形狀參數(shù)和尺寸參數(shù)，假設(shè)數(shù)據(jù)為x:

得到形狀參數(shù)shape與尺度參數(shù)scale

然后利用ks.test進(jìn)行檢驗：

上邊的jitter用來做小擾動，因為如果x中有重復(fù)數(shù)據(jù)的話ks.test會報錯，如果x中沒有重復(fù)數(shù)據(jù)則不需要jitter。shape是得到的形狀參數(shù)，scale是得到的尺度參數(shù)。

ks.test得到兩個結(jié)果，一個是D,越小越好，一個是p-value，這個值要大于0.05