R-matrix calculations for few-quark bound states