horton: Horton infiltration parameter optmisation in R

Description Usage Arguments Value Author(s) References Examples

Description

This function optimises Horton (1940) infiltration parameters: fo,fc, and k. It also predicts infiltration.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
horton(data = NULL, time, I, fc = 0.1, fo = 0.1, k = 0.1,
  group = NULL)

## Default S3 method:
horton(data = NULL, time, I, fc = 0.1, fo = 0.1,
  k = 0.1, group = NULL)

## S3 method for class 'horton'
predict(object, time = NULL, ...)

## S3 method for class 'horton'
plot(x, xlab = "Time(Minutes)", ylab = "Cumulative (cm)",
  main = NULL, layout = NULL, ...)

## S3 method for class 'horton'
summary(object)

## S3 method for class 'horton'
print(x, ...)

## S3 method for class 'horton'
coef(object, ...)

Arguments

data

dataframe. It can contain data with column names of "time" and "I"

time

character or numeric. The name of time variable in the dataframe. If the "data" parameter contains "time", this will be ignored. The unit must be in seconds.

I

cumulative infiltration (I) or infiltration rate (i) [mm]

fc

final constant rate of infiltration at saturation

fo

parameter for initial infiltration capacity

k

a constant depending primarily upon soil and vegetation

group

character. The name of the group variables if the data is from different areas.

object

Model output object

...

Any other graphical parameter

x

a return object of the function.

xlab

x label of the plot

ylab

y label of the plot

main

Title of the plot

layout

plot layout

Value

Author(s)

George Owusu

References

Horton, R. E. (1940). An approach toward a physical interpretation of infiltration capacity Proc. Sci. Soc. Amer., 5, 399-417.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
data=read.csv(system.file("ext","sys","exampleBEST.csv",package="vadose"))
horton1<-horton(data=data,time="time",I="I")
#print(gof(horton1))
plot(horton1)
predict(horton1)
coef(horton1)

#group simulation
data=read.csv(system.file("ext","sys","infiltration2.csv",package="vadose"))
horton1g<-horton(data=data,time="minutes",I="CumInfil",group="ID")
coef(horton1g)
#generic function ######################################

gowusu/vadose documentation built on May 17, 2019, 7:59 a.m.