We measured the $^{8}\mathrm{Li}(p,p)^{8}\mathrm{Li}$, $^{8}\mathrm{Li}(p,d)^{7}\mathrm{Li}$, and $^{8}\mathrm{Li}(p,\ensuremath{\alpha})^{5}\mathrm{He}$ reactions at low energies using the thick target inverse kinematics method, with a polyethylene ${[{\mathrm{CH}}_{2}]}_{n}$ target and a radioactive $^{8}\mathrm{Li}$ beam available at the Radioactive Ion Beams in Brazil facility of S\~ao Paulo. By measuring simultaneously several reaction channels $(p,p)$, $(p,d)$, and $(p,\ensuremath{\alpha})$, the still uncertain high-lying resonances of $^{9}\mathrm{Be}$, close to the proton threshold, can be studied and their parameters, such as energy, width, and spin parity can be more reliably determined. The experimental excitation functions of the reactions $^{8}\mathrm{Li}(p,p)^{8}\mathrm{Li}$, $^{8}\mathrm{Li}(p,d)^{7}\mathrm{Li}$, and $^{8}\mathrm{Li}(p,\ensuremath{\alpha})^{5}\mathrm{He}$ were analyzed using the $R$-matrix theory, which allows us to infer the properties of the resonances. Multichannel $R$-matrix analysis provides evidence for a significant clustering in the $(p,d)$ channel. The experimental data and the multichannel $R$-matrix analysis will be presented.