1、實驗（實訓）報告項目名稱 R 語言初步所屬課程名稱抽樣調查項目類型綜合性實驗實驗（實訓）日期 2013年4月1日班 級學 號姓 名指導教師浙江財經學院教務處制實驗名稱:R語言初步實驗目的:熟悉R語言實驗要求:會以下操作:1. 安裝R并導入程序包2. 查詢函數幫助3. 數據讀取和t檢驗4. 方差分析和箱線圖繪制5. 下標和條件篩選6. 了解工作路徑7. R腳本運行8. 繪圖練習9. 編寫函數作業:畫出三大抽樣分布2分布、t分布、F分布的概率密度函數，并比較參數不同時，概率密度函數的不同形態。作業說明： 作業用 word 編輯，以唯一附件的方式上傳。作業中要包含程 序及圖形第 2 頁實驗過程:1.

2、 2分布概率密度函數(p.d.f.):1nfn(x)22 (0,xnx21 2x2 e 2,x 0 ;)0程序> x<-c(0,10)> yv-c(0,0.8)> par(lwd=2)> plot(fu ncti on(x)dchisq(x,1),xlim=x,ylim=: y,col=hsv(0),xlab="n ",ylab二"de nsityof t(n)",lty=2,main="theChisqDe nsity Distributio n")>plot(fu ncti on(x)dchis

3、q(x,1),xlim=x,ylim=y,col= ” red ”，add二TRUE)>plot(fu ncti on(x)dchisq(x,2),xlim=x,ylim=: y,II II1力力yaxs="i",col=gree n,add=TRUE)>plot(fu ncti on(x)col= ” blue ” ,add=TRUE)dchisq(x,3),xlim=x,ylim=: y,>plot(fu ncti on(x)dchisq(x,5),xlim=x,ylim=y,col= ” orange ” ,add=TRUE)lege nd(&quo

4、t;topright",lege nd=paste("df=",c(1,2,3,5), lwd=3,+col=c("red", ” green ”>,"blue", ” orange ” ),bty="n")p.d.f.圖像第7頁2. t分布概率密度函數（p.d.f.）:> x<-c(-4,4)> yv-c(0,0.4)> par(lwd=2)>plot(fu nctio n(x)dn orm(x,0,1), xlim=x, ylim = y,col二hsv(0),xl

5、ab二"n ",ylab二"de nsity of t(n )",lty=2,ma in二"Ce ntral t-Den sity")>plot(fu ncti on(x)dt(x,df =20),xlim=x,ylim=y,col=hsv(0.2),add=TRUE)>plot(fu ncti on(x)dt(x,df =10),xlim=x,ylim=y,yaxs="i",col=hsv(0.4),add=TRUE)>plot(fu ncti on(x)dt(x,df= 5),xlim=x,y

6、lim=y,col=hsv(0.6),add=TRUE)>plot(fu ncti on(x)dt(x,df= 2),xlim=x,ylim=y,col二hsv(0.8),add二TRUE)>lege nd("topright",c("N(0,1)","t(20)","t(10)","t(5)","t (2) "),col二hsv(seq(0,1,0.2),bty=" n",cex=1.2,pch=15)> abli ne(v=0)p.d

7、.f.圖像3. F分布概率密度函數(p.d.f.):nin22n_,n 1n2L 1n 晉f (x)()( x) (1x)x 0fmdx)(QL)(卑 n2 n2n2,x 00, x 0程序> xv-c(0,2)> yv-c(0,0.8)> par(lwd=2)> plot(fu ncti on(x)df(x,1,1),xlim=x, ylim = y,col二hsv(0),xlab二"n ”,ylab二"de nsityof t(n )",lty=2,ma in二"TheFDen sity Distributio n"

8、)>plot(fu ncti on(x)df(x,3,1),xlim=x,ylim= y,col= ”green ” ,add=TRUE)>plot(fu ncti on(x)df(x,6,1),xlim=x,ylim= y,yaxs="i",col= ” blue ” ,add=TRUE)>plot(fu ncti on(x)df(x,3,3),xlim=x,ylim= y,col= ”orange ” ,add=TRUE)>plot(fu ncti on(x)df(x,3,6),xlim=x,ylim= y,col 二” pink ” ,add二

9、TRUE)> lege nd("topright",c(1,1),(3,1),(6,1),(3,3),(3,6), lwd=3,lege nd("topright",lege nd=paste("df1=",c(1,1),"df2=",c(3,1),"df3=",c(6,1),"df4=",c(3,3),"df5=",c(3,6),),lwd=3,col=c("red", ” green ” ,"blue", ” orange ” , ” pink ” ),bty="n")p.d.f.圖像希望以上資料對你有所幫助，附勵志名言 3 條：1、常自