Accurate joint kinematics plays an important role in estimating joint kinetics in musculoskeletal simulations. Biplanar fluoroscopic (BPF) systems have been introduced to measure skeletal kinematics with six degrees-of-freedom. The purpose of this study was to model knee kinematic coupling using knee kinematics during walking, as measured by the BPF system. Seven healthy individuals (mean age, 23 +/- 2 yr) performed treadmill walking trials at 1.2 m/s. Knee kinematics was regressed separately for the swing and stance phases using a generalized mixed effects model. Tibial anterior translation function was y = 0.20x - 3.09 for the swing phase and y = 0.31x -0.54 for the stance phase, where x was the flexion angle and y was the tibial anterior translation. Tibial lateral and inferior translation were also regressed separately for the stance phase and the swing phase. Tibial external rotation was y = - 0.002x(2) + 0.19x -0.64 for the swing phase and y = -0.19x -1.22 for the stance phase. The tibial adduction rotation function was also calculated separately for the stance and swing phase. The study presented three-dimensional coupled motion in the knee during the stance and swing phases of walking, and demonstrated the lateral pivoting motion found in previous studies. This expanded understanding of secondary knee motion functions will benefit musculoskeletal simulation and help improve the accuracy of calculated kinetics.