Function to assess the distribution of maximum genotype probabilities (maxP
), if these are available. The function
plots a violin graph showing the distribution of the samples' maxP
.
check_maxP(probgeno_df)
This function does not return any value, is simply a visualisation tool to help assess data quality.
A data frame as read from the scores file produced by function
saveMarkerModels
of R package fitPoly
, or alternatively, a data frame containing the following columns:
Name of the sample (individual)
Name of the marker
Probabilities of dosage score '0'
Probabilities of dosage score '1' etc. (up to max dosage, e.g. P4 for tetraploid population)
Maximum genotype probability identified for a particular individual and marker combination
Most probable dosage for a particular individual and marker combination
Most probable dosage for a particular individual and marker combination, if maxP
exceeds a user-defined threshold (e.g. 0.9), otherwise NA
data("gp_df")
check_maxP(gp_df)
Run the code above in your browser using DataLab