Statistical analysis
There were an estimated 2,000 community pharmacies that were operational in Malaysia. Assuming that there were a total of 4,000 community pharmacists working in community pharmacies, we estimated that we needed 187 responses in order to achieve a 95% confidence level and 7% confidence that these responses were representative of all pharmacists. Results were presented in frequencies, proportions, means, and standard deviations. Since not all participants completed the full questionnaire, data were presented based upon total responses received for each question. No data manipulation was performed. To investigate the associations between socio-demographic profile of the participants and outcomes of level of agreement, we performed a bivariate analysis after excluding those with incomplete responses. Multivariate regression analysis was conducted to explore the relationships between socio-demographic characteristics and the levels of agreement of community pharmacist within the three categories, namely safety at workplace, resilience, and support. A p value of less than 0.05 was considered to indicate statistical significance. All analyses were performed with the use of SPSS Statistics version 26 for Windows (IBM Corp., Armonk, NY, USA).