- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 2069 个
- 通用积分
- 0.0275
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 116 点
- 帖子
- 3
- 精华
- 0
- 在线时间
- 115 小时
- 注册时间
- 2018-1-2
- 最后登录
- 2022-7-15
大专生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
- 2069 个
- 通用积分
- 0.0275
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 116 点
- 帖子
- 3
- 精华
- 0
- 在线时间
- 115 小时
- 注册时间
- 2018-1-2
- 最后登录
- 2022-7-15
| 悲催 2021-4-6 19:05:21 |
---|
签到天数: 1 天 连续签到: 1 天 [LV.1]初来乍到
|
50论坛币
有大佬可以帮我检查一下是为什么吗
自己根据公式写的代码计算莫兰指数
- >w <- matrix(c(0,1,1,1,1,0,1,0,1,1,0,1,1,0,1,0),nrow=4,ncol=4,dimnames=list(c("千","昭","天","名"),c("千","昭","天","名")))
- >a <- c(80,95,67,78)
- >mean(a)
- >b <- a-mean(a)
- >i <- (length(a)/sum(w))*(t(b)%*%w%*%b)/(t(b)%*%b)
- >i
复制代码 [,1]
[1,] -0.3396985
用R语言ape包 计算莫兰指数
- >install.packages("ape")
- >library(ape)
- >Moran.I(a,w)
复制代码 $observed
[1] -0.3538526
$expected
[1] -0.3333333
$sd
[1] 0.1521452
$p.value
[1] 0.8927175
顺便问一下,有大佬自己写了moran i 的计算和检定公式吗?谢谢! |
|