# HG changeset patch
# User dereeper
# Date 1707041239 0
# Node ID ad5e3d3c809e4b3ecd3d6ac60e901347e1ac8259
# Parent  e44eb2ba80ecbb275ab92acb33e05e81d823210a
Uploaded

diff -r e44eb2ba80ec -r ad5e3d3c809e Snmf.pl
--- a/Snmf.pl	Tue Jan 08 11:22:09 2019 -0500
+++ b/Snmf.pl	Sun Feb 04 10:07:19 2024 +0000
@@ -68,7 +68,7 @@
 	print E "K=$k	$ent\n";
 
 	print $O2 "Indiv";
-	print $O3 "Indiv;Group\n";
+	print $O3 "Indiv\tGroup\n";
 	for (my $j = 0; $j <$k; $j++){
 		print $O2 "	Q$j";
 	}
@@ -89,7 +89,7 @@
 		my $ind = $individuals[$i];
 		for (my $j = 0; $j <$k; $j++){
 			my $val = $infos[$j];
-			if ($val > 0.5){$group = "Q$j";}
+			if ($val > ($threshold/100)){$group = "Q$j";}
 		}
 		if ($ind){
 			$hash_indv{$ind} = join("	",@infos);
@@ -103,7 +103,7 @@
 		my @inds = split(",",$hash_groupes{$group}{"ind"});
 		foreach my $ind(@inds){
 			if ($ind =~/\w+/){
-				print $O3 "$ind;$group\n";
+				print $O3 "$ind\t$group\n";
 				print $O2 $ind."	".$hash_indv{$ind}. "\n";
 			}
 		}