Gait analysis is a branch of biomechanics where its purpose is the study of mechanical laws relating to the way the body moves from one place to another. In most cases, the data sets for human gait analysis consist of continuous recordings of multiple physical activities, including kinematics and muscle performance. Despite the registered data being functions, the most common practice to detect any anomalies among experimental conditions consists of analyzing the vector of discrete observations or even summary measures of the curves. This fact causes an important information loss since the continuous nature of the data is being ignored. A suitable solution is to apply functional data analysis for analyzing continuous biomechanical data as functions, revealing the true nature of movement and allowing us to model and forecast the data with more precision. In the current paper, a new functional methodology for the analysis of variance with repeated measures was introduced. In particular, since functional data variability can be summarized by their first principal component scores, we proposed to turn the functional model into a multivariate one for the response of the most explicative principal components, and then, considered a semi-parametric approach to overcome the restrictive assumptions required in the classic repeated measures design. The motivation of this research was to contrast the differences in gait patterns of elementary school students when walking to school, depending on the type of bag they use to carry their school materials. The analysis reveals that gait joint movement is influenced by sex and the type of schoolbag, regardless of the load carried.
Citation: Helena Ortiz, Christian Acal, Manuel Escabias, Ana M. Aguilera. Characterizing the functional ANOVA model for repeated measures via PCA application to biomechanical data[J]. AIMS Mathematics, 2025, 10(4): 8468-8494. doi: 10.3934/math.2025390
[1] | Meshari Alesemi . Innovative approaches of a time-fractional system of Boussinesq equations within a Mohand transform. AIMS Mathematics, 2024, 9(10): 29269-29295. doi: 10.3934/math.20241419 |
[2] | Azzh Saad Alshehry, Humaira Yasmin, Ali M. Mahnashi . Analyzing fractional PDE system with the Caputo operator and Mohand transform techniques. AIMS Mathematics, 2024, 9(11): 32157-32181. doi: 10.3934/math.20241544 |
[3] | Aisha Abdullah Alderremy, Rasool Shah, Nehad Ali Shah, Shaban Aly, Kamsing Nonlaopon . Comparison of two modified analytical approaches for the systems of time fractional partial differential equations. AIMS Mathematics, 2023, 8(3): 7142-7162. doi: 10.3934/math.2023360 |
[4] | Musawa Yahya Almusawa, Hassan Almusawa . Numerical analysis of the fractional nonlinear waves of fifth-order KdV and Kawahara equations under Caputo operator. AIMS Mathematics, 2024, 9(11): 31898-31925. doi: 10.3934/math.20241533 |
[5] | M. Mossa Al-Sawalha, Khalil Hadi Hakami, Mohammad Alqudah, Qasem M. Tawhari, Hussain Gissy . Novel Laplace-integrated least square methods for solving the fractional nonlinear damped Burgers' equation. AIMS Mathematics, 2025, 10(3): 7099-7126. doi: 10.3934/math.2025324 |
[6] | Aslı Alkan, Halil Anaç . The novel numerical solutions for time-fractional Fornberg-Whitham equation by using fractional natural transform decomposition method. AIMS Mathematics, 2024, 9(9): 25333-25359. doi: 10.3934/math.20241237 |
[7] | Reetika Chawla, Komal Deswal, Devendra Kumar, Dumitru Baleanu . A novel finite difference based numerical approach for Modified Atangana- Baleanu Caputo derivative. AIMS Mathematics, 2022, 7(9): 17252-17268. doi: 10.3934/math.2022950 |
[8] | Emad Salah, Ahmad Qazza, Rania Saadeh, Ahmad El-Ajou . A hybrid analytical technique for solving multi-dimensional time-fractional Navier-Stokes system. AIMS Mathematics, 2023, 8(1): 1713-1736. doi: 10.3934/math.2023088 |
[9] | Qasem M. Tawhari . Advanced analytical techniques for fractional Schrödinger and Korteweg-de Vries equations. AIMS Mathematics, 2025, 10(5): 11708-11731. doi: 10.3934/math.2025530 |
[10] | Ritu Agarwal, Mahaveer Prasad Yadav, Dumitru Baleanu, S. D. Purohit . Existence and uniqueness of miscible flow equation through porous media with a non singular fractional derivative. AIMS Mathematics, 2020, 5(2): 1062-1073. doi: 10.3934/math.2020074 |
Gait analysis is a branch of biomechanics where its purpose is the study of mechanical laws relating to the way the body moves from one place to another. In most cases, the data sets for human gait analysis consist of continuous recordings of multiple physical activities, including kinematics and muscle performance. Despite the registered data being functions, the most common practice to detect any anomalies among experimental conditions consists of analyzing the vector of discrete observations or even summary measures of the curves. This fact causes an important information loss since the continuous nature of the data is being ignored. A suitable solution is to apply functional data analysis for analyzing continuous biomechanical data as functions, revealing the true nature of movement and allowing us to model and forecast the data with more precision. In the current paper, a new functional methodology for the analysis of variance with repeated measures was introduced. In particular, since functional data variability can be summarized by their first principal component scores, we proposed to turn the functional model into a multivariate one for the response of the most explicative principal components, and then, considered a semi-parametric approach to overcome the restrictive assumptions required in the classic repeated measures design. The motivation of this research was to contrast the differences in gait patterns of elementary school students when walking to school, depending on the type of bag they use to carry their school materials. The analysis reveals that gait joint movement is influenced by sex and the type of schoolbag, regardless of the load carried.
The purpose of this paper is to study the global behavior of the following max-type system of difference equations of the second order with four variables and period-two parameters
{xn=max{An,zn−1yn−2},yn=max{Bn,wn−1xn−2},zn=max{Cn,xn−1wn−2},wn=max{Dn,yn−1zn−2}, n∈N0≡{0,1,2,⋯}, | (1.1) |
where An,Bn,Cn,Dn∈R+≡(0,+∞) are periodic sequences with period 2 and the initial values x−i,y−i,z−i,w−i∈R+ (1≤i≤2). To do this we will use some methods and ideas which stems from [1,2]. For a more complex variant of the method, see [3]. A solution {(xn,yn,zn,wn)}+∞n=−2 of (1.1) is called an eventually periodic solution with period T if there exists m∈N such that (xn,yn,zn,wn)=(xn+T,yn+T,zn+T,wn+T) holds for all n≥m.
When xn=yn and zn=wn and A0=A1=B0=B1=α and C0=C1=D0=D1=β, (1.1) reduces to following max-type system of difference equations
{xn=max{α,zn−1xn−2},zn=max{β,xn−1zn−2}, n∈N0. | (1.2) |
Fotiades and Papaschinopoulos in [4] investigated the global behavior of (1.2) and showed that every positive solution of (1.2) is eventually periodic.
When xn=zn and yn=wn and An=Cn and Bn=Dn, (1.1) reduces to following max-type system of difference equations
{xn=max{An,yn−1xn−2},yn=max{Bn,xn−1yn−2}, n∈N0. | (1.3) |
Su et al. in [5] investigated the periodicity of (1.3) and showed that every solution of (1.3) is eventually periodic.
In 2020, Su et al. [6] studied the global behavior of positive solutions of the following max-type system of difference equations
{xn=max{A,yn−txn−s},yn=max{B,xn−tyn−s}, n∈N0, |
where A,B∈R+.
In 2015, Yazlik et al. [7] studied the periodicity of positive solutions of the max-type system of difference equations
{xn=max{1xn−1,min{1,pyn−1}},yn=max{1yn−1,min{1,pxn−1}}, n∈N0, | (1.4) |
where p∈R+ and obtained in an elegant way the general solution of (1.4).
In 2016, Sun and Xi [8], inspired by the research in [5], studied the following more general system
{xn=max{1xn−m,min{1,pyn−r}},yn=max{1yn−m,min{1,qxn−t}}, n∈N0, | (1.5) |
where p,q∈R+, m,r,t∈N≡{1,2,⋯} and the initial conditions x−i,y−i∈R+ (1≤i≤s) with s=max{m,r,t} and showed that every positive solution of (1.5) is eventually periodic with period 2m.
In [9], Stević studied the boundedness character and global attractivity of the following symmetric max-type system of difference equations
{xn=max{B,ypn−1xpn−2},yn=max{B,xpn−1ypn−2}, n∈N0, |
where B,p∈R+ and the initial conditions x−i,y−i∈R+ (1≤i≤2).
In 2014, motivated by results in [9], Stević [10] further study the behavior of the following max-type system of difference equations
{xn=max{B,ypn−1zpn−2},yn=max{B,zpn−1xpn−2},zn=max{B,xpn−1ypn−2}. n∈N0, | (1.6) |
where B,p∈R+ and the initial conditions x−i,y−i,z−i∈R+ (1≤i≤2), and showed that system (1.6) is permanent when p∈(0,4).
For more many results for global behavior, eventual periodicity and the boundedness character of positive solutions of max-type difference equations and systems, please readers refer to [11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30] and the related references therein.
In this section, we study the global behavior of system (1.1). For any n≥−1, write
{x2n=A2nXn,y2n=B2nYn,z2n=C2nZn,w2n=D2nWn,x2n+1=A2n+1X′n,y2n+1=B2n+1Y′n,z2n+1=C2n+1Z′n,w2n+1=D2n+1W′n. |
Then, (1.1) reduces to the following system
{Xn=max{1,C2n−1Z′n−1A2nB2nYn−1},Yn=max{1,D2n−1W′n−1B2nA2nXn−1},Z′n=max{1,A2nXnC2n+1D2n+1W′n−1},W′n=max{1,B2nYnD2n+1C2n+1Z′n−1},Zn=max{1,A2n−1X′n−1C2nD2nWn−1},Wn=max{1,B2n−1Y′n−1D2nC2nZn−1},X′n=max{1,C2nZnA2n+1B2n+1Y′n−1},Y′n=max{1,D2nWnB2n+1A2n+1X′n−1}, n∈N0. | (2.1) |
From (2.1) we see that it suffices to consider the global behavior of positive solutions of the following system
{un=max{1,bvn−1aAUn−1},Un=max{1,BVn−1aAun−1},vn=max{1,aunbBVn−1},Vn=max{1,AUnbBvn−1}, n∈N0, | (2.2) |
where a,b,A,B∈R+, the initial conditions u−1,U−1,v−1,V−1∈R+. If (un,Un,vn,Vn,a,A,b,B)=(Xn,Yn,Z′n,W′n,A2n,B2n,C2n−1,D2n−1), then (2.2) is the first four equations of (2.1). If (un,Un,vn,Vn,a,A,b,B)=(Zn,Wn,X′n,Y′n,C2n,D2n,A2n−1,B2n−1), then (2.2) is the next four equations of (2.1). In the following without loss of generality we assume a≤A and b≤B. Let {(un,Un,vn,Vn)}∞n=−1 be a positive solution of (2.2).
Proposition 2.1. If ab<1, then there exists a solution {(un,Un,vn,Vn)}∞n=−1 of (2.2) such that un=vn=1 for any n≥−1 and limn⟶∞Un=limn⟶∞Vn=∞.
Proof. Let u−1=v−1=1 and U−1=V−1=max{baA,aAB,abB}+1. Then, from (2.2) we have
{u0=max{1,bv−1aAU−1}=1,U0=max{1,BV−1aAu−1}=BV−1aA,v0=max{1,au0bBV−1}=1,V0=max{1,AU0bBv−1}=V−1ab, |
and
{u1=max{1,bv0aAU0}=max{1,bBV−1}=1,U1=max{1,BV0aAu0}=max{1,BV−1aAab}=BV−1aAab,v1=max{1,au1bBV0}=max{1,aabbBV−1}=1,V1=max{1,AU1bBv0}=max{1,V−1(ab)2}=V−1(ab)2. |
Suppose that for some k∈N, we have
{uk=1,Uk=BV−1aA(ab)k,vk=1,Vk=V−1(ab)k+1. |
Then,
{uk+1=max{1,bvkaAUk}=max{1,b(ab)kBV−1}=1,Uk+1=max{1,BVkaAuk}=max{1,BV−1aA(ab)k+1}=BV−1aA(ab)k+1,vk+1=max{1,auk+1bBVk}=max{1,a(ab)k+1bBV−1}=1,Vk+1=max{1,AUk+1bBvk}=max{1,V−1(ab)k+2}=V−1(ab)k+2. |
By mathematical induction, we can obtain the conclusion of Proposition 2.1. The proof is complete.
Now, we assume that ab≥1. Then, from (2.2) it follows that
{un=max{1,bvn−1aAUn−1},Un=max{1,BVn−1aAun−1},vn=max{1,abBVn−1,vn−1ABUn−1Vn−1},Vn=max{1,AbBvn−1,Vn−1abun−1vn−1}, n∈N0. | (2.3) |
Lemma 2.1. The following statements hold:
(1) For any n∈N0,
un, Un, vn, Vn∈[1,+∞). | (2.4) |
(2) If ab≥1, then for any k∈N and n≥k+2,
{un=max{1,baAUn−1,bvkaA(AB)n−k−1Un−1Un−2Vn−2⋯UkVk},Un=max{1,BaAun−1,BVkaA(ab)n−k−1un−1un−2vn−2⋯ukvk},vn=max{1,abBVn−1,vk(AB)n−kUn−1Vn−1⋯UkVk},Vn=max{1,AbBvn−1,Vk(ab)n−kun−1vn−1⋯ukvk}. | (2.5) |
(3) If ab≥1, then for any k∈N and n≥k+4,
{1≤vn≤vn−2,1≤Vn≤AaVn−2,1≤un≤max{1,bBun−2,bvkaA(AB)n−k−1},1≤Un≤max{1,BbUn−2,BVkaA(ab)n−k−1}. | (2.6) |
Proof. (1) It follows from (2.2).
(2) Since AB≥ab≥1, it follows from (2.2) and (2.3) that for any k∈N and n≥k+2,
un=max{1,bvn−1aAUn−1}=max{1,baAUn−1max{1,abBVn−2,vn−2ABUn−2Vn−2}}=max{1,baAUn−1,bvn−2ABaAUn−1Un−2Vn−2}=max{1,baAUn−1,bABaAUn−1Un−2Vn−2max{1,abBVn−1,vn−3ABUn−3Vn−3}}=max{1,baAUn−1,bvn−3(AB)2aAUn−1Un−2Vn−2Un−3Vn−3}⋯=max{1,baAUn−1,bvkaA(AB)n−k−1Un−1Un−2Vn−2⋯UkVk}. |
In a similar way, also we can obtain the other three formulas.
(3) By (2.5) one has that for any k∈N and n≥k+2,
{un≥baAUn−1,Un≥BaAun−1,vn≥abBVn−1,Vn≥AbBvn−1, |
from which and (2.4) it follows that for any n≥k+4,
{1≤un≤max{1,bBun−2,bvkaA(AB)n−k−1},1≤Un≤max{1,BbUn−2,BVkaA(ab)n−k−1},1≤vn≤max{1,avn−2A,vn−2}=vn−2,1≤Vn≤max{1,AVn−2a,Vn−2}=AVn−2a. |
The proof is complete.
Proposition 2.2. If ab=AB=1, then {(un,Un,vn,Vn)}+∞n=−1 is eventually periodic with period 2.
Proof. By the assumption we see a=A and b=B. By (2.5) we see that for any k∈N and n≥k+2,
{un=max{1,b3Un−1,b3vkUn−1Un−2Vn−2⋯UkVk},Un=max{1,b3un−1,b3Vkun−1un−2vn−2⋯ukvk},vn=max{1,a3Vn−1,vkUn−1Vn−1⋯UkVk},Vn=max{1,a3vn−1,Vkun−1vn−1⋯ukvk}. | (2.7) |
(1) If a=b=1, then it follows from (2.7) and (2.4) that for any n≥k+4,
{un=max{1,vkUn−1Un−2Vn−2⋯UkVk}≤max{1,vkUn−2Un−3Vn−3⋯UkVk}=un−1,Un=max{1,Vkun−1un−2vn−2⋯ukvk}≤Un−1,vn=max{1,vkUn−1Vn−1⋯UkVk}≤vn−1,Vn=max{1,Vkun−1vn−1⋯ukvk}≤Vn−1. | (2.8) |
We claim that vn=1 for any n≥6 or Vn=1 for any n≥6. Indeed, if vn>1 for some n≥6 and Vm>1 for some m≥6, then
vn=v1Un−1Vn−1⋯U1V1>1, Vm=V1um−1vm−1⋯u1v1>1, |
which implies
1≥v1Un−1Vn−1⋯U1V1V1um−1vm−1⋯u1v1=Vmvn>1. |
A contradiction.
If vn=1 for any n≥6, then by (2.8) we see un=1 for any n≥10, which implies Un=Vn=V10.
If Vn=1 for any n≥6, then by (2.8) we see Un=1 for any n≥10, which implies vn=un=v10.
Then, {(un,Un,vn,Vn)}+∞n=−1 is eventually periodic with period 2.
(2) If a<1<b, then it follows from (2.7) that for any n≥k+4,
{un=max{1,b3Un−1,b3vkUn−1Un−2Vn−2⋯UkVk},Un=max{1,b3un−1,b3Vkun−1un−2vn−2⋯ukvk},vn=max{1,vkUn−1Vn−1⋯UkVk}≤vn−1,Vn=max{1,Vkun−1vn−1⋯ukvk}≤Vn−1. | (2.9) |
It is easy to verify vn=1 for any n≥6 or Vn=1 for any n≥6.
If Vn=vn=1 eventually, then by (2.9) we have
{1≥vkUn−1Vn−1⋯UkVk eventually,1≥Vkun−1vn−1⋯ukvk eventually. |
Since Un≥b3un−1 and un≥b3Un−1, we see
{un=max{1,b3Un−1,b3vkUn−1Un−2Vn−2⋯UkVk}=max{1,b3Un−1}≤un−2 eventually,Un=max{1,b3un−1,b3Vkun−1un−2vn−2⋯ukvk}=max{1,b3un−1}≤Un−2 eventually, |
which implies
{un−2≥un=max{1,b3Un−1}≥max{1,b3Un−3}=un−2 eventually,Un−2≥Un=max{1,b3un−1}≥max{1,b3un−3}=Un−2 eventually. |
If Vn>1=vn eventually, then by (2.9) we have
{1≥vkUn−1Vn−1⋯UkVk eventually,Vn=Vkun−1vn−1⋯ukvk>1 eventually. |
Thus,
{un=max{1,b3Un−1,b3vkUn−1Un−2Vn−2⋯UkVk}=max{1,b3Un−1}≤un−2 eventually,Un=max{1,b3un−1,b3Vkun−1un−2vn−2⋯ukvk}=max{1,b3Vkun−1un−2vn−2⋯ukvk}≤Un−2 eventually, |
which implies
{un−2≥un=max{1,b3Un−1}≥max{1,b3Un−3}=un−2 eventually,Un=1 eventually or b3Vk eventually. |
If Vn=1<vn eventually, then by (2.9) we have Un−2=Un eventually and un=un−1 eventually. By the above we see that {(un,Un,vn,Vn)}+∞n=−1 is eventually periodic with period 2.
(3) If b<1<a, then for any k∈N and n≥k+2,
{un=max{1,b3vkUn−1Un−2Vn−2⋯UkVk}≤un−1,Un=max{1,b3Vkun−1un−2vn−2⋯ukvk}≤Un−1,vn=max{1,a3Vn−1,vkUn−1Vn−1⋯UkVk},Vn=max{1,a3vn−1,Vkun−1vn−1⋯ukvk}. | (2.10) |
It is easy to verify un=1 for any n≥3 or Un=1 for any n≥3.
If un=Un=1 eventually, then
{1≥b3vkUn−1Un−2Vn−2⋯UkVk eventually,1≥b3Vkun−1un−2vn−2⋯ukvk eventually. |
Thus, by (2.6) we have
{vn−2≥vn=max{1,a3Vn−1,vkUn−1Vn−1⋯UkVk}=max{1,a3Vn−1}≥vn−2 eventually,Vn−2≥Vn=max{1,a3vn−1,Vkun−1vn−1⋯ukvk}=max{1,a3vn−1}≥Vn−2 eventually. |
If un=1<Un eventually, then
{1≥b3vkUn−1Un−2Vn−2⋯UkVk eventually,1<b3Vkun−1un−2vn−2⋯ukvk=Un eventually. |
Thus,
{vn−2≥vn=max{1,a3Vn−1,vkUn−1Vn−1⋯UkVk}=max{1,a3Vn−1}≥vn−2 eventually,Vn=max{1,a3vn−1,Vkun−1vn−1⋯ukvk}=max{1,Vkun−1vn−1⋯ukvk}=1 eventually or Vk eventually. |
If un>1=Un eventually, then we have Vn=Vn−2 eventually and vn=1 eventually or vn=vk eventually.
By the above we see that {(un,Un,vn,Vn)}+∞n=−1 is eventually periodic with period 2.
Proposition 2.3. If ab=1<AB, then {(un,Un,vn,Vn)}+∞n=−1 is eventually periodic with period 2.
Proof. Note that Un≥BaAun−1 and Vn≥AbBvn−1. By (2.5) we see that there exists N∈N such that for any n≥N,
{un=max{1,b2AUn−1}≤un−2,Un=max{1,BaAun−1,BVkaAun−1un−2vn−2⋯ukvk},vn=max{1,a2BVn−1}≤vn−2,Vn=max{1,AbBvn−1,Vkun−1vn−1⋯ukvk}. | (2.11) |
It is easy to verify that un=1 for any n≥N+1 or vn=1 for any n≥N+1.
If un=vn=1 eventually, then by (2.11) we see that Un=Un−1 eventually and Vn=Vn−1 eventually.
If uM+2n>1=vn eventually for some M∈N, then by (2.11) and (2.4) we see that
{uM+2n=b2AUM+2n−1>1 eventually,UM+2n+1=max{1,BbUM+2n−1,BVkaAuM+2nuM+2n−1vM+2n−1⋯ukvk}≥BbUM+2n−1 eventually,vn=max{1,a2BVn−1}=1 eventually,Vn=max{1,AbBvn−1,Vkun−1vn−1⋯ukvk}≤Vn−1 eventually. |
By (2.11) we see that Un is bounded, which implies B=b.
If UM+2n−1≤BVkaAuM+2nuM+2n−1vM+2n−1⋯ukvk eventually, then
UM+2n+1=BVkaAuM+2nuM+2n−1vM+2n−1⋯ukvk≤UM+2n−1 eventually. |
Thus, UM+2n+1=UM+2n−1 eventually and uM+2n=uM+2n−2 eventually. Otherwise, we have UM+2n+1=UM+2n−1 eventually and uM+2n=uM+2n−2 eventually. Thus, Vn=Vn−1=max{1,AbB} eventually since limn⟶∞Vkun−1vn−1⋯ukvk=0. By (2.2) it follows UM+2n=UM+2n−2 eventually and uM+2n+1=uM+2n−1 eventually.
If vM+2n>1=un eventually for some M∈N, then we may show that {(un,Un,vn,Vn)}+∞n=−1 is eventually periodic with period 2. The proof is complete.
Proposition 2.4. If ab>1, then {(un,Un,vn,Vn)}+∞n=−1 is eventually periodic with period 2.
Proof. By (2.5) we see that there exists N∈N such that for any n≥N,
{un=max{1,baAUn−1},Un=max{1,BaAun−1},vn=max{1,abBVn−1},Vn=max{1,AbBvn−1}. | (2.12) |
If a<A, then for n≥2k+N with k∈N,
vn=max{1,abBVn−1}≤max{1,aAvn−2}≤⋯≤max{1,(aA)kvn−2k}, |
which implies vn=1 eventually and Vn=max{1,AbB} eventually.
If a=A, then
{vn=max{1,abBVn−1}≤vn−2 eventually,Vn=max{1,AbBvn−1}≤Vn−2 eventually. |
Which implies
{vn−2≥vn=max{1,abBVn−1}≥max{1,abBVn−3}=vn−2 eventually,Vn−2≥Vn=max{1,AbBvn−1}≥max{1,AbBvn−3}=Vn−2 eventually. |
Thus, Vn,vn are eventually periodic with period 2. In a similar way, we also may show that Un,un are eventually periodic with period 2. The proof is complete.
From (2.1), (2.2), Proposition 2.1, Proposition 2.2, Proposition 2.3 and Proposition 2.4 one has the following theorem.
Theorem 2.1. (1) If min{A0C1,B0D1,A1C0,B1D0}<1, then system (1.1) has unbounded solutions.
(2) If min{A0C1,B0D1,A1C0,B1D0}≥1, then every solution of system (1.1) is eventually periodic with period 4.
In this paper, we study the eventual periodicity of max-type system of difference equations of the second order with four variables and period-two parameters (1.1) and obtain characteristic conditions of the coefficients under which every positive solution of (1.1) is eventually periodic or not. For further research, we plan to study the eventual periodicity of more general max-type system of difference equations by the proof methods used in this paper.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Project supported by NSF of Guangxi (2022GXNSFAA035552) and Guangxi First-class Discipline SCPF(2022SXZD01, 2022SXYB07) and Guangxi Key Laboratory BDFE(FED2204) and Guangxi University of Finance and Economics LSEICIC(2022YB12).
There are no conflict of interest in this article.
[1] |
R. K. Fukuchi, C. A. Fukuchi, M. Duarte, A public dataset of running biomechanics and the effects of running speed on lower extremity kinematics and kinetics, PeerJ, 5 (2017), e3298. http://doi.org/10.7717/peerj.3298 doi: 10.7717/peerj.3298
![]() |
[2] |
E. Orantes-González, J. Heredia-Jiménez, V. M. Soto-Hermoso, The effect of school trolley load on spatiotemporal gait parameters of children, Gait Posture, 42 (2015), 390–393. https://doi.org/10.1016/j.gaitpost.2015.06.004 doi: 10.1016/j.gaitpost.2015.06.004
![]() |
[3] |
D. A. Santos, M. Duarte, A public data set of human balance evaluations, PeerJ, 4 (2016), e2648. http://doi.org/10.7717/peerj.2648 doi: 10.7717/peerj.2648
![]() |
[4] | J. O. Ramsay, B. W. Silverman, Applied functional data analysis: Methods and case studies, 1 Eds., New York, NY: Springer, 2002. https://doi.org/10.1007/b98886 |
[5] | J. O. Ramsay, B. W. Silverman, Functional data analysis, 2 Eds., New York, NY: Springer, 2005. https://doi.org/10.1007/b98888 |
[6] | J. O. Ramsay, B. W. Silverman, G. Graves, Functional data analysis with R and MATLAB, New York, NY: Springer, 2009. https://doi.org/10.1007/978-0-387-98185-7 |
[7] | F. Ferraty, P. Vieu, Nonparametric functional data analysis: Theory and practice, New York, NY: Springer, 2006. https://doi.org/10.1007/0-387-36620-2 |
[8] | L. Horváth, P. Kokoszka, Inference for functional data with applications, New York, NY: Springer, 2012. https://doi.org/10.1007/978-1-4614-3655-3 |
[9] | D. G. E. Robertson, G. E. Caldwell, J. Hamill, G. Kamen, S. Whittlesey, Research methods in biomechanics, 2Eds., Human Kinetics, 2013. |
[10] |
J. Park, M. K. Seeley, D. Francom, C. S. Reese, J. T. Hopkins, Functional vs. traditional analysis in biomechanical gait data: An alternative statistical approach, J. Hum. Kinet., 60 (2017), 39–49. https://doi.org/10.1515/hukin-2017-0114 doi: 10.1515/hukin-2017-0114
![]() |
[11] |
J. Warmenhoven, A. Harrison, M. A. Robinson, J. Vanrenterghem, N. Bargary, R. Smith, et al., A force profile analysis comparison between functional data analysis, statistical parametric mapping and statistical non-parametric mapping in on-water single sculling, J. Sci. Med. Sport, 21 (2018), 1100–1105. https://doi.org/10.1016/j.jsams.2018.03.009 doi: 10.1016/j.jsams.2018.03.009
![]() |
[12] |
L. V. de Paula, G. T. de Oliveira, R. M. Ferreira, E. R. Soares, Functional data analysis of the force unsteadiness, J. Appl. Physiol., 136 (2024), 1270. https://doi.org/10.1152/japplphysiol.00200.2024 doi: 10.1152/japplphysiol.00200.2024
![]() |
[13] | E. Crane, D. Childers, G. Gerstner, E. Rothman, Functional data analysis for biomechanics, In: Theoretical biomechanics, 2011, 77–92. http://doi.org/10.5772/22382 |
[14] |
J. Dannenmaier, C. Kaltenbach, T. Kölle, G. Krischak, Application of functional data analysis to explore movements: walking, running and jumping–A systematic review, Gait Posture, 77 (2020), 182–189. https://doi.org/10.1016/j.gaitpost.2020.02.002 doi: 10.1016/j.gaitpost.2020.02.002
![]() |
[15] | J. T. Zhang, Analysis of variance for functional data, Chapman and Hall/CRC, 2013. https://doi.org/10.1201/b15005 |
[16] |
T. Górecki, L. Smaga, A Comparison of tests for the one-way ANOVA problem for functional data, Comput. Stat., 30 (2015), 987–1010. https://doi.org/10.1007/s00180-015-0555-0 doi: 10.1007/s00180-015-0555-0
![]() |
[17] |
P. Martinez-Camblor, N. Corral, Repeated measures analysis for functional data, Comput. Stat. Data Anal., 55 (2011), 3244–3256. https://doi.org/10.1016/j.csda.2011.06.007 doi: 10.1016/j.csda.2011.06.007
![]() |
[18] |
L. Smaga, A note on repeated measures analysis for functional data, AStA Adv. Stat. Anal., 104 (2020), 117–139. https://doi.org/10.1007/s10182-018-00348-8 doi: 10.1007/s10182-018-00348-8
![]() |
[19] |
C. Acal, A. M. Aguilera, A. Sarra, A. Evangelista, T. Di Battista, S. Palermi, Functional ANOVA approaches for detecting changes in air pollution during the COVID-19 pandemic, Stoch. Environ. Res. Risk Assess., 36 (2022), 1083–1101. https://doi.org/10.1007/s00477-021-02071-4 doi: 10.1007/s00477-021-02071-4
![]() |
[20] |
C. Acal, A. M. Aguilera, Basis expansion approaches for functional analysis of variance with repeated measures, Adv. Data Anal. Classif., 17 (2023), 291–321. https://doi.org/10.1007/s11634-022-00500-y doi: 10.1007/s11634-022-00500-y
![]() |
[21] | M. Escabias, A. M. Aguilera, J. M. Heredia-Jiménez, E. Orantes-González, Functional data analysis in kinematics of children going to school, In: Functional statistics and related fields, Springer, 2017, 95–103. https://doi.org/10.1007/978-3-319-55846-2_13 |
[22] |
J. Warmenhoven, N. Bargary, D. Liebl, A. Harrison, M. A. Robinson, E. Gunning, G. Hooker, PCA of waveforms and functional PCA: A primer for biomechanics, J. Biomech., 116 (2021), 110106. https://doi.org/10.1016/j.jbiomech.2020.110106 doi: 10.1016/j.jbiomech.2020.110106
![]() |
[23] |
K. E. Roach, V. Pedoia, J. J. Lee, T. Popovic, T. M. Link, S. Majumdar, et al., Multivariate functional principal component analysis identifies waveform features of gait biomechanics related to early-to-moderate hip osteoarthritis, J. Orthop. Res., 39 (2021), 1722–1731. https://doi.org/10.1002/jor.24901 doi: 10.1002/jor.24901
![]() |
[24] |
K. Yoshida, D. Commandeur, S. Hundza, M. Klimstra, Detecting differences in gait initiation between older adult fallers and non-fallers through multivariate functional principal component analysis, J. Biomech., 144 (2022), 111342. https://doi.org/10.1016/j.jbiomech.2022.111342 doi: 10.1016/j.jbiomech.2022.111342
![]() |
[25] |
A. M. Aguilera, R. Gutiérrez, F. A. Ocaña, M. J. Valderrama, Computational approaches to estimation in the principal component analysis of a stochastic process, Appl. Stoch. Models Data Anal., 11 (1995), 279–299. https://doi.org/10.1002/asm.3150110402 doi: 10.1002/asm.3150110402
![]() |
[26] |
P. Hall, M. Hosseini-Nasab, On properties of functional principal components analysis, J. R. Stat. Soc. B, 68 (2006), 109–126. https://doi.org/10.1111/j.1467-9868.2005.00535.x doi: 10.1111/j.1467-9868.2005.00535.x
![]() |
[27] |
C. Acal, A. M. Aguilera, M. Escabias, New modeling approaches based on Varimax rotation of functional principal components, Mathematics, 8 (2020), 2085. https://doi.org/10.3390/math8112085 doi: 10.3390/math8112085
![]() |
[28] |
B. Shi, P. Wei, X. Huang, Functional principal component based landmark analysis for the effects of longitudinal cholesterol profiles on the risk of coronary heart disease, Stat. Med., 40 (2020), 650–667. https://doi.org/10.1002/sim.8794 doi: 10.1002/sim.8794
![]() |
[29] |
J. Jacques, C. Preda, Functional data clustering: A survey, Adv. Data Anal. Classif., 8 (2014), 231–255. https://doi.org/10.1007/s11634-013-0158-y doi: 10.1007/s11634-013-0158-y
![]() |
[30] |
F. Fortuna, F. Maturo, K-means clustering of item characteristic curves and item information curves via functional principal component analysis, Qual. Quant., 53 (2019), 2291–2304. https://doi.org/10.1007/s11135-018-0724-7 doi: 10.1007/s11135-018-0724-7
![]() |
[31] |
F. Maturo, Unsupervised classification of ecological communities ranked according to their biodiversity patterns via a functional principal component decomposition of Hill's numbers integral functions, Ecol. Indic., 90 (2018), 305–315. https://doi.org/10.1016/j.ecolind.2018.03.013 doi: 10.1016/j.ecolind.2018.03.013
![]() |
[32] |
F. Maturo, R. Verde, Pooling random forest and functional data analysis for biomedical signals supervised classification: Theory and application to electrocardiogram data, Stat. Med., 41 (2022), 2247–2275. https://doi.org/10.1002/sim.9353 doi: 10.1002/sim.9353
![]() |
[33] |
M. C. Aguilera-Morillo, A. M. Aguilera, M. Escabias, M. J. Valderrama, Penalized spline approaches for functional logit regression, Test, 22 (2013), 251–277. https://doi.org/10.1007/s11749-012-0307-1 doi: 10.1007/s11749-012-0307-1
![]() |
[34] |
M. Escabias, A. M. Aguilera, M. C. Aguilera-Morillo, Functional PCA and base-line logit models, J. Classif., 31 (2014), 296–324. https://doi.org/10.1007/s00357-014-9162-y doi: 10.1007/s00357-014-9162-y
![]() |
[35] |
A. Evangelista, C. Acal, A. M. Aguilera, A. Sarra, T. Di Battista, S. Palermi, High-dimensional variable selection through group Lasso for multiple function-on-function linear regression: A case study in PM10 monitoring, Environmetrics, 36 (2024), e2852. https://doi.org/10.1002/env.2852 doi: 10.1002/env.2852
![]() |
[36] |
B. Shi, R. T. Ogden, Inference in functional mixed regression models with applications to Positron Emission Tomography imaging data, Stat. Med., 40 (2021), 4640–4659. https://doi.org/10.1002/sim.9087 doi: 10.1002/sim.9087
![]() |
[37] |
A. M. Aguilera, C. Acal, M. C. Aguilera-Morillo, F. Jiménez-Molinos, J. B. Roldán, Homogeneity problem for basis expansion of functional data with applications to resistive memories, Math. Comput. Simulat., 186 (2021), 41–51. https://doi.org/10.1016/j.matcom.2020.05.018 doi: 10.1016/j.matcom.2020.05.018
![]() |
[38] | N. H. Timm, 2 Multivariate analysis of variance of repeated measurements, In: Handbook of Statistics, 1, Elsevier, 1980, 41–87. https://doi.org/10.1016/S0169-7161(80)01004-8 |
[39] |
R. J. Boik, The mixed model for multivariate repeated measures: Validity conditions and an approximate test, Psychometrika, 53 (1988), 469–486. https://doi.org/10.1007/BF02294401 doi: 10.1007/BF02294401
![]() |
[40] |
R. J. Boik, Scheffés mixed model for multivariate repeated measures: A relative efficiency evaluation, Commun. Stat.-Theory M., 20 (1991), 1233–1255. https://doi.org/10.1080/03610929108830562 doi: 10.1080/03610929108830562
![]() |
[41] |
F. Konietschke, A. C. Bathke, S. W. Harrar, M. Pauly, Parametric and nonparametric bootstrap methods for general MANOVA, J. Multivariate Anal., 140 (2015), 291–301. https://doi.org/10.1016/j.jmva.2015.05.001 doi: 10.1016/j.jmva.2015.05.001
![]() |
[42] |
A. C. Bathke, S. Friedrich, M. Pauly, F. Konietschke, W. Staffen, N. Strobl, et al., Testing mean differences among groups: Multivariate and repeated measures analysis with minimal assumptions, Multivar. Behav. Res., 53 (2018), 348–359. https://doi.org/10.1080/00273171.2018.1446320 doi: 10.1080/00273171.2018.1446320
![]() |
[43] |
E. Orantes-González, J. Heredia-Jiménez, G. J. Beneck, Children require less gait kinematic adaptations to pull a trolley than to carry a backpack, Gait Posture, 52 (2017), 189–193. https://doi.org/10.1016/j.gaitpost.2016.11.041 doi: 10.1016/j.gaitpost.2016.11.041
![]() |
[44] |
E. Orantes-González, J. Heredia-Jiménez, M. A. Robinson, A kinematic comparison of gait with a backpack versus a trolley for load carriage in children, Appl. Ergon., 80 (2019), 28–34. https://doi.org/10.1016/j.apergo.2019.05.003 doi: 10.1016/j.apergo.2019.05.003
![]() |
[45] |
N. E. Helwig, K. A. Shorter, P. Ma, E. T. Hsiao-Wecksler, Smoothing spline analysis of variance models: A new tool for the analysis of cyclic biomechanical data, J. biomech., 49 (2016), 3216–3222. https://doi.org/10.1016/j.jbiomech.2016.07.035 doi: 10.1016/j.jbiomech.2016.07.035
![]() |
[46] | C. Gu, Smoothing spline ANOVA models, 2 Eds., New York, NY: Springer, 2013. https://doi.org/10.1007/978-1-4614-5369-7 |
[47] | H. Cardot, F. Ferraty, P. Sarda, Functional linear model, Stat. Probabil. Lett., 45 (1999), 11–22. https://doi.org/10.1016/S0167-7152(99)00036-X |
[48] |
J. Jiang, H. Lin, Q. Zhong, Y. Li, Analysis of multivariate non-gaussian functional data: A semiparametric latent process approach, J. Multivariate Anal., 189 (2022), 104888. https://doi.org/10.1016/j.jmva.2021.104888 doi: 10.1016/j.jmva.2021.104888
![]() |
[49] |
N. Ling, P. Vieu, Nonparametric modelling for functional data: Selected survey and tracks for future, Statistics, 52 (2018), 934–949. https://doi.org/10.1080/02331888.2018.1487120 doi: 10.1080/02331888.2018.1487120
![]() |
[50] |
R. S. Zoh, Y. Luan, L. Xue, D. B. Allison, C. D. Tekwe, A Bayesian semi-parametric scalar-on-function regression with measurement error using instrumental variables, Stat. Med., 43 (2024), 4043–4054. https://doi.org/10.1002/sim.10165 doi: 10.1002/sim.10165
![]() |
[51] |
J. Dauxois, A. Pousse, Y. Romain, Asymptotic theory for the principal component analysis of a vector random function: Some applications to statistical inference, J. Multivariate Anal., 12 (1982), 136–154. https://doi.org/10.1016/0047-259X(82)90088-4 doi: 10.1016/0047-259X(82)90088-4
![]() |
[52] | I. Jolliffe, Principal component analysis, Springer, 2005. https://doi.org/10.1002/0470013192.bsa501 |
[53] |
S. Friedrich, E. Brunner, M. Pauly, Permuting longitudinal data in spite of the dependencies, J. Multivariate Anal., 153 (2017), 255–265. https://doi.org/10.1016/j.jmva.2016.10.004 doi: 10.1016/j.jmva.2016.10.004
![]() |
[54] |
H. Guillermou, C. Abraham, I. Annesi-Maesano, N. Molinari, Multiple hypothesis testing in allergy and hypersensitivity diseases investigation: A pedagogical perspective, JAHD, 3-4 (2024), 100014. https://doi.org/10.1016/j.jahd.2024.100014 doi: 10.1016/j.jahd.2024.100014
![]() |
[55] |
S. K. Sarkar, Z. Zhao, Local false discovery rate based methods for multiple testing of one-way classified hypotheses, Electron. J. Stat., 16 (2022), 6043–6085. https://doi.org/10.1214/22-EJS2080 doi: 10.1214/22-EJS2080
![]() |
[56] | S. Friedrich, F. Konietschke, M. Pauly, Resampling-based analysis of multivariate data and repeated measures designs with the R package MANOVA.RM, R J., 11 (2019), 380–400. |
[57] |
E. Orantes-González, J. Heredia-Jiménez, Does schoolbag carriage equally affect obese/overweight and healthy-weight children? Appl. Ergon., 90 (2021), 103236. https://doi.org/10.1016/j.apergo.2020.103236 doi: 10.1016/j.apergo.2020.103236
![]() |