|||
共有7个变量,试求变量之间的灰色关 联系数。
## 环境设置
rm(list=ls());setwd(“E:/My Documents/R/data/”);options(digits=8)
##装载数据
mydata <- read.csv(“greysystem2.csv”)
##原始数据变换–均值化变换
mydata.mean <- sapply(mydata[-1],mean) #各指标均值
mydata.transform <- mydata[-1]/matrix(rep(mydata.mean,each=5),nrow=5,ncol=8) #均值变换
##求差序 列
mydata.difference <- abs(mydata.transform[,-1] – matrix(rep(mydata.transform[,1],7),nrow=5,ncol=7))
mydata.difference.max <- max(mydata.difference) #最大值
mydata.difference.min <- min(mydata.difference) #最小值
## 求灰关联度系数
mydata.R <- lapply(mydata.difference,function(x)(mydata.difference.min+0.5*mydata.difference.max)/(x+0.5*mydata.difference.max))
sort(sapply(mydata.R,mean))
注:0.5为分辨系数
参考文献:
夏军, Excel 2000在灰关联分析中的运用, 中国卫生统计, 2004, 21(2):117-121
唐启义,冯明光. DPS数据处理系统, 第39章灰色系统分析
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-5-1 22:10
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社