private void initSampleAlleles(DuoBaumLevel level, int sampleA, int sampleB) { int m = level.marker(); for (int copy = 0; copy < nCopies; ++copy) { int state = initialRandomState(level); nodeAB1[copy] = level.parentNodeAB1(state); nodeA2[copy] = level.parentNodeA2(state); nodeB2[copy] = level.parentNodeB2(state); nodeValue[copy] = parentSum(level, sampleA, sampleB, state); allelesAB1[copy][m] = level.symbolAB1(state); allelesA2[copy][m] = level.symbolA2(state); allelesB2[copy][m] = level.symbolB2(state); } }