umx (version 1.4.9)

umx_make_TwinData: umx_make_TwinData: Simulate twin data with control over A, C, E, and moderation

Description

Makes MZ and DZ twin data, optionally moderated. Note, if you want a power calculator, see here: http://www.people.vcu.edu/~bverhulst/power/power.html

Usage

umx_make_TwinData(nMZpairs, nDZpairs = nMZpairs, a = c(avg = 0.5, min = 0, max = 1), c = NULL, e = NULL)

Arguments

nMZpairs

Number of MZ pairs to simulate

nDZpairs

Number of DZ pairs to simulate (if omitted defaults to nMZpairs)

a

value for a, defaults to an example of moderation: c(avg=.5, min=0, max=1)

c

value for c

e

value for e

Value

- list of mzData and dzData data.frames

Details

You supply the number of pairs of each zygosity that wish to simulate (nMZpairs, nDZpairs), along with the values of a, c,and e. a can take a list c(avg = .5, min = 0, max = 1). If specified will act like a moderated heritability, with average value avg, and swinging down to min and up to max across 4-SDs of the moderator.

References

- https://github.com/tbates/umx, https://tbates.github.io

See Also

Other Twin Modeling Functions: plot.MxModel, umxACESexLim, umxACEcov, umxACE, umxCF_SexLim, umxCP, umxGxE_window, umxGxE, umxIP, umxPlotACEcov, umxPlotCP, umxPlotGxE, umxPlotIP, umxSummaryACEcov, umxSummaryACE, umxSummaryCP, umxSummaryGxE, umxSummaryIP, umx

Examples

str(umx_make_TwinData(nMZpairs = 100, nDZpairs = 100, a = .5, c = .3, e = .4))
str(umx_make_TwinData(nMZpairs = 100, nDZpairs = 100, a = .5, c = .3, e = .4))
str(umx_make_TwinData(nMZpairs = 100, a = c(avg = .5, min = 0, max = 1), c = .3, e = .4))