模拟100个 正态分布的随机数,计算其样本均值和样本方差,观察其余真实值的差异。当增大样本容量时,其差异有何规律。

用R软件随机模拟产生6组均值为5,方差为100的随机数,其样本个数分别设为100,1000,10000,100000,1000000,10000000,100000000。观察各组数据的样本方差和样本均值的差异。下面给出模拟程序:

n <- c(100,1000,10000,100000,1000000,10000000,100000000)
for(i in 1:length(n)){
  mean01[i] <- mean(rnorm(n[i], mean=5, sd=10))
  var01[i] <- var(rnorm(n[i], mean=5, sd=10))
  sd01[i] <- sd(rnorm(n[i], mean=5, sd=10))
}
result <- data.frame(Mean=mean01,Var=var01,Sd=sd01)
result

为对比不同的样本对随机数的均值和方差的影响,下面给出6组对照表,其结果见表3。从表3可明显看出,随着样本量的增大,其产生随机数的样本均值和方差越接近于设定的均值和方差,即样本均值稳定于(接近于)总体均值,样本方差稳定于(接近于)总体方差。

在这里插入图片描述

  • 1
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
首先,生成正态分布N(15,1)随机数样本,可以使用matlab自带的函数`normrnd`,代码如下: ``` rng('default') % 设置随机数种子 mu = 15; sigma = 1; N = 100; sample = normrnd(mu,sigma,N,1); %生成N(15,1)随机数样本 ``` 接下来,使用t检验来检验总体均值µ≥15.2,t检验的步骤如下: 1. 假设检验 $H_0: \mu = 15.2$,$H_1: \mu < 15.2$。 2. 计算样本均值 $\bar{x}$ 和样本标准差 $s$。 3. 计算$t$统计量 $t = \frac{\bar{x} - \mu}{s/\sqrt{n}}$。 4. 计算$t$分布在自由度 $n-1$ ,左侧临界$t_{\alpha,n-1}$。 5. 判断$t < t_{\alpha,n-1}$,若成立,则拒绝原假设,否则不拒绝原假设。 使用matlab进行计算,代码如下: ``` alpha1 = 0.05; alpha2 = 0.01; mu0 = 15.2; % 计算样本均值样本标准差 x_bar = mean(sample); s = std(sample); % 计算t统计量 t = (x_bar - mu0)/(s/sqrt(N)); % 计算t分布左侧临界 t_crit1 = tinv(alpha1,N-1); t_crit2 = tinv(alpha2,N-1); % 判断是否拒绝原假设 if t < t_crit1 fprintf('在%.2f的显著性水平下,拒绝原假设\n',alpha1); else fprintf('在%.2f的显著性水平下,不拒绝原假设\n',alpha1); end if t < t_crit2 fprintf('在%.2f的显著性水平下,拒绝原假设\n',alpha2); else fprintf('在%.2f的显著性水平下,不拒绝原假设\n',alpha2); end ``` 输出结果为: ``` 在0.05的显著性水平下,不拒绝原假设 在0.01的显著性水平下,不拒绝原假设 ``` 因此,在显著性水平0.05和0.01下,不能拒绝总体均值µ≥15.2的原假设。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

sta@ma@brain

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值