I am working on the AGC circuit below and was recently made aware of the fact that Vgs(off) varies widely for most JFETs including the JFET I was planning to use in my circuit J113 which has a minimum and maximum Vgs(off) and -0.5V and -3V source. The JFET is one of the critical components of my circuit as I use it as a voltage controlled resistor.
For reference my AGC circuit is below

I would like to get a better understanding of how the Drain-Source resistance will vary with Vgs(off) so I have created the circuit below. This circuit consists only of a singular jfet and a variable input DC signal. Note that I am using the .model spice directive to vary the Vgs(off) of the jfet.
Once I simulate the circuit I use a python script and PyLTSpice to load the data and calculate the resistance across the jfet as a function of the gate's voltage.
Using the .model directive I can change Vto (Vgs) of the jfet to whatever value I want. For example if I change Vto to 0.5V I get a plot like the one below for resistance vs gate voltage.
My first question: Is varying Vto an accurate way to model the JFET with differing values of Vgs(off) or are there other parameters which would also change if Vto changes? Per the lt spice model description it looks as though beta and lambda are the only other two parameters likely to impact the resistance of the jfet.
My second question: How can I negate the impacts of variable Vgs(off) in my design? Should I design based on the performance for minimum Vgs(off)? Are there other discrete that could act like a controllable resistor that have lower variability?
Thank You


