In this paper the author proposes a method of implementing a numerical model for threshold voltage (V_th) shift in organic thin-film transistors (QTFTs) into SPICE tools. V_th shift is first numerically modeled by dividing the shift into sequentially ordered groups. The model is then used to derive a simulations model which takes into simulation parameters and calculation complexity. Finally, the numerical and simulation models are implemented in AIM-SPICE. The SPICE simulation results agree well with the V_th shift obtained from an OTFT fabricated without any optimization. The proposed method is also used to implement the stretched-exponential time dependent V_th shift in AIM-SPICE and the results show the proposed method is applicable to various types of V_th shifts.