实验报告三答.docVIP

实验报告三

课程统计软件分析与应用

学生姓名

学号

学院数学与统计学院

专业统计学

指导教师

二Ｏ一四年四月三日

1.编写一个R程序（函数）。输入一个整数n，如果n=0，则中止运算，并输出一句话：“要求输入一个正整数”；否则，如果n是偶数，则将n除2，并赋给n；否则，将3n+1赋给n。不断循环，直到n=1,才停止计算，并输出一句话：“运算成功”。这个例子是为了检验数论中的一个简单的定理。

程序：

myfun-function(n){

if(n=0){print(要求输入一个正整数)}

else{

repeat{

if(n%%2==0)n=n/2

elseif(n%%2!=0){

if(n==1){

print(运算成功);break

}

elsen=3*n+1

}

}

}

}

myfun(2)

[1]运算成功

myfun(-2)

[1]要求输入一个正整数

2.随机变量X~B（10,0.4），求分布率和分布函数，求P（X=3）。

程序：

x-0:10

y-dbinom(x,10,0.4)

z-pbinom(x,10,0.4)

y

z

结果：

y

[1]0.00604661760.04031078400.12093235200.21499084800.2508226560

[6]0.20065812480.11147673600.04246732800.01061683200.0015728640

[11]0.0001048576

z

[1]0.0060466180.0463574020.1672897540.3822806020.6331032580.833761382

[7]0.9452381180.9877054460.9983222780.9998951421.000000000

程序：

dbinom(3,10,0.4)

结果：

[1]0.2149908

3.随机变量X服从1到5的均匀分布,求P（1X3）

程序：

result-punif(3,1,5)-punif(1,1,5);result

结果：

[1]0.5

3*.随机变量X服从1到5的离散均匀分布,求P（1X3）

4.随机变量X~U（2,4），画出密度函数和分布函数，求P（2.5X3.5）。

程序：

x-seq(0,8,0.01)

y-dunif(x,2,4)

z-punif(x,2,4)

plot(x,z,type=l,col=Blue)

lines(x,y,type=l,col=Red)

result-punif(3.5,2,4)-punif(2.5,2,4);result

结果：

[1]0.5

5.随机变量X~E（4），画出密度函数和分布函数，求P（1X5）。

程序：

x-seq(0,5,0.01)

y-dexp(x,4)

z-pexp(x,4)

plot(x,z,type=l,col=Blue)

lines(x,y,type=l,col=Red)

result-pexp(5,4)-pexp(1,4)

result

结果：

result

[1]06.随机变量X~N（1,4），画出密度函数和分布函数，求P（-2X3）。

程序：

x-seq(-10,10,0.01)

y-dnorm(x,1,2)

z-pnorm(x,1,2)

plot(x,z,type=l,col=Blue)

lines(x,y,type=l,col=Red)

图像：

result-pnorm(3,1,2)-pnorm(-2,1,2)

result

结果：

result

[1]0.7745375

7.总体服从标准正态分布，用蒙特卡罗的方法模拟样本均值和样本中位数的均方误差。

程序：

set.seed(100)

options()$digi

op1-options(digits=3)

n-50###样本量

m-1000####重复抽样本的次数

tmean_median-matrix(0,m,2)

for(iin1:m){

x-rnorm(n)####从标准正态分布总体中抽出样本容量为50的一组样本。

tmean_median[i,1]-mean(x