Produce a summary of wide-format twin data, showing the number of individuals, the mean and SD for each trait, and the correlation for each twin-type.
Set MZ and DZ to summarize the two-group case.
umxSummarizeTwinData(
data = NULL,
selVars = NULL,
sep = "_T",
zyg = "zygosity",
MZ = NULL,
DZ = NULL,
MZFF = "MZFF",
DZFF = "DZFF",
MZMM = "MZMM",
DZMM = "DZMM",
DZOS = "DZOS",
digits = 2,
report = c("markdown", "html")
)
The twin data.
Collection of variables to report on, e.g. c("wt", "ht").
The separator string that will turn a variable name into a twin variable name, default= "_T" for wt_T1 and wt_T2.
The zygosity variable in the dataset, default = "zygosity".
Set level in zyg corresponding to MZ for two group case (defaults to using 5-group case).
Set level in zyg corresponding to DZ for two group case (defaults to using 5-group case).
The level of zyg corresponding to MZ FF pairs: default= "MZFF".
The level of zyg corresponding to DZ FF pairs: default= "DZFF".
The level of zyg corresponding to MZ MM pairs: default= "MZMM".
The level of zyg corresponding to DZ MM pairs: default= "DZMM".
The level of zyg corresponding to DZ OS pairs: default= "DZOS".
Rounding precision of the report (default 2).
What to return (default = 'markdown'). Use 'html' to open a web table.
formatted table, e.g. in markdown.
Other Twin Modeling Functions:
power.ACE.test()
,
umxACEcov()
,
umxACEv()
,
umxACE()
,
umxCP()
,
umxDoCp()
,
umxDoC()
,
umxGxE_window()
,
umxGxEbiv()
,
umxGxE()
,
umxIP()
,
umxReduceACE()
,
umxReduceGxE()
,
umxReduce()
,
umxRotate.MxModelCP()
,
umxSexLim()
,
umxSimplex()
,
umxSummaryACEv()
,
umxSummaryACE()
,
umxSummaryDoC()
,
umxSummaryGxEbiv()
,
umxSummarySexLim()
,
umxSummarySimplex()
,
umxTwinMaker()
,
umx
# NOT RUN {
data(twinData)
umxSummarizeTwinData(twinData, sep = "", selVars = c("wt", "ht"))
MZs = c("MZMM", "MZFF"); DZs = c("DZFF","DZMM", "DZOS")
umxSummarizeTwinData(twinData, sep = "", selVars = c("wt", "ht"), MZ = MZs, DZ = DZs)
# }
Run the code above in your browser using DataLab