Participant characteristics and data availability
All 240 participants from six sites completed baseline assessment and were allocated equally to four groups (n = 60 each). Every site contributed 10 participants to each group. Overall comfort for week 4 was available for all participants, who therefore entered the allocated-condition primary analysis. Mean age was 31.5 ± 7.4 years, 125 participants (52.1%) were women, mean body mass index was 23.0 ± 3.2 kg/m2, and mean weekly work duration was 41.5 ± 5.9 h. Table 4 presents baseline characteristics without post-randomization significance testing. The primary outcome was complete for all 240 participants. Available case sample sizes for week 4 secondary outcomes ranged from 235 to 240; detailed group-specific sample sizes for these measures are provided in the respective efficacy summaries below.
Intervention delivery and exposure separation
The adjustable-furniture package separated the groups in adjustability, anthropometric fit, and desk-height behavior. Mean anthropometric fit was 4.95 ± 0.70 with adjustable furniture and 2.46 ± 0.73 with fixed furniture. Adjustable groups recorded 5.27 ± 1.13 desk-height transitions/day and 75.84 ± 20.38 min/day above the participant-specific seated threshold; fixed groups recorded 2.07 ± 0.99 transitions/day and 19.78 ± 18.86 min/day. The latter values can reflect device repositioning, reference-height variation, or other desk-height records and are not interpreted as fixed-desk standing adherence. Raw height records, the seated-threshold definition, device allocation, and group codes were rechecked against the deposited data.
The optimized-IEQ package produced lower temperature (24.15 ± 0.75°C versus 26.30 ± 0.75°C), relative humidity (54.57% ± 4.99% versus 66.65% ± 6.11%), CO₂ (750.15 ± 120.47 ppm versus 1080.35 ± 123.98 ppm), noise (43.39 ± 3.92 dBA versus 50.94 ± 4.14 dBA), and PM2.5 (12.33 ± 4.84 µg/m3 versus 18.16 ± 4.49 µg/m3), and higher illuminance (439.30 ± 68.70 lux versus 293.03 ± 72.82 lux), than basic IEQ. Table 5 reports these measures along with adherence, deviations, per-protocol counts, and daily mean target attainment.
Thus, furniture assignment primarily changed workstation fit and desk-height behavior, whereas IEQ assignment changed the monitored thermal, ventilation-related, lighting, acoustic, and particulate conditions. These are package-level contrasts and do not isolate any single component.
Overall comfort
Week 4 overall comfort was 3.99 ± 0.60, 4.84 ± 0.65, 4.33 ± 0.55, and 5.62 ± 0.53 in the fixed-basic, adjustable-basic, fixed-optimized, and adjustable-optimized groups, respectively. Adjustable furniture increased comfort by 0.86 points (95% CI, 0.62–1.09; P < 0.001), optimized IEQ increased it by 0.34 points (95% CI, 0.12–0.55; P = 0.002), and their interaction was 0.44 points (95% CI, 0.14–0.74; P = 0.005) in the site-adjusted HC3 model. Figure 3 shows the four groups and allocated-condition population, raw week 4 comfort distributions, baseline-to-week 4 group means with 95% CIs, and the site-adjusted factorial effects.
Posture comfort and musculoskeletal outcomes
Week 4 posture comfort scores were 2.72 ± 0.92, 4.33 ± 0.93, 2.82 ± 0.92, and 4.99 ± 0.85 across the four groups. The adjustable-furniture effect was 1.61 points (95% CI, 1.28–1.95; P < 0.001; q < 0.001), the optimized-IEQ effect was 0.11 points (95% CI, −0.23 to 0.44; P = 0.534; q = 0.712), and the interaction effect was 0.55 points (95% CI, 0.09–1.01; P = 0.020; q = 0.044).
Week 4 neck discomfort scores were 4.51 ± 1.66, 3.27 ± 1.71, 4.75 ± 1.66, and 2.57 ± 1.32 across the four groups. The adjustable-furniture effect was −1.24 points (95% CI, −1.85 to −0.62; P < 0.001; q < 0.001), the optimized-IEQ effect was 0.26 points (95% CI, −0.34 to 0.87; P = 0.391; q = 0.712), and the interaction effect was −0.97 points (95% CI, −1.80 to −0.14; P = 0.022; q = 0.044).
Week 4 lower-back discomfort scores were 5.16 ± 1.78, 3.46 ± 1.83, 5.04 ± 1.59, and 3.02 ± 1.55 across the four groups. Adjustable furniture reduced discomfort by 1.70 points (effect estimate, −1.70; 95% CI, −2.36 to −1.04; P < 0.001; q < 0.001). Neither the optimized-IEQ effect (−0.12; 95% CI, −0.74 to 0.51) nor the interaction effect (−0.32; 95% CI, −1.19 to 0.56) was significant. Figure 4 presents the raw distributions of anthropometric fit and posture comfort, mean desk-height transitions per day, and above-threshold time per day with 95% CIs, and mean neck and lower-back discomfort scores with 95% CIs.
Environmental comfort and work-related outcomes
Environmental comfort scores were 5.28 ± 0.57, 5.29 ± 0.60, 5.83 ± 0.49, and 5.88 ± 0.54 across the four groups. Optimized IEQ increased environmental comfort by 0.55 points (95% CI, 0.35–0.74; P < 0.001; q < 0.001). Neither the furniture effect (−0.00; 95% CI, −0.22 to 0.21) nor the interaction effect (0.04; 95% CI, −0.24 to 0.33) was significant.
Task accuracy was 83.49% ± 6.44%, 82.99% ± 6.75%, 82.15% ± 5.96%, and 85.24% ± 6.27% across the four groups. The furniture and optimized-IEQ main effects were not significant. The interaction effect was 3.72 percentage points (95% CI, 0.44–7.00; P = 0.026), but it did not remain significant after within-family false-discovery-rate correction (q = 0.078). It was therefore interpreted as exploratory rather than confirmatory.
Task-speed data were available for 235 participants. Neither furniture (0.10 correct units/35 min; 95% CI, −2.95 to 3.16), optimized IEQ (1.51 correct units/35 min; 95% CI, −1.43 to 4.44), nor their interaction (0.24 correct units/35 min; 95% CI, −4.02 to 4.51) had a significant effect. The earlier claim of a significantly optimized IEQ effect on task speed was not reproducible and was therefore removed.
Adjustable furniture and optimized IEQ increased perceived productivity by 0.70 points (95% CI, 0.40–1.00) and 0.78 points (95% CI, 0.47–1.10), respectively (both P < 0.001). Their interaction was not significant. The corresponding effects on workspace satisfaction were 0.88 points for adjustable furniture (95% CI, 0.60–1.15), 0.46 points for optimized IEQ (95% CI, 0.21–0.72), and 0.77 points for their interaction (95% CI, 0.40–1.14); all q values were < 0.001.
Mental fatigue decreased with adjustable furniture (−0.56 points; 95% CI, −0.97 to −0.16; P = 0.006; q = 0.006) and optimized IEQ (−0.81 points; 95% CI, −1.23 to −0.39; P < 0.001; q < 0.001). The interaction effect was not significant. No significant furniture, optimized-IEQ, or interaction effects were detected for creativity. Table 6 presents group-specific week 4 values, site-adjusted HC3 factorial effects, false discovery rate q values, and the results of the per-protocol and major-deviation-exclusion analyses.
Figure 5 presents mean objective IEQ conditions, perceived environmental comfort, task accuracy and speed, perceived productivity, workspace satisfaction, mental fatigue, and creativity. All error bars represent 95% CIs, and available sample sizes are shown for outcomes with incomplete data.
Sensitivity analyses
Furniture-package adherence was 59/60 participants (98.3%) in each adjustable-furniture group. IEQ-package adherence was 57/60 participants (95.0%) in the fixed-furniture plus optimized-IEQ group and 54/60 participants (90.0%) in the adjustable-furniture plus optimized-IEQ group. Ten major deviations were recorded: two in the fixed-furniture plus basic-IEQ group, two in the adjustable-furniture plus basic-IEQ group, four in the fixed-furniture plus optimized-IEQ group, and two in the adjustable-furniture plus optimized-IEQ group.
The per-protocol population included 221 participants. The corresponding overall-comfort estimates—0.79 points for furniture (95% CI, 0.56–1.02), 0.32 points for optimized IEQ (95% CI, 0.09–0.55), and 0.50 points for the interaction (95% CI, 0.18–0.81)—were consistent with the allocated-condition analysis. Estimates from the major-deviation-exclusion analysis, which included 230 participants, were also consistent with the primary findings. Based on the verified safety and protocol deviation logs compiled from all participating sites, there were no intervention suspensions, safety-triggered workstation reassessments, intervention-related adverse events, or discontinuations due to harm during the four-week period.
Participant-clustered repeated-measures models showed furniture-by-week and IEQ-by-week changes in overall comfort (P < 0.001 and P = 0.009, respectively), with no significant three-way interaction (P = 0.274). Furniture-by-week changes were also detected for posture comfort (P < 0.001), neck discomfort (P = 0.032), lower-back discomfort (P = 0.002), and mental fatigue (P = 0.023). An IEQ-by-week change was detected for mental fatigue (P = 0.001).
Exploratory mediation analyses did not support either an ergonomic indirect effect (0.06; bootstrap 95% CI, −0.08 to 0.19) or an environmental indirect effect (0.04; bootstrap 95% CI, −0.03 to 0.12). Restricted cubic-spline analyses showed no evidence of nonlinearity for temperature (P = 0.187), CO₂ (P = 0.947), illuminance (P = 0.852), or noise (P = 0.142).
In summary, the reproducible analyses support distinct but complementary package effects. Adjustable furniture most consistently improved overall comfort and posture comfort and reduced musculoskeletal discomfort, whereas optimized IEQ improved overall comfort and environmental comfort. Their combination provided additional benefits for overall comfort and selected posture-comfort and workspace-satisfaction outcomes. Objective task speed and creativity did not differ significantly among groups, and the task-accuracy interaction was not retained after correction for multiple comparisons.
DATA AVAILABILITY:
De-identified participant-level, weekly-assessment, and workday-exposure datasets are publicly available from Zenodo (doi: https://zenodo.org/records/20657821).

Figure 1: Study flow, factorial structure, and assessment schedule. (A) Provisional screening flow inferred from the randomized dataset: employees assessed for eligibility (n = 240), eligible and consented (n = 240), excluded (n = 0), randomly allocated (n = 240), assigned to four groups (n = 60 per group), and included in the allocated-condition primary-outcome analysis at week 4. (B) The 2 × 2 furniture-by-IEQ factorial structure. (C) Baseline assessment, installation and commissioning, continuous monitoring, weekly checks, and week 4 assessment. Abbreviations: IEQ = indoor environmental quality. Please click here to view a larger version of this figure.

Figure 2: Multicomponent furniture and IEQ intervention packages. (A) Fixed and adjustable workstation configurations and desk-height ranges. (B) Participant-specific targets for chair, desk, monitor, keyboard, mouse, reach, wrist position, viewing distance, and sit-to-stand clearance. (C) Monitor placement, sampling intervals, operational targets, and the five-variable IEQ composite; PM2.5 was assessed separately. Abbreviations: IEQ = indoor environmental quality; PM2.5 = particulate matter ≤2.5 µm in aerodynamic diameter. Please click here to view a larger version of this figure.

Figure 3: Week 4 overall comfort and site-adjusted factorial effects. (A) Allocated-condition groups (n = 60 per group; total n = 240). (B) Participant-level observations, box plots showing the median, interquartile range, and whiskers extending to 1.5 × the interquartile range, and group means. (C) Baseline and week 4 means; error bars represent 95% CIs. (D) Site-adjusted HC3 factorial effects with 95% CIs. Abbreviations: CI = confidence interval; IEQ = indoor environmental quality. Please click here to view a larger version of this figure.

Figure 4: Furniture fit, desk-height behavior, posture comfort, and musculoskeletal discomfort.
(A) Week 4 anthropometric-fit scores. Participant-level observations are overlaid on box plots showing the median, interquartile range, and whiskers extending to 1.5 × the interquartile range. Group means and 95% confidence intervals (CIs) are also shown. (B) Mean desk-height transitions per day and above-threshold time per day. Error bars represent 95% CIs. Above-threshold time is defined as the number of minutes per day during which desk height exceeded the participant-specific seated reference height and does not confirm standing posture. (C) Week 4 posture-comfort scores. Participant-level observations are overlaid on box plots showing the median, interquartile range, and whiskers extending to 1.5 × the interquartile range. Group means, 95% CIs, and site-adjusted factorial effect estimates are shown. (D) Week 4 neck and lower-back discomfort scores. Group means and 95% CIs are shown together with the adjusted furniture effects for both outcomes. Abbreviations: AB = adjustable furniture plus basic IEQ; AO = adjustable furniture plus optimized IEQ; FB = fixed furniture plus basic IEQ; FO = fixed furniture plus optimized IEQ; CI = confidence interval; IEQ = indoor environmental quality. Please click here to view a larger version of this figure.

Figure 5: Objective IEQ conditions, environmental comfort, and work-related outcomes. (A) Mean objective indoor environmental quality (IEQ) conditions under the basic and optimized IEQ packages for temperature, relative humidity, carbon dioxide (CO₂), task illuminance, equivalent continuous noise level, and particulate matter ≤2.5 µm in aerodynamic diameter (PM2.5). Error bars represent 95% confidence intervals (CIs). (B) Group means and 95% CIs for thermal comfort, air-quality comfort, visual comfort, acoustic comfort, and the environmental-comfort composite. (C) Week 4 objective task performance, including task accuracy (%) and task speed (correct units per 35 min). Box plots show the median, interquartile range, and whiskers extending to 1.5 × the interquartile range; points represent participant observations, and available sample sizes are shown for outcomes with incomplete data. (D) Week 4 perceived productivity, workspace satisfaction, mental fatigue, and creativity. Points indicate group means, and error bars represent 95% CIs. Available sample sizes are shown where outcomes are incomplete. Abbreviations: CO₂ = carbon dioxide; CI = confidence interval; IEQ = indoor environmental quality; PM2.5 = particulate matter ≤2.5 µm in aerodynamic diameter. Please click here to view a larger version of this figure.
Table 1: Multicomponent intervention packages, operational targets, delivery procedures, and adherence criteria. Abbreviations: CO₂ = carbon dioxide; IEQ = indoor environmental quality; PM2.5 = particulate matter ≤2.5 µm in aerodynamic diameter. Please click here to download this file.
Table 2: Study questions and tasks, scoring and derivation rules, measurement-property limitations, practical interpretation, and assessment schedule. Abbreviations: CO₂ = carbon dioxide; IEQ = indoor environmental quality; PM2.5 = particulate matter ≤2.5 µm in aerodynamic diameter. Please click here to download this file.
Table 3: Statistical analysis framework, including primary, secondary, repeated-measures, sensitivity, and explicitly labeled exploratory analyses. Abbreviations: CI = confidence interval; CO₂ = carbon dioxide; IEQ = indoor environmental quality; PM2.5 = particulate matter ≤2.5 µm in aerodynamic diameter. Please click here to download this file.
Table 4: Baseline participant, occupational, workspace, and outcome characteristics, with site-by-group counts. Values are presented as mean ± standard deviation unless otherwise indicated. No significance testing of baseline characteristics was performed. Abbreviation: IEQ = indoor environmental quality. Please click here to download this file.
Table 5: Furniture fit and desk-height behavior, objective IEQ conditions, perceived environmental comfort, adherence, protocol deviations, per-protocol counts, and target attainment. Environmental-comfort values and objective conditions are descriptive. Target-attainment percentages represent the proportion of workday daily means within the prespecified range rather than the proportion of all one-minute intervals. PM2.5 was assessed separately from the five-variable composite. The second worksheet reports target attainment by site. Abbreviations: CO₂ = carbon dioxide; dBA = A-weighted decibels; IEQ = indoor environmental quality; PM2.5 = particulate matter ≤2.5 µm in aerodynamic diameter. Please click here to download this file.
Table 6: Week 4 outcomes, site-adjusted factorial effects, and sensitivity analyses. Group values are presented as mean ± standard deviation with the available n. Effects were estimated using factorial linear models with six site indicators, HC3 robust standard errors, and baseline adjustment when the corresponding baseline variable was available. Positive estimates indicate higher scores, whereas negative estimates for discomfort and fatigue indicate improvement. The primary outcome was not adjusted for multiplicity. For secondary outcomes, q values were calculated using the Benjamini–Hochberg procedure within outcome families. The second worksheet reports the per-protocol and major-deviation-exclusion analyses. Abbreviations: CI = confidence interval; IEQ = indoor environmental quality. Please click here to download this file.