Parameter Estimation Research Articles

Knowing what to know: Implications of the choice of prior distribution on the behavior of adaptive design optimization.

Adaptive design optimization (ADO) is a state-of-the-art technique for experimental design (Cavagnaro et al., 2010). ADO dynamically identifies stimuli that, in expectation, yield the most information about a hypothetical construct of interest (e.g., parameters of a cognitive model). To calculate this expectation, ADO leverages the modeler's existing knowledge, specified in the form of a prior distribution. Informative priors align with the distribution of the focal construct in the participant population. This alignment is assumed by ADO's internal assessment of expected information gain. If the prior is instead misinformative, i.e., does not align with the participant population, ADO's estimates of expected information gain could be inaccurate. In many cases, the true distribution that characterizes the participant population is unknown, and experimenters rely on heuristics in their choice of prior and without an understanding of how this choice affects ADO's behavior. Our work introduces a mathematical framework that facilitates investigation of the consequences of the choice of prior distribution on the efficiency of experiments designed using ADO. Through theoretical and empirical results, we show that, in the context of prior misinformation, measures of expected information gain are distinct from the correctness of the corresponding inference. Through a series of simulation experiments, we show that, in the case of parameter estimation, ADO nevertheless outperforms other design methods. Conversely, in the case of model selection, misinformative priors can lead inference to favor the wrong model, and rather than mitigating this pitfall, ADO exacerbates it.

Knowing what to know: Implications of the choice of prior distribution on the behavior of adaptive design optimization.

Adaptive design optimization (ADO) is a state-of-the-art technique for experimental design (Cavagnaro et al., 2010). ADO dynamically identifies stimuli that, in expectation, yield the most information about a hypothetical construct of interest (e.g., parameters of a cognitive model). To calculate this expectation, ADO leverages the modeler's existing knowledge, specified in the form of a prior distribution. Informative priors align with the distribution of the focal construct in the participant population. This alignment is assumed by ADO's internal assessment of expected information gain. If the prior is instead misinformative, i.e., does not align with the participant population, ADO's estimates of expected information gain could be inaccurate. In many cases, the true distribution that characterizes the participant population is unknown, and experimenters rely on heuristics in their choice of prior and without an understanding of how this choice affects ADO's behavior. Our work introduces a mathematical framework that facilitates investigation of the consequences of the choice of prior distribution on the efficiency of experiments designed using ADO. Through theoretical and empirical results, we show that, in the context of prior misinformation, measures of expected information gain are distinct from the correctness of the corresponding inference. Through a series of simulation experiments, we show that, in the case of parameter estimation, ADO nevertheless outperforms other design methods. Conversely, in the case of model selection, misinformative priors can lead inference to favor the wrong model, and rather than mitigating this pitfall, ADO exacerbates it.

The effect of surface on knee landing mechanics and muscle activity during a single-leg landing task in recreationally active females

ObjectiveInvestigate the effect of surface on frontal plane knee angle, knee moment and muscle activity. DesignRandomised cross over. SettingUniversity Laboratory. MethodsTwenty females performed single-leg hop-landings onto sand, grass and firm surfaces. Kinematic, kinetic and muscle activity data were obtained. Compatibility curves were used to visualise parameter estimates alongside P- values, and S-value transforms. ResultsKnee angle for firm-sand (mean difference (d)‾ = −2.2°; 95% compatibility interval (CI): −4.6 to 0.28, p = 0.083, s = 3.6) and firm-grass (d‾ = −1.9; 95% CI: −4.3 to 0.5, p = 0.125, S = 3) yielded <4 bits of reputational information against the null hypothesis (H). 5 bits (p = 0.025) of information against H were observed for knee moment between firm-sand (d‾ = 0.17 N m/kg-1. m-1; 95% CI: 0.02 to 0.31) with similar effects for firm-grass (d‾ = 0.14 N m/kg-1. m-1; 95% CI: −0.02 to 0.29, p = 0.055, S = 4). Muscle activity across surfaces ranged from almost no (S = 1) reputational evidence against H (Quadriceps and Hamstrings) to 10–13 ‘bits’ against H for lateral gastrocnemius (lower on sand). ConclusionsOur study provides valuable information for practitioners of the observed effect sizes for lower-limb landing mechanics across surfaces in asymptomatic females.

Uncertain statistics models for characterizing annual total energy production

