{'fraction_non_alphanumeric': 0.027219352993624144, 'fraction_numerical': 0.010160994566223104, 'mean_word_length': 4.914825870646766, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Gaussian Tree Constraints Applied to Acoustic Linguistic Functional Data

Evolutionary models of languages are usually considered to take the form of trees. With the development of so-called tree constraints the plausibility of the tree model assumptions can be addressed by checking whether the moments of observed variables lie within regions consistent with trees. In our linguistic application, the data set comprises acoustic samples (audio recordings) from speakers of five Romance languages or dialects. We wish to assess these functional data for compatibility with a hereditary tree model at the language level. A novel combination of canonical function analysis (CFA) with a separable covariance structure provides a method for generating a representative basis for the data. This resulting basis is formed of components which emphasize language differences whilst maintaining the integrity of the observational language-groupings. A previously unexploited Gaussian tree constraint is then applied to component-by-component projections of the data to investigate adherence to an evolutionary tree. The results indicate that while a tree model is unlikely to be suitable for modeling all aspects of the acoustic linguistic data, certain features of the spoken Romance languages highlighted by the separable-CFA basis may indeed be suitably modeled as a tree.


INTRODUCTION
As interest in functional data analysis (FDA) has increased, so has the development of theory (Ramsay and Silverman, 2005;Horváth and Kokoszka, 2012) and the range of applications (e.g. brain imaging (Sørensen et al. (2013)), climatology (Besse et al. (2000)), and medical research (Ratcliffe et al. (2002))). This has been facilitated in part through better access to functional data and through greater availability of computational power for analyses. The use of FDA in statistical phonetics has recently attracted attention (e.g. Koenig et al. (2008); Mooshammer (2010)). Such analyses, which involve acoustic functional data, have provided particularly promising and interesting results in a diverse range of settings. For example, Grabe et al. (2007) uses a polynomial basis expansion to examine pitch variation in English, while Aston et al. (2010) investigates Qiang, a Sino-Tibetan language, and finds previously unidentified gender differences amongst speakers via a functional principal components based modeling approach.
The acoustic structure of spoken words can be used to investigate areas of linguistic interest in a similar way that discrete (alphabetic) representations of speech have been utilized to make cross-language comparisons, and more recently inferences regarding proto languages (Hock (1991), Bouchard-Côté et al. (2013)). The differences and similarities between spoken languages suggest that any meaningful functional observations taken across languages are unlikely to be independently, identically distributed. As such it is probable that the language relationships form a tree or network structure, which may be informative about possible historical developments of these languages. If this alternative (acoustic) approach can be used to corroborate known and uncontroversial language relationships, then our methods offer great potential for less certain language relationships. For instance, this would be useful for languages where there are few historical records but in which inference of a family tree is reasonably supported by the contemporary data (e.g. African language families), or alternatively, in cases where reconstruction of a family tree is disputed, such as Greenberg's classification of native American languages (Bolnick et al. (2004)).
Relationships between languages have long been described as phylogenetic trees constructed using linguistic factors (e.g. Schleicher (1860)) where all non-leaf variables are unobserved and represent features of the past languages before their divergence. Greenberg (1954) developed some of the first quantitative methods which were used to investigate evolutionary relationships between languages. More recently there have been large scale attempts to reconstruct trees or networks of languages (e.g. Nakhleh et al. (2005) for the Indo-European language family, Nicholls and Ryder (2011) for the Semitic language family). Some researchers have shifted away from describing the evolutionary language relationships via trees toward using networks (for example, Forster and Toth (2003), Nelson-Sathi et al. (2011)). However, trees have a somewhat more natural interpretation in terms of evolutionary structure, and assessing the suitability of a tree model for language data would therefore be of interest to researchers in linguistics. The focus of this paper is to examine functional acoustic data from speakers of five Romance languages (French, Italian, Portuguese, American Spanish, and Iberian Spanish) to provide insight at an exploratory level as to whether a tree may be adequate for describing certain features of these language relationships.
To address questions of tree-amenability we appeal to the notion of tree constraints. The theory of tree constraints is embedded in the area of algebraic statistics, a field which has a significant recent literature related to phylogenetics (e.g. Sturmfels and Sullivant (2005), Allman and Rhodes (2008)). It has been known for some time that covariance functions of data on observed variables respecting an evolutionary tree must obey particular algebraic and semi-algebraic constraints, e.g. Settimi and Smith (2000). Recently these have become much better understood (for example Drton and Sullivant (2007), Allman et al. (2009), Allman et al. (2012)) and fully characterized in some cases (e.g. the binary case Zwiernik and Smith (2011), Zwiernik and Smith (2012)). In this paper, building on developments in Shiers and Smith (2012), a Gaussian analogue of the binary tree constraints is applied to the covariances of component-by-component projections of the data.
By considering the data component-wise, a more realistic and nuanced analysis can be performed which permits some observed features of linguistic data to be tree-amenable and indicates that that for others that any evolutionary tree is unlikely to provide a good explanation.
Section 2 describes the data set and the preprocessing in preparing an audio recording for FDA. In a similar spirit to the work on object data by Wang and Marron (2007), for the application to Romance languages presented here, the observation units of interest are two-dimensional functional data objects known as spectrograms (time-frequency descriptions of the data). These are formed from transformations of audio recordings of people speaking single words. When regarded as functional data, observations are in fact stored as high dimensional objects. Therefore, in Section 3 tools from FDA are employed to transform high dimensional speaker data to a lower dimension. This is achieved through the novel approach of using between-language covariance as described in canonical function analysis (CFA) and combining it with a tensor decomposable covariance structure.
Having achieved the required reduction in dimension, Section 4 provides a brief summary of tree constraints. A fundamental but yet to be exploited constraint for use with Gaussian data is then introduced. Statistics associated with the violation of or adherence to the constraint are then constructed from the acoustic language data to answer the question of tree suitability. In Section 4.3, in preparation for use with this Gaussian tree constraint, we describe the construction of a between-language cross-covariance matrix using the scores (projections) of the acoustic data. In Section 4.4, the general effectiveness of the Gaussian tree constraint is investigated via simulations to assess its ability to correctly accept or reject tree-amenability. Section 4.5 entails further simulation, tailored to mimic aspects of the acoustic data, so the results are more immediately relevant to the setting. Section 5 addresses the ultimate aim of the study: to determine whether any evolutionary relationships between the Romance words represented in the acoustic data set can be described by a Gaussian evolutionary tree model. These tree constraints are then used to explore tree-amenability for each of the components of the chosen basis. It is found that a subset of the components adhere to the tree constraint. This suggests that some features of the acoustic linguistic data which distinguish between languages could have evolved in a tree-like manner whilst others have not. This is consistent with the current understanding that the development of Romance languages has been complex, involving much cross-language interaction (Harris and Vincent (1988)) in addition to a historically well-documented common origin from Latin. Thus, attempting to fit a tree model to the entire data set is likely to be misguided based on the empirical evidence presented here. More appropriately, a different model can be fitted for each component, using the tree constraints to indicate whether to restrict the space of models to trees.

ACOUSTIC FUNCTIONAL DATA SET
The data set of interest comprises audio recordings originating from speakers of one of five different Romance languages: French, Italian, Portuguese, Spanish (American), and Spanish (Iberian)while two dialects of Spanish are being used in this study, they are treated as different spoken languages in this analysis as we are interested in pronunciation rather than textual representation, the difference between "dialect" and "language" being a matter of degree of difference rather than an absolute quantitative difference. Each recording is of some individual saying an integer from 'one' to 'ten' in their particular language. Recordings were made at a rate of 16000 samples per second and a resolution of 16 bits. In total there are 219 word recordings and each can be classified by the language, the gender of the speaker and the number being spoken. Observations of the same word being spoken in different languages are treated as sharing the same word attribute. For example the word 'four' includes recordings of 'quatre' (French) and 'quattro' (Italian) as well the word 'four' in other languages. Integers were chosen because these have no ambiguity in terms of translation making comparison of their use across languages straightforward. Furthermore, the cardinals 'one' to 'ten' of Romance languages (among many other words) stem from shared Latin forms (Price (1992)). This suggests that these words might also be suitable when comparing languages acoustically.
In this paper the observations are modeled as functional data as is becoming increasingly common in studies involving sound recordings (e.g. Holan et al. (2010)). Such models make the reasonable assumption that the data have been obtained by observing an underlying function at finitely many discrete points along a continuum, and that this underlying function is smooth (i.e. a certain number of derivatives exist).
The overall duration of a word can vary significantly per speaker as can the timings of intraword elements (for instance syllables). Thus, to adjust for these differences all observations within a word grouping undergo registration (also known as alignment or warping, see Ramsay and Silverman (2005); Lucero and Koenig (2000)). A short-time Fourier transform is taken of each audio recording to produce a spectrogram. A spectrogram is a two-dimensional representation of audio signal energy intensity in frequency-time space (Fulop (2011)). Spectrograms are a natural choice for representing power with functional data, though approaches such as Mel-frequency cepstra (Davis and Mermelstein (1980)) can provide possible alternative representations. As part of the preprocessing, the standardization of word duration results in the time dimension being measured in generic time units. The value stored at a frequency-time point is a function of the power (or amplitude), with frequencies binned every 100Hz up to the Nyquist frequency of 8000Hz. The resulting spectrograms are stored as matrices of 81 frequency by 100 time points. Figure 1 is the spectrogram of a female French speaker saying the word 'quatre'. Broadly, this interpolated plot indicates that there is greater power in the lower frequencies, and that the beginning and the end portions of the standardized time period are quieter.

Notation
The underlying function of each spectrogram is denoted x d,g l,m (f, t) with the two dimensions f and t referring to frequency and time respectively. Recall that each spectrogram is derived from a spoken word -the subscripts and superscripts encode observational information: l = 1, . . . , n l denotes the language being spoken; d = 1, . . . , n d indicates the word being spoken; m = 1, . . . , m ld is a counter where m ld is the number of observations of word d from language l; g refers to the gender of the speaker.
It is well documented that there are differences in the acoustics of male and female speakers which go beyond a simple shift in the spoken frequencies. For instance Nittrouer et al. (1990); Pépiot (2013). Parris and Carey (1996) present a statistical method for discriminating between speaker gender of short acoustic recordings. In their analysis of seven Indo-European languages (of which Romance is a subset), gender was correctly identified on average 98% of the time.
This suggests that there are commonalities in acoustic gender differences across Indo-European languages. In light of this result, it is judged that gender should be adjusted for at the macro level: wherex g is the difference between the mean of all samples with gender g and the mean of all samples. Henceforth it will be the gender adjusted function that will be the object of interest in this paper.
The mean spectrograms for language l, word d are defined in (1), for language l in (2), and the where m l· = n d d=1 m ld , m ·· = n = n l l=1 m l· , and for t ∈ T , f ∈ F. The parameters m ·· and n will be used interchangeably depending on whether summation is being emphasized.

GROUP BASED PROJECTIONS OF FUNCTIONAL DATA
Dimension reduction is a well-studied area of statistics with tools such as principal component analysis (PCA) and multidimensional scaling (e.g. see Jolliffe (2002); Cox and Cox (2010)) having widespread use. Functional counterparts of such techniques have also been formulated, for example functional principal component analysis (FPCA) (Castro et al. (1986); Rice and Silverman (1991); Yao and Lee (2006)), and functional multidimensional scaling (Mizuta (2006)).
The acoustic data presented in this paper benefit from a dimension reduction in two main ways.
First and foremost, dimension reduction provides a route to feature extraction whilst also reducing unwanted noise. Second, if subsequent to the reduction it is found that n ≥ p then techniques which make use of inverse covariances can be implemented straightforwardly. If this is not so, standard estimates will produce singular sample covariance matrices. Of course these benefits must be balanced against potential information loss from the data reduction. One approach to feature extraction which mitigates against this loss is to find an ordered basis which prioritizes one or more characteristics of interest. Thus by projecting data onto the first few components of such a basis the most prominent aspects of the data are retained whilst what remains is treated as noise. In the cases of PCA and FPCA, the dimension reduction is optimized so as to efficiently capture modes of variation. Such techniques are often used in linguistic and semantic analyses, for example Lee et al. (2001) and Wenyin et al. (2010). However, as our focus is on macro-language comparisons, we argue that the feature of interest is the between-to within-language covariance, and it is this which should directly inform the method selected to construct a basis. When data is known a priori to be grouped then canonical function analysis (CFA) and its multivariate analogue canonical variate analysis (CVA) are standard techniques implemented to select variables to discriminate between groups. These tools are therefore the starting points for our analysis.

Canonical Function Analysis
Here we present CFA as a tool for FDA to produce a basis which maximizes between-to withingroup variation (subject to the basis component functions being uncorrelated) with the intention of achieving an efficient dimension reduction. The finer details of CFA can be found in Kiiveri (1992). To illustrate CFA, consider a set of one-dimensional functional data denoting the betweengroup covariance function as B(u 1 , u 2 ) and within-group covariance function as W (u 1 , u 2 ) (with u 1 , u 2 ∈ U, U ⊂ R) where both functions are considered to be bounded and piecewise continuous.
The aim is to identify canonical functions h q (u) such that between-group variation is maximized relative to within-group variation under the restriction that each canonical function is uncorrelated to every other. This is expressed as a generalized eigen-equation which simplifies to: with eigenvalues λ q and eigenfunctions h q . There may be countably infinite solutions to equation (4) but in discretized estimation, only a maximum of s (say) will have non-zero λ q . Pairs of canonical functions and real numbers (h 1 (u), λ 1 ), . . . , (h s (u), λ s ) can be found by solving (4) numerically, where λ 1 , . . . , λ s is a monotone decreasing sequence. An r-dimensional projection of the data is obtained using the first r canonical functions, and this projection is such that the between-to within-group covariance is maximally retained.

Separable Covariance Functions
Recall that the Romance data set comprises spectrograms which have time and frequency directions. A straightforward model for describing how these directions interact is that of separable covariance.
The assumption underpinning this model can be encapsulated as there being no dependency between the (standardized) time and frequency of the data. This is, of course, a significant simplification of the likely underlying model. However, particularly for computational considerations, this assumption can very useful. Recall that a covariance function C is said to be separable if: where C f and C t are functions only of their arguments. The factored covariances provide an understanding of how frequency or time dimensions of the spectrograms vary when the other has been averaged out. However, the main purpose of the assumption becomes apparent for the Romance data set subsequently (as described in Section 3.3.1) when use of the separable model of covariance overcomes the challenge of covariance rank deficiency.
Under the separable covariance assumption for two-dimensional data the CFA optimality equation equivalent to (4) is: It can be easily shown that the solutions to this equation can be obtained as the product of the solutions to two CFAs performed on the frequency and time covariances separately. Thus given any canonical function pairs (h q f (f 2 ), λ q f ) and (h qt (t 2 ), λ qt ) from a frequency and time CFA respectively, the products provide a solution to (5): Moreover, any solution to (5) can be obtained from such products. This result is useful when proceeding to obtain numerical solutions to (5).

Canonical Variate Analysis
Although less frequently implemented than PCA, the theory of CVA as a multivariate tool has been well developed in Krzanowski (1990, Chapter 11). However, beyond being a purely multivariate technique, CVA can also be used with functional data as an approximation to CFA as is presented in Kiiveri (1992). The technicalities of implementing CVA do not differ whether in a functional or multivariate setting, although it is sometimes necessary to interpret their outputs differently, as is encountered with other tools (e.g. Fervaha and Remington (2012)).
In practice spectrograms are often discretized representations of underlying functions, and so each function x d l,m is instead given by a matrix X d l,m with time-frequency dimensions n f × n t (i.e. the number of sample points of the frequency and time). These finite approximations tend to be high dimensional and so the question of dimension reduction is pertinent. By concatenating the rows of these matrix representations of spectrograms the data corresponds to the vector description of CVA. As CVA considers each covariance entry independently of its adjacent values, this does not affect the implementation of CVA. The only notable downside of concatenation is that it can obscure visual representation and description of the data.
Recall that the aim of CVA is to find successive uncorrelated vectors a that form linear combi- that maximize the ratio of the between-groups covariance (B) to the within-groups covariance (W ). In the context of the paper, B describes the variation between the per-language mean spectrograms and the grand mean spectrogram, whereas the W describes the variation between individual observations and the associated per-language mean spectrograms.
Finding the optimal a is equivalent to solving (W −1 B − λI)a T = 0 where λ ∈ R. This reduces to performing an eigenanalysis on W −1 B -the eigenvector corresponding to the largest eigenvalue is the optimal a. As with CFA, canonical pairs (a T r , λ r ) are sought. These are found through a full eigenanalysis of W −1 B such that λ 1 > λ 2 > . . . > λ s > 0 where s = min(p, n l − 1) is the number of non-zero eigenvalues of W −1 B. Thus (a T r , λ r ) produces the rth greatest ratio of between-to within-language variability. Hence, the optimal projection to r dimensions requires only (a T 1 , λ 1 ), . . . , (a T r , λ r ). Theoretically, CVA is designed for use with Gaussian data and assumes that within-group covariance is equal across groups. If these assumptions hold true then CVA is optimal for identifying modes of variability which distinguish between groups.

Separable-CVA
As mentioned in Section 3.2, the overall solutions to a CFA optimality problem with a separable covariance structure can be found as the product of solutions to CFAs of the decomposed covari-ance functions. We propose combining a tensor decomposable covariance structure with CVA in order to obtain numerical solutions to the decomposition of the separable-CFAs. This, when taking products, also gives solutions to the overall CFA. While separable covariance structures have been adopted elsewhere in the literature (e.g. Jones and Moriarty (2012); Aston and Kirch (2012)), this is a novel approach for both CVA and CFA. Even though the assumption behind separable covariance is strong, the accuracy of the assumption for CVA and CFA only impacts on basis efficiency not basis validity. If the data is far from separable, then simply a higher number of dimensions will be needed to retain the same amount of information.
The main purpose of assuming a tensor-decomposable covariance structure is to overcome the obstacle of rank-deficient sample covariance matrices caused by the length of the observations exceeding the number of observations (i.e. p > n). This is not just a problem with the Romance speaker data set but is commonly encountered with functional data sets due to their often highdimensionality (e.g. Long et al. (2005)). Rank deficiency obstructs using CVA to obtain numerical solutions to CFA. Theoretically in CFA an inverse function W −1 is neither required nor is usually bounded, whereas in CVA W −1 is needed for the eigenanalysis of W −1 B but cannot be obtained because in this case W is singular.
In the observational matrix setting, C is separable if: where ⊗ is the standard Kronecker product. Using known results of the Kronecker product (see Lancaster and Tismenetsky (1985) for example), the separability assumption in the multivariate setting implies: where the estimates of separate within-and between-language covariance matrices in the frequency for the time direction. Treating each frequency and time sample as a separate observation leads to the product covariance matrices W and B having higher ranks than previously. Explicitly, for to be nonsingular, we need that nn f ≥ n t and nn t ≥ n f . This is equivalent nt) . This contrasts to the previous condition n ≥ p = n f n t . So the new requirement is usually significantly more relaxed, and CVA can often then be implemented.
and the Kronecker product of the corresponding eigenvectors results in matrices denoted (c 1 , c 2 , . . . , c n f nt ) of size n f × n t , solving the overall CVA. It should be noted that while this basis defined is based on an assumption of separability, it nevertheless provides a complete basis of the space. So although when separability does not hold the basis is less efficient and is rather longer than it needs to be, the basis is still valid. For further details, see an analogous argument for separable PCA in Aston and Kirch (2012).

