R語言上機(jī)實(shí)驗(yàn)三

理學(xué)院實(shí)驗(yàn)報(bào)告班級(jí)：學(xué)號(hào)班級(jí)：學(xué)號(hào):姓名:實(shí)驗(yàn)編號(hào):實(shí)驗(yàn)三：概率和分布的R實(shí)現(xiàn)一、 實(shí)驗(yàn)?zāi)康呐c要求：1、 會(huì)用R給出常見分布的概率密度、概率、分位數(shù)和隨機(jī)數(shù)。2、 會(huì)利用sample命令進(jìn)行隨機(jī)抽樣，prod,choose命令計(jì)算概率。3、 會(huì)利用R繪制各類分布的圖形。4、 會(huì)利用choose,prod命令計(jì)算古典概率。二、 實(shí)驗(yàn)內(nèi)容：1.從一副撲克牌(52張)中隨機(jī)抽5張，求下列概率(1)抽到的是10,J,Q,K,A；4/choose(52,5)1.539077e-06抽到的是同花順。9*choose(4,1)/choose(52,5)[1]1.385169e-05注：同花順是指5張同一色牌能按從小到大連續(xù)排序，如2<3<4<5<6,3<4<5<6<7,…，10<J<Q<K<A都為同花順。2.模擬隨機(jī)游動(dòng)：(1)從兩點(diǎn)分布中產(chǎn)生1000個(gè)隨機(jī)數(shù)；x<-rbinom(1000,1,0.5)x—l-l.一/ 5 — —fl]00100000111000□0001111101110111110101000012[50]llllOllllClDlUOlOCilllCllOUllOOOJOlllllliaDlll[93]<]110110G000001D01010Cllll0(j010<]ia0001Ollia0J[14S]lOlllOiOlOOlOllllllLlOiLllliJlOOlOlLllOOillOLlOJ[1S7]0011111G1101001001001101001100OIOJOOOIGOlie[Z-56]OOOlOllOiaOOOlOOliaOCIllDCiOOOLOLlDllOllOLlOC[2SS]lllOlllGlOlllOlOlOOlOlOllCjlllOllOlOOOOGOllC[3-5^]llQDQOiCLIOlDOlOOOOllUlOllCiilOCILIOlllQlQlOlOJ[3S3]looiiiciGiiioiiaiooaoooiiocjiiiiiaioociaciGDiiJ[4^2]1011010iiioiooaoiioooiooiooiiiiiiioooioociijrdeiininnnnnriiiiiinniiinnniiniiinnniinnnnniiniir⑵用函數(shù)ifelse()將上面隨機(jī)數(shù)中的0替換成-1;>ifelse(x==0廣1,1)

