This function determines the most limiting factor based on von Liebig law of the minimum given results of the predicted boundary line values for the different factors of interest. Boundary lines for various factors are fitted and the factor that predicts the minimum response for a particular point is considered as the most limiting factor (Casanova et al. 1995).
Examples
N<-rnorm(10,50,5)#assuming these are predicted responses using the fitted BL for N,P,K
K<-rnorm(10,50,4)
P<-rnorm(10,50,6)
limfactor(N,K,P)
#> [[1]]
#> Rs Lim_factor
#> 1 43.82230 P
#> 2 51.09964 K
#> 3 40.12164 K
#> 4 45.05137 P
#> 5 39.34444 P
#> 6 41.56472 P
#> 7 43.68471 N
#> 8 48.82520 N
#> 9 44.86946 K
#> 10 52.25433 K
#>
#> [[2]]
#> [1] 61.92102
#>