Home > Software design >  Generate square symmetric matrix in R
Generate square symmetric matrix in R

Time:01-19

I have a matrix in R as follows:

dat <- matrix(c(2,3,5,7,8,4), ncol = 6)
colnames(dat) <- c("A(1,1)", "A(1,2)", "A(1,3)", "A(2,2)", "A(2,3)", "A(3,3)")

How can I create a square symmetric matrix based on an apply function that has the following form:

A(1,1) A(1,2) A(1,3)
A(2,1) A(2,2) A(2,3)
A(3,1) A(3,2) A(3,3)

Note that A(1,2)=A(2,1)

CodePudding user response:

This isn't based on apply but it shows how you can use lower.tri and upper.tri to create a symmetric matrix given a vector of elements. It shows the 3x3 case, but will easily generalize to larger n:

dat <- matrix(rep(0,9),nrow = 3)
dat[lower.tri(dat,diag = TRUE)] <- c(2,3,5,7,8,4)
dat[upper.tri(dat)] <- t(dat)[upper.tri(dat)]
print(dat)

Result:

     [,1] [,2] [,3]
[1,]    2    3    5
[2,]    3    7    8
[3,]    5    8    4
  • Related