r - Doing a matrix with absolute frequencies -
i have 2 matrices: matrix 1, vector :
matrix(c(0.16, 0.24, 0.16, 0.08, 0.09, 0.12, 0.06, 0.04, 0.04, 0.01))
actually, created using function "for", utilizing punnett square frequencies 0.4, 0.3, 0.2, 0.1.
and second, have matrix 2:
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] "a 1 1" "b 1 b 1" "c 1 c 1" "d 1 d 1" "e 1 e 1" "f 1 f 1" "g 1 g 1" "h 1 h 1" "i 1 1" "j 1 j 1" [2,] "a 1 2" "b 1 b 2" "c 1 c 2" "d 1 d 2" "e 1 e 2" "f 1 f 2" "g 1 g 2" "h 1 h 2" "i 1 2" "j 1 j 2" [3,] "a 1 3" "b 1 b 3" "c 1 c 3" "d 1 d 3" "e 1 e 3" "f 1 f 3" "g 1 g 3" "h 1 h 3" "i 1 3" "j 1 j 3" [4,] "a 1 4" "b 1 b 4" "c 1 c 4" "d 1 d 4" "e 1 e 4" "f 1 f 4" "g 1 g 4" "h 1 h 4" "i 1 4" "j 1 j 4" [5,] "a 2 2" "b 2 b 2" "c 2 c 2" "d 2 d 2" "e 2 e 2" "f 2 f 2" "g 2 g 2" "h 2 h 2" "i 2 2" "j 2 j 2" [6,] "a 2 3" "b 2 b 3" "c 2 c 3" "d 2 d 3" "e 2 e 3" "f 2 f 3" "g 2 g 3" "h 2 h 3" "i 2 3" "j 2 j 3" [7,] "a 2 4" "b 2 b 4" "c 2 c 4" "d 2 d 4" "e 2 e 4" "f 2 f 4" "g 2 g 4" "h 2 h 4" "i 2 4" "j 2 j 4" [8,] "a 3 3" "b 3 b 3" "c 3 c 3" "d 3 d 3" "e 3 e 3" "f 3 f 3" "g 3 g 3" "h 3 h 3" "i 3 3" "j 3 j 3" [9,] "a 3 4" "b 3 b 4" "c 3 c 4" "d 3 d 4" "e 3 e 4" "f 3 f 4" "g 3 g 4" "h 3 h 4" "i 3 4" "j 3 j 4" [10,] "a 4 4" "b 4 b 4" "c 4 c 4" "d 4 d 4" "e 4 e 4" "f 4 f 4" "g 4 g 4" "h 4 h 4" "i 4 4" "j 4 j 4"
i want new table (lines=100, columns=10) each column of matrix 2, having frequencies of matrix 1. have 16 times "a1a1", 24 times "a1a2" , on. same others columns. want frequencies (matrix 1) , number of columns lines of matrix 2.
try
out <- matrix(rep(c(m2),m1[row(m2)]*100),100,10) dim(out) #[1] 100 10
data
m1 <- matrix(c(0.16, 0.24, 0.16, 0.08, 0.09, 0.12, 0.06, 0.04, 0.04, 0.01)) m2 <- structure(c("a 1 1", "a 1 2", "a 1 3", "a 1 4", "a 2 2", "a 2 3", "a 2 4", "a 3 3", "a 3 4", "a 4 4", "b 1 b 1", "b 1 b 2", "b 1 b 3", "b 1 b 4", "b 2 b 2", "b 2 b 3", "b 2 b 4", "b 3 b 3", "b 3 b 4", "b 4 b 4", "c 1 c 1", "c 1 c 2", "c 1 c 3", "c 1 c 4", "c 2 c 2", "c 2 c 3", "c 2 c 4", "c 3 c 3", "c 3 c 4", "c 4 c 4", "d 1 d 1", "d 1 d 2", "d 1 d 3", "d 1 d 4", "d 2 d 2", "d 2 d 3", "d 2 d 4", "d 3 d 3", "d 3 d 4", "d 4 d 4", "e 1 e 1", "e 1 e 2", "e 1 e 3", "e 1 e 4", "e 2 e 2", "e 2 e 3", "e 2 e 4", "e 3 e 3", "e 3 e 4", "e 4 e 4", "f 1 f 1", "f 1 f 2", "f 1 f 3", "f 1 f 4", "f 2 f 2", "f 2 f 3", "f 2 f 4", "f 3 f 3", "f 3 f 4", "f 4 f 4", "g 1 g 1", "g 1 g 2", "g 1 g 3", "g 1 g 4", "g 2 g 2", "g 2 g 3", "g 2 g 4", "g 3 g 3", "g 3 g 4", "g 4 g 4", "h 1 h 1", "h 1 h 2", "h 1 h 3", "h 1 h 4", "h 2 h 2", "h 2 h 3", "h 2 h 4", "h 3 h 3", "h 3 h 4", "h 4 h 4", "i 1 1", "i 1 2", "i 1 3", "i 1 4", "i 2 2", "i 2 3", "i 2 4", "i 3 3", "i 3 4", "i 4 4", "j 1 j 1", "j 1 j 2", "j 1 j 3", "j 1 j 4", "j 2 j 2", "j 2 j 3", "j 2 j 4", "j 3 j 3", "j 3 j 4", "j 4 j 4"), .dim = c(10l, 10l))
Comments
Post a Comment