EQXYresid(Rquake)
EQXYresid()所属R语言包：Rquake
Calculate Residuals
计算残差
译者：生物统计家园网 机器人LoveR
描述----------Description----------
given an earthquake hypocenter and a list of station information, retrieve the station residuals.
由于地震震源和电台信息的列表，检索站残差。
用法----------Usage----------
EQXYresid(XY, vel = list(), h1 = c(0, 0, 0, 0), PLOT = FALSE)
参数----------Arguments----------
参数：XY
matrix of station location and arrival times.
矩阵站的位置和到达时间。
参数：vel
list, RSEIS velocity model
列表，RSEIS速度模型
参数：h1
hypocenter location, c(x,y,z,t)
震源位置（X，Y，Z，T，C）
参数：PLOT
logical, TRUE=plot the residuals
逻辑，TRUE =图的残差
Details
详细信息----------Details----------
The XY mtrix is in cartesian coordinates, i.e. it has been projected into units of km. Only 1D velocity models are used at this time. Only residuals of P and S wave arrivals are estimated.
XY mtrix是在直角坐标系中，也就是说，它已被投影到单位为公里。在这个时候只有一维的速度模型。只有残差的P波和S波到达的估计。
值----------Value----------
vector, right hand side of the least squares problem.
矢量，右手侧的最小二乘问题。
（作者）----------Author(s)----------
Jonathan M. Lees<jonathan.lees@unc.edu>
参见----------See Also----------
travel.time1D,UPdateEQLOC
travel.time1D，UPdateEQLOC
实例----------Examples----------
## Not run: [＃不运行：]
XY = GLOB.XY(pstas$lat, pstas$lon, proj)
elcor = rep(0, length(pstas$lat))
DZ = pstas$z - mean(stas$z)
elcor[pstas$phase=="P"] = DZ[pstas$phase=="P"]/v$vp[1]
elcor[pstas$phase=="S"] = DZ[pstas$phase=="S"]/v$vs[1]
XY$cor = elcor
XY$phase = pstas$phase
XY$sec = pstas$sec
eqXY = GLOB.XY(sol[1], sol[2], proj)
res = EQXYresid(XY, vel=v , h1=c(eqXY$x, eqXY$y, sol[2],sol[4]) , PLOT=FALSE)
## End(Not run)[＃（不执行）]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-28 20:26:23
