2.3. Participants and Procedure
The target population comprised adult residents of India with sufficient English proficiency to complete an online questionnaire. A non-probability convenience sampling strategy, supplemented by elements of snowball sampling, was adopted. The questionnaire was distributed electronically in December 2025 through social media platforms — including Instagram, WhatsApp, LinkedIn, and Facebook — in addition to personalised email invitations. Recipients were invited to forward the questionnaire to other potential respondents. No restrictions were imposed on age, gender, dietary pattern, or geographic location within India.
Prior to the commencement of the questionnaire, respondents were informed of the purpose of the study, the voluntary and anonymous nature of participation, and their right to withdraw at any time without consequence. Proceeding to the first question constituted informed consent.
A total of 259 responses were obtained. Following the application of listwise deletion to address missing values on the key analytic variables, the final analytic sample comprised
N = 255 complete cases. The level of missingness was minimal and non-systematic (< 1% per variable); the missing values occurred exclusively on the four single-item perceptual measures (perceived unnaturalness, disgust, ethical concern, and health-related concern), and each of the four excluded cases omitted at least one of these items. The achieved sample size satisfied both the absolute minimum of 200 observations and the 10:1 ratio criterion conventionally recommended for structural equation modelling with maximum likelihood estimation [
22].
2.4. Measures
The final instrument comprised eight sections. Section A consisted of a single screening item on prior awareness of cultured meat (binary: Yes / No). Section B presented a neutrally worded informational note describing the production process and characteristics of cultured meat, in line with the informed-choice paradigm widely adopted in this literature [
3]; the framing of such information has itself been shown to affect cultured-meat acceptance [
23]. Sections C through F contained the core TPB constructs. Attitude toward trying cultured meat was measured by four semantic differential items rated on 7-point scales (negative–positive, unpleasant–pleasant, unattractive–attractive, not worthwhile–worthwhile), following Ajzen [
24]; all four items were oriented so that higher scores denote a more positive attitude, and no reverse-coding was required. Subjective norms were measured by two 7-point Likert items (1 = strongly disagree, 7 = strongly agree) capturing the injunctive component of perceived social approval. Perceived behavioural control was measured by two 7-point Likert items assessing perceived ease of access and personal autonomy, in line with Ajzen’s [
25] treatment of perceived behavioural control and self-efficacy. Behavioural intention was measured by two 7-point Likert items addressing general intention and intention to purchase contingent on availability and reasonable price.
Section G assessed four psychological antecedents of attitude through single-item 7-point Likert measures: perceived unnaturalness, disgust, ethical concern for animal welfare, and health-related concern. The use of single-item measures for these constructs is consistent with established practice in cultured meat acceptance research, where each construct is theoretically concrete, unambiguous, and narrow in scope [
8,
9]. Section H collected demographic and background information, including age, gender, dietary pattern (omnivore, flexitarian, vegetarian, vegan), country of residence, education level, place of residence, and self-rated financial situation. The dietary-pattern item presented the general label “vegetarian” without distinguishing lacto- from lacto-ovo-vegetarian diets; this lack of granularity, which carries particular cultural significance in the Indian context, is acknowledged as a limitation.
The items were adapted from established TPB instruments for cultured meat and novel food acceptance research [
3,
7,
8,
9]. The adapted items are reproduced in line with standard scholarly use; the full instrument is available as
Supplementary Materials.
2.5. Statistical Analysis
All statistical analyses were conducted in Python (version 3.13), using the scipy (1.14), statsmodels (0.14), semopy (2.4), numpy (2.2), and pandas (2.2) libraries. The significance threshold was set at α = 0.05 for all inferential tests.
Composite scores for the four multi-item TPB constructs were computed as the arithmetic mean of the constituent items. Internal consistency was assessed by Cronbach’s α and, given the documented limitations of α under unequal factor loadings, by McDonald’s ω estimated from confirmatory factor analysis loadings [
26]. Descriptive statistics were computed for all study variables, and distributional properties were evaluated using the Shapiro-Wilk test. Pearson and Spearman correlations were computed in parallel in order to assess the robustness of bivariate associations under potential non-normality and the ordinal scaling of Likert-type items. Variance inflation factors were inspected to evaluate multicollinearity, with the conventional threshold of
VIF < 5 adopted [
27].
The hypothesised model was tested using latent-variable structural equation modelling with maximum likelihood estimation. The measurement model specified four latent constructs (Attitude, Subjective Norms, Perceived Behavioural Control, Intention), each indicated by its respective items, while the four psychological antecedents of attitude were retained as observed exogenous predictors. Global model fit was assessed using the χ² statistic, the Comparative Fit Index (
CFI), the Tucker-Lewis Index (
TLI), the Root Mean Square Error of Approximation (
RMSEA), and the Standardised Root Mean Square Residual (
SRMR), with the conventional thresholds recommended by Bentler and Bonett [
28], Hu and Bentler [
29], and Kline [
22] adopted.
In view of the single-source self-report design, common method variance was assessed as a potential threat to construct validity [
30]. The Unmeasured Latent Method Construct (ULMC) approach was employed as a post hoc robustness check [
31]. A latent method factor was added to the baseline model and constrained to load on all multi-item indicators while remaining orthogonal to the substantive constructs. Indirect effects of the four perceptual predictors on intention through attitude were estimated using the percentile bootstrap method with 5,000 resamples; an indirect effect was deemed significant when its 95% confidence interval excluded zero [
32,
33].
Group differences in attitude and intention across dietary patterns (H7) were tested using a multi-method approach. A multivariate analysis of variance (MANOVA) with Pillai’s Trace was used to test the joint vector of attitude and intention, selected for its robustness to violations of multivariate normality and unequal group sizes [
34]. One-way analysis of variance with Tukey HSD post hoc comparisons was conducted separately for attitude and intention, with effect sizes reported as η², ω², and Cohen’s f [
35]. Non-parametric Kruskal-Wallis tests, including all four dietary groups, served as sensitivity checks. The Vegan subgroup (
n = 4) was excluded from the parametric tests owing to insufficient size but was retained in the non-parametric analyses.
2.6. Robustness and Triangulation Analyses
The primary confirmatory analysis was the latent structural equation model with Holm-Bonferroni multiplicity correction and the ULMC robustness check (
Section 2.5). In order to evaluate the robustness of those primary results and to triangulate the findings across statistical paradigms, a pre-specified set of secondary analyses was undertaken. These secondary analyses were planned in advance of estimation of the main structural model but were not formally preregistered. All analyses were conducted on the same analytic sample (
N = 255) and were implemented in Python (version 3.13). The complete reproducibility package, including data, scripts, and a pinned requirements file, is available as
Supplementary Materials.
Confirmatory versus exploratory tests were pre-classified. Hypotheses H1 to H6 and H7 were treated as confirmatory; the exploratory hypothesis H8 (health-related concern) was analysed separately and is reported with this label throughout.
Item-level psychometric diagnostics. For each multi-item scale, Cronbach’s α-if-item-deleted and the corrected item-total correlations were computed in addition to the scale-level α and McDonald’s ω reported above.
Multiplicity correction. Within the confirmatory family H1–H6, raw p-values from the latent structural paths were adjusted using the Holm–Bonferroni procedure at a familywise α of 0.05. The Benjamini–Hochberg false discovery rate procedure (q = 0.05) was applied to the exploratory hypothesis H8.
Bootstrap robust standard errors. The latent structural model was re-estimated in 1,000 non-parametric bootstrap resamples (resampled with replacement; case-level resampling). For each structural path the bootstrap mean, bootstrap standard error, percentile and bias-corrected 95% confidence intervals, and a two-sided bootstrap p-value were computed. Replicates yielding standardised coefficients with absolute value above 1.5 were treated as numerical failures (arising from near-singular factor loadings in the resampled covariance structure) and were excluded from the summary. The number of stable replicates is reported for every path.
Sensitivity power analysis. For every structural path, a post-hoc Wald-type sensitivity power analysis was conducted using the bootstrap-derived standard error as a realistic estimate of the path uncertainty. Achieved power was computed at the uncorrected α of 0.05 and at the most conservative Holm-corrected α of 0.05 / 6. The minimum detectable standardised effect at power = 0.80 was reported for each path under both α thresholds. This analysis was used to distinguish substantive null findings from findings whose non-significance may reflect inadequate power.
Binary-outcome sensitivity. The continuous intention measure was dichotomised at the scale midpoint (intention > 4 indicating willingness to try cultured meat). A logistic regression with heteroskedasticity-consistent (HC3) robust standard errors was estimated to test the core TPB predictors against the dichotomised outcome. The same procedure was applied to a dichotomised attitude outcome for the H4–H6 / H8 sensitivity model. Odds ratios with 95% confidence intervals, McFadden’s pseudo-R², and the area under the receiver operating characteristic curve (AUC) are reported.
Dominance analysis. The relative importance of the predictors in two ordinary least-squares models — intention regressed on the three core TPB constructs, and attitude regressed on the four perceptual predictors — was quantified by Budescu’s general-dominance decomposition [
36]. For every predictor the average incremental
R² across all 2ᵖ⁻¹ subsets of co-predictors was computed and normalised to the full-model
R². This procedure provides an order-independent decomposition of explained variance and is reported alongside the structural-model coefficients to confirm the relative-importance ranking.
Triangulation across estimation paradigms. Three independent paradigms were used to confirm the patterns identified by covariance-based SEM. First, a Bayesian path model with weakly informative priors (b ~ Normal(0, 1) on standardised coefficients) was estimated in PyMC with four chains of 2,000 tuning and 2,000 sampling iterations. Posterior 95% credible intervals, posterior probability of direction, and Savage–Dickey approximations of the Bayes factor BF₁₀ against a point-null at zero are reported. Second, an elastic net regression (l1-ratio and penalty tuned by 10-fold cross-validation) and a random-forest regressor with permutation importance were estimated for both outcomes (intention and attitude). Predictors retained by the elastic net and ranked highly by random-forest permutation importance were compared with the set of predictors identified by the structural model. Third, the same path model was re-estimated under the variance-based partial least squares structural equation modelling (PLS-SEM) paradigm using a transparent in-house implementation of the canonical Wold–Tenenhaus algorithm with Mode A outer estimation and the path-weighting inner scheme. Five hundred non-parametric bootstrap resamples were used for inference on the PLS path coefficients. Convergence of the path estimates across the covariance-based SEM, the Bayesian path model, the regularised and non-parametric machine-learning analyses, and the PLS-SEM analysis was treated as triangulation evidence for the structural conclusions.
Two further sensitivity analyses were conducted as additional robustness checks on the Perceived Behavioural Control findings. To probe the borderline reliability of the two-item Perceived Behavioural Control scale, the latent structural model was re-estimated twice, with Perceived Behavioural Control represented by each of its two indicators alone (single-indicator specifications, with the item treated as a perfect measure of the construct). To address the high collinearity between Subjective Norms and Perceived Behavioural Control, an alternative hierarchical model specifying Subjective Norms as an antecedent of Perceived Behavioural Control (Subjective Norms → Perceived Behavioural Control → Intention) was estimated and compared against the parallel baseline model using the Akaike and Bayesian information criteria.
Generative AI statement. The authors used a generative AI assistant (Claude, Anthropic) to assist with the restructuring and formatting of the Python analysis code and with the code used to generate the manuscript figures. The study design, data collection, statistical analyses, results, and their interpretation were carried out by the authors; all AI-assisted code was independently executed and verified against the reported outputs, and the authors take full responsibility for the content.