Révision dadb6a4d src/R/nucleominer.R
b/src/R/nucleominer.R | ||
---|---|---|
596 | 596 |
tmp_nuc[[paste("index_nuc_", strain_ref2, sep="")]] = index_nuc_strain_ref2 |
597 | 597 |
# common |
598 | 598 |
tmp_nuc[["llr_score"]] = signif(llr_score,5) |
599 |
tmp_nuc[["common_wp_pval"]] = signif(1-pchisq(2*tmp_nuc[["llr_score"]], df=2),5) |
|
600 |
tmp_nuc[["diff"]] = abs(abs((roi_strain_ref2$begin - roi_strain_ref2$end)) / 2 - abs((nuc_strain_ref2$lower_bound - nuc_strain_ref2$upper_bound)) / 2) |
|
599 | 601 |
common_nuc = dfadd(common_nuc, tmp_nuc) |
600 | 602 |
} |
601 | 603 |
} |
Formats disponibles : Unified diff