|
先试小样本,看看是不是你想要的.
母体A=15
子体n=10
simulations:s=6
A=c(seq(1:15))
data=matrix(0,nrow = n,ncol=s)
> data
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 8 12 3 13 7 6
[2,] 3 14 1 2 11 8
[3,] 9 1 11 10 6 12
[4,] 14 4 14 6 4 7
[5,] 11 2 9 8 13 13
[6,] 13 15 5 3 2 9
[7,] 10 11 10 1 15 4
[8,] 2 3 12 4 10 15
[9,] 6 9 4 15 8 11
[10,] 5 10 2 11 5 10
[,1] [,2]
[1,] 10 0
[2,] 0 5
[1] "a"
[1] 10
[1] "b"
[1] 0
[1] "c"
[1] 0
[1] "d"
[1] 5
[1] "Yule'Q"
[1] 1
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 0 0 0 0 0
[2,] 0 0 0 0 0 0
[3,] 0 0 0 0 0 0
[4,] 0 0 0 0 0 0
[5,] 0 0 0 0 0 0
[6,] 0 0 0 0 0 0
[,1] [,2]
[1,] 6 4
[2,] 4 1
[1] "a"
[1] 6
[1] "b"
[1] 4
[1] "c"
[1] 4
[1] "d"
[1] 1
[1] "Yule'Q"
[1] -0.4545455
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 -0.4545455 0 0 0 0
[2,] 0 0.0000000 0 0 0 0
[3,] 0 0.0000000 0 0 0 0
[4,] 0 0.0000000 0 0 0 0
[5,] 0 0.0000000 0 0 0 0
[6,] 0 0.0000000 0 0 0 0
[,1] [,2]
[1,] 7 3
[2,] 3 2
[1] "a"
[1] 7
[1] "b"
[1] 3
[1] "c"
[1] 3
[1] "d"
[1] 2
[1] "Yule'Q"
[1] 0.2173913
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1 -0.4545455 0.2173913 0 0 0
[2,] 0 0.0000000 0.0000000 0 0 0
[3,] 0 0.0000000 0.0000000 0 0 0
[4,] 0 0.0000000 0.0000000 0 0 0
[5,] 0 0.0000000 0.0000000 0 0 0
[6,] 0 0.0000000 0.0000000 0 0 0
..
..
..
> matrix
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1.0000000 -0.4545455 0.2173913 0.2173913 0.2173913 -0.4545455
[2,] -0.4545455 1.0000000 0.9459459 0.2173913 -1.0000000 -0.4545455
[3,] 0.2173913 0.9459459 1.0000000 -0.4545455 -1.0000000 -1.0000000
[4,] 0.2173913 0.2173913 -0.4545455 1.0000000 0.7142857 0.2173913
[5,] 0.2173913 -1.0000000 -1.0000000 0.7142857 1.0000000 0.7142857
[6,] -0.4545455 -0.4545455 -1.0000000 0.2173913 0.7142857 1.0000000
|