Clustered tensegrity structures integrated with continuous cables are lightweight, foldable, and deployable. Thus, they can be used as flexible manipulators or soft robots. The actuation process of such soft structure has high probabilistic sensitivity. It is essential to quantify the uncertainty of actuated responses of the tensegrity structures and to modulate their deformation accurately. In this work, we propose a comprehensive data-driven computational approach to study the uncertainty quantification (UQ) and probability propagation in clustered tensegrity structures, and we have developed a surrogate optimization model to control the flexible structure deformation. An example of clustered tensegrity beam subjected to a clustered actuation is presented to demonstrate the validity of the approach and its potential application. The three main novelties of the data-driven framework are: (1) The proposed model can avoid the difficulty of convergence in nonlinear Finite Element Analysis (FEA), by two machine learning methods, the Gauss Process Regression (GPR) and Neutral Network (NN). (2) A fast real-time prediction on uncertainty propagation can be achieved by the surrogate model, and (3) Optimization of the actuated deformation comes true by using both Sequence Quadratic Programming (SQP) and Bayesian optimization methods. The results have shown that the proposed data-driven computational approach is powerful and can be extended to other UQ models or alternative optimization objectives.