插值和空间分析（二）_变异函数分析（R语言）

方法1、散点图

hscat(log(zinc)~1, meuse, (0:9)*100)

方法2、变异函数云图

library(gstat)
cld <- variogram(log(zinc) ~ 1, meuse, cloud = TRUE)
svgm <- variogram(log(zinc) ~ 1, meuse)
d <- data.frame(gamma = c(cld$gamma, svgm$gamma),dist = c(cld$dist, svgm$dist),id = c(rep("cloud", nrow(cld)), rep("sample variogram", nrow(svgm))))
xyplot(gamma ~ dist | id, d,scales = list(y = list(relation = "free", #ylim = list(NULL, c(-.005,0.7)))),limits = list(NULL, c(-.005,0.7)))),layout = c(1, 2), as.table = TRUE,panel = function(x,y, ...) {if (panel.number() == 2)ltext(x+10, y, svgm$np, adj = c(0,0.5)) #$panel.xyplot(x,y,...)},xlim = c(0, 1590),cex = .5, pch = 3
)

plot(variogram(log(zinc) ~ 1, meuse)) // 对每一个距离去平均

sel <- plot(variogram(zinc ~ 1, meuse, cloud = TRUE), digitize = TRUE)
plot(sel, meuse)

v <- variogram(log(zinc) ~ 1, meuse)
print(xyplot(gamma ~ dist, v, pch = 3, type = 'b', lwd = 2, col = 'darkblue',panel = function(x, y, ...) {for (i in 1:100) {meuse$random = sample(meuse$zinc)v = variogram(log(random) ~ 1, meuse)llines(v$dist, v$gamma, col = 'grey')}panel.xyplot(x, y, ...)},ylim = c(0, 0.75), xlab = 'distance', ylab = 'semivariance'
))

转载于:https://www.cnblogs.com/takeaction/p/4123296.html

插值和空间分析（二）_变异函数分析（R语言）相关推荐

【嵌入式】Libmodbus源码分析(二)-常用接口函数分析
00. 目录文章目录 00. 目录 01. modbus常用接口函数概述 02. modbus辅助接口函数 03. modbus功能接口函数 04. modbus数据处理 05. 附录 01. mo ...
白话空间统计二十一：密度分析（五）带宽与核表面曲率的关系
白话空间统计二十一:密度分析(五) 上次讲密度分析的时候,有同学问道带宽的问题,实际上我翻了一下以前写的文章,在密度分析一.二里面,都对这个有过描述,详细的可以回去翻一下(很老的文章了--可能要翻到 ...
白话空间统计二十一：密度分析（一）
白话空间统计二十一:密度分析(一) 密度分析这个概念其实很早就想写了,也有无数同学都问过我,虾神你能不能讲讲那些漂亮的热度图是怎么做的啊?比如下面这种: 如果说,这是互联网地图里面,最让人喜闻乐见的一 ...
在r中弄方差分析表_医学统计与R语言： qvalue
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发 (FalseDiscoveryRate(FDR)=Expected(FalsePositive/(FalsePositive+TruePos ...
多元有序logistic回归_医学统计与R语言：多分类logistic回归HosmerLemeshow拟合优度检验...
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发输入1:multinominal logistic regression install.packages("nnet" ...
二元置信椭圆r语言_医学统计与R语言：圆形树状图（circular dendrogram）
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发输入1: "ggraph") 结果1: name 输入2: <- graph_from_data_frame(my ...
二元置信椭圆r语言_医学统计与R语言：多分类logistic回归HosmerLemeshow拟合优度检验...
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发输入1:multinominal logistic regression "nnet") 结果1: test (mult ...
r语言library什么意思_医学统计与R语言：百分条图与雷达图
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发百分条图-输入1: library(ggplot2) 结果1: year 输入2: percentbar <- gather(perc ...
语言nomogram校准曲线图_医学统计与R语言：Meta 回归作图（Meta regression Plot）
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发输入1: install.packages("metafor") library(metafor) dat.bcg 结果 ...
oracle数据泵导入导出_【软件】R语言数据导入与导出
"R语言导入文本和xlsx文件数据的方法,以及数据与图片的输出" 许多数据往往保存在TXT文件或Excel文件中,该如何将这些文件导入R语言进行分析呢?另外,使用R语言处理完数据之 ...

插值和空间分析（二）_变异函数分析（R语言）

插值和空间分析（二）_变异函数分析（R语言）相关推荐

最新文章

热门文章