Objetivou-se com este trabalho estimar parâmetros genéticos para características de produção de leite (produção de leite em até 305 dias – PL305), peso (peso aos 120 – P120, aos 365 – P365 e aos 550 dias – P550) e reprodutiva (idade ao primeiro parto - IPP) em bovinos da raça Guzerá utilizando modelos multicaracterísticos e comparar metodologias que se baseiam na utilização de dados censurados de idade ao primeiro parto. Os dados foram cedidos pela Associação Brasileira dos Criadores de Zebu (ABCZ) e Embrapa Gado de Leite em parceria com o Centro Brasileiro de Melhoramento do Guzerá (CBMG2). O arquivo de pedigree (120.599 animais) incluiu animais com registros fenotípicos e todos os ancestrais conhecidos. Registros censurados foram definidos como dados de IPP que extrapolaram os limites de 740 e 1860 dias. Os registros de IPP (69.157 informações) foram analisados por meio de sete diferentes metodologias: método linear convencional (LM); de simulação (SM); de penalidade (PM); modelos bicaracterística limiar-linear em que se considerou (TLcens) ou não (TLmiss) informações prévias sobre os registros censurados; e a metodologia de análise de sobrevivência por meio do modelo de riscos proporcionais de Weibull segmentado em que se considerou (PWPHcens) ou não (PWPH) os registros censurados. Para as análises de sobrevivência, os valores para o critério de informação da deviance (DIC) sugerem o uso de 0 e 2 nós na função base para os métodos PWPH e PWPHcens, respectivamente. Os componentes de variância genética aditiva estimados para os métodos LM, PM e TLmiss foram similares. As estimativas de herdabilidade para IPP variaram de 0,19 (TLcens) à 0,28 (SM) e 0,40 (PWPH) e 0,46 (PWPHcens). De forma geral, as correlações entre os valores genéticos obtidos por meio das diferentes metodologias e a porcentagem de touros selecionados em comum variaram de -0,28 (SM x PWPH) à 0,99 (TLmiss x LM), indicando reordenamento moderado entre os animais. As comparações baseadas na metodologia de validação cruzada, indicam a metodologia TLmiss como a melhor opção para predição dos valores genéticos dos animais para a característica IPP na população Guzerá. Para estimação dos parâmetros genéticos utilizando modelos multicaracterísticos, foram considerados os efeitos
sistemáticos de sexo e idade ao parto. Os efeitos aleatórios genético aditivo e de grupo de contemporâneo (rebanho, ano e estação de parto) foram considerados. Adicionalmente, os efeitos aleatórios genético aditivo materno e de ambiente permanente materno foram considerados para a característica peso à desmama. As estimativas de herdabilidade foram 0,29 (PL305), 0,42 (P120), 0,49 (P365), 0,56 (P550) e 0,25 (IPP). As correlações genéticas entre as características de peso foram maiores que 0,83 e entre PL305 e as demais foram de 0,25 (PL305 x P205), 0,32 (PL305 x P365) e 0,36 (PL305 x P550). A característica IPP foi negativamente correlacionada com as características de leite e de peso. Os resultados sugerem que a seleção para a produção de leite não compromete a seleção para características de peso e reprodutivas, bem como é passível a inclusão de dados censurados de IPP na avaliação genética por meio do uso de um modelo de limiar-linear em bovinos da raça Guzerá.