Description Usage Arguments Details Value Author(s) See Also Examples

Mean and variance of tensors again tensors.

1 2 | ```
mean.tensor(x,along,...,na.rm=FALSE)
var.tensor(x,y=NULL,...,along,by=NULL,na.rm=FALSE,mark="'")
``` |

`x` |
(set of) dataset(s) of tensors represented by a tensor |

`y` |
a second dataset of connected tensors represented by a tensor |

`along` |
the indices indexing the datasets |

`...` |
here for generic compatibility with the compositions package |

`by` |
the indices indexing the set of datasets |

`na.rm` |
a boolean, if FALSE and missings are in the dataset a error is given. If TRUE pairwise exclusion is used. |

`mark` |
the to mark the second instance of indices in var(x,...) |

Let denote *a* the along dimension,
*i_1,...,i_k* and *j_1,...,j_l*
the data dimension, and
b the by dimension, then the mean is given by:

*
M^x_{bi_1,...,i_k}=\frac1{n}∑_{a}x_{abi_1,...,i_k}*

the covariance by

*
C_{ab i_1,...,i_kj_1,...,j_l}=\frac1{n-1}∑_{a}
(x_{abi_1,...,i_k}-x_{bi_1,...,i_k})(y_{abj_1,...,j_l}-M^y_{bj_1,...,j_l})*

and the variance by

*
V_{ab i_1,...,i_kj_1,...,j_l}=\frac1{n-1}∑_{a}
(x_{abi_1,...,i_k}-M^x_{bi_1,...,i_k})(x_{abi'_1,...,i'_k}-M^x_{bi'_1,...,i'_l})*

`mean` |
gives a tensor like x without the along dimensions representing the a mean over all tensors in the dataset. It is not necessary to have a by dimension since everything not in along is automatically treated parallel |

`var(x,...) ` |
Gives the covariate tensor representing the covariance of x and y. The data tensor indices of x any y should be different, since otherwise duplicated names exist in the result. |

`var(x,...) ` |
Gives the covariate representation of the variance of x. All data indices (i.e. all indices neither in by nor in along are duplicated. One with and one without the given mark. |

K.Gerald van den Boogaart

tensorA

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ```
d1 <- c(a=2,b=2)
d2 <- c("a'"=2,"b'"=2)
# a mean tensor:
m <- to.tensor(1:4,d1)
# a positive definite variance tensor:
V <- delta.tensor(d1)+one.tensor(c(d1,d2))
V
# Simulate Normally distributed tensors with these moments:
X <- (power.tensor(V,c("a","b"),c("a'","b'"),p=1/2) %e%
to.tensor(rnorm(1000*2*2),c(i=1000,d2))) + m
# The mean
mean.tensor(X,along="i")
# Full tensorial covariance:
var.tensor(X,along="i")
# Variance of the slices X[[b=1]] and X[[b=2]] :
var.tensor(X,along="i",by="b")
# Covariance of the slices X[[b=1]] and X[[b=2]] :
var.tensor(X[[b=1]],X[[a=~"a'",b=2]],along="i")
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.