From c217b644228685c061d3f3b1650efada1e9fb5f7 Mon Sep 17 00:00:00 2001
From: sandywang <sandywang.rest@gmail.com>
Date: Tue, 9 May 2017 14:42:19 +0800
Subject: [PATCH] Fixed a bug when transmit the paramater of thresholding
 methods

---
 RunFun/gretna_RUN_ThresMat.m | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/RunFun/gretna_RUN_ThresMat.m b/RunFun/gretna_RUN_ThresMat.m
index cb17f9e..a4d2304 100644
--- a/RunFun/gretna_RUN_ThresMat.m
+++ b/RunFun/gretna_RUN_ThresMat.m
@@ -38,8 +38,16 @@ function gretna_RUN_ThresMat(Mat, OutputFile, SType, TType, Thres, NType)
 else
     error('Error: Invalid Matrix Sign');
 end
-
-A=arrayfun(@(t) logical(gretna_R2b(Mat, TType, t)), Thres,...
+if TType==1
+    TFlag='s';
+elseif TType==2
+    TFlag='r';
+else
+    error('Invalid Threshold Method');
+    return
+end
+    
+A=arrayfun(@(t) logical(gretna_R2b(Mat, TFlag, t)), Thres,...
     'UniformOutput', false);
 if NType==2
     A=cellfun(@(bin) bin.*Mat, A, 'UniformOutput', false);