In the field of annual total energy production, it is hard to determine whether the experimental data is linear or nonlinear growth characteristics. Therefore, different regression models are used to test the properties of the observed data. This article uses actual data to perform parameter estimation and residual analysis of linear and nonlinear regression models. In addition, uncertain nonlinear time series model is used to fit the experimental data in order to find a more suitable statistics model after obtaining the growth properties of the observed data. A test is introduced to prove that the disturbance terms presented in the regression and time series models are actually uncertain variables, not random variables. Moreover, the uncertain hypothesis tests show that the uncertain time series model is in better agreement with the real data of the annual total energy production.

Enhanced Remote Sensing Image Compression Method Using Large Network with Sparse Extracting Strategy

Deep neural networks based on hyper-encoders play a critical role in estimating prior distributions in remote sensing image compression issues. However, most of the existing encoding methods suffer from a problem on the hyper-encoding side, namely the mismatch of extraction ability with the encoder. This ability bias results in likelihood features that fail to extract sufficient information from latent representations. To solve this problem, the feature extraction capabilities of the hyper-encoder are enhanced to better estimate the Gaussian likelihood of the latent representation in end-to-end network optimization. Specifically, residual blocks and a parameter estimation module are incorporated to balance the performance of the encoder and the hyper-encoder. Furthermore, it is observed that the well-trained compression model tends to generate a fixed pattern of latent representations. Therefore, we incorporate a nonlocal cross-channel graph (NCG) on the backside of the encoder. Specifically, it aggregates features between similar latent representations in a graphical manner to further enhance the side information extraction capability of the hyper-encoder. Considering the computational cost, a sparse graph strategy is further developed to dynamically select the most relevant latent representations for aggregation operations, which greatly reduces the computational effort. The proposed algorithm is named nonlocal cross-channel efficient graph (NCEG). A long-dependent residual network is selected as the backbone, and a sparse attention module is inserted into the encoder/decoder side to enhance the perceptual field of the network. The experimental results on two evaluation datasets demonstrate that the proposed method achieves satisfactory results compared to other learning-based methods.

Robust Design Based on Cost-Quality Model in Micro-Manufacturing

This paper proposes a novel total cost model for the micro‐products' entire life cycle that takes into account the uncertainty of the model parameters. The total cost includes pre-sale manufacturing and post-sale warranty costs. Additionally, different marketing strategies are also given based on the weight of internal and external costs. Furthermore, limited data and unknown effects in experiments may cause large errors in parameter estimates. This could prevent the achievement of reliable designs. To address this, robust optimization and interval estimation are used. This approach reduces the impact of uncertainty on parameter estimates. It ensures optimality and robustness in micro-manufacturing parameters. Example analysis and numerical simulation results show that the proposed method assists companies in selecting the optimal manufacturing parameter level that aligns with their marketing strategies. Besides, considering uncertainty factors can ensure that the optimization results remain guaranteed, even under the worst-case scenarios.

Different Statistical Inference Algorithms for the New Pareto Distribution Based on Type-II Progressively Censored Competing Risk Data with Applications

In this research, the statistical inference of unknown lifetime parameters is proposed in the presence of independent competing risks using a progressive Type-II censored dataset. The lifetime distribution associated with a failure mode is assumed to follow the new Pareto distribution, with consideration given to two distinct competing failure reasons. Maximum likelihood estimators (MLEs) for the unknown model parameters, as well as reliability and hazard functions, are derived, noting that they are not expressible in closed form. The Newton–Raphson, expectation maximization (EM), and stochastic expectation maximization (SEM) methods are employed to generate maximum likelihood (ML) estimations. Approximate confidence intervals for the unknown parameters, reliability, and hazard rate functions are constructed using the normal approximation of the MLEs and the normal approximation of the log-transformed MLEs. Additionally, the missing information principle is utilized to derive the closed form of the Fisher information matrix, which, in turn, is used with the delta approach to calculate confidence intervals for reliability and hazards. Bayes estimators are derived under both symmetric and asymmetric loss functions, with informative and non-informative priors considered, including independent gamma distributions for informative priors. The Monte Carlo Markov Chain sampling approach is employed to obtain the highest posterior density credible intervals and Bayesian point estimates for unknown parameters and reliability characteristics. A Monte Carlo simulation is conducted to assess the effectiveness of the proposed techniques, with the performances of the Bayes and maximum likelihood estimations examined using average values and mean squared errors as benchmarks. Interval estimations are compared in terms of average lengths and coverage probabilities. Real datasets are considered and examined for each topic to provide illustrative examples.

