Commits

opticall committed 2a219e3

revised X chromosome inference to allow male hets

  • Participants
  • Parent commits 76a3850

Comments (0)

Files changed (1)

File opticall/opticall.cpp

 			mat_tijs(i,j) = mvt_pdf(data.block(i,0,1,dim), mus.block(j,0,1,dim), staticcovars[j], vs(j,0))*((1-p_proportion)*js(j,0)+p_proportion*probs_j(i,j));
 		}
         
-        //gender specific fixes
-        if (chrom_type == "X" && snp_genders(i,0) == 1)
-        {
-            mat_tijs(i,2) = 0.0;
-        }
+        //gender specific fixes - THIS DIDN'T GIVE GOOD RESULTS, so commented out
+        //if (chrom_type == "X" && snp_genders(i,0) == 1)
+        //{
+        //    mat_tijs(i,2) = 0.0;
+        //}
         
 		MatrixXd unnorm_block = mat_tijs.block(i,0,1,js.rows());
 		mat_tijs.block(i,0,1,js.rows()) = unnorm_block*(1.0/unnorm_block.sum());
 	//char * output_prob_ext = ".probs";
 	//char * output_call_ext = ".calls";
 
-	string version = "0.6.3";
+	string version = "0.6.4";
 	cout << "opticall version " << version << endl;
 	cout << "thank you for choosing opticall for your genotyping needs" << endl;