抽样调查实验一

1实 验（实训）报 告项 目 名 称 R 语言初步 所属课程名称 抽样调查 项 目 类 型 综合性实验 实验(实训)日期 2013 年 4 月 1 日 班 级 学 号 姓 名 指导教师 浙江财经学院教务处制实验名称: R 语言初步实验目的: 熟悉 R 语言实验要求:会以下操作:21. 安装 R 并导入程序包2. 查询函数帮助3. 数据读取和 t 检验4. 方差分析和箱线图绘制5. 下标和条件筛选6. 了解工作路径7. R 脚本运行8. 绘图练习9. 编写函数作业: 画出三大抽样分布 分布、t 分布、F 分布的概率密度函数，2并比较参数不同时，概率密度函数的不同形态。作业说明：作业用 word 编辑，以唯一附件的方式上传。作业中要包含程序及图形。3实验过程:1 分布2概率密度函数(p.d.f.): )(xfn0,0,)2(121xen程序 x y par(lwd=2) plot(function(x) dchisq(x,1), xlim=x, ylim = y, col=hsv(0),xlab=“n“,ylab=“density of t(n)“,lty=2,main= “the Chisq Density Distribution“) plot(function(x) dchisq(x,1),xlim=x, ylim = y, col=”red”,add=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), lwd=3,+col=c(“red“,”green”,“blue“,”orange”),bty=“n“)p.d.f.图像42t 分布概率密度函数(p.d.f.): 21)(21)(ntntft程序 x y par(lwd=2) plot(function(x) dnorm(x,0,1), xlim=x, ylim = y, col=hsv(0),xlab=“n“,ylab=“density of t(n)“,lty=2,main=“Central t - Density“) plot(function(x) dt(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(“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.图像53F 分布概率密度函数(p.d.f.): 0,0, )1()()2()(22112, 121 xxnxnf nn程序 x y par(lwd=2) plot(function(x) df(x,1,1), xlim=x, ylim = y, col=hsv(0),xlab=“n“,ylab=“density of t(n)“,lty=2,main=“The F Density Distribution“) 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,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),“df