Different Statistical Inference Algorithms for the New Pareto Distribution Based on Type-II Progressively Censored Competing Risk Data with Applications

In this research, the statistical inference of unknown lifetime parameters is proposed in the presence of independent competing risks using a progressive Type-II censored dataset. The lifetime distribution associated with a failure mode is assumed to follow the new Pareto distribution, with consideration given to two distinct competing failure reasons. Maximum likelihood estimators (MLEs) for the unknown model parameters, as well as reliability and hazard functions, are derived, noting that they are not expressible in closed form. The Newton–Raphson, expectation maximization (EM), and stochastic expectation maximization (SEM) methods are employed to generate maximum likelihood (ML) estimations. Approximate confidence intervals for the unknown parameters, reliability, and hazard rate functions are constructed using the normal approximation of the MLEs and the normal approximation of the log-transformed MLEs. Additionally, the missing information principle is utilized to derive the closed form of the Fisher information matrix, which, in turn, is used with the delta approach to calculate confidence intervals for reliability and hazards. Bayes estimators are derived under both symmetric and asymmetric loss functions, with informative and non-informative priors considered, including independent gamma distributions for informative priors. The Monte Carlo Markov Chain sampling approach is employed to obtain the highest posterior density credible intervals and Bayesian point estimates for unknown parameters and reliability characteristics. A Monte Carlo simulation is conducted to assess the effectiveness of the proposed techniques, with the performances of the Bayes and maximum likelihood estimations examined using average values and mean squared errors as benchmarks. Interval estimations are compared in terms of average lengths and coverage probabilities. Real datasets are considered and examined for each topic to provide illustrative examples.

Modeling and Optimal Control of Infectious Diseases

We propose a stochastic model of infectious disease transmission that is more realistic than those found in the literature. The model is based on jump-diffusion processes. However, it is defined in such a way that the number of people susceptible to be infected decreases over time, which is the case for a population of fixed size. Next, we consider the problem of finding the optimal control of the proposed model. The dynamic programming equation satisfied by the value function is derived. Estimators of the various model parameters are obtained.

Modeling and Optimal Control of Infectious Diseases

We propose a stochastic model of infectious disease transmission that is more realistic than those found in the literature. The model is based on jump-diffusion processes. However, it is defined in such a way that the number of people susceptible to be infected decreases over time, which is the case for a population of fixed size. Next, we consider the problem of finding the optimal control of the proposed model. The dynamic programming equation satisfied by the value function is derived. Estimators of the various model parameters are obtained.

Adaptive field-oriented control of current-fed induction motors: a discrete-time implementation

Indirect field oriented control (IFOC) is the industry standard for high performance current-fed induction motors. It is well-known that this controller preserves stability in the face of large rotor resistance variations – which are inevitable in practical scenarios – but its performance is significantly degraded. To overcome this drawback many adaptation schemes have been reported, both in drives and control-oriented publications, with the overwhelming majority of these developments done in continuous-time. This, in spite of the fact that the IFOC, as well as their proposed adaptive versions, involve the solution of highly nonlinear differential equations, whose discretization is far from obvious. In this paper we try to reverse this trend and present a discrete-time estimator of the motor parameters and an adaptive IFOC based on it. It is shown that the proposed estimator and adaptive controller are globally convergent under very weak prior knowledge and excitation assumptions. Simulations, that validate the theoretical claims and illustrate the robustness of the schemes, are also presented.

A Study on the Object-Based High-Resolution Remote Sensing Image Classification of Crop Planting Structures in the Loess Plateau of Eastern Gansu Province