Application of CVA as an Approximation to CFA
As motivated, the separable-CVA is used to approximate the separable-CFA of the Romance languages data to achieve a dimension reduction based on components which maximize between-to within-language variability. This is a suitable approximation to make as the functional spectrograms have been sufficiently densely sampled during discretization.
[ Figure 2 about here.] Interpolated plots of the between-and within-language covariances for both the time and frequency directions of the spectrograms are given in Figure 2. The B f plot displays positive covariance which increases as the frequency increases. The time covariance B t is positive throughout with an internal minimum and the highest covariances corresponding to beginning and end time points. These particularly high time corner covariances could indicate genuine similarities at the beginning and ends of spoken words or could be an artefact of the spectrogram registration. The covariances W t and W f both have a ridge along the diagonal which reassuringly suggests that similar time and frequency points have strong covariances within languages.
When selecting a dimension r to project to, it is unusual to have anything but an arbitrary albeit sensible method for selecting r. However, in some acoustic contexts (e.g. Hadjipantelis et al. (2012)) thresholds can be proposed based on sounds which are audible to humans. Otherwise, equivalent techniques to those employed with PCA (Jolliffe (2002)) can be used. For this linguistic study Figure 3 shows the cumulative variation explained by selecting particular numbers of components. This indicates that almost 95% of the between-to within-language variance can be explained by a single component, and an additional two components take this figure to over 96%.
[ Figure 3 about here.] Once a dimension r has been selected, each observation from the language data can be pro- . . , c r and wherẽ X is 1 × p and formed by concatenating the n f rows of length n t of the observation X. The sub and superscripts of the observations are omitted solely for notational clarity.
[ Figure 4 about here.] To demonstrate the effectiveness of projection of the spectrograms in even two dimensions, the projections of the means of the word observations are plotted in Figure 4. The results of the separable-CVA are encouraging as there are clear groupings of the projected word means from the same languages. Furthermore, the first two dimensions appear to reflect aspects known about the languages. The second dimension seems to distinguish the four distinct languages, and then the first dimension is able to separate the Spanish dialects (while also providing clear distinctions between some of the other languages). Given that the acoustic data are undoubtedly noisy, this indicates the effectiveness of separable-CVA at selecting components which discriminate on a group basis. Note that whilst CVA operates on all languages simultaneously rather than in a pairwise manner, this does not necessarily imply languages in close proximity post-projection share particular acoustic features. However, this can be examined in more detail, for example through studying the Hadamand matrices of each projection as is illustrated in Figure 8.

