R之应用及相关资源.pptVIP

R之應用及相關資源 演講者:簡孝竑 2009/11/23 指令 功能 dname( ) 密度函數 pname( ) 累積機率 qname( ) 分位數 rname( ) 產生亂數 分配 指令 Beta beta Binomial binom Cauchy cauchy Chi-square chisq Exponential exp F f Gamma gamma Geometric geom Hypergeometric hyper Negative binomial nbinom Normal norm Poisson pois Student’s t t Uniform unif 操作分配函數的配合指定 常用機率分配的指令 機率分配函數 離散型機率分配(以BINOMIAL分配為例)-1 eb-rbinom(50,20,0.22) eb [1] 4 1 4 5 4 6 5 4 5 2 4 4 5 6 3 6 3 9 3 3 3 4 3 3 4 6 4 2 2 5 3 3 6 4 1 6 [37] 1 5 5 4 1 6 6 2 6 3 4 8 4 5 eb1-table(eb) eb1 eb 1 2 3 4 5 6 8 9 4 4 10 13 8 9 1 1 barplot(eb1) 離散型機率分配(以BINOMIAL分配為例)-2 plot(ecdf(eb),main=CDF of eb) 連續型機率分配(以GAMMA分配為例)-1 eg-rgamma(100,5,3) hist(eg,prob=T) lines(density(eg),col=2) 連續型機率分配(以GAMMA分配為例)-2 plot(ecdf(eg)) eg1sort(eg) egy-ppoints(100) plot(egy~eg1,type=l,col=2,main=CDF of eg) 連續型機率分配(以GAMMA分配為例)-3 plot(ecdf(eg)) lines(egy~eg1,type=l,col=2) abline(h=0.5,lty=3,col=gray50) 中央極限定理 藉由軟體來模擬中央極限定理。一母體平均值為 ，變異數為 ，則從該母體隨機抽取n個樣本，當n越大，不論何種分配， ,with large n。 par(mfcol=c(2,5)) x-matrix(,100,5) for(n in 1:5){ + for(i in 1:100){ + x[i,n]=((mean(rbinom(50*(2^n),1,0.6))-0.6)*sqrt(50*(2^n)))/sqrt(0.6*0.4)} + ;hist(x[,n],main=paste(n=,50*(2^n))) + ;qqnorm(x[,n],main=paste(qqplot with n=,50*(2^n))) + ;qqline(x[,n],main=paste(qqplot with n=,50*(2^n)),col=red)} 生成binomal然後利用迴圈(for())使得n越來越大 (n=100,400,800,1600) 觀察各自的直方圖以及qqplot 中央極限定理-動畫 網路資源-1 1. .tw/huangwj/book/R.pdf R?基本使用手冊- 許湘伶/國立中山大學應用數學研究所 博士生 黃文璋/國立高雄大學應用數學系 教授 羅夢娜/國立中山大學應用數學系 教授 編著九十三年八月 2 ..tw/chinese/epaper/00099005-1.htm 專題報導_R language簡介-范素瑋 / 臺灣大學生態學與演化生物學所博士候選人 3. http://bm2.genes.nig.ac.jp/RGM2/index.php?clear=all http://addictedtor.free.fr/graphiques/ R的圖庫 * *