Usage
gaRoth (
w_1h = range(SFM_metric[,-1][,1]),
w_10h = range(SFM_metric[,-1][,2]) ,
w_100h = range(SFM_metric[,-1][,3]) ,
w_Live_Herb = range(SFM_metric[,-1][,4]) ,
w_Live_Woody = range(SFM_metric[,-1][,5]) ,
s_1h = range(SFM_metric[,-1][,6]) ,
s_10h = range(SFM_metric[,-1][,7]) ,
s_100h = range(SFM_metric[,-1][,8]) ,
s_Live_Herb = range(SFM_metric[,-1][,9]) ,
s_Live_Woody = range(SFM_metric[,-1][,10]) ,
delta = range(SFM_metric[,-1][,11]) ,
mx.dead = range(SFM_metric[,-1][,12]) ,
h_1h = range(SFM_metric[,-1][,13]) ,
h_10h = range(SFM_metric[,-1][,14]) ,
h_100h = range(SFM_metric[,-1][,15]) ,
h_Live_Herb = range(SFM_metric[,-1][,16]) ,
h_Live_Woody = range(SFM_metric[,-1][,17]) ,
m, u, slope, modeltype, obs,
method="rmse", maxiter=50, popSize = 20, pcrossover = 0.8,
pmutation = 0.1, elitism = base::max(1, round(popSize * 0.05)),
...)
Arguments
w_1h
a vector of minimum and maximum fuel load [t/ha] for fuel class 1-hour; (0,0) if fuel class is absent.
w_10h
a vector of minimum and maximum fuel load [t/ha] for fuel class 10-hour; (0,0) if fuel class is absent.
w_100h
a vector of minimum and maximum fuel load [t/ha] for fuel class 100-hour; (0,0) if fuel class is absent.
w_Live_Herb
a vector of minimum and maximum fuel load [t/ha] for live herbs; (0,0) if fuel class is absent.
w_Live_Woody
a vector of minimum and maximum fuel load [t/ha] for live woody fuels; (0,0) if fuel class is absent.
s_1h
a vector of minimum and maximum surface-to-volume ratio [m2/m3] for fuel class 1-hour; (0,0) if fuel class is absent.
s_10h
a vector of minimum and maximum surface-to-volume ratio [m2/m3] for fuel class 10-hour; (0,0) if fuel class is absent.
s_100h
a vector of minimum and maximum surface-to-volume ratio [m2/m3] for fuel class 100-hour; (0,0) if fuel class is absent.
s_Live_Herb
a vector of minimum and maximum surface-to-volume ratio [m2/m3] for live herbs; (0,0) if fuel class is absent.
s_Live_Woody
a vector of minimum and maximum surface-to-volume ratio [m2/m3] for live woody fuels; (0,0) if fuel class is absent.
delta
a vector of minimum and maximum fuel bed depth [cm].
mx.dead
a vector of minimum and maximum dead fuel moisture of extinction [percent].
h_1h
a vector of minimum and maximum heat content [kJ/kg] for fuel class 1-hour; (0,0) if fuel class is absent.
h_10h
a vector of minimum and maximum heat content [kJ/kg] for fuel class 10-hour; (0,0) if fuel class is absent.
h_100h
a vector of minimum and maximum heat content [kJ/kg] for fuel class 100-hour; (0,0) if fuel class is absent.
h_Live_Herb
a vector of minimum and maximum heat content [kJ/kg] for live herbs; (0,0) if fuel class is absent.
h_Live_Woody
a vector of minimum and maximum heat content [kJ/kg] for live woody fuels; (0,0) if fuel class is absent.
modeltype
select S(tatic) or D(ynamic) Rothermel model.
obs
a vector of observed ROS [m/min].
m
a data frame of percent moisture on a dry weight basis [percent] for fuel classes 1-hour, 10-hour, 100-hour, live herbs and live woody, respectively (5 values or columns; 0 if fuel class is absent).
u
a vector of midflame windspeed [km/h].
slope
a vector of percent site slope [percent].
method
method of forecast error measure, see error. maxiter
the maximum number of iterations to run before the GA search is halted.
popSize
the population size.
pcrossover
the probability of crossover between pairs of chromosomes. Typically this is a large value and by default is set to 0.8.
pmutation
the probability of mutation in a parent chromosome. Usually mutation occurs with a small probability, and by default is set to 0.1.
elitism
the number of best fitness individuals to survive at each generation. By default the top 5/100 individuals will survive at each iteration.
...
additional arguments to be passed to ga.