symb2I("A",c("A","B","C"))
#> [,1] [,2] [,3]
#> [1,] 1 0 0
symb2I("A+B",c("A","B","C"))
#> [,1] [,2] [,3]
#> [1,] 1 1 0
symb2I("NoA",c("A","B","C"))
#> [,1] [,2] [,3]
#> [1,] 0 0 0
symb2I(c("A","B+C"),c("A","B","C")) # gives a 2 by 3 matrix.
#> [,1] [,2] [,3]
#> [1,] 1 0 0
#> [2,] 0 1 1