A quantidade de informações disponíveis representa um desafio para a engenharia geotécnica, pois, quanto maior o volume de dados, maior o potencial de redução das incertezas e de aprimoramento da tomada de decisão. Assim, é essencial compreender o comportamento da distribuição dos parâmetros, com intuito de representar a anisotropia geotécnica. Para isso, é necessário obter a função de distribuição acumulada, com intuito de subsidiar indicadores de desempenho mais assertivos, como a probabilidade de falha e o índice de confiabilidade. O presente trabalho traz uma metodologia de aderência da distribuição acumulada, para 18 parâmetros, utilizando o software R e o pacote fitdistrplus. Os resultados obtidos foram validados a partir de critérios de informação de Akaike (AIC) e Bayesiano (BIC), para diferentes parâmetros, e a metodologia mostrou-se uma ferramenta fundamental para a realização de simulações numéricas probabilísticas mais precisas, principalmente para análises de variáveis aleatórias (Random Variable – RV).