TREE-AMENABILITY
The aim of this study is to examine the suitability of a tree for acoustic functional data. It is of particular interest to know whether there are certain features of these spoken Romance languages which could have developed over time in the manner of an evolutionary tree. CFA (or CVA) is compatible with this aim as it effectively identifies components with features which distinguish between languages. Thus a projection to r dimensions provides r different combinations of characteristics which could potentially be modeled as a tree. In order to help assess the plausibility of these r hypotheses we turn to the topic of tree constraints. These are algebraic and semi-algebraic constraints which mathematically must be obeyed if a data set corresponds to the leaves of a Bayesian network which is a tree. Note here that the non-leaf vertices correspond to various unobserved ancestral languages. While there has been considerable recent interest in understanding graphical structures, including trees, in very general settings (Loh and Wainwright, 2013), these approaches are design for settings when all nodes (both internal and leaves) are observed. Here we are concerned with determining constraints when the internal nodes are unobserved. Examining whether these constraints are respected for a given data set determines whether that data set adheres to a tree structure (i.e. whether the data set is tree-amenable). Provided with r component-by-component projections, a diagnostic for tree-amenability can be applied to each (see e.g. Shiers and Smith (2012) for a binary example). This can then be used as an exploratory tool to give insight as to whether the distinguishing characteristics captured by a component could have an evolutionary tree structure. We begin by giving a brief overview of the concept of tree constraints and how they form a natural choice for use with data sets such as the acoustic recordings in this study.

