如何通过python实现R矩阵功能

Question

当我使用'R'代码时 要实现此功能， 如下：

mir=   #data.frame
  s S
1 a A
2 b B
3 c C
4 d D
5 a B
6 b D
7 d C
class(mir)
"data.frame"

table(mir)
   S
s   A B C D
  a 1 1 0 0
  b 0 1 0 1
  c 0 0 1 0
  d 0 0 1 1

mir %*% t(mir)
   s
s   a b c d
  a 2 1 0 0
  b 1 2 0 1
  c 0 0 1 1
  d 0 1 1 2

结果表示a-a共享（2：A + B）a-b共享（1：B）...所以我想实现该矩阵，如何通过“ python”实现此功能？