The timely and accurate acquisition of information on the distribution of the crop planting structure in the Loess Plateau of eastern Gansu Province, one of the most important agricultural areas in Western China, is crucial for promoting fine management of agriculture and ensuring food security. This study uses multi-temporal high-resolution remote sensing images to determine optimal segmentation scales for various crops, employing the estimation of scale parameter 2 (ESP2) tool and the Ratio of Mean Absolute Deviation to Standard Deviation (RMAS) model. The Canny edge detection algorithm is then applied for multi-scale image segmentation. By incorporating crop phenological factors and using the L1-regularized logistic regression model, we optimized 39 spatial feature factors—including spectral, textural, geometric, and index features. Within a multi-level classification framework, the Random Forest (RF) classifier and Convolutional Neural Network (CNN) model are used to classify the cropping patterns in four test areas based on the multi-scale segmented images. The results indicate that integrating the Canny edge detection algorithm with the optimal segmentation scales calculated using the ESP2 tool and RMAS model produces crop parcels with more complete boundaries and better separability. Additionally, optimizing spatial features using the L1-regularized logistic regression model, combined with phenological information, enhances classification accuracy. Within the OBIC framework, the RF classifier achieves higher accuracy in classifying cropping patterns. The overall classification accuracies for the four test areas are 91.93%, 94.92%, 89.37%, and 90.68%, respectively. This paper introduced crop phenological factors, effectively improving the extraction precision of the shattered agricultural planting structure in the Loess Plateau of eastern Gansu Province. Its findings have important application value in crop monitoring, management, food security and other related fields.

A noncommutative Bianchi I model with radiation

In this work, we study the dynamical evolution of a homogeneous and anisotropic, noncommutative (NC) Bianchi I (BI) model coupled to a radiation perfect fluid. Our first motivation is determining if the present model tends to a homogeneous and isotropic NC Friedmann–Robertson–Walker (FRW) model, during its evolution. In order to simplify our task, we use the Misner parametrization of the BI metric. In terms of that parametrization the BI metric has three metric functions: the scale factor [Formula: see text] and the two parameters [Formula: see text], which measure the spatial anisotropy of the model. Our second motivation is trying to describe the present accelerated expansion of the universe using noncommutativity (NCTY). The NCTY is introduced by two nontrivial Poisson brackets between some geometrical as well as matter variables of the model. We recover the description in terms of commutative variables by introducing some variables transformations that depend on the NC parameter. Using those variables’ transformations, we rewrite the total NC Hamiltonian of the model in terms of commutative variables. From the resulting Hamiltonian, we obtain the dynamical equations for a generic perfect fluid. In order to solve these equations, we restrict our attention to a model in which the perfect fluid is radiation. We solve, numerically, these equations and compare the NC solutions to the corresponding commutative ones. The comparison shows that the NC model may be considered as a possible candidate for describing the accelerated expansion of the universe. Finally, we obtain estimates for the NC parameter and compare the main results of the NC BI model coupled to radiation with the same NC BI model coupled to other perfect fluids. As our main result, we show that the solutions, after some time, produce an isotropic universe. Based on that result, we can speculate that the solutions may represent an initial anisotropic stage of our Universe, that later, due to the expansion, became isotropic.

Population dynamics of spotted scat, Scatophagus argus (Linnaeus 1766) from the Sundarbans of Bangladesh

In Bangladesh, the mangrove areas of the Sundarbans are a hotspot for the spotted scat (Scatophagus argus), as the estuarine tributeries are significant zones for fishing, local marketing and foreign trading. Due to the widespread capture of spotted scats in the wild for a variety of uses, parameters of population were estimated using length-frequency and body weight data gathered from various sites in the southwestern forest areas of the country to comprehend the stock's status. Between May 2022 and April 2023, 1280 S. argus individuals were sampled from Khulna, Bagerhat and Satkhira. To gain insight into the population dynamics of S. argus, estimations of the growth parameters, recruitment pattern, mortalities, and level of exploitation (E) were made. For S. argus, a bit of negetive allomentric growth (b = 2.95) was observed. Moreover, the growth performance index (φ′), the asymptotic length (TL∝) and the growth coefficient (K) were 2.83, 34.65, cm and 0.56 yr−1, respectively. Natural mortality (1.08 yr−1) was lower than fishing mortality (1.12 yr−1). In the wild, this species recruits throughout the year, reaching its peak between October and January. By fishing gear, 50 % of the S. argus were caught between 7.71 cm and 9.43 cm. Moreover, high fishing mortality was assumed in both the small group (11–15 cm) and the large group (29–33 cm) of S. argus by virtual population analysis (VPA). It was discovered that the E for S. argus was 0.51, over the maximum permissible limit (Emax = 0.48). In summary, it is discovered that the S. argus population is overexploited. Therefore, in order to preserve this popular fish species and avoid local extinction, spotted scat resources should be taken into account under appropriate fisheries management guidelines. These findings will guide policymakers to manage the species accordingly.