Constraints on Observed Covariances Respecting Trees with Hidden
Variables Tree constraints are useful for studies investigating evolutionary relationships between both observable and unobservable variables. These evolutionary relationships can be expressed graphically as phylogenetic trees where traditionally the objects of interest have been biological species (e.g. Felsenstein (1983)), however, this idea naturally extends to other fields such as linguistics (e.g. Dunn et al. (2005)). Phylogenetic trees are of particular statistical interest when formally considered as Bayesian networks describing conditional independence relationships (e.g. Dutkowski and Tiuryn (2007)). In the case of discrete graphs in which all variables are assumed to be observed, Loh and Wainwright (2013) make use of precision matrices to assess graph structure. However, typically Bayesian networks have both observed and hidden variables and these are much more complex to analyze. Furthermore, the data set in question is functional and thus the analyses are embedded in Gaussian assumptions and observations treated as Gaussian processes as opposed to individual points. This is where Gaussian tree constraints can provide useful insight into the difficult question of tree-amenability.
A graph T is a tree if there exists a unique path between any two connected vertices. The trees of interest are those with interior hidden nodes H ∈ H and manifest leaf nodes X ∈ X (denoted by white and black circles respectively). The observed leaf variables can be thought of as contemporary languages, and the latent interior variables as past versions of languages. Attention is restricted to strictly trivalent trees, i.e. trees where no interior node has more than three adjacent nodes. Any tree from this class with more than three leaves can be expressed as a bifurcating tree, a common model in phylogenetics (e.g. Felsenstein (1978)). The tripod tree (shown in Figure 5) describes any such relationship between three observed variables. To see this, simply note that any three connected leaves of a strictly trivalent tree share a unique interior node, and that the conditional independence of these four nodes is minimally described by the tripod tree. Hence, the tripod tree is a fundamental component for analyses involving tree constraints.
[ Figure 5 about here.] The distributions associated with moments of the observed variables of such tree models are known for some settings. The case of univariate binary random variables was expanded in Settimi and Smith (2000) and has recently attracted considerable interest (for example Allman et al. (2012)). These advances have encouraged some authors to proceed to use the known geometry of these spaces to support inference and learning over the space of tree models (see Drton and Sullivant (2007)). These focus on the polynomial constraints that are implicit in these models. Zwiernik and Smith (2011) fully characterized the space for binary trees by extending the understanding beyond algebraic relationships to also include the active semialgebraic constraints, whilst Allman et al. (2012) broadened the results to a k-state n-pod tree. Our interest is in the use of the constraints which define these geometries for assessing whether data could have originated from a tree model, in particular whether the spoken languages in the acoustic data set could be modeled as a phylogenetic tree.

