Commits

opticall committed acd2d68

opticall 0.6.5 - fixed snpcount when using inblock

Comments (0)

Files changed (1)

opticall/opticall.cpp

 	//char * output_prob_ext = ".probs";
 	//char * output_call_ext = ".calls";
 
-	string version = "0.6.4";
+	string version = "0.6.5";
 	cout << "opticall version " << version << endl;
 	cout << "thank you for choosing opticall for your genotyping needs" << endl;
 
 	individuals = readHeaderInfoFromIntFile(fname.c_str());
 	cout << "number of inds is " << individuals.size() << endl;
 
+        int snpcount_block;
+        snpcount_block = snpCount(fname_block.c_str(),true);
 	
 	vector < vector< string > > snpinfo2;
 
 	/********************************/
 	int blockattempts = 0;
 	do {
-
-		MatrixXd sample  =  fetchRandomIntensities(fname_block, sample_block_size,  snpcount, individuals, combineExclusionLists(initial_outliers,excluded_samples)    );
+                cout << "file being used for block is " << fname_block << endl;
+		MatrixXd sample  =  fetchRandomIntensities(fname_block, sample_block_size,  snpcount_block, individuals, combineExclusionLists(initial_outliers,excluded_samples)    );
 		
 		//only transform the sample during the rescue stage
 		//sample = transform_data(sample);
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.