Locally sparse and robust partial least squares in scalar-on-function regression

We present a novel approach for estimating a scalar-on-function regression model, leveraging a functional partial least squares methodology. Our proposed method involves computing the functional partial least squares components through sparse partial robust M regression, facilitating robust and locally sparse estimations of the regression coefficient function. This strategy delivers a robust decomposition for the functional predictor and regression coefficient functions. After the decomposition, model parameters are estimated using a weighted loss function, incorporating robustness through iterative reweighting of the partial least squares components. The robust decomposition feature of our proposed method enables the robust estimation of model parameters in the scalar-on-function regression model, ensuring reliable predictions in the presence of outliers and leverage points. Moreover, it accurately identifies zero and nonzero sub-regions where the slope function is estimated, even in the presence of outliers and leverage points. We assess our proposed method’s estimation and predictive performance through a series of Monte Carlo experiments and an empirical dataset—that is, data collected in relation to oriented strand board. Compared to existing methods our proposed method performs favorably. Notably, our robust procedure exhibits superior performance in the presence of outliers while maintaining competitiveness in their absence. Our method has been implemented in the robsfplsr package in .

The filtering-based recursive least squares identification and convergence analysis for nonlinear feedback control systems with coloured noises

The coloured noise is ubiquitous in industrial processes. This paper addresses the identification problem for the nonlinear feedback systems with coloured noise. Firstly, a direct identification scheme based on the least squares principle is developed to estimate the whole parameters of the nonlinear feedback systems and the convergence analysis is carried out through the stochastic stability theory. Secondly, for the purpose of improving the estimation accuracy, a filtering-based identification framework is proposed by constructing a linear filter for filtering the input data, output data and the coloured noise, and the coloured noise is transformed into a white noise. This identification scheme based on the filtering technique can effectively reduce the adverse effects caused by coloured noise and parameter estimation accuracy is enhanced compared with the direct least squares algorithm. Meanwhile, the convergence analysis of the filtering-based identification algorithm is given to provide a theoretical analysis. Finally, the simulation example is carried out by performance test and comparison analysis and simulation results show the effectiveness of the proposed identification methods.

Pharmacokinetic-pharmacodynamic modelling and simulation of methotrexate dosing in patients with rheumatoid arthritis.

To develop a non-linear mixed-effects population pharmacokinetic and pharmacodynamic (PK-PD) model describing the change in the concentration of methotrexate polyglutamates in erythrocytes (ery-MTX-PGn with "n" number of glutamate, representing PK component) and how this relates to modified 28-joint Disease Activity Score incorporating erythrocyte sedimentation rate (DAS-28-3) for rheumatoid arthritis (RA), representing PD component. An existing PK model was fitted to data from a study consisting of 117 RA patients. The estimation of population PK-PD parameters was performed using stochastic approximation expectation maximisation algorithm in Monolix 2021R2. The model was used to perform Monte Carlo simulations of a loading dose regimen (50mg subcutaneous methotrexate as loading doses, then 20mg weekly oral methotrexate) compared to a standard dosing regimen (10mg weekly oral methotrexate for 2 weeks, then 20mg weekly oral methotrexate). Every 40 nmol/L increase in ery-MTX-PG3-5 total concentration correlated with 1-unit reduction in DAS-28-3. Significant covariate effects on the therapeutic response of methotrexate included the use of prednisolone in the first 4 weeks (positive use correlated with 25% reduction in DAS-28-3 when other variables were constant) and patient age (every 10-year increase in age correlated with 3.4% increase in DAS-28-3 when other variables were constant). 4 methotrexate loading doses led to a higher percentage of patients achieving a good/moderate response compared to the standard regimen (Week 4: 87.6% vs. 39.8%; Week 10: 64.7% vs. 57.0%). A loading dose regimen was more likely to achieve higher ery-MTX-PG concentration and better therapeutic response after 4 weeks of methotrexate treatment.

Estimating error rates for single molecule protein sequencing experiments.

