我的问题:
如何在Javascript中创建具有给定均值和标准差(sd)的随机数列表?
示例:
我想创建一个包含5个随机数的列表,范围在1到10之间。结果均值应为5,标准差应为2.
到目前为止我做了什么:
我的想法是(http://jsfiddle.net/QASDG/3/):
注意:我还没有添加SD,因为代码并不优雅。
我发现该站点的另一个脚本声明:"此代码的目标是生成一个集合,其中包含5个正态分布的随机数,平均值为1.0,标准差为0.5。 " http://rosettacode.org/wiki/Random_numbers#JavaScript。我已经尝试通过添加" Math.floor"来根据我的需要调整它。并且改变了" for" -loop(i< 5)中的条件。
function randomNormal() {
return Math.floor(Math.cos(2 * Math.PI * Math.random()) * Math.sqrt(-2 * Math.log(Math.random())))
}
var a = []
for (var i=0; i < 5; i++){
a[i] = randomNormal() / 2 + 1
}
我想如何修改代码:
我不太确定我是否理解了此代码的数学部分。根据网页,它将创建一个正常的分布(这是伟大的,但没有必要)。所以我需要的是:
非常感谢!
答案 0 :(得分:11)
如何在JavaScript中创建具有给定均值和标准差(sd)的随机数列表?
这似乎是一个关于随机创建具有完全指定mean和精确指定standard deviation的数字列表的问题(而不是关于从特定概率分布中绘制随机数的问题)给定的均值和sd)。
一个简单的解决方案是绘制一个随机数列表,然后将此列表移位并缩放以获得所需的均值和sd,如 stats.stackexchange 中的described in this answer。
说,我们在1到10之间生成以下5个随机数:
4.527991433628388
6.3254986488276055
5.123502737960912
7.3331068522336125
9.069573681037484
此列表的平均值为6.475934670737601
,sd为1.8102412442104023
。
然后,我们将列表中的每个数字转换为:
newNum = oldSD * (oldNum - oldMean) / oldSD + newMean
通过将新均值设置为5并将新sd设置为2,我们得到以下转换列表:
2.847863379160965
4.83379450402964
3.505799227476338
5.947025358346529
7.865517530986525
计算此列表的平均值和sd确认它们确实是5和2。
下面是在JavaScript中演示此方法的代码:
// create a list of 5 random numbers between 1 and 10
var list = randomList(5, 1, 10);
// transform the list to have an exact mean of 5 and sd of 2
var newList = forceDescriptives(list, 5, 2);
// display the transformed list and descriptive statistics (mean and sd)
console.log('Transformed random list', newList, descriptives(newList));
// display the original list and descriptive statistics (mean and sd)
console.log('Original random list', list, descriptives(list));
/* demo functions */
function randomList(n, a, b) {
// create a list of n numbers between a and b
var list = [],
i;
for (i = 0; i < n; i++) {
list[i] = Math.random() * (b - a) + a;
}
return list;
}
function descriptives(list) {
// compute mean, sd and the interval range: [min, max]
var mean,
sd,
i,
len = list.length,
sum,
a = Infinity,
b = -a;
for (sum = i = 0; i < len; i++) {
sum += list[i];
a = Math.min(a, list[i]);
b = Math.max(b, list[i]);
}
mean = sum / len;
for (sum = i = 0; i < len; i++) {
sum += (list[i] - mean) * (list[i] - mean);
}
sd = Math.sqrt(sum / (len - 1));
return {
mean: mean,
sd: sd,
range: [a, b]
};
}
function forceDescriptives(list, mean, sd) {
// transfom a list to have an exact mean and sd
var oldDescriptives = descriptives(list),
oldMean = oldDescriptives.mean,
oldSD = oldDescriptives.sd,
newList = [],
len = list.length,
i;
for (i = 0; i < len; i++) {
newList[i] = sd * (list[i] - oldMean) / oldSD + mean;
}
return newList;
}
请注意,由于浮点运算的精度限制,演示有时会计算精确均值和sd略有偏差(如2.0000000000000004
而不是2
的sd)。另请注意,根据所需的平均值和sd,可能无法生成所需范围内的数字列表,如 math.stackexchange 中this answer中所述。这意味着转换后的列表可能与原始列表的范围非常不同。另请注意,可能无法使用给定的整数均值和整数sd生成给定范围内的随机整数列表。