由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Statistics版 - 真心请教各位统计高手一个关联序列生成问题
相关主题
怎样生成离散分布的随机数啊?有直接的方法吗,比如像randn()这样的函数[合集] 请教:如何在iml里生成有条件的随机数?
R里面的函数set.seed() 的作用讨论关于‘’Probability‘’ 的定义的问题,外行请教
突然对直线拟合的R不明白起来了[R] Generate a set of random number summing up to N
关于R的Simplex的错误信息在R里怎么跳到指定的地方?
如何把model fitting statistics 读出来(R)求助一个问题关于随机数的生成
为什么不能把ABLINE加到散点图中了请问SAS怎么用do loop生成随机数字?
R code 问题【包子】生成RAW SAS DATASET问题
[合集] 关于proc reg 一问, 急~~~~~~~Re: reliability coefficient and its CI
相关话题的讨论汇总
话题: var话题: 随机数话题: rnorm话题: sample话题: sqrt
进入Statistics版参与讨论
1 (共1页)
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
3
y <- rnorm()
y <- x+y
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 的大作中提到】
:
: 可以解释一下为什么这么做吗?谢谢!

1 (共1页)
进入Statistics版参与讨论
相关主题
Re: reliability coefficient and its CI如何把model fitting statistics 读出来(R)
网上看到一道题为什么不能把ABLINE加到散点图中了
一个关于multicollinearity的问题R code 问题
再问multinomial logit[合集] 关于proc reg 一问, 急~~~~~~~
怎样生成离散分布的随机数啊?有直接的方法吗,比如像randn()这样的函数[合集] 请教:如何在iml里生成有条件的随机数?
R里面的函数set.seed() 的作用讨论关于‘’Probability‘’ 的定义的问题,外行请教
突然对直线拟合的R不明白起来了[R] Generate a set of random number summing up to N
关于R的Simplex的错误信息在R里怎么跳到指定的地方?
相关话题的讨论汇总
话题: var话题: 随机数话题: rnorm话题: sample话题: sqrt