非参数统计实验一
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实验报告
课程名称: ___________
实验项目名称: _______
姓名:
学号:
班级:
指导教师:
2016 年3 月10 日
新疆财经大学实验报告
实验要求与数据:
1、产生50个标准正态分布的随机数并画出它们的正态分布概率图形.
2、画出参数变化的概率密度分布图,正态分布的u, c变化,产生100个均值为5, 标准差
为1的正态分布的随机数,再产生100个均值为6,标准差为1的正态分布的随机数,画概率密度图
3、设随机变量X取区间[-5,5]上步长为0.1的各值,计算X的服从参数为5的t分布的概率,并画出概率密度函数图形,同时画出标准正态概率密度曲线,观察二者的区别.
4、根据已知数据,求秩和正态得分
2012
1. set.seed(1)
x<-seq(-5,5,le ngth.out=50)
yv-d norm(x,0,1)
plot(x,y,col="red",xlim=c(-5,5),ylim=c(0,1),type='l',xaxs="i", yaxs="i",ylab='de nsity',xlab二”,ma in二"The Normal Den sity Distributio n")
lin es(x,d no rm(x,0,0.5),col二"gree n")
lin es(x,d no rm(x,0,2),col二"blue")
lin es(x,d no rm(x,-2,1),col="ora nge")
lege nd("topright",lege nd=paste("m=",c(0,0,0,-2),"sd=",c(1,0.5,
2,1)),lwd=1,col=c("red","gree n","blue","ora nge"))
The Normal Density Distribution
2.画出参数变化的概率密度分布图,正态分布的u, 变化,产生100个
均值为5,标准差为1的正态分布的随机数,再产生100个均值为6,标准差为1的正态分布的随机数,画概率密度图
命令:
set.seed(1)
x <- seq(0,10,le ngth.out=100)
y <- dno rm(x,5,1)
plot(x,y,col二"red",xlim二c(0,10),ylim二c(0,1),type二T,
xaxs="i", yaxs="i",ylab='de nsity',xlab二”,
mai n="The Normal Den sity Distributi on")
lin es(x,d norm(x,6,1),col="gree n")
legend("topright",legend=paste("m=",c(5,6)," sd=", c(1,1)), lwd=1, col=c("red", "gree n"))
o
3. set.seed(1)
x<-seq(-5,5,le ngth.out=1000)
y<-dt(x,1,0)
plot(x,y,col="red",xlim=c(-5,5),ylim=c(0,0.5),type="l",xaxs="i ",yaxs="i",ylab='de nsity',xlab二”,ma in二"The T Den sity Distributio n")
lin es(x,dt(x,5,0),col 二"gree n") lin es(x,dt(x,15,0),col 二"blue") lin es(x,dt(x,100,0),col="ora nge")
legen d("topleft",lege nd=paste("df=",c(1,5,15,100)," ncp=",c(0,0 ,0,0)),lwd=1,col=c("red","gree n","ora nge") set.seed(1)
x<-seq(-5,5,len gth.out=1000) yv-duni f(x,0,1)
plot(x,y,col="red",xlim=c(0,10),ylim=c(0,2),type="l",xaxs="i", yaxs 二"i",ylab 二'den sity',xlab='',ma in 二"The Uniform Den sity Distributio n") lin es(x,d un if(x,0,0.5),col="gree n")
in
p d
4 -2 0 2 4
The T Density Distribution
寸
o
Eo CMC
lin es(x,d un if(x,0,2),col 二"blue")
lin es(x,d un if(x,1,6),col="ora nge") lin es(x,du ni f(x,2,4),col 二"purple")
lege nd("topright",lege nd=paste("m=",c(0,0,0,1,2),"sd=",c(1,0.5 ,2,6,4)),lwd=1,col=c("red","gree n","ora nge","purple")
The Uniform Density Distribution
4.根据已知数据用SPS 统计软件得到的结果为
畫 RATJOOI
P
q
R*畠
Pp Qq
ppp
m
20352 4QW W -L53 651bB 1 WO -153
657 » 1000
-1.53
:
1G1C5 <6DDD 13 -1 15 4573® 10BW 粘 417 »
sm
64 00 石:旬于弄
1E&L5
11 ODO 1-9 -■K
费兔00 3 000 19 砸M 3 000 19 •船 a 县
3QW
普
展啊
5WQ 再 碌W TWO 44 —16 和mt 昭阴 1 DM .31 -4-9 2 000 13 -1 15
5WW 4m 25 -.67 羅糕玮 U522 帖.DOO 背 ■32 4$4J8® TWO
44 4-12 « 4 600
務
24 乌■衣齐丽
I37BB 12 DM
44 -陀
6297 M
31
-49 毁CO GOOA
-J2
13*M
輿
MWW
4WD 2& ■6F 6&1 W 7WQ*
13
-115
12501 14 DDO 1C
牺押曲 em W 8 3T&W
um 11S
常巳;町县
12175 5QDD B3 32 3743 g 15000 S4 153 刚M 1S0® 的 153 11^43 15 tiM
的
羁
竝騁那
11
的
4$ 占吃M
&0C4 主1 -.4曲
1H 阪
12CVD 67 W 17 004 7$ 57 跆具 115TO ^DDO B1 曲 ?W8® 14BCK 1 16 412 » 5 5(n
&9
24 晋t 县
11505 2Q£H) 曲
1 15 轴&3備 13000 S1 S9 37S W
勢
105T9
TOM
1^3
4J&3 00
9 MO
*
1€ 地M
11 000
4$
in
10