dissdomassoc  R Documentation 
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 (Rsquare) of the dissimilarities in one domain that can be reproduced from the dissimilarities in the other domains.
dissdomassoc(domdiss, jointdiss = NULL, what = c("pearson","R2"),
dnames=names(domdiss), weights=NULL, w.rank=FALSE)
domdiss 
List of symmetrical matrices or dist objects: the pairwise dissimilarities per domain. 
jointdiss 

what 
String or vector of strings: requested association measures among 
dnames 
String vector of length equal to number of domains. Names of domains. Default is 
weights 
Vector of nonnegative weights. If 
w.rank 
Logical. If weights are provided and 
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 nonweighted ranks.
When "pearson"
and/or "spearman"
correlations are requested, pvalues (probability to get a stronger correlation under the zero correlation assumption) are automatically computed.
An object of class ddomassoc
, which is a list of tables of the requested association measures and tables of pvalues of the Pearson and Spearman correlations when applicable. The summary
method organizes the correlations and their pvalues in table form.
Gilbert Ritschard
Piccarreta (2017). Joint Sequence Analysis: Association and Clustering. Sociological Methods and Research, 46(2), 252–287. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1177/0049124115591013")}.
seqdomassoc
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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.