etasap(SAPP)
etasap()所属R语言包：SAPP
Maximum Likelihood Estimates of The ETAS Model
最大似然估计的ETAS模型
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Compute the maximum likelihood estimates of five parameters of ETAS model. This function consists of two (exact and approximated) versions of the calculation algorithm for the maximization of likelihood.
ETAS模型参数计算的最大似然估计。此功能包括两个版本（精确和近似）的计算算法的最大化的可能性。
用法----------Usage----------
etasap(time, mag, threshold=0.0, reference=0.0, parami,
zts=0.0, tstart, zte, approx=2, tmpfile=NULL, nlmax=1000, plot=TRUE)
参数----------Arguments----------
参数：time
the time measured from the main shock(t=0).
测定的时间，从主震（t = 0时）。
参数：mag
magnitude.
级。
参数：threshold
threshold magnitude.
阈值的大小。
参数：reference
reference magnitude.
基准幅度。
参数：parami
initial estimates of five parameters mu, K, c, alpha and p.
五个参数的初步估计mu，K，c，alpha和p。
参数：zts
the start of the precursory period.
开始前兆期间。
参数：tstart
the start of the target period.
的对象期间的开始。
参数：zte
the end of the target period.
的对象期间的结束。
参数：approx
>0 : the level for approximation version, which is one of the five levels 1, 2, 4, 8 and 16. The higher level means faster processing but lower accuracy. <br> =0 : the exact version.
> 0：近似版本，这是一个五个级别1，2，4，8和16的电平。水平越高，意味着更快的处理速度，但精度较低。参考= 0：确切的版本。
参数：tmpfile
write the process of maximum likelihood procedure to tmpfile.<br> If "" print the process to the standard output and if NULL (default) no report.
写的过程中，最大似然法tmpfile。<BR>如果“打印输出到标准输出的过程中，如果NULL（默认）没有报告。
参数：nlmax
the maximum number of steps in the process of minimizing.
最小化的过程中的步骤的最大数目。
参数：plot
logical. If TRUE (default) the graph of cumulative number and magnitude of earthquakes against the ordinary time is plotted.
逻辑。如果是TRUE（默认值）的累积数量和规模的地震对普通的图形绘制。
Details
详细信息----------Details----------
The ETAS model is a point-process model representing the activity of earthquakes of magnitude M_z and larger occurring in a certain region during a certain interval of time. The total number of such earthquakes is denoted by N. The seismic activity includes primary activity of constant occurrence rate mu in time (Poisson process). Each earthquake ( including aftershock of another earthquake) is followed by its aftershock activity, though only aftershocks of magnitude M_z and larger are included in the data. The aftershock activity is represented by the Omori-Usu formula in the time domain. The rate of aftershock occurrence at time t following the ith earthquake (time: t_i, magnitude: M_i) is given by
ETAS模型是一个点过程模型级M_z和发生在某一个区域在一定的时间间隔的地震活动。这样的地震总数的N表示。地震活动主要活动不断发生率mu的时间（泊松过程）。每次地震及其余震活动（包括另一次地震的余震），其次是，虽然只是余震的震级M_z和更大的包含在数据中。余震活动表示在时域由大森乌苏式的。余震发生率在时间ti日地震（时间：t_i，大小：M_i）
for t>t_i where K, alpha, c, and p are constants, which are common to all aftershock sequences in the region. The rate of occurrence of the whole earthquake series at time t becomes
为 t>t_i 其中K，alpha，c和p是常量，它们共同在该区域的所有余震序列。整个地震系列的发生率在时间t成为
The summation is done for all i satisfying t_i < t. Five parameters mu, K, c, alpha and p represent characteristics of seismic activity of the region.
的总和是为所有的i满意t_i < t的。五参数mu，K，c，alpha和p代表该区域的地震活动特征。
值----------Value----------
参数：ngmle
negative max log-likelihood.
负最大似然。
参数：param
list of maximum likelihood estimates of five parameters mu, K, c, alpha and p.
列表中的最大似然估计的五个参数mu，K，c，alpha和p。
参数：aic2
AIC/2.
AIC / 2。
参考文献----------References----------
Computer Science Monographs, No.33, Statistical Analysis of Seismicity - updated version (SASeies2006). The Institute of Statistical Mathematics.
实例----------Examples----------
data(main2003JUL26) # The aftershock data of 26th July 2003 earthquake of M6.2 [2003年7月26日的6.2级大地震的余震数据]
x <- main2003JUL26
etasap(x$time, x$magnitude, 2.5, 6.2,
c(0, 0.63348E+02, 0.38209E-01, 0.26423E+01, 0.10169E+01),, 0.01, 18.68)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-29 22:15:37
