1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
args = commandArgs(trailingOnly=T)
N1 = as.numeric(args[1])
N2 = as.numeric(args[2])
N3 = as.numeric(args[3])
infile = args[4]
outfile = args[5]
svg(file=outfile,height=8,width=8)
M = as.matrix(read.table(infile))
L = dim(M)[1]
library(stats)
points = cmdscale(M)
args = commandArgs(trailingOnly=T)
N1 = as.numeric(args[1])
N2 = as.numeric(args[2])
N3 = as.numeric(args[3])
print(c(N1,N2,N3))
#library(MASS)
#mds = isoMDS(M)
#points = mds$points
minx = min(points[,1])
maxx = max(points[,1])
miny = min(points[,2])
maxy = max(points[,2])
xw = (maxx-minx)
yw = (maxy-miny)
plot(points,type="n",ylab="",xlab="",axes=F)
colorseq1 = hsv(seq(from=0.83333, to=1,length=N1),
seq(from=0.2,to=1, length=N1),
1)
colorseq2 = hsv(seq(from=0.5,to=0.666666,length=N2),
seq(from=0.2,to=1, length=N2),
1)
colorseq3 = hsv(seq(from=0.166666,to=0.333333,length=N3),
seq(from=0.2,to=1, length=N3),
1)
lines(points[1:N1,], pch=20,col=hsv(1, 0.05, 0.9))
lines(points[(N1+1):(N1+N2),], pch=20,col=hsv(0.66666, 0.05, 0.85))
lines(points[(N1+N2+1):(N1+N2+N3),],pch=20,col=hsv(0.33333, 0.05, 0.85))
points(points,pch=20,col=c(colorseq1,colorseq2,colorseq3))
xwidth = (maxx-minx)*0.025
ywidth = (maxy-miny)*0.025
nbox = 12
i = (0:(nbox-1))*xwidth
start.x = maxx-nbox*xwidth*1.1
start.y = miny+2.1*ywidth
height = ywidth
rect(start.x+i,start.y,start.x+xwidth+i,start.y+height, col=colorseq1[floor(seq(from=0,to=N1,length=nbox+1))])
start.x = start.x
start.y = start.y - ywidth
height = ywidth
rect(start.x+i,start.y,start.x+xwidth+i,start.y+height, col=colorseq2[floor(seq(from=0,to=N2,length=nbox+1))])
start.x = start.x
start.y = start.y - ywidth
height = ywidth
rect(start.x+i,start.y,start.x+xwidth+i,start.y+height, col=colorseq3[floor(seq(from=0,to=N3,length=nbox+1))])
# dev.copy2eps(file="NNI-SPR.eps")
text(points[1,1],points[1,2]-yw*0.001,pos=1,"S1")
text(points[N1,1],points[N1,2]-yw*0.001,pos=1,"E1")
text(points[N1+1,1],points[N1+1,2]-yw*0.001,pos=1,"S2")
text(points[N1+N2,1],points[N1+N2,2]-yw*0.001,pos=1,"E2")
text(points[N1+N2+1,1],points[N1+N2+1,2]-yw*0.001,pos=1,"S3")
text(points[N1+N2+N3,1],points[N1+N2+N3,2]-yw*0.001,pos=1,"E3")
|
