N4212cop(0.4,0.6,para=1)
PSP(0.4,0.6) # same two values
N4212cop(0.4,0.6,para=10)
D <- simCOP(n=400, cop=N4212cop, para=2)
D <- simCOP(n=400, cop=N4212cop, para=10)
D <- simCOP(n=400, cop=N4212cop, para=100)
Run the code above in your browser using DataLab