-
1 votesanswersviews
使用gstat进行局部变异函数的局部块克里格法
我一直无法使用R中的gstat包找到本地变异函数的局部块克里金特有关的信息 . 澳大利亚精确农业中心有一个名为VESPER的免费软件可以做到这一点，并且从我读过的内容应该可以在R中，我可以使用一些帮助将一个for循环放在一起，使gstat函数在本地工作 . 使用meuse数据集作为示例，我已经能够计算并将全局变异函数拟合到数据集： library(gstat) data(meuse) ... -
1 votesanswersviews
来自克里金数据的ggmap上的热图
我想在ggmap上绘制热图 . library(ggmap) turku<-get_map('turku', zoom=13) turkumap<-ggmap(turku, extent="device", legend="topleft") turkumap turkumap+geom_density2d(mapping=aes... -
1 votesanswersviews
如何在R中使用gstat拟合半变异函数模型？
我有一个csv文件文件包含3月1日下午12点的大气PM10浓度数据 . please, download.我想在R中使用gstat包绘制一个半变异函数 . 我试着在R中编写这些代码 . 但是对于这些数据，我不能适应模型 . library(sp) library(gstat) seoul3112<-read.csv("seoul3112.csv") ... -
0 votesanswersviews
R：差异变异函数拟合，克里金结果差
我试着在雅加达湾做一个克里金 . 我有一套测量点，有适当的坐标和属性（pH，盐度，......） 为了做克里金法，我首先需要为变差函数找到一个模型 . 当我使用"variogram"函数输出不完美但是应该没问题，但是当我尝试拟合变异函数时，我得到一个警告信息：在fit.variogram中（ph.vgm，model = vgm（0.12，"Sph"，0.1）... -
6 votesanswersviews
如何反向转换正态分数转换后的数据
我每天有61个测量站的降雨量，为期12年（8000平方公里） . 目标是创造5公里和25公里分辨率网格日降雨量产品 . 由于车站的数量很少，并且即使在雨季也不是所有车站都有雨，我选择使用气候变异图 . 典型日（irain）的雨量计测量如下，具有由NA表示的少量缺失值 . 7.8 4.4 15.4 19.1 5.8 0 42 6.4 21 21 0 0 ... -
1 votesanswersviews
GeoRglm R包中二项式数据的回归克里金法
我正在使用R包 geoRglm 的 binom.krige() 函数来确定具有多个连续和离散协变量的二元（0,1）响应变量的空间预测 . 使用带有二项式logit链接函数的 glm() ，我发现响应变量显示出对几个协变量的显着依赖性 . 我使用 krige.glm.control() 将趋势包含到 binom.krige() 中，其中我将两个趋势模型指定为 > trend.d=trend.s...