andersson: R optmisation of Andersson four-parameter particle size...

Description Usage Arguments Value Author(s) References See Also Examples

Description

The function optimises 4 parameters of Andersson (1990) particle size distribution function

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
andersson(data = NULL, D = NULL, fr = NULL, p = NULL, group = NULL,
  a = 1, b = 1, df1 = 1, do = 1)

## Default S3 method:
andersson(data = NULL, D = NULL, fr = NULL, p = NULL,
  group = NULL, a = 1, b = 1, df1 = 1, do = 1)

## S3 method for class 'andersson'
plot(x, main = "Particle Size Distribution Function",
  xlab = "Particle Diameter", ylab = "Fraction", ...)

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

## S3 method for class 'andersson'
predict(object, D, ...)

Arguments

data

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

D

character. The corresponding diameter in the data.

fr

character. The corresponding fraction of the particles.

p

porosity[g/m^3]

group

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

a

an initial value of 'a' parameter to be optimised

b

an initial value of 'b' parameter to be optimised

df1

an initial value of 'df1' parameter to be optimised

do

an initial value of 'do' parameter to be optimised

x

a return object of the function.

main

Title of the plot

xlab

x label of the plot

ylab

y label of the plot

...

Any other graphical parameter

object

Model output object

Value

Author(s)

George Owusu

References

Andersson, S. (1990). Markfysikaliska undersokningar i odlad jord XXVI. Om mineraljordens och mullens rumsutfyllande egenskaper. En teoretisk studie. (In Swedish) (pp. 7). Uppsala: Swedish University of agricultural sciences.

See Also

zhuang3, zhuang4,lass3,logarithmic, logistic2, logistic3,fredlund3,fredlund4, jaky,andersson,gompertz, haverkamp

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
data=read.csv(system.file("ext","sys","soil2.csv",package="vadose"))
single<- subset(data, ID=="30B20_1")
mod=andersson (data=single,p="sand",D="D",fr="Sand.")
plot(mod)
## Not run: 
#group simulation
mod2=andersson (data=data,D="D",fr="FractionSand",group="ID")
mod2$coef

## End(Not run)
#generic function

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