A Fundamental Constraint on the Covariance of Gaussian Tree Models
Most of the results in the literature of tree constraints are only applicable in settings with binary random variables. Alternative constraints specific to Gaussian tree models are required in order to perform equivalent tree-amenability diagnostics for Gaussian functional data such as the Romance data set. Here we describe a fundamental univariate Gaussian tree constraint which has not yet been exploited for the purpose of investigating tree-amenability.
Consider the precision matrix Σ −1 related to the tripod tree in Figure 5 with one latent and three manifest Gaussian random variables. It is well known in the Gaussian setting that if X i is independent of X j conditional on all other observed variables then the corresponding entry Σ −1 ij = 0 (see Lauritzen (1996, Chapter 5) for example). So the precision matrix for the univariate Gaussian tripod tree has the form: where the final row/column relates to the interior hidden variable. The covariance matrix resulting from taking the inverse of the precision matrix can be expressed algebraically in terms of entries of Σ −1 . Using the resulting entries of the covariance Σ = [σ ij ] it is then straightforward to calculate that a necessary condition for the leaves of this tree to be the margin of the tripod tree given in Figure 5 is that We shall refer to this constraint as the positivity constraint. For a strictly trivalent tree with n l observed leaf nodes, there are n l 3 tripod trees that must be valid -one for each triple -and thus n l 3 such positivity constraints. Whilst there are further constraints imposed on the observed moment space of trees, the derivations are considerably more involved and not the emphasis of this applied paper. Given that applications of Gaussian tree-constraints are not apparent in the literature, there is still much to explore focusing solely on this fundamental Gaussian tree constraint.
These constraints can be used with the linguistic data set as a diagnostic for assessing treeamenability. More precisely, for the r component-by-component projections, each can be assessed for tree-amenability. Thus some components may be found to violate the tree constraints whereas others may satisfy them. This may suggest which attributes of the spoken languages have evolved in a tree-like manner and which have not. Note that although we have functional data, the formulation of this positivity constraint does not assume this, and thus the constraint is equally valid for multivariate Gaussian data.

