Resumo. As avaliações rotineiras do caráter tempo de cozimento em feijão (Phaseolus vulgaris L.) podem ser efetuadas de distintas maneiras resultando em diferentes variáveis. Por vez, a análise estatística univariada não considera as interdependêcias entre as variáveis, podendo omitir importantes informações a respeito dos genótipos. Com isso, o objetivo do trabalho foi dispor uma proposta alternativa para análise do tempo de cozimento em feijão, permitindo a discriminação entre genótipos. O experimento utilizado para esta abordagem foi conduzido em condições de campo na safra agrícola do ano 2017/18 em Lages, Santa Catarina, Brasil. Os tratamentos foram compostos por doze genótipos, sendo quatro genitores, estruturados em dois cruzamentos BAF50 x BAF07 e BAF09 x IPR 88 Uirapuru, com suas gerações F2, F3, F8 e F9. O delineamento utilizado foi blocos casualizados, com dois blocos e duas observações em cada unidade experimental. Posteriormente a colheita, a variável resposta tempo de cocção dos grãos de feijão foi mensurada com o cozedor Mattson, sendo considerado o tempo de cocção das 13 hastes iniciais. Na análise multivariada, as variáveis tempo de cocção da segunda (TCH2), décima segunda (TCH12) e décima terceira haste (TCH13) foram utilizadas com base em sua significância pelo método de seleção de variáveis passo a passo (stepwise). A análise de variância multivariada demonstrou diferença entre os genótipos (P<0,05). A partir da matriz de dissimilaridade com as distâncias de Mahalanobis e o dendrograma de agrupamento, foi possível verificar as distâncias dos genótipos derivados dos cruzamentos BAF50 x BAF07 e BAF09 x IPR 88 Uirapuru. Com isso, a análise multivariada possibilitou a discriminação dos genótipos, adicionalmente o cruzamento BAF50 x BAF07 demonstrou maiores estimativas de dissimilaridade nas progênies.