「全局溢出」当一个区域的特征变化影响到所有区域的结果时,就会产生全局溢出效应.这甚至适用于区域本身,因为影响可以传递到邻居并返回到自己的区域(反馈).具体来说,全球溢出效应影响到邻居.邻居到邻居.邻居到邻居等等. 「局部溢出」是指影响只落在附近或近邻的情况,在它们影响邻邻区域之前就消失了. 对应全局与局部溢出,存在全局与局部自相关检验.全局自相关检验指标主要有 moran'I 指数.Geary 指数 C 统计量以及 Getis-Ord global G 统计量:局部自相检验指标主要有局部mora…

R data analysis examples 功效分析 power analysis for one-sample t-test单样本t检验 例1.一批电灯泡,标准寿命850小时,标准偏差50,40小时的差值是巨大的,此研究设定效应值d= (850-810)/50,希望有90%的可能检测到,即功效值为0.9,还希望有95%的把握不误报显著差异, 问需要多少支电灯泡. H0=850,HA=810 library('pwr') pwr.t.test(d=(850-810)/50,power=0.…

正态检验与R语言 1.Kolmogorov–Smirnov test 统计学里, Kolmogorov–Smirnov 检验(亦称:K–S 检验)是用来检验数据是否符合某种分布的一种非参数检验,通过比较一个频率分布f(x)与理论分布g(x)或者两个观测值分布来判断是否符合检验假设.其原假设H0:两个数据分布一致或者数据符合理论分布.拒绝域构造为:D=max| f(x)- g(x)|,当实际观测值D>D(n,α)则拒绝H0,否则则接受H0假设.由于KS检验不需要知道数据的分布情况,在小样本的统计分…

twoway 系列命令绘制地图 ❝ 在谈到用stata绘制地图时,很多人首先想到的是spmap命令.其实,最常见的twoway系列命令就可以完成我们的大多数绘图需求. 对于spmap命令的使用,在线资源已经有很多了,这主要介绍使用最常用的twoway系列命令如何进行地图可视化.spmap命令的使用只做简单介绍.(后台回复「20200412」获取相关数据) ❞ 首先,使用shp2dta命令将 shp 文件转为 stata 能够识别的 .dta 文件.shp2dta语法格式如下: shp2dta u…

R语言进行DW检验: library(lmtest) dw = dwtest(fm1) > dw Durbin-Watson test data: fm1 DW = 2.4994, p-value = 0.8706 DW检验的原假设为:误差不相关! 因为dw>0.05所以不拒绝原假设,即认为误差是不相关的. 误差自相关会产生的后果: 1．参数估计量仍然是线性的.无偏的,但非有效. 2．OLS估计量的被估方差是有偏的且会被低估,因而会使相应的t值变大. 3．模型的t和F统计检验失效. 4．用通常…

前言 很多朋友说在R里没法使用高德地图,这里给出一个基于leaflet包的解决方法. library(leaflet) # 添加高德地图 m <- leaflet() %>% addTiles( 'http://webrd02.is.autonavi.com/appmaptile?lang=zh_cn&size=1&scale=1&style=8&x={x}&y={y}&z={z}', tileOptions(tileSize=256, minZ…

做线性回归的时候,检验回归方程和各变量对因变量的解释参数很容易搞混乱,下面对这些参数进行一下说明: 1.t检验:t检验是对单个变量系数的显著性检验   一般看p值:    如果p值小于0.05表示该自变量对因变量解释性很强. 2.F检验:F检验是对整体回归方程显著性的检验,即所有变量对被解释变量的显著性检验 3.P值:P值就是t检验用于检测效果的一个衡量度,t检验值大于或者p值小于0.05就说明该变量前面的系数显著,选的这个变量是有效的. 4.R方:拟合优度检验 5.调整后的R方: 小结: t检…

A IMA模型是一种著名的时间序列预测方法,主要是指将非平稳时间序列转化为平稳时间序列,然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型.ARIMA模型根据原序列是否平稳以及回归中所含部分的不同,包括移动平均过程(MA).自回归过程(AR).自回归移动平均过程(ARMA)以及ARIMA过程.其中ARIMA(p,d,q)称为差分自回归移动平均模型,AR是自回归, p为自回归项: MA为移动平均,q为移动平均项数,d为时间序列成为平稳时所做的差分次数. 通常的建立ARIMA…

基于空间自相关,R语言克里金插值 library(gstat) Warning message: In scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : EOF within quoted string library(sp) data(meuse) head(meuse) x y cadmium copper lead zinc elev dist om ffreq soil lime landuse 1…

本章開始学习<Forecasting: principles and practice> 1 getting started 1.1 事件的可预言性 一个时间能不能被预言主要取决于以下三点 1. 对事件的影响因素的了解程度.比方彩票号码.没有内在的影响因素不能被预測 2. 可用数据量的多少,数据量太少没法预測 3. 预測结果本身的影响,比方预測汇率,可能大家知道预測的会长,那么人们就会採取对应的措施使预測结果不准. 1.2 经常使用预測模型 解释性模型,如其模型内包括其影响因素,通过影响因素来…