w*****a 发帖数: 218 | 1 给定一随机数列, 例如 X = (0.1,0.3,0.02,0.6,0.22,0.9,...)
需要生成另一随机数列 Y = (y1,y2,y3,y4,...)
使 Y 与 X 的关联等于 c
怎么作? | w*****a 发帖数: 218 | 2 版上的高人都在哪啊?
帮在下一把!
【在 w*****a 的大作中提到】 : 给定一随机数列, 例如 X = (0.1,0.3,0.02,0.6,0.22,0.9,...) : 需要生成另一随机数列 Y = (y1,y2,y3,y4,...) : 使 Y 与 X 的关联等于 c : 怎么作?
| s*r 发帖数: 2757 | | D******n 发帖数: 2836 | 4 baozi please.
if your sequence is long enough, do this
Y <- a*X+E
where a is a constant
E is a random sample from normal(0,B).(B=sigma^2)
let A=var(X);
to get the desire c. u can manipulate a and B, s.t
a^2 *A /(a^2*A+B) = c^2;
for example, if u sample E from N(0,1)
u should use a = sqrt(p^2 * B /(1-p^2) /A );
【在 w*****a 的大作中提到】 : 给定一随机数列, 例如 X = (0.1,0.3,0.02,0.6,0.22,0.9,...) : 需要生成另一随机数列 Y = (y1,y2,y3,y4,...) : 使 Y 与 X 的关联等于 c : 怎么作?
| w******a 发帖数: 25 | 5
可以解释一下为什么这么做吗?谢谢!
【在 D******n 的大作中提到】 : baozi please. : if your sequence is long enough, do this : Y <- a*X+E : where a is a constant : E is a random sample from normal(0,B).(B=sigma^2) : let A=var(X); : to get the desire c. u can manipulate a and B, s.t : a^2 *A /(a^2*A+B) = c^2; : for example, if u sample E from N(0,1) : u should use a = sqrt(p^2 * B /(1-p^2) /A );
| D******n 发帖数: 2836 | 6 from definition of correlation coefficient
c= cov(x,y)/sqrt(var(x)*var(y))
可以解释一下为什么这么做吗?谢谢!
【在 w******a 的大作中提到】 : : 可以解释一下为什么这么做吗?谢谢!
|
|