REVIEW 5 major objections 5 minor 47 references
Analysis of in-vivo skin anisotropy using elastic wave measurements and Bayesian modelling
T0 review · 5 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper reports that in-vivo skin anisotropy increases logarithmically with age and that stiffness along Langer lines rises with age, based on elastic-wave measurements from 78 volunteers aged 3 to 93, and proposes ellipse eccentricity…
desk verdict A useful new anisotropy metric and a plausible age trend, but the functional-form claims and unquantified ellipse-fit failures need referee attention before the headline holds up. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is an ellipse fitted by direct least squares to polar-plot RRT data measured at 10-degree increments around the skin site. From the fitted ellipse, the semi-major axis a and semi-minor axis b give the eccentricity e = $\sqrt$(1 - $b^{2}$/$a^{2}$) as the anisotropy measure, the area A = pi*a*b as the average all-direction arrival time, and the semi-minor axis as the arrival time along the fastest direction, i.e. the Langer line. The conversion from arrival time to stiffness uses the Rayleigh-wave relation E = rho*$v^{2}$*(3.284). The statistical machinery is a Bayesian multivariate outcome regression on log eccentricity, area/1000, and semi-minor axis, with a correlated error term and Markov chain Monte Carlo inference, which accounts for the dependence among the ellipse-derived outcomes; a projected-normal circular regression tests whether stretching changes the ellipse angle.
What would settle it
Compute a per-subject ellipse goodness-of-fit measure, such as residual standard deviation or an asymmetry index, regress it on age, and re-run the Bayesian model excluding poor fits; if the 0.95 HPDI for the age coefficient on log eccentricity still excludes zero after exclusion, the age trend survives, but if the excluded cases cluster in children or the elderly, the reported logarithmic increase could be an artifact of fit quality.
Extended reading notes
Core claim
The discovery is that in-vivo human skin anisotropy, measured as the eccentricity of an ellipse fitted to angular Reviscometer arrival-time data, increases logarithmically with age: the posterior mean age coefficient on log eccentricity is 0.007 with 95% HPDI [0.005, 0.009]. In the same multivariate model, stiffness in the Langer-line direction, represented by the ellipse's semi-minor axis, increases with age (coefficient -0.391, HPDI [-0.696, -0.085]), while the ellipse area, representing average stiffness across directions, does not show a significant age effect. Gender does not significantly change anisotropy, but males show significantly stiffer skin both on average and along Langer lines. Applied stretch in the Langer-line direction significantly increases eccentricity and decreases the semi-minor axis while conserving the ellipse tilt angle, evidence that the level of skin tension strongly affects both measurements. The paper further claims that the eccentricity measure is more robust than the anisotropic ratio, because in simulations the ratio drifts upward and develops extreme outliers under noise while the eccentricity remains near its true value.
Load-bearing premise
The entire age trend rests on the assumption that an ellipse adequately fits every subject's angular arrival-time data, so that its eccentricity truly summarizes the anisotropy; the authors note that some subjects' raw data lacked the required symmetry and produced poor ellipse fits, but they do not quantify how many or how those cases were handled.
Editorial extensions
If this is right
- Surgical incision orientation and wound-closure planning could be tailored by age, since older skin is more direction-dependent and therefore more sensitive to the orientation of cuts.
- Elastic-wave measurements may become a practical way to infer in-vivo skin tension, because an applied stretch changes eccentricity and semi-minor axis in a consistent direction while conserving the Langer-line angle.
- Universal cosmetic surgery practices applied to very young or elderly patients are called into question, since the level of anisotropy varies strongly across the lifespan.
- The eccentricity metric, being less sensitive to outliers than the anisotropic ratio, could replace the ratio in future angular RRT studies.
- The results imply that directional stiffening and overall stiffness are partly decoupled: with age, stiffness rises only along Langer lines, not on average across directions.
Reading between the lines
- The authors do not test whether the age trend generalizes beyond the volar forearm site; a natural extension would be to repeat the protocol on additional body locations and see whether the logarithmic anisotropy trajectory is site-specific.
- The ellipse-fit assumption is the main unresolved risk; a systematic check of per-subject fit residuals by age group, not reported in the paper, would show whether the steep childhood-to-adulthood rise is partly an artifact of poorer fits in younger subjects.
- The calibration relation E = rho*v^2*(3.284) comes from linear-elastic, isotropic, unstressed materials; extending the wave-speed model to include pre-tension could convert the observed stretch effect into an actual numerical estimate of in-vivo skin tension, which the paper stops short of doing.
- The finding that gender affects stiffness but not anisotropy suggests that overall connective-tissue stiffness and directional organization are governed by at least partly separate biological mechanisms, so interventions targeting one may not alter the other.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript proposes using the eccentricity of a least-squares ellipse fitted to angular Reviscometer Resonance Running Time (RRT) measurements as a measure of in-vivo skin anisotropy, compares this metric with the classic anisotropy ratio through simulations, and then fits a Bayesian multivariate regression of log eccentricity, ellipse area, and semi-minor axis on age, gender, and stretch configuration for 78 subjects aged 3–93. The authors report that anisotropy increases with age, that stiffness increases along Langer lines with age, that gender does not affect anisotropy but affects stiffness, and that applied stretch changes both anisotropy and stiffness while preserving the direction of Langer lines.
Significance. The paper has two potentially useful contributions. First, the eccentricity metric is simple, and the simulation study (10,000 replicates per condition) provides credible evidence that it is more robust to noise and outliers than the max/min anisotropy ratio; this is a practical contribution for skin-mechanics measurement. Second, the dataset spans a wider age range than most comparable studies, and the code is publicly available, which makes the results reproducible. If the statistical issues listed below are fixed, the study could be a useful reference for age-related skin anisotropy; the Bayesian framing and explicit reporting of posterior intervals are strengths.
major comments (5)
- [Section 3.2, Eq. (5)] The statement that "the eccentricity increases logarithmically with age" is the opposite of what Eq. (5) implies. With log(Eccentricity) as the outcome and Age as a linear predictor, the model is log(e) = a + 0.007·Age + ..., so e increases exponentially with age (e ∝ exp(0.007·Age)), not logarithmically. The interpretation in the abstract, Section 3.2, and the Discussion should be corrected; if the authors intend to claim a logarithmic dependence, they should instead fit eccentricity against log(Age). The statement that the steepest increase occurs "from childhood into adulthood" also does not follow from the fitted model on the untransformed scale.
- [Section 3.2, Table 3] The claim that "the skin stiffness increases linearly with age" in the Langer-line direction is not supported by the fitted model. The outcome is the semi-minor axis length in RRT units (an arrival time), not stiffness. Combining Eq. (1) with the linear model gives E ∝ (semi-minor axis)⁻² under the authors' calibration, which is not a linear function of age. The manuscript should either model stiffness directly after converting RRT to wave speed, or explicitly state that the linear relation is for arrival time / semi-minor axis, not for stiffness.
- [Section 4, Section 2.2] The Discussion concedes that "the raw measurements from some subjects did not exhibit this symmetry and thus the ellipse fit was poor," yet the paper reports no count of such subjects, no measure of per-subject fit quality, and no analysis of whether poor fits are related to age. Because eccentricity is the outcome variable in Eq. (5) and is treated as error-free, poor ellipse fits introduce unmodeled measurement error that can bias the age coefficient. A sensitivity analysis that quantifies per-subject fit residuals, tests their correlation with age, and re-estimates the model after excluding poor fits is needed before the headline age trend can be accepted.
- [Section 2.4, Eq. (5)] Each subject contributes two rows (natural and stretched configurations), but the model has no subject-level random effect and therefore treats these repeated measurements as independent. This ignores within-subject correlation and can make the reported 0.95 HPDIs too narrow. The model should include a subject-level random intercept or otherwise account for the repeated-measures structure.
- [Section 2.1, Table 2, Section 1] The sample-size reporting is internally inconsistent: the Introduction states 72 subjects, Section 2.1 states 78 subjects (37 female, 41 male), and Table 2, as printed, does not sum to 78 and contains rows where the number of females exceeds the total (e.g., "11-20 3 10"). The age range also varies between the abstract (3–93) and Section 2.1 (3–92). These discrepancies must be corrected and the demographic table verified, as the "sizeable dataset" claim depends on this information.
minor comments (5)
- [Section 2.4] Specify the actual prior distributions used in rstanarm rather than only "default non-informative priors," and report MCMC diagnostics such as R-hat and effective sample size.
- [Table 1] Table 1 has formatting and unit issues: ρ is presumably in kg/m³, not kg/m², and the "RRT Conversion" values should be defined more clearly.
- [Figure 7] Figure 7 reports R² = 0.5105 from a simple log-linear model, but the Bayesian regression in Eq. (5) includes additional covariates; the relationship between this figure and the main model should be clarified.
- [Section 3.2] The gender coefficient is described with female as the baseline, but the coding of the gender covariate is not stated explicitly; it should be defined in the model description.
- [Section 2.2] The phrase "axisymmetric measurements" in Section 4 may be misleading; the data are expected to be symmetric under 180° rotation, not axisymmetric in the strict sense. Consider using "two-fold rotational symmetry" or similar.
Circularity Check
No significant circularity: the age-anisotropy finding is inferred from subject data with an external calibration constant; the acknowledged ellipse assumption is a stated modeling premise, not a concealed input.
full rationale
The derivation chain is self-contained with respect to its empirical claims. The central result (posterior mean age coefficient 0.007 on log-eccentricity, and -0.391 on semi-minor axis RRT; Section 3.2, Table 3) comes from fitting the multivariate regression in Eq. (5) directly to per-subject eccentricity, area, and semi-minor-axis measurements derived from Reviscometer data; no fitted parameter from that regression is recycled as a prediction. The stiffness calibration in Eq. (1) uses Young moduli measured by tensile tests on three elastomers (Table 1), an external input independent of the in-vivo age regressions. The simulation study (Section 2.3) is an in-sample performance comparison: data are generated from ellipses and the proposed eccentricity metric is benchmarked against that same parametric family, so it is a self-consistency check rather than independent validation of the ellipse assumption; however, the paper states this assumption explicitly ('Assuming an ellipse is an appropriate model for our data', Section 2.2) and concedes the limitation that 'the raw measurements from some subjects did not exhibit this symmetry and thus the ellipse fit was poor' (Section 4). This is a transparent modeling premise and robustness concern, not circularity. The only overlapping-author citation used in the reasoning (Deroy et al. 2017 for associating the fast-wave direction with Langer lines) is not load-bearing for the magnitude/anisotropy-age claim, since eccentricity is direction-independent. Separate correctness risks exist but are not circular: the phrase 'eccentricity increases logarithmically with age' inverts the functional form implied by Eq. (5), where log(e) is linear in age so e grows exponentially in age; and the number/systematic age-dependence of poor ellipse fits is unquantified. These are interpretation and robustness issues, not self-referential derivations.
Assumptions & free parameters
free parameters (1)
- RRT-to-time conversion factor =
0.284 microseconds per RRT unit
assumptions (5)
- domain assumption RRT readings are proportional to Rayleigh wave travel time, and E = rho * v^2 * 3.284 applies to in vivo skin
- domain assumption Angular RRT data are adequately modeled by an ellipse
- domain assumption The fastest wave direction corresponds to the Langer line
- domain assumption Natural and stretched measurements are independent conditional on covariates
- ad hoc to paper Age effects are linear on the log-eccentricity scale and on the semi-minor axis
Cite this review
Pith. "Pith review of Analysis of in-vivo skin anisotropy using elastic wave measurements and Bayesian modelling." pith.science (2026). https://pith.science/paper/AA4MAKQA
@misc{pith2026250602248,
author = {Pith},
title = {Pith review of: Analysis of in-vivo skin anisotropy using elastic wave measurements and Bayesian modelling},
year = {2026},
howpublished = {\url{https://pith.science/paper/AA4MAKQA}},
note = {Machine review of arXiv:2506.02248}
}
read the original abstract
In vivo skin exhibits viscoelastic, hyper-elastic and non-linear characteristics. It is under a constant non-equibiaxial tension in its natural configuration and is reinforced with oriented collagen fibers, giving rise to anisotropic behaviour. Understanding the complex mechanical behaviour of skin has relevance across many sectors including pharmaceuticals, cosmetics and surgery. However, there is a dearth of quality data characterizing human skin anisotropy in vivo. The available data is usually confined to limited population groups and/or limited angular resolution. Here, we use elastic waves travelling through the skin to obtain measurements from 78 volunteers from 3 to 93 years old. Using a Bayesian framework, we analyse the effect that age, gender and level of skin tension have on the skin anisotropy and stiffness. First, we propose a new measurement of anisotropy based on the eccentricity of angular data and conclude that it is a more robust measurement compared to the classic ``anisotropic ratio". We then find that in vivo skin anisotropy increases logarithmically with age, while the skin stiffness increases linearly along the direction of Langer Lines. We also conclude that gender does not significantly affect the skin anisotropy level, but does affect the overall stiffness, with males having stiffer skin on average. Finally, we find that skin tension significantly affects both the anisotropy and stiffness measurements, indicating that elastic wave measurements have promising applications in determining in vivo skin tension. In contrast to earlier studies, these results represent a comprehensive assessment of the variation of skin anisotropy with age and gender using a sizeable dataset and robust modern statistical analysis. This data has implications for the planning of surgical procedures and the adoption of universal cosmetic surgery practices for young or elderly patients.
Reference graph
Works this paper leans on
-
[1]
The Journal of the Acoustical Society of America 117:3469--3477
Bayón A, Gascón F, Nieves FJ (2005) Estimation of dynamic elastic constants from the amplitude and velocity of rayleigh waves. The Journal of the Acoustical Society of America 117:3469--3477
work page 2005
-
[2]
Skin Research and Technology 15:55--67
Boyer G, Laquièze L, Le Bot A, et al (2009) Dynamic indentation on human skin in vivo: ageing effects. Skin Research and Technology 15:55--67
work page 2009
-
[3]
Cremers J (2021) bpnreg: Bayesian projected normal regression models for circular data. R package version 2.0.0
work page 2021
-
[4]
Frontiers in Psychology 9:2040
Cremers J, Klugkist I (2018) One direction? a tutorial for circular data analysis using r with examples in cognitive psychology. Frontiers in Psychology 9:2040
work page 2018
-
[5]
Dahlgren RM, Elsnau WH (1984) Measurement of skin condition by sonic velocity
work page 1984
-
[6]
Deroy C, Destrade M, Alinden AM, et al (2017) Non-invasive evaluation of skin tension lines with elastic waves. Skin Res Technol
work page 2017
-
[7]
International Journal of Cosmetic Science 22:421--435
Diridollou S, Black D, Lagarde J, et al (2000) Sex- and site-dependent variations in the thickness and mechanical properties of human skin in vivo. International Journal of Cosmetic Science 22:421--435
work page 2000
-
[8]
The Journal of Surgical Research 54:21--28
Dombi G, Haut R, Sullivan W (1993) Correlation of high-speed tensile strength with collagen content in control and lathyritic rat skin. The Journal of Surgical Research 54:21--28
work page 1993
Show all 47 references
-
[9]
International Journal of Peptide Research and Therapeutics
Dulińska-Molak I, Pasikowska M, Pogoda K, et al (2014) Age-related changes in the mechanical properties of human fibroblasts and its prospective reversal after anti-wrinkle tripeptide treatment. International Journal of Peptide Research and Therapeutics
2014
-
[10]
American Journal of Epidemiology 153:1222--1226
Dunson DB (2001) Commentary: Practical advantages of bayesian analysis of epidemiologic data. American Journal of Epidemiology 153:1222--1226
2001
-
[11]
IEEE Transactions on Pattern Analysis and Machine Intelligence 21:476--480
Fitzgibbon A, Pilu M, Fisher R (1999) Direct least square fitting of ellipses. IEEE Transactions on Pattern Analysis and Machine Intelligence 21:476--480
1999
-
[12]
nikolai chernov
Gama J, Chernov N (2015) conicfit: Algorithms for fitting circles, ellipses and conics based on the work by prof. nikolai chernov. R package version 1.0.4
2015
-
[13]
Cambridge University Press, Cambridge, UK, (Ch
Gelman A, Hill J (2007) Data Analysis Using Regression and Multilevel/Hierarchical Models. Cambridge University Press, Cambridge, UK, (Ch. 11-15)
2007
-
[14]
Chapman and Hall/CRC
Gelman A, Carlin J, Stern H, et al (1995) Bayesian Data Analysis (1st ed.). Chapman and Hall/CRC
1995
-
[15]
Goodrich B, Gabry J, Ali I, et al (2020) rstanarm: Bayesian applied regression modeling via stan
2020
-
[16]
Medical Engineering & Physics 28:259--266
Hendriks F, Brokken D, Oomens C, et al (2006) The relative contributions of different skin layers to the mechanical behavior of human skin in vivo using suction experiments. Medical Engineering & Physics 28:259--266
2006
-
[17]
Exp Gerontol 36:363--372
Hermanns-Lê T, Jonlet F, Scheen A, et al (2001) Age- and body mass index-related changes in cutaneous shear wave velocity. Exp Gerontol 36:363--372
2001
-
[18]
Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine 232:323--343
Joodaki H, Panzer MB (2018) Skin mechanical properties and modeling: A review. Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine 232:323--343
2018
-
[19]
Skin Research and Technology 10:96--103
Khatyr F, Imberdis C, Varchon D, et al (2004) Model of the viscoelastic behaviour of skin in vivo and study of anisotropy. Skin Research and Technology 10:96--103
2004
-
[20]
Trends in Cognitive Sciences 14:293--300
Kruschke JK (2010) What to believe: Bayesian methods for data analysis. Trends in Cognitive Sciences 14:293--300
2010
-
[21]
Acta Biomaterialia 88:141--148
Laiacona D, Cohen J, Coulon K, et al (2019) Non-invasive in vivo quantification of human skin tension lines. Acta Biomaterialia 88:141--148
2019
-
[22]
Langer K (1978) On the anatomy and physiology of the skin. i. the cleavability of the cutis. (translated from langer, k. (1861). zur anatomie und physiologie der haut. i. Über die spaltbarkeit der cutis. sitzungsbericht der mathematisch-naturwissenschaftlichen classe der kaise...
1978
-
[23]
IEEE Transactions on Bio-Medical Engineering 57:953--959
Liang X, Boppart S (2009) Biomechanical properties of in vivo human skin from dynamic optical coherence elastography. IEEE Transactions on Bio-Medical Engineering 57:953--959
2009
-
[24]
In: Esprit Basic Research Series
Maurel W, Thalmann PD, Wu Y, et al (1998) Biomechanical models for soft tissue simulation. In: Esprit Basic Research Series. Springer
1998
-
[25]
Statistics in Medicine 38:2074--2102
Morris TP, White IR, Crowther MJ (2019) Using simulation studies to evaluate statistical methods. Statistics in Medicine 38:2074--2102
2019
-
[26]
The Quantitative Methods for Psychology 14:99--119
Muth C, Oravecz Z, Gabry J (2018) User-friendly bayesian regression modeling: A tutorial with rstanarm and shinystan. The Quantitative Methods for Psychology 14:99--119
2018
-
[27]
In: Biomechanics of Skin
Ní Annaidh A, Destrade M (2019) Tension lines of the skin. In: Biomechanics of Skin. Springer International Publishing, Cham, p 265--280
2019
-
[28]
Journal of the Mechanical Behavior of Biomedical Materials 5:139--148
Ní Annaidh A, Bruyère K, Destrade M, et al (2012) Characterization of the anisotropic mechanical properties of excised human skin. Journal of the Mechanical Behavior of Biomedical Materials 5:139--148
2012
-
[29]
Skin Research and Technology 17:101--107
Ohshima H, Tada A, Kanamaru A, et al (2011) Relevance of the directionality of skin elasticity to aging and sagging of the face. Skin Research and Technology 17:101--107
2011
-
[30]
Journal of the Mechanical Behavior of Biomedical Materials 41:241--250
Ottenio M, Tran D, Annaidh AN, et al (2015) Strain rate and anisotropy effects on the tensile failure characteristics of human skin. Journal of the Mechanical Behavior of Biomedical Materials 41:241--250
2015
-
[31]
Medical Engineering & Physics 30:599--606
Pailler-Mattei C, Bec S, Zahouani H (2008) In vivo measurements of the elastic mechanical properties of human skin by indentation tests. Medical Engineering & Physics 30:599--606
2008
-
[32]
Skin Research and Technology 20:463--472
Pamplona DC, Weber HI, Leta FR (2014) Optimization of the use of skin expanders. Skin Research and Technology 20:463--472
2014
-
[33]
Plast Reconstr Surg Glob Open 5:e1614
Paul SP (2017) Biodynamic excisional skin tension lines for excisional surgery of the lower limb and the technique of using parallel relaxing incisions to further reduce wound tension. Plast Reconstr Surg Glob Open 5:e1614
2017
-
[34]
Ann R Coll Surg Engl 100:330--337
Paul SP (2018) Biodynamic excisional skin tension lines for surgical excisions: untangling the science. Ann R Coll Surg Engl 100:330--337
2018
-
[35]
R Core Team (2020) R: A language and environment for statistical computing
2020
-
[36]
Skin Pharmacol Physiol
Ruvolo EJ, Stamatas G, Kollias N (2007) Skin viscoelasticity displays site- and age-dependent angular anisotropy. Skin Pharmacol Physiol
2007
-
[37]
Nature Reviews Methods Primers 1
Schoot R, Depaoli S, King R, et al (2021) Bayesian statistics and modelling. Nature Reviews Methods Primers 1
2021
-
[38]
Computer-Aided Design 45:551--555
Seo H, Kim Sj, Cordier F, et al (2013) Estimating dynamic skin tension lines in vivo using 3d scans. Computer-Aided Design 45:551--555. Solid and Physical Modeling 2012
2013
-
[39]
J Korean Med Sci 29:751--757
Son D, Harijan A (2014) Overview of surgical scar prevention and management. J Korean Med Sci 29:751--757
2014
-
[40]
J Mech Behav Biomed Mater 118:104340
Stowers C, Lee T, Bilionis I, et al (2021) Improving reconstructive surgery design using gaussian process surrogates to capture material behavior uncertainty. J Mech Behav Biomed Mater 118:104340
2021
-
[41]
J Mech Behav Biomed Mater 103:103551
Thieulin C, Pailler-Mattei C, Abdouni A, et al (2020) Mechanical and topographical anisotropy for human skin: Ageing effect. J Mech Behav Biomed Mater 103:103551
2020
-
[42]
Journal of Investigative Dermatology 113:732--739
Vexler A, Polyansky I, Gorodetsky R (1999) Evaluation of skin viscoelasticity and anisotropy by measurement of speed of shear wave propagation with viscoelasticity skin analyzer. Journal of Investigative Dermatology 113:732--739
1999
-
[43]
The American Statistician 70:129--133
Wasserstein RL, Lazar NA (2016) The asa statement on p-values: Context, process, and purpose. The American Statistician 70:129--133
2016
-
[44]
Skin Research and Technology 16:413--421
Xin S, Man W, Fluhr JW, et al (2010) Cutaneous resonance running time varies with age, body site and gender in a normal chinese population. Skin Research and Technology 16:413--421
2010
-
[45]
Wear 271:2364--2369
Zahouani H, Boyer G, Pailler-Mattei C, et al (2011) Effect of human ageing on skin rheology and tribology. Wear 271:2364--2369. 18th International Conference on Wear of Materials
2011
-
[46]
, " * write output.state after.block = add.period write newline
ENTRY address archive author booktitle chapter doi edition editor eid eprint howpublished institution journal key keywords month note number organization pages publisher school series title type url volume year archivePrefix primaryClass adsurl adsnote version label extra.labe...
-
[47]
write newline
" write newline "" before.all 'output.state := FUNCTION add.period duplicate empty 'skip "." * add.blank if FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap dupl...
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.