Displaying 1 result from an estimated 1 matches for "potyield3".
2002 Feb 11
0
profile
...ta1 <- 0.3791
eta2 <- 0.6332
PopStd <- 90468
beta <- Beta
DIs <- Dls
MnmN <- Nmin
OptN <- Nopt1
# define model function
Y.model <- function(gN, MnmN, OptN, DIs, beta, eta1, eta2,
Popn, Dmax, AWC, SumEp, PotYield3, Nsupply)
{
Ymax<- 1-ifelse(Popn<=PopStd, eta1, eta2)*log(Popn/PopStd)
Ymax <- Ymax*PotYield3*Popn/1000
Ymax <- Ymax*ifelse(Dmax<=DIs*AWC, 1, 1 - beta*(Dmax
-DIs*AWC)/SumEp)
Nstar <- (Nsupply- MnmN*Ymax)...