disease_vs_normal_KL: Identify perturbed ligand-receptor interactions between two...
In yuliangwang/talklr: Uncover Perturbed Ligand-Receptor Pairs

This function calculates Kullback-Leibler divergence of ligand-receptor interaction scores between two conditions (e.g., disease vs. normal)

disease_vs_normal_KL(
  ligand_mat_disease,
  receptor_mat_disease,
  ligand_mat_normal,
  receptor_mat_normal,
  pseudo_count,
  method
)

`ligand_mat_disease`	ligand expression matrix in the disease/perturbed condition
`receptor_mat_disease`	receptor expression matrix in the disease/perturbed condition. The same ligand-receptor pair should be in the same order as ligand_mat_disease
`ligand_mat_normal`	ligand expression matrix in the normal/baseline condition. Same order as ligand_mat_disease
`receptor_mat_normal`	receptor expression matrix in the normal/baseline coniditon. The same ligand-receptor pair should be in the same order as ligand_mat_normal
`pseudo_count`	a small value added to the expression matrices to avoid dividing by 0.
`method`	"product" or "min": ligand-receptor interaction score is either defined by multiplying ligand and receptor expression values or take the minimum