-
Notifications
You must be signed in to change notification settings - Fork 4
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
🔖 release v0.9.7 #28
base: master
Are you sure you want to change the base?
🔖 release v0.9.7 #28
Conversation
a new output variable has been added to keep track of the average percentage of samples in which SNPs are "in runs" the new output variable is named avg_pct and is a column added at the end of the output data.frame
snpInsideRuns(Cpp) is now internal
check runs where SNP are major or equal the threshold specified
rewrite code to calculate percentages correctly
correctly skip empty subset
no more copy values from dataframe
…Runs Set guides(fill="none") since ggplot2 deprecated guides(fill=FALSE)
👽 set guides(fill="none") in plot_InbreedingChr and plot_DistributionRuns
Codecov ReportAll modified and coverable lines are covered by tests ✅
Additional details and impacted files@@ Coverage Diff @@
## master #28 +/- ##
==========================================
+ Coverage 51.11% 59.01% +7.89%
==========================================
Files 6 6
Lines 1254 1359 +105
==========================================
+ Hits 641 802 +161
+ Misses 613 557 -56 ☔ View full report in Codecov by Sentry. |
⚡ refactor tableRuns to better manage memory
- added a comment on the need for chromosomes to be numerically coded - in the vignette, the description of the `summary_ROH_mean_chr` dataset from `summaryRuns()` has been corrected - the documentation of the function `summaryRuns` has been expanded
Ready to merge for me |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think that issue #41 should be addressed before merging
@@ -243,7 +243,16 @@ Froh_inbreedingClass <- function(runs, mapFile, Class=2){ | |||
#' | |||
#' @return A list of dataframes containing the most relevant descriptives | |||
#' statistics on detected runs. The list conveniently contains 9 dataframes that can | |||
#' be used for further processing and visualization, or can be written out to text files | |||
#' be used for further processing and visualization, or can be written out to text files: | |||
#' 1) summary_ROH_count_chr: n. of runs per chromosome and breed/group; |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This list will be rendered as a plain text. Consider using list with roxygen, for example:
#' \enumerate{
#' \item{summary_ROH_count_chr: }{n. of runs per chromosome and breed/group}
#' \item{summary_ROH_percentage_chr: }{percent distribution of runs per chromosome in each breed/group (sum to 1)}
#' \item{summary_ROH_count: }{n. of runs per size-class (Mb) in each breed/group}
#' <add all the other items>
#' }
#' 6) summary_ROH_mean_class: average size of runs (Mb) per size-class (Mb) in each breed/group | ||
#' 7) result_Froh_genome_wide: genome-wide inbreeding ($F_{ROH}$) for each individual in the dataset | ||
#' 8) result_Froh_chromosome_wide inbreeding ($F_{ROH}$) for each individual (and chromosome) in the dataset | ||
#' 9) result_Froh_class: genome-wide inbreeding ($F_{ROH}$) for each individual in the dataset per size-class (Mb) of runs |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
There could be a 10th item when providing snpInRuns=TRUE
as parameter
#' | ||
#' @param genotypeFile genotype (.ped) file path | ||
#' @param mapFile map file (.map) file path | ||
#' @param mapFile map file (.map) file path (chromosomes should be numeric) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
will be better to check chromosomes in readMapFile
directly and throw an error here?
Improving the output of the function
tableRuns
Description
On the
devel
branch, the code of the functiontableRuns
has been modified to add a column namedavg_pct
to the output dataframeRelated Issue
close #43
close #39
close #31
close #36
close #17
Motivation and Context
It adds a piece of relevant information to the output of the function
tableRuns
(ROH/HRR islands)How Has This Been Tested?
Types of changes
Checklist: