This function calculates the suitability scores and class of the land units.
suit(
crop,
terrain = NULL,
water = NULL,
temp = NULL,
mf = "triangular",
sow_month = NULL,
minimum = NULL,
maximum = "average",
interval = NULL,
sigma = NULL
)a string for the name of the crop;
a data frame for the terrain characteristics of the input land units;
a data frame for the water characteristics of the input land units;
a data frame for the temperature characteristics of the input land units;
membership function with default assigned to "triangular"
fuzzy model. Other fuzzy models included are "trapezoidal" and
"gaussian".
sowing month of the crop. Takes integers from 1 to 12 (inclusive), representing the twelve months of the year. So if sets to 1, the function assumes sowing month to be January.
factor's minimum value. If NULL (default), minimum is
set to 0. But if numeric of length one, say 0.5, then minimum
is set to 0.5, for all factors. To set multiple minimums for multiple factors,
simply concatenate these into a numeric vector, the length of this vector should be equal
to the number of factors in input land units parameters. However, it can also be set to
"average", please refer to the online documentation for more, link in the "See Also" section below.
maximum value for factors. To set multiple maximums for multiple factors,
simply concatenate these into a numeric vector, the length of this vector should be equal
to the number of factors in input land units parameters. However, it can also be set to
"average", please refer to the online documentation for more, link in the "See Also" section below.
domains for every suitability class (S1, S2, S3). If fixed (NULL), the
interval would be 0 to 25% for N (Not Suitable), 25% to 50% for S3 (Marginally Suitable),
50% to 75% for S2 (Moderately Suitable), and 75% to 100% for (Highly Suitable). If "unbias",
the package will take into account the shape of the membership function, and provide the
appropriate suitability class intervals. However, it can also be customized by specifying the
limits of the suitability classes. Please refer to the online documentation for more, link in the "See Also" section below.
If mf = "gaussian", then sigma represents the constant sigma in the
Gaussian formula.
A list of outputs of target characteristics, with the following components:
"terrain" - a list of outputs for terrain characteristics
"soil" - a list of outputs for soil characteristics
"water" - a list of outputs for water characteristics
"temp" - a list of outputs for temperature characteristics
These components are only available when specified as the target characteristics in either
of the arguments above, that is, if terrain argument is specified above, then the "terrain"
and "soil" components will be available in the output list. This is also true if water and temp
are specified in the arguments above.
Each of the components returned above contains a list of outputs as well with the following components:
"Factors Evaluated" - a character of factors that matched between the input land units factor and the targetted crop requirement factor
"Suitability Score" - a data frame of suitability scores for each of the matched factors
"Suitability Class" - a data frame of suitability classes for each of the matched factors
"Factors' Minimum Values" - a numeric of minimum values used in the membership function for computing the suitability scores
"Factors' Minimum Values" - a numeric of maximum values used in the membership function for computing the suitability scores
"Factors' Weights" - a numeric of weights of the factors specified in the input crop requirements
"Crop Evaluated" - a character of the name of the targetted crop requirement dataset
https://alstat.github.io/ALUES/
# NOT RUN {
library(ALUES)
rice_suit <- suit("ricebr", water=MarinduqueWater, temp=MarinduqueTemp, sow_month = 1)
lapply(rice_suit[["water"]], function(x) head(x)) # access results for water suitability
lapply(rice_suit[["temp"]], function(x) head(x)) # access results for temperature suitability
rice_suit <- suit("ricebr", terrain=MarinduqueLT)
lapply(rice_suit[["terrain"]], function(x) head(x))
lapply(rice_suit[["soil"]], function(x) head(x))
# }
Run the code above in your browser using DataLab