| {0.9, 0, 0.05, 0.05},
{0.3, 0.3, 0.3, 0.1},
{0.25, 0.25, 0.25, 0.25},
{0, 0.9, 0, 0.1},
{0.2, 0.2, 0.2, 0.4} };
PPM ppm = new PPM(ppm_matrix, "cAVNCt PPM (slightly different from another cAVNCt matrix)");
PWM pwm_from_ppm = ppm.to_pwm(background, ppm_effective_count);
pwmCollection.add(pwm_from_ppm);
//////////////////////////////////////////////////////////////////////////////////////////////////
// Now we should get collection of SNPs with sequence around them. Each SNP should have only two variants (if one need 3 or 4 variants, it's possible to add several SNPs: 1-2, 2-3, 1-3 for instance)
|