Check the alignment of harvester and applicator/planter for mixed treatment problems where multiple input rates are associated with yield monitor data
Examples
#--- load trial design ---#
data(td_single_input)
#--- check the alignment of harvester and applicator/planter ---#
machine_alignment <- check_alignment(td_single_input)
#--- check the degree of mixed treatment problem ---#
machine_alignment$overlap_data
#> [[1]]
#> ha_strip_id strip_id area ha_area total_intersecting_ha_area
#> <int> <int> <num> <num> <num>
#> 1: 7 19 4802.257 4914.659 4802.257
#> 2: 8 19 4802.257 4914.096 4802.257
#> 3: 9 19 4802.257 4913.533 4802.257
#> 4: 10 18 4800.567 4912.969 4800.567
#> 5: 11 18 4800.567 4912.406 4800.567
#> 6: 12 18 4800.567 4911.843 4800.567
#> 7: 13 17 4798.877 4911.279 4798.877
#> 8: 14 17 4798.877 4910.716 4798.877
#> 9: 15 17 4798.877 4910.153 4798.877
#> 10: 16 16 4797.187 4909.590 4797.187
#> 11: 17 16 4797.187 4909.026 4797.187
#> 12: 18 16 4797.187 4908.463 4797.187
#> 13: 19 15 4795.497 4907.900 4795.497
#> 14: 20 15 4795.497 4907.336 4795.497
#> 15: 21 15 4795.497 4906.773 4795.497
#> 16: 22 14 4793.807 4906.210 4793.807
#> 17: 23 14 4793.807 4905.647 4793.807
#> 18: 24 14 4793.807 4905.083 4793.807
#> 19: 25 13 4792.117 4904.520 4792.117
#> 20: 26 13 4792.117 4903.957 4792.117
#> 21: 27 13 4792.117 4903.393 4792.117
#> 22: 28 12 4790.427 4902.830 4790.427
#> 23: 29 12 4790.427 4902.267 4790.427
#> 24: 30 12 4790.427 4901.703 4790.427
#> 25: 31 11 4788.737 4901.140 4788.737
#> 26: 32 11 4788.737 4900.577 4788.737
#> 27: 33 11 4788.737 4900.014 4788.737
#> 28: 34 10 4779.625 4897.279 4779.625
#> 29: 35 10 4779.625 4892.004 4779.625
#> 30: 36 10 4779.625 4886.726 4779.625
#> 31: 37 9 4763.788 4881.447 4763.788
#> 32: 38 9 4763.788 4876.168 4763.788
#> 33: 39 9 4763.788 4871.851 4763.788
#> 34: 40 8 4756.631 4870.038 4756.631
#> 35: 41 8 4756.631 4868.389 4756.631
#> 36: 42 8 4756.631 4867.833 4756.631
#> 37: 43 7 4755.824 4867.802 4755.824
#> 38: 44 7 4755.824 4867.772 4755.824
#> 39: 45 7 4755.824 4867.741 4755.824
#> 40: 46 6 4755.733 4867.711 4755.733
#> 41: 47 6 4755.733 4867.680 4755.733
#> 42: 48 6 4755.733 4867.650 4755.733
#> 43: 49 5 4755.641 4867.619 4755.641
#> 44: 50 5 4755.641 4867.589 4755.641
#> 45: 51 5 4755.641 4867.558 4755.641
#> 46: 52 4 4755.550 4867.528 4755.550
#> 47: 53 4 4755.550 4867.497 4755.550
#> 48: 54 4 4755.550 4867.467 4755.550
#> 49: 55 3 4755.458 4867.436 4755.458
#> 50: 56 3 4755.458 4867.406 4755.458
#> 51: 57 3 4755.458 4867.375 4755.458
#> 52: 58 2 4755.367 4867.345 4755.367
#> 53: 59 2 4755.367 4867.314 4755.367
#> 54: 60 2 4755.367 4867.284 4755.367
#> 55: 61 1 4755.275 4867.253 4755.275
#> 56: 62 1 4755.275 4867.222 4755.275
#> 57: 63 1 4755.275 4867.192 4755.275
#> ha_strip_id strip_id area ha_area total_intersecting_ha_area
#> intersecting_pct dominant_pct
#> <num> <num>
#> 1: 0.9771291 1
#> 2: 0.9772411 1
#> 3: 0.9773531 1
#> 4: 0.9771212 1
#> 5: 0.9772333 1
#> 6: 0.9773454 1
#> 7: 0.9771134 1
#> 8: 0.9772255 1
#> 9: 0.9773376 1
#> 10: 0.9771055 1
#> 11: 0.9772176 1
#> 12: 0.9773298 1
#> 13: 0.9770976 1
#> 14: 0.9772098 1
#> 15: 0.9773219 1
#> 16: 0.9770897 1
#> 17: 0.9772019 1
#> 18: 0.9773141 1
#> 19: 0.9770818 1
#> 20: 0.9771941 1
#> 21: 0.9773063 1
#> 22: 0.9770739 1
#> 23: 0.9771862 1
#> 24: 0.9772985 1
#> 25: 0.9770660 1
#> 26: 0.9771783 1
#> 27: 0.9772907 1
#> 28: 0.9759755 1
#> 29: 0.9770279 1
#> 30: 0.9780833 1
#> 31: 0.9758968 1
#> 32: 0.9769533 1
#> 33: 0.9778189 1
#> 34: 0.9767132 1
#> 35: 0.9770441 1
#> 36: 0.9771558 1
#> 37: 0.9769962 1
#> 38: 0.9770023 1
#> 39: 0.9770085 1
#> 40: 0.9769958 1
#> 41: 0.9770019 1
#> 42: 0.9770080 1
#> 43: 0.9769954 1
#> 44: 0.9770015 1
#> 45: 0.9770076 1
#> 46: 0.9769949 1
#> 47: 0.9770010 1
#> 48: 0.9770072 1
#> 49: 0.9769945 1
#> 50: 0.9770006 1
#> 51: 0.9770067 1
#> 52: 0.9769941 1
#> 53: 0.9770002 1
#> 54: 0.9770063 1
#> 55: 0.9769936 1
#> 56: 0.9769997 1
#> 57: 0.9770059 1
#> intersecting_pct dominant_pct
#>
#--- visualize the degree of mixed treatment problem ---#
machine_alignment$g_overlap[[1]]
#> `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.
#> Warning: Removed 1 row containing missing values or values outside the scale range
#> (`geom_bar()`).
