DynCorr: Dynamical Correlation

In functionaldata/tPACE: Functional Data Analysis and Empirical Dynamics

Dynamical Correlation

Description

Calculates the Dynamical Correlation for 2 paired dense regular functional data observed on the same grid.

Usage

DynCorr(x, y, t)

Arguments

x	a n by m matrix where rows representing subjects and columns representing measurements, missings are allowed.
y	a n by m matrix where rows representing subjects and columns representing measurements, missings are allowed.
t	a length m vector of time points where x,y are observed.

Value

A length n vector of individual dynamic correlations. The dynamic correlation can be obtained by taking average of this vector.

References

Dubin J A, Müller H G. Dynamical correlation for multivariate longitudinal data (2005). Journal of the American Statistical Association 100(471): 872-881. Liu S, Zhou Y, Palumbo R, Wang, J.L. (2016). Dynamical correlation: A new method for quantifying synchrony with multivariate intensive longitudinal data. Psychological methods 21(3): 291.

Examples

set.seed(10)
n=200             # sample size
t=seq(0,1,length.out=100)       # length of data
mu_quad_x=8*t^2-4*t+5
mu_quad_y=8*t^2-12*t+6
fun=rbind(rep(1,length(t)),-t,t^2)
z1=matrix(0,n,3)
z1[,1]=rnorm(n,0,2)
z1[,2]=rnorm(n,0,16/3)
z1[,3]=rnorm(n,0,4)
x1_quad_error=y1_quad_error=matrix(0,nrow=n,ncol=length(t))
for (i in 1:n){
  x1_quad_error[i,]=mu_quad_x+z1[i,]%*%fun+rnorm(length(t),0,0.01)
  y1_quad_error[i,]=mu_quad_y+2*z1[i,]%*%fun +rnorm(length(t),0,0.01)
}
dyn1_quad=DynCorr(x1_quad_error,y1_quad_error,t) 

functionaldata/tPACE documentation built on Aug. 16, 2022, 8:27 a.m.