DFanalysis: DFanalysis

In Omisc: Univariate Bootstrapping and Other Things

DFanalysis

Description

There are three possible models to be fit. The default is the Rodgers and Kohler formulation of the DF model (Rodgers & Kohler, 2005). The non-default (if RK=F), is to fit the original DeFries-Fulker model. The third option is only used when dominance coefficients are provided, and is based on the formulation by Waller (Waller 1994).

Usage

DFanalysis(
  data = NULL,
  proband,
  sibling,
  Rs,
  Ds = NULL,
  RK = T,
  robust = T,
  DE = T,
  betasonly = F,
  typicalSE = F
)

Arguments

data	A dataframe. This is not necessary as the variables can be passed directly via the other arguments.
proband	Called "proband" for historical reasons this is the variable on the left hand side of the regression.
sibling	The right hand side version of proband. This would be the matched sibling scores.
Rs	This is the vector of relatedness coefficients
Ds	A vector of dominance coefficients. 1 for MZ twins, .25 for DZ twins and full siblings. The default is null, and no value should be provided if using the ACE model. This should only have a non-null value when fitting an ADE model. There is an RK version of this model, however it is not based on published work. The RK version uses double entered (and mean centered) data in order to drop the intercept term and the extraneous regression coefficient (both of which can be constrained to 0 when the phenotypic mean is 0). Initial simulations suggest that this formulation provides accurate parameter estimates, however the original formulation can be used by simply setting RK=F. It is assumed that, if RK=T, that DE=T (i.e., do NOT double enter data prior to analysis if using the ADE model).
RK	Use the Rodgers and Kohler simplified version of the DF model (recommended). Data should not be double entered prior to analysis.
robust	Use the Kohler and Rodgers robust standard errors (recommeneded when using double entered data)
DE	Will the data need to be double entered?
betasonly	If TRUE only the beta weights from the regression analysis will be returned.
typicalSE	Should the typical regression standard errors be used? Default is false.

Value

The results from MyLM

Examples

TwinData<-DFSimulated(2000,2000,.3,.3)
p<-TwinData[,1]
s<-TwinData[,2]
r<-TwinData[,3]
DFanalysis(data=NULL, p,s,r)

Omisc documentation built on Aug. 9, 2022, 5:09 p.m.