How can I calculate degree centrality correctly when there are multiple edges within two nodes?

I am a novice who is working on a data to calculate directed degree centrality by editing R code my senior already made.