For symmetric *C* and *D* (with zero diagonals), we implement the Siemiatycki moment computation.

1 |

`Px` |
list of linear combinations |

`Py` |
list of linear combinations |

`n` |
sample size |

`mycoef` |
global coefficients we need for the function. |

`first` |
first permutation moment of quadratic form |

`second` |
second permutation moment |

`third` |
third permutation moment |

`fourth` |
fourth permutation moment |

YH Zhou, G Mayhew, Z Sun, X Xu, F Zou, FA Wright, 2013 *Space-time clustering and the permutation moments of quadratic forms*, Stat 2(1), 292-302

`quadp`

,`lincombfun`

.

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.