代码如下：
library("diagram")
#a <- read.table(file="clipboard",header=TRUE)
write.table(a,file="a.txt",sep=";")
aa <- a
c <- matrix(nrow=ncol(aa),ncol=ncol(aa),0)
i <- 1
while(i<=ncol(aa)){
j <- 1
while(j<=ncol(aa)){
if(2>3&&aa[i,j]<10){
aa[i,j] <- 0
}
if(aa[i,j]>0 && aa[j,i]>0){
c[i,j] <- c[j,i] <- 0.06
}
j <- j+1
}
i <- i+1
}
c[1,10] <- 0.45
c[6,1] <- 0.45
c[15,1] <- 0.7
c[1,15] <- 0.1
c[1,6] <- -0.65
c[1,7] <- -0.95
c[7,1] <- 0.45
c[7,5] <- 0.75
c[5,7] <- -0.95
c[11,5] <- 0.65
c[12,5] <- -0.85
c[12,10] <- -0.38
c[10,12] <- 0.45
c[10,7] <- -1.25
c[10,15] <- 0.45
c[14,12] <- 0.35
c[14,11] <- 0.4
c[11,10] <- -0.45
c[5,14] <- 0.35
c[5,6] <- 0.8
c[3,5] <- 0.8
c[5,2] <- -0.65
col <- matrix(nrow=ncol(aa),ncol=1,"black")
col[1]="red"
col[5]="darkorchid"
col[10]="darkorange"
col2 <- matrix(nrow=ncol(aa),ncol=ncol(aa),"black")
col2[,1]="red"
col2[,5]="darkorchid"
col2[,10]="darkorange"
plotmat(aa, pos = c(1,4,5,4,1), curve = c,
box.lcol = col,arr.lcol=col2,
name = colnames(aa), lwd = 1.5, box.lwd = 2,
cex.txt = 0.8, box.cex = 0.7, box.size = 0.05,
arr.length = 0.5, box.type = "circle", box.prop = 1,
shadow.size = 0.01, self.cex = 0.6, my = -0.075, mx = -0.01,
relsize = 0.8, self.shiftx = c(0, 0, 0.125, -0.12, 0.125, 0),
self.shifty = 0, main = "上海一周订单状态流程图")
