0

我需要生成许多(100k 阶)可变大小的随机向量,例如:

sum(v_i)=0sum(v_i^2)=1

对于每个向量,总和为零,平方和为一。

这是我写的代码:

c1 <- rnorm(ki[i])
c1 <- c1-mean(c1)
e1 <- c1/sqrt(sum(c1^2)) 

e1我需要的最终向量在哪里:sum(e1)=0sum(e1^2)=1

是否有可能加快这一代?

4

0 回答 0