mmollina / MAPpoly

Compare a5d949e ... +0 ... f82defb

No flags found

Use flags to group coverage reports by test type, project and/or folders.
Then setup custom commit statuses and notifications for each flag.

e.g., #unittest #integration

#production #enterprise

#frontend #backend

Learn more about Codecov Flags here.

Showing 1 of 1 files from the diff.

@@ -514,7 +514,7 @@
Loading
514 514
        seq.num <- tail.temp$maps[[i]]$seq.num 
515 515
      if(length(seq.num) < extend.tail)
516 516
        seq.num <- tail(cur.map$maps[[i]]$seq.num, extend.tail)
517 -
      ## If the tail do not contain the marker responsable for carrying 
517 +
      ## If the tail do not contain the marker responsible for carrying 
518 518
      ## multiple linkage phases through the rounds of insertion,
519 519
      ## extend the tail so it contains that marker
520 520
      if(max(check_ls_phase(all.ph)) >= length(seq.num)){
@@ -537,15 +537,15 @@
Loading
537 537
    ## information obtained in the current round
538 538
    twopt.phase.number<-length(input.ph$config.to.test)
539 539
    ## If this number is higher than phase.number.limit,
540 -
    ## procede to the next iteration
540 +
    ## proceed to the next iteration
541 541
    if(length(input.ph$config.to.test) > phase.number.limit) {
542 542
      if(verbose)
543 543
        cat(crayon::italic$yellow(paste0(ct ,": not included (linkage phases)\n", sep = "")))
544 544
      ct <- ct + 1
545 545
      next()
546 546
    }
547 547
    ## Appending the marker to the numeric 
548 -
    ## sequence and makeing a new sequence
548 +
    ## sequence and making a new sequence
549 549
    seq.num <- c(seq.num, input.seq$seq.num[ct])
550 550
    seq.cur <- make_seq_mappoly(input.obj = get(input.seq$data.name, pos=1),
551 551
                                arg = seq.num,
@@ -613,12 +613,14 @@
Loading
613 613
    } else { selected.map <- LOD < thres.hmm }
614 614
    all.ph.temp <- update_ph_list_at_hmm_thres(cur.map.temp, Inf)
615 615
    cur.map.temp$maps <- cur.map.temp$maps[selected.map]
616 +
    id <- !sapply(cur.map.temp$maps, is.null)
617 +
    cur.map.temp$maps <- cur.map.temp$maps[id]
616 618
    cur.map <- cur.map.temp
617 -
    all.ph <- add_mrk_at_tail_ph_list(all.ph, all.ph.temp, M[selected.map,,drop=FALSE])
619 +
    all.ph <- add_mrk_at_tail_ph_list(all.ph, all.ph.temp, M[id,,drop=FALSE])
618 620
    ct <- ct + 1
619 621
  }
620 -
  ##### Reestimating final map ####
621 -
  ## Reestimating final map with higher tolerance
622 +
  ##### Re-estimating final map ####
623 +
  ## Re-estimating final map with higher tolerance
622 624
  seq.final <- make_seq_mappoly(input.obj = get(input.seq$data.name, pos=1), 
623 625
                                arg = as.numeric(names(all.ph$config.to.test[[1]]$P)), 
624 626
                                data.name = input.seq$data.name)

Learn more Showing 1 files with coverage changes found.

Changes in R/filters.R
-1
+1
Loading file...
Files Coverage
calc_genoprob.cpp 100.00%
calc_genoprob_based_on_phased_marker_blocks.cpp 47.95%
calc_loglike_given_map.cpp 88.00%
combinatorial.cpp 100.00%
est_hmm_map_based_on_phased_mrk_blocks.cpp 42.51%
est_map_hmm_given_dose.cpp 83.55%
est_map_hmm_given_prior.cpp 94.21%
genotypic_counts.cpp 100.00%
hmm_elements.cpp 80.00%
mappoly_init.c 0.00%
pairwise_estimation.cpp 100.00%
read_mappoly_vcf.cpp 55.83%
Folder Totals (12 files) 76.71%
Project Totals (59 files) 71.66%
Loading