模拟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
    评论
好的,我可以回答你关于正态分布随机数的问题。 首先,使用MATLAB中的normrnd函数可以产生正态分布随机数。代码如下: x = normrnd(20, sqrt(41), [1, n]); 其中,20为期望,sqrt(41)为标准差,n为样本容量。这个代码会返回一个大小为n的随机数向量x。 接下来,我们可以用MATLAB中的函数进行各种统计量的计算,代码如下: % 计算均值、方差、标准差、变异系数、偏度、峰度 mean_val = mean(x); var_val = var(x); std_val = std(x); cv_val = std_val / mean_val; skew_val = skewness(x); kurt_val = kurtosis(x); % 计算上下四分位数、四分位极差 [q1, q3] = prctile(x, [25, 75]); iqr_val = q3 - q1; % 计算均值 m1 = mean_val; m2 = median(x); m3 = (q1 + 2 * m2 + q3) / 4; 这些统计量的含义如下: - 均值:表示随机数的平均值。 - 方差:表示随机数的离散程度。 - 标准差:表示随机数的离散程度,但受均值影响较小。 - 变异系数:表示标准差与均值之比,用于比较两个不同均值和标准差的随机数的离散程度。 - 偏度:表示随机数分布的偏斜程度,正表示右偏,负表示左偏。 - 峰度:表示随机数分布的峰度程度,正表示高峰,负表示低峰。 - 上下四分位数:表示随机数中位数左右各25%的数据。 - 四分位极差:表示上下四分位数之差,用于衡量数据的变异程度。 - 三均值:表示随机数分布的偏斜程度,分别以均值、中位数和上下四分位数表示。 以上统计量的计算方法均可以利用MATLAB中的函数进行计算

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

打赏作者

sta@ma@brain

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

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

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

打赏作者

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

抵扣说明:

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

余额充值