我有以下代码:
1. A=rnorm(100,mean=0, sd=30)
2. B1=A+rnorm(100,mean=0,sd=10); B2=A+rnorm(100,mean=0,sd=30); B3=A+rnorm(100,mean=0,sd=50)
3. d=data.frame(A, B1,B2,B3)
4. c=cor(d)^2
5. r=c[2,1]/c[4,1]
我想重述几句1-4
。直到满足条件1.6 <=r <=3.5
。
即我想生成A, B1, B2, B3
,以便它们满足r
值之间的标准1.6
,3.5
如果值r
落在该范围之外,我想重复1-4.
,直到它落在该范围内。