非参数统计实验一

.z.实验报告课程名称：实验工程名称：姓名：学号：班级：指导教师：2016年3月10日新疆财经大学实验报告课程名称非参数统计实验类型综合实验实验工程名称非参数统计根底实验时数2学生姓名专业统计专业实验地点统计实验室实验日期年月日主要仪器设备〔实验软件〕R软件，SPSS与SAS软件实验目的通过实验掌握非参数统计与参数统计的区别，学习求秩的各种方法和正态得分计算法，掌握R统计软件的应用实验内容原理1、实验用参数改变分布说明参数的意义。2、非参数统计主要用秩，学习用秩的方法和正态得分3、学习R语言，并学会用R绘出统计图讨论心得成绩评定根据实验情况和实验报告质量写出真实评价：综合评分指导教师签名：时间：年月日实验要求与数据：1、产生50个标准正态分布的随机数并画出它们的正态分布概率图形.2、画出参数变化的概率密度分布图，正态分布的u,σ变化，产生100个均值为5,标准差为1的正态分布的随机数,再产生100个均值为6,标准差为1的正态分布的随机数,画概率密度图3、设随机变量*取区间[-5,5]上步长为0.1的各值,计算*的服从参数为5的t分布的概率,并画出概率密度函数图形，同时画出标准正态概率密度曲线，观察二者的区别.4、根据数据，求秩和正态得分2012年新疆城镇居民人均消费支出〔元/年〕城市总消费支出食品肉类蔬菜类乌鲁木齐市137855297863500克拉玛依市2039265111075657木垒自治县161654273627417精河县115793998782412焉耆县115054063898378库车县1495650861122436喀什市1057943831253394和田市1464960671647596伊宁市1250146981021376奎屯市134805400844651塔城市145224948807412沙湾县119434236713512福海县118624182686391哈巴河县121753743438302石河子市1554559599186061.set.seed(1)*<-seq(-5,5,length.out=50)y<-dnorm(*,0,1)plot(*,y,col="red",*lim=c(-5,5),ylim=c(0,1),type='l',*a*s="i",ya*s="i",ylab='density',*lab='',main="TheNormalDensityDistribution")lines(*,dnorm(*,0,0.5),col="green")lines(*,dnorm(*,0,2),col="blue")lines(*,dnorm(*,-2,1),col="orange")legend("topright",legend=paste("m=",c(0,0,0,-2),"sd=",c(1,0.5,2,1)),lwd=1,col=c("red","green","blue","orange"))2.画出参数变化的概率密度分布图，正态分布的u,σ变化，产生100个均值为5,标准差为1的正态分布的随机数,再产生100个均值为6,标准差为1的正态分布的随机数,画概率密度图命令：set.seed(1)*<-seq(0,10,length.out=100)y<-dnorm(*,5,1)plot(*,y,col="red",*lim=c(0,10),ylim=c(0,1),type='l',*a*s="i",ya*s="i",ylab='density',*lab='',main="TheNormalDensityDistribution")lines(*,dnorm(*,6,1),col="green")legend("topright",legend=paste("m=",c(5,6),"sd=",c(1,1)),lwd=1,col=c("red","green"))3.set.seed(1)*<-seq(-5,5,length.out=1000)y<-dt(*,1,0)plot(*,y,col="red",*lim=c(-5,5),ylim=c(0,0.5),type="l",*a*s="i",ya*s="i",ylab='density',*lab='',main="TheTDensityDistribution")lines(*,dt(*,5,0),col="green")lines(*,dt(*,15,0),col="blue")lines(*,dt(*,100,0),col="orange")legend("topleft",legend=paste("df=",c(1,5,15,100),"ncp=",c(0,0,0,0)),lwd=1,col=c("red","green","orange")set.seed(1)*<-seq(-5,5,length.out=1000)y<-dunif(*,0,1)plot(*,y,col="red",*lim=c(0,10),ylim=c(0,2),type="l",*a*s="i",ya*s="i",ylab='density',*lab='',main="TheUniformDensityDistribution")lines(*,dunif(*,0,0.5),col="green")lines(*,dunif(*,0,2),col="blue")lines(*,dunif(*,1,6),col="orange")lines(*,dunif(*,2,4),col="purple")le