資源描述:
《《抽樣調(diào)查實(shí)驗(yàn)》word版》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。
1、實(shí)驗(yàn)(實(shí)訓(xùn))報(bào)告項(xiàng)目名稱(chēng)R語(yǔ)言初步所屬課程名稱(chēng)抽樣調(diào)查項(xiàng)目類(lèi)型綜合性實(shí)驗(yàn)實(shí)驗(yàn)(實(shí)訓(xùn))日期2013年4月1日班級(jí)學(xué)號(hào)姓名指導(dǎo)教師浙江財(cái)經(jīng)學(xué)院教務(wù)處制實(shí)驗(yàn)名稱(chēng):R語(yǔ)言初步實(shí)驗(yàn)?zāi)康?熟悉R語(yǔ)言實(shí)驗(yàn)要求:會(huì)以下操作:61.安裝R并導(dǎo)入程序包2.查詢函數(shù)幫助3.數(shù)據(jù)讀取和t檢驗(yàn)4.方差分析和箱線圖繪制5.下標(biāo)和條件篩選6.了解工作路徑7.R腳本運(yùn)行8.繪圖練習(xí)9.編寫(xiě)函數(shù)作業(yè):畫(huà)出三大抽樣分布分布、t分布、F分布的概率密度函數(shù),并比較參數(shù)不同時(shí),概率密度函數(shù)的不同形態(tài)。作業(yè)說(shuō)明:作業(yè)用word編輯,以唯一附件的方式
2、上傳。作業(yè)中要包含程序及圖形。6實(shí)驗(yàn)過(guò)程:1.分布概率密度函數(shù)(p.d.f.):程序>x<-c(0,10)>y<-c(0,0.8)>par(lwd=2)>plot(function(x)dchisq(x,1),xlim=x,ylim=y,col=hsv(0),xlab="n",ylab="densityoft(n)",lty=2,main="theChisqDensityDistribution")>plot(function(x)dchisq(x,1),xlim=x,ylim=y,col=”red”,a
3、dd=TRUE)>plot(function(x)dchisq(x,2),xlim=x,ylim=y,yaxs="i",col=”green”,add=TRUE)>plot(function(x)dchisq(x,3),xlim=x,ylim=y,col=”blue”,add=TRUE)>plot(function(x)dchisq(x,5),xlim=x,ylim=y,col=”orange”,add=TRUE)legend("topright",legend=paste("df=",c(1,2,3,5
4、)),lwd=3,+col=c("red",”green”,"blue",”orange”),bty="n")>p.d.f.圖像62.t分布概率密度函數(shù)(p.d.f.):程序>x<-c(-4,4)>y<-c(0,0.4)>par(lwd=2)>plot(function(x)dnorm(x,0,1),xlim=x,ylim=y,col=hsv(0),xlab="n",ylab="densityoft(n)",lty=2,main="Centralt-Density")>plot(function(x)dt
5、(x,df=20),xlim=x,ylim=y,col=hsv(0.2),add=TRUE)>plot(function(x)dt(x,df=10),xlim=x,ylim=y,yaxs="i",col=hsv(0.4),add=TRUE)>plot(function(x)dt(x,df=5),xlim=x,ylim=y,col=hsv(0.6),add=TRUE)>plot(function(x)dt(x,df=2),xlim=x,ylim=y,col=hsv(0.8),add=TRUE)>legend
6、("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)>abline(v=0)p.d.f.圖像63.F分布概率密度函數(shù)(p.d.f.):程序>x<-c(0,2)>y<-c(0,0.8)>par(lwd=2)>plot(function(x)df(x,1,1),xlim=x,ylim=y,col=hsv(0),xlab="n",ylab="densityoft(n)
7、",lty=2,main="TheFDensityDistribution")>plot(function(x)df(x,3,1),xlim=x,ylim=y,col=”green”,add=TRUE)>plot(function(x)df(x,6,1),xlim=x,ylim=y,yaxs="i",col=”blue”,add=TRUE)>plot(function(x)df(x,3,3),xlim=x,ylim=y,col=”orange”,add=TRUE)>plot(function(x)df(x
8、,3,6),xlim=x,ylim=y,col=”pink”,add=TRUE)>legend("topright",c((1,1),(3,1),(6,1),(3,3),(3,6)),lwd=3,legend("topright",legend=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("r