R之應用及相關(guān)資源

R之應用及相關(guān)資源演講者:簡孝竑2009/11/23指令功能dname()密度函數(shù)pname()累積機率qname()分位數(shù)rname()產(chǎn)生亂數(shù)分配指令BetabetaBinomialbinomCauchycauchyChi-squarechisqExponentialexpFfGammagammaGeometricgeomHypergeometrichyperNegativebinomialnbinomNormalnormPoissonpoisStudent’sttUniformunif操作分配函數(shù)的配合指定常用機率分配的指令機率分配函數(shù)離散型機率分配(以BINOMIAL分配為例)-1>eb<-rbinom(50,20,0.22)>eb[1]414546545244563639333433464225336416[37]15541662634845>eb1<-table(eb)>eb1eb123456894410138911>barplot(eb1)離散型機率分配(以BINOMIAL分配為例)-2>plot(ecdf(eb),main="CDFofeb")連續(xù)型機率分配(以GAMMA分配為例)-1>eg<-rgamma(100,5,3)>hist(eg,prob=T)>lines(density(eg),col=2)連續(xù)型機率分配(以GAMMA分配為例)-2>plot(ecdf(eg))>eg1<sort(eg)>egy<-ppoints(100)>plot(egy~eg1,type="l",col=2,main="CDFofeg")連續(xù)型機率分配(以GAMMA分配為例)-3>plot(ecdf(eg))>lines(egy~eg1,type="l",col=2)>abline(h=0.5,lty=3,col="gray50")中央極限定理藉由軟體來模擬中央極限定理。一母體平均值為

，變異數(shù)為，則從該母體隨機抽取n個樣本，當n越大，不論何種分配，,withlargen。>par(mfcol=c(2,5))>x<-matrix(,100,5)>for(nin1:5){+for(iin1: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("qqplotwithn=",50*(2^n)))+;qqline(x[,n],main=paste("qqplotwithn=",50*(2^n)),col="red")}生成binomal然後利用迴圈(for())使得n越來越大(n=100,400,800,1600)觀察各自的直方圖以及qqplot中央極限定理-動畫網(wǎng)路資源-11.

http://.tw/huangwj/book/R.pdf R

基本使用手冊- 許湘伶/國立中山大學應用數(shù)學研究所博士生

黃文璋/國立高雄大學應用數(shù)學系教授

羅夢娜/國立中山大學應用數(shù)學系教授

編著九十三年八月2

..tw/chinese/epaper/0009/20090620_9005-1.htm

專題報導_Rlanguage簡介-范素瑋

/臺灣大學生態(tài)學與演化生物學所博士候選人3.http:/