CODING IN R STUDIO!

How can I write a code that computes the number of genes in common for each pair of phenotypes and then store it in a matrix with dimensions N x N where N is the # of phenoypes in the matrix?

Example of data I am working with:

	A	B	C	D	E	F	G	H	I	J	K	L	M	N
Gene1	0	0	0	0	1	0	1	1	1	0	0	0	0	0
Gene2	1	0	0	0	0	0	0	0	0	1	0	0	1	0
Gene3	0	0	1	1	0	0	0	0	1	0	0	1	0	0
Gene4	0	1	1	0	0	0	1	0	0	0	1	1	0	1
Gene5	1	0	0	0	1	1	0	0	0	0	0	0	1	1
Gene6	0	0	1	1	0	0	0	1	1	0	0	0	0	0
Gene7	0	0	0	0	1	0	0	0	1	1	0	0	1	0
Gene8	0	0	0	0	0	0	1	0	0	0	0	1	0	0
Gene9	1	0	1	0	0	0	0	1	0	0	0	0	0	1

'1' means that the gene has the phenotype, '0' means that the gene does not have it.