Motivation:The branching processes model yields unevenly stochastically distributed data that consists of sparse and dense regions. The work tries to solve the problem of an precise evaluation of a parameter for this type of a model. The application of branching processes model to cancer cell evolution has many difficulties like high dimensionality and rare appearance of a result of interest. Moreover we would like to solve the ambitious task of obtaining the coefficients of the model reflecting the relationship of driver genes mutations and cancer hallmarks on the basis of personal data of variant allele frequencies. Results: The Approximate Bayesian computation method based on Isolation kernel is designed. The method includes a transformation row data to a Hilbert space (mapping) and measure the similarity between simulation points and maxima weighted Isolation kernel mapping related to observation point. Also, we designed heuristic algorithm to find parameter estimation without gradient calculation and dimension-independent. The advantage of the proposed machine learning method is shown for multidimensional test data as well as for an example of cancer cell evolution.