Quote: 
That GWINC link is more than a year old. You're best off just updating your CVS checkout of the code, or getting a new zip file from someone else if your CVS is broken. When I run gwinc with nomm.m, I get R_BNS = 189.5 Mpc.

I just saw you comment. I'll find an update version for GWINC.
Anyway, I have a code to plot the result. I will use it on an updated code.
Some material parameters in the calculation are:
 The concentration of TiO2 is 25 % based on Harry 2010 paper (Class. Quantum Grav. 27 (2010) 084006).
 Young's Modulus of Ta2O5 doped with TiO2 (25%) is ~ 140 +/ 6 GPa (one sigma). (Abernathy et al. 2014).
 Loss angle of TiO2 doped Ta2O5 is (2.3 +/ 0.4) x 10^4 [Martin 2008, Martin2009] (@300K, heat treated)
 Loss angle of SiO2 is (6.3 +/ 2.6) x10^5 [LMA2014, Crooks2006, Crooks2004]
In the IFOModel_rnd.m file which is a copy of IFOModel.m for material params, I use normrnd(mean,sigma) to generate the random value of the material parameters.
Note, for loss angle of SiO2, I have to use abs command to make sure that all the generated values are greater than zero.
This is because the mean is comparable to the standard deviation, and sometime it gives negative values.
Note: I have not taken the coherent between the loss and Young's modulus of Ta2O5 into account yet. I have to read how they measure this more carefully.
Here are prelim results from the above numbers.
above: a histogram of BNS range, due to uncertainties in loss angles/young's moduli of the coatings.
The mean of the histogram is slightly less than the nominal value from GWINC because the mean values of loss angles for fused silica (6e5) is slightly higher than the original value (4e5) used in the code.
above: histograms and Gaussian fits for BR noise (blue/green) and total noise and its Gaussian fit (red/cyan) at 100 Hz in the strain unit.
