That's what you *should* do if you want it clean. See LTspice help file: '.PARAM -- User-Defined Parameters' for an example of passing parameters to a subcircuit. You'd probably also want to parameterize the gain.
Its already been half-way set up for it with:
.MODEL Rb RES
+ R=1
.MODEL Re RES
+ R=1
The '1' in those two model statements is a unity multiplier for the values in the R1 and R2 component lines respectively.
Try:
* DTC114EE DIGITAL NPN BJT model
* Date: 2007/01/11 mod. 2020/09/04 Ian.M
* BJT with resistors, parameterized for Monte Carlo
*****************C B E
.SUBCKT DTC114EE 1 2 3
.params Krb=1 Kre=1 Kbeta=1 ; default to unit multiplier
R1 2 Base Rb 10k
R2 Base 3 Re 10k
Q1 1 Base 3 QDTC1
.MODEL Rb RES
+ R={Krb}
.MODEL Re RES
+ R={Kre}
.MODEL QDTC1 NPN
+ IS=20.000E-15
+ BF={198.80*Kbeta}
+ VAF=15
+ IKF=.28787
+ ISE=20.001E-15
+ NE=1.6289
+ BR=8.4556
+ VAR=100
+ IKR=.2645
+ ISC=8.0686E-12
+ NC=2.3678
+ NK=.75224
+ RE=.2
+ RB=10.670
+ RC=1.1274
+ CJE=11.342E-12
+ MJE=.38289
+ CJC=4.0230E-12
+ MJC=.34629
+ TF=500.00E-12
+ XTF=580.91
+ VTF=195.13
+ ITF=19.351
+ TR=200.00E-9
+ XTB=1.5000
.ENDS DTC114EE
Then you can pass in per-instance parameters to your DTC114EE model, set in the generic Component Attribute Editor:
e.g.
Krb=1.1 Kre=0.9 Kbeta=0.8
which gives you a DTC114EE with its series base resistor 10% high, its shunt base resistor 10% low and 20% below nominal gain, or for your specific application you'd need to use multipliers derived from your Monte Carlo or Worst Case functions.