# Domain association measures Measures of association between domains are computed as the association between the pairwise dissimilarities in the domains. Measures are: Pearson correlation, Spearman correlation, global Cronbach alpha, and Cronbach alpha for each subset of the domains. The function can also return the share of variance (R-square) of the dissimilarities in one domain that can be reproduced from the dissimilarities in the other domains. ```r dissdomassoc(domdiss, jointdiss = NULL, what = c("pearson","R2"), dnames=names(domdiss), weights=NULL, w.rank=FALSE) ``` ## Arguments - `domdiss`: List of symmetrical matrices or dist objects: the pairwise dissimilarities per domain. - `jointdiss`: `NULL` (default), matrix or dist object: pairwise dissimilarities for joint dimensions. - `what`: String or vector of strings: requested association measures among `'pearson'`, `'spearman'`, `'R2'` `'cronbach'`, `'cron.subsets'`. Can also be `'all'` for all measures but `'spearman'`. - `dnames`: String vector of length equal to number of domains. Names of domains. Default is `names(domdiss)`. - `weights`: Vector of non-negative weights. If `NULL`, no weights are applied. (Currently Cronbach measures ignore weights!). - `w.rank`: Logical. If weights are provided and `'spearman'` is selected, should we use weighted ranks. Caution: computation of weighted ranks considerably increases computation time! ## Details The `dissdomassoc` function computes the domain association measures proposed by Piccarreta (2017) . These are for each pair of domains the correlation (Pearson or Spearman) between the domain specific pairwise dissimilarities. When `'R2'` is requested, the function computes the share of variance of the dissimilarities in one domain that is reproduced by the dissimilarities in all other domains. The Cronbach alpha measures the coherence between the domains. With `'cron.subsets'`, Cronbach alpha is computed for the entire set of domains as well as for each possible subset of two or more domains. When a `jointdiss` distance matrix or object is provided, correlations of each individual with this joint domain are also computed. The R2 of the joint domain takes account of all domains. However, R2's for the domains remain unchanged, i.e. they ignore the joint domain. Likewise, `jointdiss` does not affect the Cronbach alpha. Spearman correlations are based on rank values of the distances. The computation of weighted ranks may take a while for large number of sequences (`> 500`, i.e. 124750 distances per domain). Therefore, weighted ranks are only used when explicitly requested by setting `w.rank = TRUE`. When `w.rank = FALSE` (default), Spearman is computed as the weighted Pearson correlation between non-weighted ranks. When `"pearson"` and/or `"spearman"` correlations are requested, p-values (probability to get a stronger correlation under the zero correlation assumption) are automatically computed. ## Returns An object of class `ddomassoc`, which is a list of tables of the requested association measures and tables of p-values of the Pearson and Spearman correlations when applicable. The `summary` method organizes the correlations and their p-values in table form. ## References Piccarreta (2017). Joint Sequence Analysis: Association and Clustering. **Sociological Methods and Research**, 46(2), 252–287. tools:::Rd_expr_doi("10.1177/0049124115591013") . ## See Also `seqdomassoc` ## Examples ```r data(biofam) ## Building one channel per type of event left, children or married ## Using only first 200 sequences bf <- as.matrix(biofam[1:200, 10:25]) children <- bf==4 | bf==5 | bf==6 married <- bf == 2 | bf== 3 | bf==6 left <- bf==1 | bf==3 | bf==5 | bf==6 ## weights weights <- biofam[1:200,"wp00tbgs"] ## Building sequence objects child.seq <- seqdef(children) marr.seq <- seqdef(married) left.seq <- seqdef(left) ## distances by channel dchild <- seqdist(child.seq, method="OM", sm="INDELSLOG") dmarr <- seqdist(marr.seq, method="OM", sm="INDELSLOG") dleft <- seqdist(left.seq, method="OM", sm="INDELSLOG") dbiofam <- list(dchild,dmarr,dleft) names(dbiofam) <- c("child","marr","left") ## Association between domains and R2 by domain rass <- dissdomassoc(dbiofam, weights=weights) ## Joint distances using additive trick with domain INDELSLOG costs. mcdist <- seqMD(channels=list(child.seq, marr.seq, left.seq), what="diss", method="OM", sm =list("INDELSLOG", "INDELSLOG", "INDELSLOG")) rassj <- dissdomassoc(dbiofam, jointdiss=mcdist, what=c("all"), weights=weights) rassj[["Pearson.Rsquare"]] rass[["Pearson.Rsquare"]] summary(rassj) ``` ## Author(s) Gilbert Ritschard

dissdomassoc function