}
}
}
- else if (matrix == "sim9") { }
+ else if (matrix == "sim9" || matrix == "sim2") { }
else {
m->mothurOut("[ERROR]: No model selected! \n");
m->control_pressed = true;
//burn-in for sim9
if(matrix == "sim9") {
- for(int i=0;i<10000;i++) trial.swap_checkerboards (co_matrix, rowtotal, columntotal, ncols, nrows);
+ for(int i=0;i<10000;i++) trial.swap_checkerboards (co_matrix, ncols, nrows);
}
//populate null matrix from probability matrix, do this a lot.
else if(matrix == "sim4") {
for(int i=0;i<nrows;i++) {
- previous = 0.0;
count = 0;
while(count < rowtotal[i]) {
previous = 0.0;
//swap_checkerboards takes the original matrix and swaps checkerboards
else if(matrix == "sim9") {
- trial.swap_checkerboards (co_matrix, rowtotal, columntotal, ncols, nrows);
+ trial.swap_checkerboards (co_matrix, ncols, nrows);
}
else {
m->mothurOut("[ERROR]: No null model selected!\n\n"); m->control_pressed = true;
}
}
/**************************************************************************************************/
-int TrialSwap2::swap_checkerboards (vector<vector<int> > &co_matrix, vector<int> rowtotal, vector<int> columntotal, int ncols, int nrows)
+int TrialSwap2::swap_checkerboards (vector<vector<int> > &co_matrix, int ncols, int nrows)
{
try {
- int ncols = co_matrix[0].size(); int nrows = co_matrix.size();
+ //int ncols = co_matrix[0].size(); int nrows = co_matrix.size();
int i, j, k, l;
i = m->getRandomIndex(nrows-1);
while((j = m->getRandomIndex(nrows-1) ) == i ) {;if (m->control_pressed) { return 0; }}