RSCHMM_ComputeSenSpec.m
1f50b9a012abmaster

RSCHMM_ComputeSenSpec.m
View Options

			function [Sens,Spec,Acc] = RSCHMM_ComputeSenSpec(G_actual,G_GT)

			TP = 0;
			TN = 0;
			FP = 0;
			FN = 0;

			Acc = 0;

			N = size(G_actual,1);

			for i = 1:N
			for j = 1:N
			if i ~= j

			% In the case of a real edge
			if G_GT(i,j) ~= 0

			if G_actual(i,j) ~= 0
			TP = TP + 1;

			if sign(G_actual(i,j)) == sign(G_GT(i,j))
			Acc = Acc + 1;
			end
			elseif G_actual(i,j) == 0
			FN = FN + 1;
			end

			% In the case of no edge
			else
			if G_actual(i,j) ~= 0
			FP = FP + 1;
			elseif G_actual(i,j) == 0
			Acc = Acc + 1;
			TN = TN + 1;
			end
			end
			end
			end
			end

			Acc = Acc/(N^2-N);
			Sens = TP/(TP+FN);
			Spec = TN/(TN+FP);


			end