lm.morantest.sad(spdep)
lm.morantest.sad()所属R语言包:spdep
Saddlepoint approximation of global Moran's I test
鞍点逼近的全局Moran的I测试
译者:生物统计家园网 机器人LoveR
描述----------Description----------
The function implements Tiefelsdorf's application of the Saddlepoint approximation to global Moran's I's reference distribution.
该的功能实现Tiefelsdorf的应用鞍点逼近全局Moran I的参考分布。
用法----------Usage----------
lm.morantest.sad(model, listw, zero.policy=NULL, alternative="greater",
spChk=NULL, resfun=weighted.residuals, tol=.Machine$double.eps^0.5,
maxiter=1000, tol.bounds=0.0001, zero.tol = 1e-07, Omega=NULL,
save.M=NULL, save.U=NULL)
## S3 method for class 'moransad'
print(x, ...)
## S3 method for class 'moransad'
summary(object, ...)
## S3 method for class 'summary.moransad'
print(x, ...)
参数----------Arguments----------
参数:model
an object of class lm returned by lm; weights may be specified in the lm fit, but offsets should not be used
类的一个对象lm返回lm;权重可以指定lm适合,但偏移量不应该使用
参数:listw
a listw object created for example by nb2listw
例如创建一个listw对象的nb2listw
参数:zero.policy
default NULL, use global option value; if TRUE assign zero to the lagged value of zones without neighbours, if FALSE assign NA
默认为空,请使用全局选项的值,如果真没有邻居的滞后值的区域分配了零,如果为FALSE分配NA
参数:alternative
a character string specifying the alternative hypothesis, must be one of greater (default), less or two.sided.
一个字符串,指定其他假设,必须有一个更大的(默认),少或two.sided的。
参数:spChk
should the data vector names be checked against the spatial objects for identity integrity, TRUE, or FALSE, default NULL to use get.spChkOption()
应的数据向量空间对象的名称进行核对身份完整性,TRUE,否则返回FALSE,默认为空,使用get.spChkOption()
参数:resfun
default: weighted.residuals; the function to be used to extract residuals from the lm object, may be residuals, weighted.residuals, rstandard, or rstudent
默认:weighted.residuals要使用的函数提取残差lm对象,可能是residuals,weighted.residuals,rstandard或rstudent
参数:tol
the desired accuracy (convergence tolerance) for uniroot
所需的精度(收敛容差)为uniroot
参数:maxiter
the maximum number of iterations for uniroot
最大迭代次数为uniroot
参数:tol.bounds
offset from bounds for uniroot
偏移从uniroot界
参数:zero.tol
tolerance used to find eigenvalues close to absolute zero
公差找到特征值接近绝对零度
参数:Omega
A SAR process matrix may be passed in to test an alternative hypothesis, for example Omega <- invIrW(listw, rho=0.1); Omega <- tcrossprod(Omega), chol() is taken internally
特别行政区过程矩阵可以传递到测试的假说,例如Omega <- invIrW(listw, rho=0.1); Omega <- tcrossprod(Omega),chol()内服
参数:save.M
return the full M matrix for use in spdep:::exactMoranAlt
返回完整的M矩阵用于在spdep:::exactMoranAlt
参数:save.U
return the full U matrix for use in spdep:::exactMoranAlt
返回完整的铀材料使用中spdep:::exactMoranAlt
参数:x
object to be printed
反对要打印
参数:object
object to be summarised
反对总结
参数:...
arguments to be passed through
通过参数
Details
详细信息----------Details----------
The function involves finding the eigenvalues of an n by n matrix, and numerically finding the root for the Saddlepoint approximation, and should therefore only be used with care when n is large.
该功能可找到一个n×n矩阵的特征值和数值模拟发现的鞍点逼近的根,因此只能当n足够大时,小心使用。
值----------Value----------
A list of class moransad with the following components:
类moransad以下组件的列表:
参数:statistic
the value of the saddlepoint approximation of the standard deviate of global Moran's I.
的鞍点近似的标准偏离值的全局Moran的一
参数:p.value
the p-value of the test.
的p值的测试。
参数:estimate
the value of the observed global Moran's I.
的全局Moran所观察到的一值
参数:alternative
a character string describing the alternative hypothesis.
一个字符串,描述了另一种假设。
参数:method
a character string giving the method used.
一个字符串提供的方法使用。
参数:data.name
a character string giving the name(s) of the data.
给予(s)的数据的名称的字符串。
参数:internal1
Saddlepoint omega, r and u
鞍欧米茄,r和u
参数:internal2
f.root, iter and estim.prec from uniroot
f.root,ITER和estim.prec,从uniroot
参数:df
degrees of freedom
自由度
参数:tau
eigenvalues (excluding zero values)
特征值(不包括零值)
(作者)----------Author(s)----------
Roger Bivand <a href="mailto:Roger.Bivand@nhh.no">Roger.Bivand@nhh.no</a>
参考文献----------References----------
and local Moran's Ii reference distributions and their numerical evaluation.
参见----------See Also----------
lm.morantest
lm.morantest
实例----------Examples----------
eire <- readShapePoly(system.file("etc/shapes/eire.shp", package="spdep")[1],
ID="names", proj4string=CRS("+proj=utm +zone=30 +units=km"))
eire.nb <- poly2nb(eire)
#data(eire)[数据(爱尔兰)]
e.lm <- lm(OWNCONS ~ ROADACC, data=eire)
lm.morantest(e.lm, nb2listw(eire.nb))
lm.morantest.sad(e.lm, nb2listw(eire.nb))
summary(lm.morantest.sad(e.lm, nb2listw(eire.nb)))
e.wlm <- lm(OWNCONS ~ ROADACC, data=eire, weights=RETSALE)
lm.morantest(e.wlm, nb2listw(eire.nb), resfun=rstudent)
lm.morantest.sad(e.wlm, nb2listw(eire.nb), resfun=rstudent)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|