A method of controlling a heating, ventilation, and air-conditioning, HVAC, system, the HVAC system comprising a plurality of HVAC devices, the method comprising: • training a plurality of comfort models to be used for controlling the heating, ventilation, and air-conditioning, HVAC, system, wherein each comfort model of the plurality of comfort models is associated with a respective user of a plurality of users, wherein the training of a respective comfort model of the plurality of comfort models comprising: a user device associated with the respective comfort model providing a plurality of user feedbacks, each user feedback of the plurality of user feedbacks describing a thermal comfort of the respective user associated with the user device, wherein each user feedback of the plurality of user feedbacks is associated with a respective point in time, wherein the plurality of user feedbacks comprises two or more cold feedbacks indicating that the user of the user device associated with the comfort model feeling cold; • determining a weight matrix, the weight matrix includes, for each trained comfort model of the plurality of trained comfort models, a respective weight value for each HVAC device of the plurality of HVAC devices; • determining respective running parameters of each HVAC device of the plurality of HVAC devices such that a total predicted thermal comfort representing all predicted thermal comforts, which are generated by the plurality of comfort models responsive to inputting the running parameters of each of the plurality of HVAC devices into each comfort model, is increased, considering the determined weight value of each HVAC device of the plurality of HVAC devices associated with each comfort model of the plurality of comfort models; and • controlling the HVAC system in accordance with the determined running parameters. • training a plurality of comfort models to be used for controlling the heating, ventilation, and air-conditioning, HVAC, system, wherein each comfort model of the plurality of comfort models is associated with a respective user of a plurality of users, wherein the training of a respective comfort model of the plurality of comfort models comprising: a user device associated with the respective comfort model providing a plurality of user feedbacks, each user feedback of the plurality of user feedbacks describing a thermal comfort of the respective user associated with the user device, wherein each user feedback of the plurality of user feedbacks is associated with a respective point in time, wherein the plurality of user feedbacks comprises two or more cold feedbacks indicating that the user of the user device associated with the comfort model feeling cold; • determining a weight matrix, the weight matrix includes, for each trained comfort model of the plurality of trained comfort models, a respective weight value for each HVAC device of the plurality of HVAC devices; • determining respective running parameters of each HVAC device of the plurality of HVAC devices such that a total predicted thermal comfort representing all predicted thermal comforts, which are generated by the plurality of comfort models responsive to inputting the running parameters of each of the plurality of HVAC devices into each comfort model, is increased, considering the determined weight value of each HVAC device of the plurality of HVAC devices associated with each comfort model of the plurality of comfort models; and • controlling the HVAC system in accordance with the determined running parameters.