Rによるデータ分析

Rによる有向比例分散分析

有向比例分散分析 のRによる実施例です。

比の歪度の行列の作成

行の項目が分子、列の項目が分母になる比についての歪度を、行列形式でまとめます。

GM1という行列は、例えば、歪度の行列の２行目１列と、１行目２列で、高い方の成分を１、低い方の成分を０として、有効グラフの作成用に作っています。

library(igraph)
library(e1071)
setwd("C:/Rtest")
Data <- read.csv("Data.csv", header=T)
n <- ncol(Data)
SKEW <- matrix(0,nrow=n,ncol=n)
rownames(SKEW)<-colnames(Data)
colnames(SKEW)<-colnames(Data)
GM1 <- SKEW
for (i in 1:n) {
for (j in 1:n) {
if( i > j){
SKEW[i,j] <- skewness(Data[,j]/Data[,i], type=2)
SKEW[j,i] <- skewness(Data[,i]/Data[,j], type=2)
if( abs(SKEW[i,j]) < abs(SKEW[j,i])){
GM1[i,j] <- 1
} else {
GM1[j,i] <- 1
}
}
}
}
SKEW # 尖度の行列


plot(graph.adjacency(GM1,weighted=T, mode = "directed"))

Rによるデータ分析