This is a five- parameter function to determine water retention hydraulic parameters contained in the Fredlund-Xing water retention model
SSfredlund(x, thr, ths, alp, nscal, mscal)A list of objects returned by standard R regression functions such as nls and lm
Suction potential/head as contained in the x-column of the xy water retention table or data
This is the residual moisture content. It's the moisture content when suction potential is very high (almost at the drying point)
This is the saturated moisture content. It's the moisture content when suction potential is very low (almost at the saturation point)
It's the inverse of air-entry potential or bubbling pressure
This is a parameter or index for the pore-size distribution
This is a scaling parameter which is also related to the index for the pore-size distribution
Christian Thine Omuto
A Fredlund-Xing model is a self-starting type of nonlinear curve fitting model for fitting water retention characteristics using experimental data. If warnings are given during the estimation process, try to increase the iteration limit using control function. Sometimes the warnings can be suppressed using (warn=-1)/ignored since they do not terminate the process
Fredlund DG and Xing A. 1994. Equations for the soil-water characteristic curve. Canadian Geotechnical Journal 31: 521-532
SSomuto, SSvgm, SSkosugi
## Attach sample data (isric data) from the package
data(isric)
pf <- subset(isric, Sample=="Benin2")
library(minpack.lm)
suppressWarnings(
fredlund.ns <- nlsLM(y ~ SSfredlund(x,thr,ths,alp,nscal,mscal), data=pf,
control = nls.lm.control(maxiter=200)))
coef(fredlund.ns)
plot(fitted(fredlund.ns)~pf$y)
Run the code above in your browser using DataLab