Constructing a Suitable Covariance Statistic
In pursuit of assessing tree-amenability of the r projections of the Romance data using the positivity constraint σ ij σ ik σ jk ≥ 0, it is clear that a sample covariance of the scores must be constructed.
Recall that the relationships of interest in this study are at the language level and thus betweenlanguage covariances (each 5 × 5) are the appropriate statistics to produce, one for each of the r components. One approach to calculating the entries of these matrices is to treat the mean score of each word in a language as an observation and then measure the distance from the overall word mean projection. Then using appropriate weights, a between-language covariance matrix can be where recall n d is the number of unique words. Note that this between-group covariance differs from that used in the CVA -this is of the projected data, with word means treated as observations. This is a valid construction in the sense that (8) is an inner product (see Istratescu (1987) for instance). The sample matricesΣ Y i will be rank deficient if n l ≥ n d . Also, observe that if for at least one word d the number of observations is unequal across languages then the weighted word meanȳ i d differs from the unweighted version. This relaxes a zero-sum condition on the the rows or columns ofΣ Y i permitting the covariance matrix to be full rank. In the alternate case of a balanced observational design, full rank can be achieved through an alternative construction (for example adding the unweighted word means back to each language-word mean).

Simulation
Before analyzing the Romance language data, a short simulation is performed. The purpose of the simulation is to demonstrate the effectiveness of the tripod positivity constraint at identifying treeamenability. Also, by varying the simulation sample size the robustness of the positivity constraint is examined.
In order to mimic the linguistic data set, a tree with five leaves is used to generate data (see Figure 6). Standard structural equations are used (Bollen (1989)) adapted for a Gaussian tree to model the dependencies between nodes. The interior node H 1 is the root and this variable is simulated from a Gaussian distribution. Data on nodes adjacent to the root are generated as linear combinations of this root simulation with the addition of Gaussian noise. In a p-dimensional scenario, the root node is simulated as follows: where V h 1 is a p × p symmetric positive matrix, constructed as A + A T + pI p with entries of A generated independently from the continuous uniform distribution U (0, 1). 0 p is a p-vector of zeros, and I p is the p × p identity matrix. Subsequent nodes are simulated from previous ones, for example: Each entry of λ h 1 x 1 ∼ U (−a, a), a ∈ R. Within a repetition, all entries of λ h 1 x 1 are fixed for all observations, although do differ depending on the node transition being considered. Additionally, x 1 ∼ N p (0 p , b 2 I p ), b ∈ R, and x 1 is pairwise independent with all other random variables.
Following the directions of the arrows, the remaining variables are simulated equivalently.
[ Figure 6 about here.] A four-dimensional data set D is generated for all eight nodes of the tree though only the observed leaf nodes are of interest and are utilized with the tree constraint. Another data set D * is also generated using a corrupted version of the tree where each transition between nodes provides an opportunity for sign-reversal of entries of the data. Given this corruption can occur on any dimension of any of the data, the simulation is no longer from a Gaussian tree model.
Although higher dimensional data could be generated, under standard CVA only a maximum of four (n l − 1) eigenvalues would be non-zero (the use of separable CVA relaxes this assumption in our analysis). Thus there would not be much information gained from the increased dimension but the computational resource required would certainly be greater. Four sample sizes are investigated: 50000, 5000, 500, and 50. These are average sample sizes as the simulation is designed such that the number of observations differs across languages allowing the use of (8) for the reasons relating to rank noted in Section 4.3. Each simulated data set undergoes the same basis change using CVA to approximate CFA as described in Section 3.3. For a full analysis, all four dimensions are retained for projecting the data, and thus for each data set four covariance matrices are calculated using the construction given in (8) [ There are three notable features highlighted by Table A

Further Simulation Based on Characteristics of Observed Data
In order to better assess the fundamental Gaussian tree constraint when applied to the data, a further simulation is performed using characteristics of the acoustic data set (e.g. sample size, eigendecomposition, sample variance). This is achieved by generating data from cross-covariance matrices Σ Y (as in (8)) for which tree-amenability status is known. Then by reversing the eigenbasis projection as performed for the Romance data set, new observations are obtained. These new samples can then be assessed using the fundamental Gaussian tree constraint, and it can be assessed whether the source Σ Y being tree-amenable or not has an affect on the new sample being deemed tree-amenable.
Three scenarios are considered: (A) cross-covariance Σ Y which is tree-amenable; (B) Σ Y which is not tree-amenable; (C) Σ Y is tree-amenable 50% of the time independently for each component and replication. For (A) and (B) the Σ Y i calculated from the Romance data set are used (for the analysis presented in Section 5). Each of these 15 Σ Y i are thus known to be treeamenable or not for the fundamental Gaussian tree constraint, and furthermore, the number of 5 3 = 10 constraints satisfied is known and can be denoted That is, it satisfies the 10 positivity constraints comprising products of covariances relating to the 10 possible choices of three languages. Note that it is not possible for T (Σ Y i ) to be exactly 1, 2, 8 or 9 due to the appearance of each off-diagonal entry of Σ Y i in three of the ten combinations. For the first 15 Σ Y i from the actual data set, the percentage of tree-amenable samples out 1000 generated is recorded against T (Σ Y i ). The four Σ Y i which are tree-amenable provide higher proportions of tree-amenable samples than the remaining covariance matrices. Considering the full results (see Table A.2 in the Appendix), a positive correlation is seen between T (Σ Y i ) and the percentage of samples simulated from Σ Y i which are then found to be tree-amenable. Considering scenario (C), a mix of the best and worst performing of the actual Σ Y i are used, having respectively 35% and 7% sample tree-amenability. The simulation results in an overall 18% tree amenability across all sampled components, suggesting that proportionately a similar level of tree-amenability is retained even when tree and non-tree components are combined.
This second simulation approach closely resembles the linguistic data set both in structure and parameter values. Thus these simulations provide comfort that even for the relatively small sample size the analysis is able to produce reasonable results. Scenarios (A) and (B) provide evidence that components which truly satisfy a low number of constraints are unlikely to give a false positive in terms of tree-amenability. Usefully they also identify that this risk grows as the number of con-straints satisfied increases, and this can thus be considered when assessing the results. Scenario (C) demonstrates that tree-amenable components can be recovered even when combined with non-tree components, and furthermore, that these appear to be retained in the output almost proportionally to the input. The sampling in all these simulation scenarios may also provide a proxy for distributional results of tree-amenability, where the higher the proportion of tree-amenability in a sample the more robust the conclusion. This hypothesis is an open problem which we look to address in a subsequent paper through derivation of explicit moments of relevant parameters.

ASSESSING TREE-AMENABILITY OF ROMANCE LAN-GUAGES
As illustrated in Figure 4, an effective projection of the Romance data can be performed in even two dimensions. However, to account for a higher proportion of between-to within-language variability further components can be included for the projection. A threshold of cumulative betweenlanguage variability is set at 97.5% which for these data equates to including all components which account for at least 0.05% explanatory power. Application of this threshold provides a dimension reduction from 8100 to 15. Each one of these 15 components c 1 , . . . , c 15 accounts for some mode of variability between languages. Although the earlier components have high explanatory power, the latter components may isolate directions of variability which are of more interest from a linguistic perspective.
Applying the positivity constraint to each of the 15 component covariance matrices results in four of the components (c 1 , c 2 , c 4 , c 6 ) adhering to the positivity constraint. Recall that the simulations in Section 4.4 hint that for less explanatory components the constraint appears more likely to identify true tree components than false positives. This suggests that the identified treeamenable components are more likely to truly be so, whereas the rejected components (e.g. c 3 , c 5 ) may or may not be truly tree-amenable.
To develop an insight into which aspects of the languages are being identified by the separable- is some symmetry in power between Italian and Portuguese. This is exhibited well in frequency ranges 300-800Hz, 1000-1500Hz, 3500-4500Hz, and 6000-6500Hz which show reflections along the line of zero power. Each of these frequency ranges can be phonetically interpreted. The 300-800Hz range likely relates to the first formant F1 being different, due to vowel differences. The 1000-1500Hz range likely relates to nasality (Portuguese is more "nasal" than Italian, and therefore has less energy in this portion of the spectrum than Italian). The 3500-4500Hz range is in the region of the third formant and could correspond to differences in lip rounding between speakers of the languages. The variation at the highest frequencies, around 6000 Hz, are likely to be due to idiosyncratic differences in speakers (since humans cannot readily control speech frequencies in that range) or in the recordings (equipment or recording location). Examining a smaller range of the data that excludes these high frequencies does not affect the main results of the analyses (data not shown).
Considering the time perspective, it is clear that the interesting time ranges are approximately 1-20 and 70-100. This suggests that the differences in the earliest and latest portions may be particularly effective at separating the languages. These results appear robust to trimming of the data set as well as to standardization of the covariance (data not shown), which suggests that the analysis is identifying more than just a feature of the data registration.
These projections are particularly effective at indicating graphically the dominant features of components which are often obscured when displayed numerically. Of course more detailed analyses are required to determine whether these are general features of the languages or simply of the particular data set studied. However, in either case, this type of exploratory data analysis is clearly a helpful one to perform before any more detailed analysis has taken place, especially if based on the firm assumption that the data originates from a tree.
[ Figure 8 about here.] This application demonstrates a method for isolating and identifying distinguishing aspects of variability in acoustic functional data which may be of evolutionary interest. It shows that it is possible to identify prominent features which render particular components effective for distinguishing the language groups. However, it also highlights the challenge of precise physical interpretation of particular components, a task which appears notably more complex due to having both a time and a frequency dimension. It would be of interest to express these differences back in the sound domain, although given the difficulties in inverting spectrograms to sound, this is not a trivial task.
However, it is the subject of ongoing work, including experiments with other parametric acoustic representations that are more easily inverted.

DISCUSSION
This paper presents a method for assessing tree-amenability of Gaussian functional data via Gaussian tree constraints. Through simulation the capability of the positivity constraint and its level of robustness to sample size have been demonstrated. Application of this fundamental tree constraint has indicated that the implementation is practically feasible. Moreover, it has shown that meaning-ful (albeit high-level) interpretations can be made making this particularly useful as an exploratory data tool.
The application to a linguistic data set comprising recordings from Romance language speakers provides several interesting preliminary results. First, it suggests that at least four components of the new basis are tree-amenable, and hence, the linguistic features these particular components represent may have a tree structure. Second, examining projections of the two tree-amenable components with highest explanatory power indicates that the second component is sufficient to distinguish languages, and including the first dimension separates dialects. Finally, closer study of the first two tree components indicates that broadly it is the beginning and end of utterances which are identified as effective language discriminants. Whilst these findings are somewhat speculative and the scope of the analysis is limited to a small set of words, it does provide a starting point for describing interesting tree and non-tree like features of grouped data which may otherwise be obscured. Furthermore, if applied on a larger scale to languages lacking established historical pathways, these techniques may offer fresh indications of the plausibility of different hypotheses concerning their historical development.
In the process of preparing data for tree constraint diagnostics, the concept of decomposable covariance structure has been combined with CVA to produce separable-CVA. Although not central to the aim of the paper, it is worth underlining its usefulness. As with standard CVA it identifies modes of variability which effectively distinguish grouped data, but furthermore, in many circumstances it overcomes the common functional data problem of variable dimension exceeding sample size. This makes separable-CVA an appealing tool to use in practice and one which worked well for separating these languages (Figure 4).
Whilst the range of techniques applied in this study can be carried across to other functional data sets, it is important to consider the assumptions which underpin the methods employed. Although no data set will truly obey all of the underlying assumptions, it is important nonetheless to understand the purpose of and the consequences of not satisfying each.
The theory of CVA assumes that the data are Gaussian and that the within-language covariances are equal. Considering the acoustic data, Gaussianity is unlikely to strictly hold. However, it may be sufficiently close to doing so and the first and second order moments may still answer questions about language relationships. If there is information available that the Gaussianity assumption is violated (as may well be evident in larger studies than the one presented here), then applying a copula transformation (see Genest et al. (1995)) will provide us with transformed variables that are marginally Gaussian (a necessary condition for joint Gaussianity). The use of such copula techniques would broaden the scope of these diagnostics and shall be described formally in future work.
In practice the within-language covariances are not likely to be equal (particularly given the sample size of the study). Nevertheless, pooling the covariances will emphasize shared commonality of covariances, and pooling may produce a more stable estimate than otherwise would be possible. Importantly, in circumstances where either assumption is violated CVA still produces a valid basis, the downside being it will not necessarily be as efficient at capturing the variability thus hindering the dimension reduction.
Covariance separability may also be unrealistic as there is often some correlation between the two separated dimensions (e.g. frequency and time), but in practice this is not a problem. CVA is an optimality tool and thus deviations from the separability assumption only decrease the efficacy of the basis obtained. That is, if covariance separability holds perfectly CVA determines an ordered basis which linearly maximizes between-to within-language variation as intended. Otherwise the ordered basis is suboptimal, the consequence being more dimensions may be required to capture a suitable proportion of the variance.
Therefore, although the three stated assumptions are not guaranteed to hold, infringement of any of them may impair performance but is unlikely to completely negate such an exploratory analysis. This permits a greater range of data sets for which a similar study could be applicable.
Aside from improving performance by selecting a data set which better meets the assumptions outlined above, the application to the Romance languages could instead be improved through a larger data set with greater breadth of words. Currently, the between-covariances of the projected data are only calculated using the means of the ten unique word observations. Furthermore, some languages have more recordings than others. By increasing the number of recordings for all groups, whilst also including a wider variety words, the outcomes of the analyses are likely to be more robust, with the caveat that inclusion of new words must be carefully considered to ensure they are linguistically suitable across all languages.
In this paper we check only some of the conditions necessary for consistency with a phylogenetic tree. However, it is possible to enhance the power of such diagnostics to include all such checks of possible violations of phylogenetic tree constraints. One method is to utilise the 4-point conditions based on the metric properties noted in Buneman (1971). If these conditions are satisfied then the tree structure can be uniquely identified. For one particular constraint based on tetrads, the work of Bollen and Ting (1993) gives a basis for possible test statistics and Drton et al. (2008) provides some useful distributional results for particular moments. An exploratory approach could use an extension of the binary graphical tree diagnostics described in Shiers and Smith (2012) to the Gaussian domain and is currently under investigation. Of course to fully and formally evaluate these methods it will be necessary to evaluate the probabilistic properties of our methods.
We are currently investigating a Bayesian approach where replicates are generated from Wishart covariance matrices, with properties derived from the projected data. These then provide posterior distributions across a graphical model space (Atay-Kayis and Massam (2005)). For graphs with high posterior probability, trees with hidden variables could then be induced. As with any formal testing framework, it will be important to consider potential effects of multiple comparisons, and to adjust appropriately for them (see Howell (2012, Chapter 12)).
Each of these methods could enrich the analysis, as the proposed trees could be compared to existing knowledge about the relationships between the languages. Contrasting results could even suggest new directions of research using traditional linguistic methods.
There is much potential for applications in fields beyond linguistics to make use of additional algebraic and semi-algebraic constraints. However, even existing constraints are currently underutilized. For example, despite the relatively simple derivation of the fundamental Gaussian tree constraint, this study marks its first use as a diagnostic for tree-amenability. The inclusion of further tree constraints, both in univariate and multivariate settings could provide further insight into phylogenetic relationships in applications. However, derivation and application of such constraints is non-trivial and is the subject of ongoing work.
APPENDIX: In general, the higher the value of T (Σ Y i ), the higher the percentage of tree-amenable components in the sample. Scenario (C) is a sample from a mixture of components 1 and 7, the resulting percentage of tree-amenability being half that of component 1 individually. Figure 1: Post-registration spectrogram of female French speaker saying 'quatre'. It can be seen that there is greater power in the lower frequencies, and that the very beginning and end of the word are unsurprisingly two of the quietest regions. Figure 2: Sample between-language and within-language covariances of speech data for frequency and time directions. There is a clear ridge along the diagonal of the within-group covariances indicating that similar times and frequencies are highly positively correlated. The higher correlations in the high frequencies of the within-group frequency covariance are associated with recordings at lower audio sampling rates capturing fewer details in these ranges of frequencies. Rerunning the analyses performed in Section 5 whilst excluding these higher frequencies produces broadly similar results. Figure 4: Two dimensional separable-CVA projection of means of word observations with ellipses representing the standard deviation around the mean of plotted points. The second dimension is effective at distinguishing at the language level. The first dimension then discriminates between the Spanish dialects. As each component can relate to a different set of speech qualities, the proximities of the languages can differ depending on the dimension projection. For example, in the second dimension, the nasality of the speech appears to contribute to the separation of Italian and Portuguese. This is further explored in Section 5.
H X 2 X 1 X 3 Figure 5: Tripod tree. A strictly trivalent tree with n l leaves contains n l 3 conditional independence relationships which can each be expressed as a tripod tree. Figure 6: Five leaf tree. Figure 7: For each sample size, interpolated plots are given of the estimated relationship between the explanatory power of components known to be tree-amenable and the chances that the components are correctly identified as being tree-amenable. This broadly suggests that the larger the sample size and/or the higher the explanatory power of a component, the higher the chances of a component being correctly identified as tree-amenable.