I rework Zhenghua's script so that each BDE score can be analyzed. I also left the values in the negative format. BDE score 7 and 8 are the best network .
##### Rework
> bde=vector()
> bde[1] = -3792.715
> bde[2] = -3681.208
> bde[3] = -3630.236
> bde[4] = -3627.093
> bde[5] = -3615.076
> bde[6]= -3627.093
> bde[7] = -3603.086
> bde[8]= -3603.086
> bde[9]= -3621.687
>
> per1 =vector()
> for(i in 1:9){
+ per1[i] = exp(-log(sum(exp(bde-bde[i]))))
+ }
>
> per1
[1] 2.208718e-83 5.902217e-35 8.088602e-13 1.874395e-11
[5] 3.102972e-06 1.874395e-11 4.999984e-01 4.999984e-01
[9] 4.175005e-09