Previous changeset 24:9955b23db68c (2016-12-02) Next changeset 26:ce25fb581ca3 (2016-12-05) |
Commit message:
Uploaded |
modified:
shm_csr.r |
b |
diff -r 9955b23db68c -r 8d1c4c75f81b shm_csr.r --- a/shm_csr.r Fri Dec 02 05:51:13 2016 -0500 +++ b/shm_csr.r Mon Dec 05 08:26:01 2016 -0500 |
[ |
@@ -462,6 +462,12 @@ frequency_bins_data_by_class = frequency_bins_data +frequency_bins_data_by_class = frequency_bins_data_by_class[order(frequency_bins_data_by_class$best_match_class, frequency_bins_data_by_class$frequency_bins),] + +frequency_bins_data_by_class$frequency_bins = gsub("-", " to ", frequency_bins_data_by_class$frequency_bins) +frequency_bins_data_by_class[frequency_bins_data_by_class$frequency_bins == "20", c("frequency_bins")] = "20 or higher" +frequency_bins_data_by_class[frequency_bins_data_by_class$frequency_bins == "0", c("frequency_bins")] = "0 or lower" + write.table(frequency_bins_data_by_class, "frequency_ranges_classes.txt", sep="\t",quote=F,row.names=F,col.names=T) frequency_bins_data = data.frame(data.table(dat)[, list(frequency_count=.N), by=c("best_match", "best_match_class", "frequency_bins")]) @@ -472,6 +478,11 @@ frequency_bins_data$frequency = round(frequency_bins_data$frequency_count / frequency_bins_data$class_sum * 100, 2) +frequency_bins_data = frequency_bins_data[order(frequency_bins_data$best_match, frequency_bins_data$frequency_bins),] +frequency_bins_data$frequency_bins = gsub("-", " to ", frequency_bins_data$frequency_bins) +frequency_bins_data[frequency_bins_data$frequency_bins == "20", c("frequency_bins")] = "20 or higher" +frequency_bins_data[frequency_bins_data$frequency_bins == "0", c("frequency_bins")] = "0 or lower" + write.table(frequency_bins_data, "frequency_ranges_subclasses.txt", sep="\t",quote=F,row.names=F,col.names=T) |