The practical application of new single molecule protein sequencing (SMPS) technologies requires accurate estimates of their associated sequencing error rates. Here, we describe the development and application of two distinct parameter estimation methods for analyzing SMPS reads produced by fluorosequencing. A Hidden Markov Model (HMM) based approach, extends whatprot, where we previously used HMMs for SMPS peptide-read matching. This extension offers a principled approach for estimating key parameters for fluorosequencing experiments, including missed amino acid cleavages, dye loss, and peptide detachment. Specifically, we adapted the Baum-Welch algorithm, a standard technique to estimate transition probabilities for an HMM using expectation maximization, but modified here to estimate a small number of parameter values directly rather than estimating every transition probability independently. We demonstrate a high degree of accuracy on simulated data, but on experimental datasets, we observed that the model needed to be augmented with an additional error type, N-terminal blocking. This, in combination with data pre-processing, results in reasonable parameterizations of experimental datasets that agree with controlled experimental perturbations. A second independent implementation using a hybrid of DIRECT and Powell's method to reduce the root mean squared error (RMSE) between simulations and the real dataset was also developed. We compare these methods on both simulated and real data, finding that our Baum-Welch based approach outperforms DIRECT and Powell's method by most, but not all, criteria. Although some discrepancies between the results exist, we also find that both approaches provide similar error rate estimates from experimental single molecule fluorosequencing datasets.

Design analysis and development of a gantry robot for multi-layer 3-D concrete printing with simulation and experimental validation

The present study investigates the design analysis and development of a framed gantry robot (FGR) for 3-D concrete printing. Finite element simulations are performed for design validation and parameter estimation to conduct detailed analysis of the concrete printer. The kinematic and dynamic models of FGR are used for workspace analysis, motion control study, and motor sizing for printer development. The motion trajectory of the printer nozzle tip is evaluated using MATLAB simulations. An integrated programmable logic controller (PLC) with an in-built human-machine interface (HMI) through control algorithms has been implemented for interactive, user-friendly, and precise motion control. A separate control system is developed to monitor the pumping power and feed rate of mortar mix during printing. The simulated control scheme is validated experimentally using trajectory tracking for various input printing paths through HMI, with a maximum tracking error of 1 mm. Successful concrete printing trials have been conducted using various mix-design, following extensive evaluation of different printer and material parameters at the laboratory scale. The proposed design concept, mechanisms, interactive control interface, and tracking accuracy have been validated through simulation and experimental results, for ensuring quality concrete prints by the developed printer. Material testing parameters have been evaluated experimentally for modeling various mix-designs. To further evaluate printing performance, simulations of layer deformation in printed structures are conducted using finite element analysis, incorporating non-linear material models and various contact layer interactions. The simulated results of printed layers of different shapes of structures show a close resemblance with the experimental concrete printing trials.

The impact of diffusion tensor imaging tractography settings on muscle fascicle architecture and diffusion parameter estimates: Tract length, completion, and curvature are most sensitive to tractography settings.

Diffusion-tensor (DT)-MRI tractography provides information about properties relevant to muscle health and function, including estimates of architectural properties such as fascicle length, pennation angle, and curvature and diffusion properties such as mean diffusivity (MD) and fractional anisotropy (FA). Tractography settings, including integration algorithms, thresholds for early tract termination, and tract smoothing approaches, impact the accuracy of the muscle property estimates. However, muscle DT-MRI tractography is performed using a variety of these settings, complicating comparisons between different studies. The effects of different tractography settings on muscle architecture estimates have not been fully explored, and optimized settings for muscle tractography have not yet been determined. We examined the influence of integration algorithm and termination check settings combined with a range of step sizes, termination criteria, and smoothing polynomial orders on tract characteristics, completion/reason for termination, and goodness of fit between fiber tracts and smoothing polynomials using 3-T DT-MR images of the lower leg muscles of seven healthy adults. We found that tract length and completion were highly sensitive to strict FA and intersegment angle thresholds (25%-69% reduction in complete fiber tracts from lowest to highest minimum FA threshold and 11%-36% reduction from highest to lowest intersegment angle threshold). Higher order polynomials (third and fourth order vs. second order) better fit the muscle fiber trajectories, but curvature estimates were highly sensitive to smoothing polynomial order (3.9-6.6 m-1 increase for second- vs. fourth-order fitting polynomials). Step size impacted curvature estimates, albeit to a lesser degree. Integration algorithm had little impact, and mean pennation angle, and tract-based FA and MD, were relatively insensitive to all parameters. The results demonstrate which muscle diffusion measures and architectural estimates are most sensitive to varying tractography settings and support the need for consistent reporting of tractography details to aid interpretation and comparison of results between studies.