ThreepSeriesModel(SoilR)
ThreepSeriesModel()所属R语言包：SoilR
Implementation of a three pool model with series structure
三池模型，该模型与串联结构的实现
译者：生物统计家园网 机器人LoveR
描述----------Description----------
This function creates a model for three pools connected in series. It is a wrapper for the more general function GeneralModel.
这个函数创建一个模型3个串联连接的。这是一个更普遍的功能GeneralModel的包装。
用法----------Usage----------
参数----------Arguments----------
参数：t
A vector containing the points in time where the solution is sought.
一种向量，包含的点的时间寻求溶液。
参数：ks
A vector of lenght 3 containing the values of the decomposition rates for pools 1, 2, and 3.
长度3的一种向量，含有池的1，2，和3中的分解速率的值。
参数：a21
A scalar with the value of the transfer rate from pool 1 to pool 2.
标量的传输速率从池1池2的值。
参数：a32
A scalar with the value of the transfer rate from pool 2 to pool 3.
标量的传输速率2至3池池的值。
参数：C0
A vector of length 3 containing the initial amount of carbon for the 3 pools.
长度为3的一种向量，含有碳的初始量为3个游泳池。
参数：In
A scalar or data.frame object specifying the amount of litter inputs by time.
一个标量或数据框对象，指定的凋落物输入量。
参数：xi
A scalar or data.frame object specifying the external (environmental and/or edaphic) effects on decomposition rates.
一个标量或数据框对象指定的外部（环境和/或土壤）的分解率的影响。
参数：solver
A function that solves the system of ODEs. This can be euler or ode or any other user provided function with the same interface.
的功能，解决了系统的微分方程。这可能是euler或ode或任何其他用户提供的函数具有相同的接口。
参数：pass
if TRUE Forces the constructor to create the model even if it is invalid
如果TRUE部队的构造函数来创建模型，即使它是无效的
值----------Value----------
A Model Object that can be further queried
一个模型对象，可以进一步查询
（作者）----------Author(s)----------
Carlos A. Sierra <csierra@bgc-jena.mpg.de>, Markus Mueller <mamueller@bgc-jena.mpg.de>
参见----------See Also----------
ThreepParallelModel, ThreepFeedbackModel
ThreepParallelModel，ThreepFeedbackModel
实例----------Examples----------
t_start=0
t_end=10
tn=50
timestep=(t_end-t_start)/tn
t=seq(t_start,t_end,timestep)
ks=c(k1=0.8,k2=0.4,k3=0.2)
C0=c(C10=100,C20=150, C30=50)
In = 50
Ex1=ThreepSeriesModel(t=t,ks=ks,a21=0.5,a32=0.2,C0=C0,In=In,xi=fT.Q10(15))
Ct=getC(Ex1)
Rt=getReleaseFlux(Ex1)
plot(t,rowSums(Ct),type="l",ylab="Carbon stocks (arbitrary units)",xlab="Time (arbitrary units)",lwd=2,ylim=c(0,sum(Ct[1,])))
lines(t,Ct[,1],col=2)
lines(t,Ct[,2],col=4)
lines(t,Ct[,3],col=3)
legend("topright",c("Total C","C in pool 1", "C in pool 2","C in pool 3"),lty=c(1,1,1,1),col=c(1,2,4,3),lwd=c(2,1,1,1),bty="n")
plot(t,rowSums(Rt),type="l",ylab="Carbon released (arbitrary units)",xlab="Time (arbitrary units)",lwd=2,ylim=c(0,sum(Rt[1,])))
lines(t,Rt[,1],col=2)
lines(t,Rt[,2],col=4)
lines(t,Rt[,3],col=3)
legend("topright",c("Total C release","C release from pool 1", "C release from pool 2","C release from pool 3"),lty=c(1,1,1,1),col=c(1,2,4,3),lwd=c(2,1,1,1),bty="n")
Inr=data.frame(t,Random.inputs=rnorm(length(t),50,10))
plot(Inr,type="l")
Ex2=ThreepSeriesModel(t=t,ks=ks,a21=0.5,a32=0.2,C0=C0,In=Inr)
Ctr=getC(Ex2)
Rtr=getReleaseFlux(Ex2)
plot(t,rowSums(Ctr),type="l",ylab="Carbon stocks (arbitrary units)",xlab="Time (arbitrary units)",lwd=2,ylim=c(0,sum(Ctr[1,])))
lines(t,Ctr[,1],col=2)
lines(t,Ctr[,2],col=4)
lines(t,Ctr[,3],col=3)
legend("topright",c("Total C","C in pool 1", "C in pool 2","C in pool 3"),lty=c(1,1,1,1),col=c(1,2,4,3),lwd=c(2,1,1,1),bty="n")
plot(t,rowSums(Rtr),type="l",ylab="Carbon released (arbitrary units)",xlab="Time (arbitrary units)",lwd=2,ylim=c(0,sum(Rtr[1,])))
lines(t,Rtr[,1],col=2)
lines(t,Rtr[,2],col=4)
lines(t,Rtr[,3],col=3)
legend("topright",c("Total C release","C release from pool 1", "C release from pool 2","C release from pool 3"),lty=c(1,1,1,1),col=c(1,2,4,3),lwd=c(2,1,1,1),bty="n")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-30 11:36:15
