我在r
中生成了一个错误u1<-rnorm(n, mean=0, sd=1);
我想生成一个随机的正常变量x~i.i.d.N(n,0,4),但这与u1无关,即corr(x, u1)
为零。
我怎样才能在R?中定义它?
答案 0 :(得分:2)
生成两个不相关的变量:
n <- 1e3
u1 <- rnorm(n, mean = 0, sd = 1)
x <- rnorm(n, mean = 0, sd = 4)
请注意,估计相关系数不会完全为0,但会接近0。
> cor(x, u1)
[1] -0.01885482
答案 1 :(得分:1)
任何两个随机高斯变量的理论相关性为0,只要n足够大。
set.seed(1)
n=1000000
u1<-rnorm(n, mean=0, sd=1)
x<-rnorm(n, mean=0, sd=4)
cor(x, u1)
输出:
0.0007637239