Abstract-We investigated the constitutive parameters in the rate-and state-dependent friction (RSF) law by conducting numerical simulations, using the friction data from large-scale biaxial rock friction experiments for Indian metagabbro. The sliding surface area was 1.5 m long and 0.5 m wide, slid for 400 s under a normal stress of 1.33 MPa at a loading velocity of either 0.1 or 1.0 mm/s. During the experiments, many stick-slips were observed and those features were as follows. (1) The friction drop and recurrence time of the stick-slip events increased with cumulative slip displacement in an experiment before which the gouges on the surface were removed, but they became almost constant throughout an experiment conducted after several experiments without gouge removal. (2) The friction drop was larger and the recurrence time was shorter in the experiments with faster loading velocity. We applied a one-degree-of-freedom spring-slider model with mass to estimate the RSF parameters by fitting the stick-slip intervals and slip-weakening curves measured based on spring force and acceleration of the specimens. We developed an efficient algorithm for the numerical time integration, and we conducted forward modeling for evolution parameters (b) and the state-evolution distances (L c ), keeping the direct effect parameter (a) constant. We then identified the confident range of b and L c values.Comparison between the results of the experiments and our simulations suggests that both b and L c increase as the cumulative slip displacement increases, and b increases and L c decreases as the loading velocity increases. Conventional RSF laws could not explain the large-scale friction data, and more complex state evolution laws are needed.