inst/tests/automaticR/BoundFcExamples.R

In SoilR: Models of Soil Organic Matter Decomposition

# A function that is valid in a time interval bounded by tstart and tend 
  

 x <- BoundFc(
	map = function(t){
				c(
					1+sin(t),
					2+sin(2*t)
				)
			},
	starttime = 1,
	endtime = 10,
	format = "Delta14C"
)
# A function that is valid in a time interval bounded by tstart and tend 
  

 x <- BoundFc(
	map = function(t){
				c(
					1+sin(t),
					2+sin(2*t)
				)
			},
	starttime = 1,
	endtime = 10,
	format = "AbsoluteFractionModern"
)
# a vector of times, a vector of scalar fractions per time step and a vector lag 
  

 x <- BoundFc(
	times = 0:99,
	data = C14Atm_NH[1:100,2],
	lag = c(1.1,1.2),
	format = "Delta14C"
)
# a vector of times, a vector of scalar fractions per time step and a vector lag 
  

 x <- BoundFc(
	times = 0:99,
	data = C14Atm_NH[1:100,2],
	lag = c(1.1,1.2),
	format = "AbsoluteFractionModern"
)
# a vector of times,a vector of scalar fractions per time step and a scalar lag 
  

 x <- BoundFc(
	times = 0:99,
	data = C14Atm_NH[1:100,2],
	lag = 1.1,
	format = "Delta14C"
)
# a vector of times,a vector of scalar fractions per time step and a scalar lag 
  

 x <- BoundFc(
	times = 0:99,
	data = C14Atm_NH[1:100,2],
	lag = 1.1,
	format = "AbsoluteFractionModern"
)
# We could also imagine time series data
# stored in an array consisting of
# many stacked vectors, one for each time step.
# and combine both to a list. 
 
times <- seq(1,10,by=0.1)
a <- array(dim=c(2,length(times)))
a[1,] <- -0.1*(sin(times)+1.1)
a[2,] <- -0.2*(sin(times)+1.2) 

 x <- BoundFc(
	map = list(times=times,data=a,format="Delta14C")
)
# We could also imagine time series data
# stored in an array consisting of
# many stacked vectors, one for each time step.
# and combine both to a list. 
 
times <- seq(1,10,by=0.1)
a <- array(dim=c(2,length(times)))
a[1,] <- -0.1*(sin(times)+1.1)
a[2,] <- -0.2*(sin(times)+1.2) 

 x <- BoundFc(
	map = list(times=times,data=a,format="AbsoluteFractionModern")
)