Skip to content

Commit

Permalink
adding allele freq distribution in the plots
Browse files Browse the repository at this point in the history
  • Loading branch information
jeffersonfparil committed Apr 17, 2024
1 parent 13b4807 commit c6daaf7
Show file tree
Hide file tree
Showing 9 changed files with 5,600 additions and 5,600 deletions.
1,248 changes: 624 additions & 624 deletions res/cocksfoot-Concordance.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,248 changes: 624 additions & 624 deletions res/cocksfoot-Concordance_high_confidence_data.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,248 changes: 624 additions & 624 deletions res/cocksfoot-Mean_absolute_error.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,248 changes: 624 additions & 624 deletions res/cocksfoot-Mean_absolute_error_high_confidence_data.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,876 changes: 938 additions & 938 deletions res/grape-Concordance.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,876 changes: 938 additions & 938 deletions res/grape-Mean_absolute_error.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
4 changes: 2 additions & 2 deletions res/perf_plot.R
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,7 @@ plot_metrics = function(df, dataset, vec_2_metrics=c("mae_frequencies", "concord
q = seq(0, 1, by=0.1)
for (algo in vec_algorithm) {
# algo = vec_algorithm[1]
par(mar=c(5,5,1,1))
par(mar=c(5,5,1,6))
idx_col_mae = which(grepl("^mae_", colnames(subdf)))
plot(x=c(0,1), y=range(subdf[, idx_col_mae], na.rm=TRUE), type="n", xlab="Observed allele frequencies (5% MAF)", ylab="Imputation error\n(mean absolute error)", main=algo, las=1)
colour = vec_colours[vec_algorithm==algo]
Expand Down Expand Up @@ -191,7 +191,7 @@ plot_metrics = function(df, dataset, vec_2_metrics=c("mae_frequencies", "concord
colour = vec_colours[vec_algorithm==algo]
idx_col = which(grepl("var_x_mae_", colnames(subdf)))
n_bins = length(idx_col) / 2
par(mar=c(5,5,1,1))
par(mar=c(5,5,1,6))
idx_temp = which(subdf$algo == algo)
y_temp = data.frame(y=unlist(subdf[idx_temp, idx_col[seq(2, length(idx_col), by=2)]]), b=rep(1:n_bins, each=length(idx_temp)), m=rep(subdf$missing_rate[idx_temp], time=n_bins))
y_temp = aggregate(y ~ b + m, data=y_temp, FUN=mean, na.rm=TRUE)[,3]
Expand Down
1,226 changes: 613 additions & 613 deletions res/soybean-Concordance.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
1,226 changes: 613 additions & 613 deletions res/soybean-Mean_absolute_error.svg
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.

0 comments on commit c6daaf7

Please sign in to comment.