Usage
bptwin(formula, data, id, zyg, DZ, OS, weight = NULL,
biweight = function(x) 1/min(x), strata = NULL,
messages = 1, control = list(trace = 0), type = "ace",
eqmean = TRUE, pairsonly = FALSE, samecens = TRUE,
allmarg = samecens & !is.null(weight), stderr = TRUE,
robustvar = TRUE, p, indiv = FALSE, constrain,
bound = FALSE, debug = FALSE, ...)Arguments
formula
Formula specifying effects of covariates
on the response.
data
data.frame with one observation pr
row. In addition a column with the zygosity (DZ or MZ
given as a factor) of each individual much be specified
as well as a twin id variable giving a unique pair of
numbers/factors to each twin pair.
id
The name of the column in the dataset
containing the twin-id variable.
zyg
The name of the column in the dataset
containing the zygosity variable.
DZ
Character defining the level in the zyg
variable corresponding to the dyzogitic twins. If this
argument is missing, the reference level (i.e. the first
level) will be interpreted as the dyzogitic twins.
OS
Optional. Character defining the level in the
zyg variable corresponding to the oppposite sex dyzogitic
twins.
weight
Weight matrix if needed by the chosen
estimator. For use with Inverse Probability Weights
biweight
Function defining the bivariate weight in
each cluster
messages
Control amount of messages shown
control
Control argument parsed on to the
optimization routine
type
Character defining the type of analysis to be
performed. Should be a subset of "aced" (additive genetic
factors, common environmental factors, unique
environmental factors, dominant genetic factors).
eqmean
Equal means (with type="u")?
pairsonly
Include complete pairs only?
stderr
Should standard errors be calculated?
robustvar
If TRUE robust (sandwich) variance
estimates of the variance are used
p
Parameter vector p in which to evaluate
log-Likelihood and score function
indiv
If TRUE the score and log-Likelihood
contribution of each twin-pair
constrain
Development argument
allmarg
Should all marginal terms be included
bound
Development argument
debug
Development argument
...
Additional arguments to lower level functions