A product quality prediction method for mass customization is provided. When a production system has a status change, data of sets of process parameters and actual measurement values of workpiece samples processed before the status change occurs, and data of sets of process parameters and actual measurement values of few workpiece samples processed after the status change occurs are used for build or retrain a prediction model, thereby predicting a metrology value of a next workpiece.