关系网络图代码

admin

1. plottest111=function(r,rpoints,mat)
   
2. {
   
3. a.col=ncol(mat)
   
4.   if(a.col%%2==0)
   
5.    {
   
6.      n=2*(1:a.col)-1
   
7.      xplace=r*cos(pi/2+2*pi/a.col/2*n)
   
8.      yplace=r*sin(pi/2+2*pi/a.col/2*n)
   
9.    }else
   
10.    {
    
11.      xplace=r*cos(pi/2+(2*pi/a.col*((1:a.col)-1)))
    
12.      yplace=r*sin(pi/2+(2*pi/a.col*((1:a.col)-1)))
    
13.     }
    
14. plot(0,0,col="white",xlim=c(-1.25*r,1.25*r),ylim=c(-1.25*r,1.25*r),frame.plot=F,axes=F,xlab="",ylab="")
    
15. for(i in 1:a.col)
    
16. {
    
17.    for(j in 1:i)
    
18.        if(mat[i,j]!=0)
    
19.          {
    
20.            a=xplace[c(i,j)];b=yplace[c(i,j)]
    
21.            lines(x=a,y=b,col="gray",lwd=1.25*mat[i,j])
    
22.          }
    
23. }
    
24. par(new=T)
    
25. plot(x=xplace,y=yplace,col=rainbow(100)[50],pch=16,cex=rpoints,xlim=c(-1.25*r,1.25*r),ylim=c(-1.25*r,1.25*r),frame.plot=F,axes=F,xlab="",ylab="")
    
26. points(xplace,yplace,cex=rpoints,pch=1,col="red")
    
27. text(x=xplace,y=yplace,names(mat))
    
28. }
    
29. 
    
30. ###参数部分：r 画面的半径
    
31. r=3;rpoints=4
    
32. 
    
33. par(mfrow=c(2,4))
    
34. 
    
35. ###点的坐标以及画图
    
36. for(i in 3:10)
    
37. {
    
38. b=matrix(sample(0:5,replace=T,size=i^2,prob=c(0.6,0.1,0.08,0.08,0.08,0.06)),ncol=i)
    
39. b=as.data.frame(b)
    
40. names(b)=letters[1:i]
    
41. rownames(b)=letters[1:i]
    
42. b[lower.tri(b)]=b[upper.tri(b)]
    
43. b[!(upper.tri(b)|lower.tri(b))]=0
    
44. plottest111(b,r=r,rpoints=rpoints)
    
45. }
    

复制代码