3. Materials and Methods
3.1. Study Design and Setting
This study employed a cross-sectional online survey design to examine associations between preferred future work arrangements and sociodemographic, experiential, occupational-perception, and attitudinal characteristics among adult respondents residing in the Roanoke Metropolitan Statistical Area (MSA), Virginia. The study was conducted after the acute disruption of the COVID-19 pandemic, during a period in which remote and hybrid work arrangements remained relevant to longer-term workforce decisions.
The Roanoke MSA is located in southwestern Virginia and comprises the independent cities of Roanoke and Salem and the counties of Botetourt, Craig, Franklin, and Roanoke, according to the metropolitan statistical area delineation issued by the U.S. Office of Management and Budget [
30]. The region had a population of approximately 315,000 residents according to the 2020 U.S. Census [
31]. As a mid-sized metropolitan region, Roanoke provides a relevant setting for examining preferred future work arrangements outside the major urban centers that have received greater attention in the remote-work literature. The location and constituent jurisdictions of the Roanoke Metropolitan Statistical Area are shown in
Figure 1.
The regional context is relevant to workforce and planning discussions because mid-sized metropolitan areas may differ from larger metropolitan regions in employment distribution, transportation infrastructure, commuting conditions, and housing characteristics. However, the present study examines stated preferences for future work arrangements rather than observed changes in commuting behavior, residential relocation, land use, or employment location.
3.2. Survey Administration and Sample Construction
Data were collected through an online survey administered using the QuestionPro platform between 15 April and 6 August 2023. The survey targeted respondents residing in the Roanoke MSA and collected information on demographic characteristics, current and preferred work arrangements, pandemic-era work-from-home experience, perceptions of remote-work possibilities within respondents’ professions, attitudes toward physical interaction in work settings, and selected residential location considerations.
The original QuestionPro export contained
1,122 response records. A sequential data-screening procedure was applied to construct the primary analytical sample. First, records not identified as completed surveys were excluded, removing 214 incomplete responses and leaving 908 completed responses. Second, 176 completed responses identified as duplicates in the QuestionPro export were removed, leaving 732 completed, nonduplicate responses. Third, because the study focuses on adult work-style preferences, 53 responses from individuals reporting an age below 18 years and 26 responses with missing age information were excluded, leaving 653 verified adult respondents. Finally, because the study is specifically situated in the Roanoke MSA, 17 adult responses without verified residence in one of the six constituent MSA jurisdictions were excluded. The final primary analytical sample therefore comprised
636 completed, nonduplicate responses from verified adult residents of the Roanoke MSA. The sequential construction of the primary analytical sample is summarized in
Table 1.
The final sample included respondents from all six jurisdictions constituting the Roanoke MSA: Roanoke City, Roanoke County, Salem City, Craig County, Franklin County, and Botetourt County.
The study employed a nonprobability online survey sample. Consequently, the findings are interpreted as associations within the analytical sample rather than as population estimates for all adults or workers in the Roanoke MSA. Participation may have been influenced by access to the online survey, willingness to complete the questionnaire, interest in work-arrangement issues, and availability during the survey period.
Item-level missing responses were retained as missing and were not imputed. Accordingly, valid sample sizes vary across analyses. Each statistical test is reported using the number of respondents with valid responses for all variables involved in that analysis.
3.3. Survey Instrument
The survey instrument collected information across five primary analytical domains: sociodemographic characteristics, current and preferred work arrangements, pandemic-era work-from-home experience, perceptions of remote-work possibilities within respondents’ professions, and attitudes toward physical interaction in work settings. The survey also included residential location-preference items, which were retained in the cleaned dataset but excluded from the primary analyses in the present manuscript.
Respondents reported gender, age group, marital status, highest educational attainment, monthly earnings, housing tenure, household size, and desired family size. Following exclusion of responses from individuals younger than 18 years, age was represented using six categories: 18–24, 25–34, 35–44, 45–54, 55–64, and above 64 years.
Educational attainment was recorded using eight categories: high school, trade/vocational/technical education, some college, associate degree, bachelor’s degree, master’s degree, professional degree, and doctorate. Monthly earnings were reported as $1,250 or less, $1,251–$3,333, greater than $3,333, or prefer not to say. Marital status was recorded as married, single, separated, divorced, widowed, or prefer not to say. Housing tenure was categorized as owner or renter.
Reported residential location was used to determine geographic eligibility for the primary analytical sample. Respondents were retained only when they reported residence in Roanoke City, Salem City, Botetourt County, Craig County, Franklin County, or Roanoke County.
Respondents identified both their current work arrangement and their preferred future work arrangement using the same four response categories:
Fully in-person, working at a physical office or location;
Hybrid, dominated by in-person work, with the majority of work time spent at a physical office or location and some remote work;
Hybrid, dominated by remote work, with the majority of work time spent working from home or another location and some in-person work; and
Fully remote, working from home or another location.
Preferred future work arrangement constituted the primary outcome variable. Current work arrangement was used only for paired analysis examining whether preferred future arrangements differed from respondents’ reported current arrangements.
Perceived efficiency while working from home during the COVID-19 period was measured using the statement: “I was efficient working from home during COVID-19.” Responses were recorded on a five-point Likert scale ranging from strongly disagree to strongly agree. This item was treated as an ordinal measure of respondents’ self-reported work-from-home effectiveness during the pandemic period.
Two items assessed respondents’ perceptions of remote-work possibilities in relation to their professions. The first asked respondents to characterize the prevalence or practicality of remote work within their profession, with response options indicating that remote work was not practical, applicable to very few activities, applicable to approximately half of activities, prevalent, applicable to almost all activities, or unknown.
The second item asked respondents whether they believed their profession was likely to shift, either partially or completely, toward remote work in the future. Response options distinguished between professions perceived as highly unlikely to shift, possibly shifting, or most likely shifting toward remote work.
These items represent respondents’ perceptions of occupational remote-work feasibility and future occupational change. They do not constitute independently verified measures of whether respondents’ occupations can in fact be performed remotely.
Respondents rated the importance of physical interaction with colleagues on a five-point scale ranging from not at all important to extremely important. This item was included to examine whether perceived importance of face-to-face interaction was associated with preferred future work arrangement. The measure was treated as attitudinal rather than as a direct indicator of preference for in-person work, because respondents may value physical interaction while still preferring remote or hybrid arrangements.
The survey additionally collected ratings concerning residential location considerations, including proximity to the workplace, proximity to the central business district, access to amenities, proximity to quality schools, access to major roads, access to public transportation, sense of community, safety, housing affordability, preference for urban or rural lifestyles, spacious living environments, and career-growth opportunities.
These residential location-preference items were excluded from the primary analyses reported in the present manuscript. This decision was made to maintain a focused analytical scope centered on preferred future work arrangements and their sociodemographic, experiential, occupational-perception, and attitudinal correlates. The residential location-preference variables may be examined separately in future research focused specifically on spatial or housing implications of remote-work preferences.
3.4. Variables and Analytical Coding
The primary outcome variable was preferred future work arrangement, measured using four ordered response categories: (1) fully in-person; (2) hybrid, dominated by in-person work; (3) hybrid, dominated by remote work; and (4) fully remote. For analyses assessing monotonic relationships, the variable was coded from 1 to 4, with higher values representing preference for increasingly remote work arrangements. For contingency-table analyses, the four categories were treated as categorical response levels without assuming equal intervals between categories.
The principal explanatory variables were organized into three domains: sociodemographic characteristics, experiential and occupational-perception variables, and an attitudinal variable.
The sociodemographic variables included gender, age group, educational attainment, monthly earnings, marital status, and housing tenure. Age group was treated as an ordered categorical variable using the six adult age categories recorded in the final dataset. Educational attainment was retained using the eight reported educational categories. Monthly earnings were retained using three substantive categories: $1,250 or less, $1,251–$3,333, and greater than $3,333. Housing tenure was represented as owner or renter.
For descriptive reporting, responses such as prefer not to say were retained and reported where present. For inferential analyses, such responses were treated as missing because they do not represent a substantive level of the associated explanatory variable. Missing values were not imputed.
Perceived work-from-home efficiency during COVID-19 was analyzed as an ordinal variable using its five-point Likert response structure, with higher values indicating greater self-reported efficiency while working from home.
Perceived prevalence or practicality of remote work within the respondent’s profession was treated as a categorical occupational-perception variable. Perceived likelihood that the respondent’s profession would shift toward remote work in the future was also treated as a categorical occupational-perception variable. These variables were interpreted as perceived professional context rather than objectively verified occupational characteristics.
Importance of physical interaction with colleagues was analyzed as an ordinal variable using its five-point response structure. Higher values indicate greater stated importance of physical interaction. No directional assumption was imposed in advance regarding whether higher importance of physical interaction would correspond to more in-person or more remote preferences.
Current work arrangement was measured using the same four categories as preferred future work arrangement. It was used only for paired comparison between current and preferred work arrangements among respondents who provided valid responses to both items. Because current work arrangement contained substantial item-level missing data, the paired analysis was conducted on a complete-case subset rather than on the full primary analytical sample.
Household size, desired family size, and residential location-preference items were retained in the cleaned dataset but excluded from the primary analyses reported in this manuscript. Their exclusion prevents expansion of the paper beyond its central purpose and reduces the risk of introducing multiple exploratory analyses not directly aligned with the primary study aims.
3.5. Statistical Analysis
Statistical analyses were conducted in Python. All statistical tests were two-sided, and statistical significance was evaluated at an alpha level of 0.05. Because the study uses cross-sectional data from a nonprobability survey sample, the results are interpreted as associations within the analytical sample rather than as causal effects or population-level predictions.
Frequency distributions and percentages were calculated for the demographic characteristics, current work arrangements, preferred future work arrangements, pandemic-era work-from-home efficiency, occupational-perception items, and importance of physical interaction. Because missing responses were not imputed, percentages for each item were calculated using valid responses to that item as the denominator. Missing responses were reported separately.
Chi-square tests of independence were used to examine associations between preferred future work arrangement and each of the following categorical sociodemographic variables: gender, age group, educational attainment, monthly earnings, marital status, and housing tenure.
Cramér’s V was calculated for each chi-square test as a measure of association strength. Effect-size values were interpreted cautiously and were not described as predictive or causal effects. Expected cell frequencies were inspected for all chi-square analyses. Where sparse response categories produced inadequate expected frequencies, category consolidation or sensitivity analysis was undertaken only where theoretically defensible and was reported transparently.
For inferential analyses involving monthly earnings and marital status, prefer not to say responses were treated as missing and excluded from the corresponding analysis.
Current and preferred work arrangements were measured using the same four ordered response categories. Because these paired outcomes comprise more than two categories, the standard McNemar test was not used.
The Stuart–Maxwell test of marginal homogeneity was used as the primary paired analysis to assess whether the distribution of preferred future work arrangements differed significantly from the distribution of current work arrangements among respondents with valid responses to both items. The Bowker test of symmetry was used as a supplementary paired analysis to examine whether movement between work-arrangement categories was directionally symmetric.
For interpretive purposes, respondent transitions were additionally summarized descriptively as:
These descriptive transition categories were used to clarify the direction of observed differences and were interpreted in conjunction with the paired categorical test results.
Spearman rank correlations were used to examine monotonic associations between preferred future work arrangement and two ordinal measures: perceived work-from-home efficiency during COVID-19 and importance of physical interaction with colleagues. Because preferred future work arrangement was coded from fully in-person to fully remote, a positive Spearman coefficient indicates that higher values on the explanatory variable are associated with preference for increasingly remote work arrangements. The direction of each association was checked against group-level descriptive statistics before interpretation.
Kruskal–Wallis H tests were used as supplementary analyses to examine whether distributions of perceived work-from-home efficiency and importance of physical interaction differed across the four preferred future work-arrangement groups. When a supplementary omnibus comparison was statistically significant and pairwise interpretation was needed, post-hoc comparisons were conducted with correction for multiple comparisons within that exploratory comparison family.
Chi-square tests of independence were used to examine associations between preferred future work arrangement and the two categorical occupational-perception variables: perceived prevalence or practicality of remote work within the respondent’s profession and perceived likelihood that the respondent’s profession would shift toward remote work in the future. Cramér’s V was reported for these associations.
An adjusted multivariable analysis was conducted to examine whether differences in preferred future work arrangement remained evident when theoretically relevant sociodemographic, experiential, occupational-perception, and attitudinal variables were considered jointly. Because preferred future work arrangement was measured using four ordered categories ranging from fully in-person to fully remote, an ordinal logistic regression model was initially evaluated.
The adjusted analysis included gender, age group, educational attainment, monthly earnings, housing tenure, marital-status grouping, perceived work-from-home efficiency during COVID-19, perceived practicality or prevalence of remote work within the respondent’s profession, perceived likelihood of future professional shift toward remote work, and importance of physical interaction with colleagues. Variables were included on theoretical grounds rather than selected solely according to significance in the bivariate analyses.
The proportional-odds assumption required for interpretation of the ordinal logistic regression model was evaluated using a global parallel-lines diagnostic. The assumption was not supported, LR χ² = 49.950, df = 20, p < 0.001. Consequently, a multinomial logistic regression model was estimated as the final adjusted robustness analysis, with fully in-person work specified as the reference outcome category.
The adjusted model was estimated using complete cases for all included variables (n = 431). For the adjusted model only, marital status was consolidated into married versus not married because modelling the original marital-status categories separately resulted in sparse-cell instability across the work-arrangement outcome categories. The original marital-status categories were retained in the primary descriptive and chi-square analyses. Responses of prefer not to say for monthly earnings and marital status, as well as I do not know responses for perceived professional remote-work practicality, were treated as missing in the adjusted model because they do not represent substantive ordered or categorical levels for model interpretation.
For ordered explanatory variables, relative risk ratios were interpreted per one-category increase in the corresponding measure. Model reporting includes the complete-case analytical sample size, outcome reference category, relative risk ratios, 95% confidence intervals, p-values, overall likelihood-ratio model test, and McFadden pseudo-R². Because this model was conducted as an adjusted robustness analysis, it was not included in the Benjamini–Hochberg correction applied to the predefined primary bivariate test family.
The predefined primary inferential family comprised eleven tests aligned with the study aims: six sociodemographic chi-square tests, one Stuart–Maxwell paired test, one Spearman correlation for perceived work-from-home efficiency, two occupational-perception chi-square tests, and one Spearman correlation for physical-interaction importance.
Because these primary analyses evaluate multiple related hypotheses concerning preferred future work arrangement, the Benjamini–Hochberg false discovery rate procedure was applied to the eleven-test primary family at
q = 0.05 [
32]. For each primary analysis, the Results section reports the valid sample size, test statistic, degrees of freedom where applicable, unadjusted p-value, Benjamini–Hochberg adjusted p-value, and relevant effect-size measure.
Cramér’s V was reported as the effect-size measure for chi-square tests of independence. Spearman’s ρ was reported for monotonic associations involving ordinal variables. For the paired current-versus-preferred work-arrangement analysis, the Stuart–Maxwell statistic was reported as the primary test of marginal homogeneity.
The Bowker test of symmetry, Kruskal–Wallis H tests, and any subsequent pairwise comparisons were treated as supplementary analyses rather than as additional primary hypothesis tests. Their purpose was to support interpretation of the primary results rather than to expand the primary test family. When pairwise post-hoc comparisons were undertaken following a supplementary omnibus test, their p-values were adjusted within the corresponding comparison family.
The adjusted multivariable model was treated as a robustness analysis and was not included in the Benjamini–Hochberg correction applied to the primary bivariate test family. Its purpose was to assess whether the observed associations remained evident after relevant variables were considered jointly.
For all analyses, item-level missing responses were excluded only from analyses requiring the corresponding variable, and no missing values were imputed. Consequently, valid sample sizes vary between analyses and are reported explicitly in the Results section.
3.6. Ethical Considerations
This study used de-identified survey data collected as part of the first author’s doctoral research at Virginia Tech. The study protocol was reviewed by the Virginia Tech Institutional Review Board and was determined to be exempt under IRB protocol #22-120. Participation in the survey was voluntary, and respondents were informed that they could withdraw from the survey at any time. Survey responses were treated confidentially and reported only in aggregate form. No personally identifiable information was included in the analysis. Responses from individuals younger than 18 years, cases without verifiable adult eligibility, and cases without verified residence in the Roanoke Metropolitan Statistical Area were excluded from the final analytical sample.