Skip to contents

This function returns the estimated motion variance calculated by the dynamic Bivariate Gaussian Bridges or dynamic Brownian Bridges

Usage

getMotionVariance(x,...)

Arguments

x

a DBBMM, DBBMMStack, DBBMMBurstStack,dBMvariance, dBMvarianceBurst, dBMvarianceStack, dynBGB or dBGBvariance object

...

Currently not implemented

Value

- a numeric vector of variances if a DBBMM, DBBMMBurstStack, dBMvariance or dBMvarianceBurst object is provided

- a list of variances per individual if a DBBMMStack or dBMvarianceStack object is provided

- a matrix of the orthogonal and parallel variances if dynBGB or dBGBvariance object is provided

Author

Bart Kranstauber & Anne Scharf

Examples

data(leroydbbmm)
data(dbbmmstack)
getMotionVariance(leroydbbmm)[1:50] ## with a DBBMM object
#>  [1]        NA        NA        NA        NA        NA        NA        NA
#>  [8]        NA        NA        NA 0.4319190 1.2668702 1.6479655 2.6416615
#> [15] 2.6559567 3.1656637 3.0886946 3.3808019 3.2778283 3.5915840 3.8014788
#> [22] 4.2408500 4.1745619 4.0447812 4.0637181 3.8881253 3.7385165 3.5424646
#> [29] 3.4088808 3.9744276 3.7983815 3.2225860 3.0644442 1.4951994 1.2049541
#> [36] 0.8849378 1.0263173 0.9058313 0.9877923 1.2483205 1.2421544 1.5202994
#> [43] 1.5037193 1.4333958 1.4333958 1.3575906 1.1031375 0.8771256 0.6833499
#> [50] 0.4228217
str(getMotionVariance(dbbmmstack)) ## with a DBBMMStack object
#> List of 2
#>  $ Leroy  : Named num [1:919] NA NA NA NA NA NA NA NA NA NA ...
#>   ..- attr(*, "names")= chr [1:919] "Leroy1" "Leroy2" "Leroy3" "Leroy4" ...
#>  $ Ricky.T: Named num [1:1000] NA NA NA NA NA NA NA NA NA NA ...
#>   ..- attr(*, "names")= chr [1:1000] "Ricky.T1" "Ricky.T2" "Ricky.T3" "Ricky.T4" ...