>ifelse(x==Dr-_f_11-1-111-1-1-1-1111-1-1111[1]1-11-1[24]11-11111-1-1-1-11-11111-11-1-1-11[47]-1-11-1-1-11-1-11111-1111-1-1-11-11【70]-1-111-11-1111-11-1-1-1-1-11-1-1-1-1-1[93]1-11-1-111-1-1-11-11-1-11-1-1111-1-1[116]111-111-111-1-1-1111-1-1-1-11-1-11[139]-1-1-111-1-1111-1-111111-1-11-1-1-1[162]-11-1-1-1-1-111-111-1-11-11-111111[135]-111-1-11-1-11-11-11-1-1-11-1-1-1-1-11[208]-11-1111-1-1111-1-1-1-1-111-1-11-11[231]-11-11-111-111-1-1-1-1-1-1-1-1-11-1-1-1[254]1111-111-11-1111-1-11-1111111[277]1-1-1-1-1-11-1-1-11-1-1-111-1-11-1-1-1-1[300]111111111-1-1111-1111-1-1111[323]11-111-111-11-11-11-11-11-1-1-111[346]1-11111-111-1-1111-1-1-111-1-111[369]1-111-1-1-111-11-11-1-11-1-1-11-1-11[392]1-1-111-1-1-111-1-111-11-11-111-1-1[415]11-11-1-1-11-1-1-11-11-1-11-1-111-1-1[438]1-1-111-11-111-1-1-11-1-1-1-1-111-1-1[461]1-11-11-11-111-1-1111-1-1-1111-11[484]-1-1-11-1-11-1-111-1-111-1-111111-1[507]1-11-11-1-1-1-1-1-1-1-111-11-111-1-11[530]-1-1-1-1-1-11111-111-111-1-111111[553]11-1-1-1-1-1-1-11111-11-1-1-1-11-1-11[576]-11111-1-1-1-11-111-1-11-11-111-1-1[599]-1-111-11-1-1-111-111-11111-1-11-1【m]1111-1-11-1111-1-11-111111111[645]1-111111111-1-1-1-11-11-1-111-11-1111-11-1-1-11-1-1-1-1-1-111-1-111-1【頃］1-1-11-111-1111-1-11111111-111[5]-1-11-1-1-1-1-1-1-111-1111-1-111-1-11[737]1-111-1-1-111-1-1-1-1-11-1-111-1-1-1-1E。]111111111111-1-1-1-11-1-11-1-1-1[783]-1-1-111-11-1111-1-11-1-11-1-11-1-11[8Q6]-1111-1-11-1111-1-11-1-1-11-1-1-1-1-1[829]-1-1-1-1-11-1-1-1-1-1111-1-1111-111-1[852]-1-1-1-1-1-1-11-111-111-11-11-1-1-1-1-1[375]-1111111-11-1-1-1-1-11-1-1-11-11-11[398]-1-1-1-1-1-11-1-111-1-11-1-111111-1-1[921]1111111-111-1111-1111111-1-1[944]1-1111-1-111-111-11-1-1-1-11-1-1-1-1[967]-1-1-11-1-1-1111-11-1-1-1-11-1-1-1-1-11[990]11-1-1-111-11-1-1(3)用函數(shù)cumsum()作出累積和;y<-ifelse(x==0,-1,1)cumsum(y)>cumsun(y)1.1]二二E5]152]Lg[8?]U03]U20][137][1S4][Ml]USS][205][222][心]1256][273][290][307]1324]2釘][358][375][M2]L"]〔霏們[E]0-11-23□1□5q70-1-11 -2 -33 2 3 2 3 ^1 -2 -3 -4 -3 -4 -3S-7 -7 -8 -7-2 -3 T -5 -4 -3S -9 -s -7 -e -3 -s冬54&2 3i 3 25 司 3-3 -47 -■€ -5 T3 T -5 -5 -5￡ -7 -S -7 -€一身-34-5-4-5-6-7-9邛-9-10-12-10-11-10一。-10-11-10-11-10-11-10-。-S-7-￡-7一6-5-€-7-6-7-B-7-B-7-0-7-S-S-10-5-ID-11-12-13-14-13-14-13-14-13-12-11-12-13-12-11-10-11-12-13-1E-14-13-14-15-14-1E-14-1E-14-IE-14~1S-14-13-14-13-12-13-14-1B-16-17-1S-1&-20-21-20-21-23-22-21-20-19-2D-19-1B-1S-1B-28-13-17-16-17-13-17-13-17-1石-15-14-13-12-11-12-13-14-15-1€-15-1€-17-IS-17-IS-15-2Q-1?-18-13-20一19-20-21-2Z-25-22-21-3D-IS-13-17-IE-1E-14-15-!￡-15-14-13-14-13-12-11-12-13-12-11-10-97777?—-1-4-5-ID67217■：-11-1F-10-10-11-12—D-7—D-7一二-3-2-3T-1-2-10-1-5一4-5-fcL—■—D一S一E-e-7-￡-7-8-u-1。-11-12-11-15-1€-15-IT-1E-1R一】匚-14一1n-10-15-10-13(4)使用命令plot(作出隨機(jī)游動(dòng)的示意圖.>plot(cumsum(y))漉RGraphics:Dnnce2(ACTIVE]D1H?1(WWE—。UXMM-UUUW-理-w2€01000008001000Index-10 12 3 4 5-10 12 3 4 53.在同一個(gè)圖形中畫出統(tǒng)計(jì)的四大分布密度曲線(dnorm,dchisq,dt,df),注意不同分布有不同的線型、顏色和寬度，還有圖形都要在同一方框中，最后用圖例說明(legend)。curve(dnorm(x,0,1),xlim=c(-1,5),ylim=c(0,0.5),col=1,lwd=1,lty=1)curve(dchisq(x,1),xlim=c(-1,5),ylim=c(0,0.5),lwd=2,lty=2,col=2,add=T)curve(dt(x,1),xlim=c(0,8),ylim=c(0,0.5),lwd=3,lty=3,col=3,add=T)curve(dt(x,1,1),xlim=c(0,8),ylim=c(0,0.5),lwd=4,lty=4,col=4,add=T)legend('topright',c("dnorm","dchisp","dt","df"),lty=c(1,2,3,4),col=c(1,2,3,4),lwd=c(1,2,3,4))curve(dnorm(x,0,1),xlim=c(-1,5),ylim=c(0,0.5),col=1,lwd=1,lty=1)curve(dchisq(x,1),xlim=c(-1,5),ylim=c(0,0.5),lwd=2,lty=8,col=2,add=T)curve(dt(x,1),xlim=c(0,8),ylim=c(0,0.5),lwd=5,lty=3,col=7,add=T)curve(dt(x,1,1),xlim=c(0,8),ylim=c(0,0.5),lwd=4,lty=4,col=4,add=T)legend('topright',c("dnorm","dchisp","dt","df"),lty=c(1,8,3,4),col=c(1,2,7,4),lwd=c(1,2,5,4))MlRGraphirE.：Device2(ACTI'-E] dnorm—dcmspdt 4.除本章給出的標(biāo)準(zhǔn)分布外，非標(biāo)準(zhǔn)的隨機(jī)變量X的抽樣可通過格式點(diǎn)離散化方法實(shí)現(xiàn).^p(x)^X的密度函數(shù)，其抽樣步驟如下在乂的取值范圍內(nèi)等間隔地選取^個(gè)點(diǎn)x,x,...,x,例如取N=1000;1 2 N計(jì)算’(了)；i=1,2,...,N;i正則化p(x);i=1,2,...,N,使其成為離散的分布律，即每一項(xiàng)除以尹p(xi);i=1按離散分布抽樣方法使用命令sample()從七i=1,2,...,N有放回地抽取〃個(gè)數(shù)，例如n=1000.注：前面4小步是用來編一個(gè)函數(shù)，功能是對(duì)給定的概率密度產(chǎn)生隨機(jī)數(shù)，形式應(yīng)與rnorm差不多。試以標(biāo)準(zhǔn)正態(tài)分布為例來說明.為與R中的正態(tài)抽樣函數(shù)rnorm()進(jìn)行比較，將作圖區(qū)域分為左右兩部分，使用rnorm()抽取n=1000個(gè)標(biāo)準(zhǔn)正態(tài)隨機(jī)數(shù)，并在左側(cè)區(qū)域畫出相應(yīng)的直方圖和核密度估計(jì)曲線；用格子點(diǎn)離散化抽樣方法完成抽樣，并在右側(cè)區(qū)域畫出相應(yīng)的直方圖和核密度估計(jì)曲線，離散化所用的N=1000，n=1000，取點(diǎn)范圍為[-4,4].N=seq(-4,4,length=1000)f<-function(x)+dnorm(x)/