clear;clc; %% Generating Random Networks [ndrugs, ntargets] = size(adjacencyR); nedges = sum(adjacencyR(:)); A=zeros(ndrugs,ntargets); t=randperm(numel(A)); A(t(1:nedges))=1; clearvars t %% Calculate NBI through drug target adjacency matrix %A=[1, 1, 0, 1;0,1, 1, 0; 0, 1, 1,1]; Ky = diag(1./sum(A)); n = size(A,1); m = size(A,2); Ky(isinf(Ky)| isnan(Ky))=0; kx = sum(A,2)'; Nx = repmat(1./kx,n,1); Nx(isinf(Nx) | isnan(Nx))=0; kx(isinf(kx) | isnan(kx))=0; W = (A*Ky)'; W1 = A*W; W2 = Nx .* W1; R = W2*A; clearvars Ky n m kx Nx W W1 W2