R/DeFriesFulker.R

In NlsyLinks: Utilities and kinship information for research with the NLSY

#roxygen2 documentation in AceUnivariate.R
DeFriesFulkerMethod1 <- function( dataSet, oName_S1, oName_S2, rName="R" ) { 
  lmDetails <- stats::lm(dataSet[, oName_S1] ~ 1 + dataSet[, oName_S2] + dataSet[, rName] + dataSet[, oName_S2]*dataSet[, rName])
  brief <- summary(lmDetails)
  coeficients <- coef(brief)
  nDouble <- length(brief$residuals) 
  #b0 <- coeficients["(Intercept)", "Estimate"]
  b1 <- coeficients["dataSet[, oName_S2]", "Estimate"]  
  #b2 <- coeficients["R", "Estimate"]
  b3 <- coeficients["dataSet[, oName_S2]:dataSet[, rName]", "Estimate"]
  eSquared <- 1 - (b1+b3)
  
  details <- list(lm=lmDetails)
  aceEstimate <- NlsyLinks::CreateAceEstimate(aSquared=b3, cSquared=b1, eSquared=eSquared, caseCount=nDouble, details=details)
  return( aceEstimate )
}

#roxygen2 documentation in AceUnivariate.R
DeFriesFulkerMethod3 <- function( dataSet, oName_S1, oName_S2, rName="R" ) { 
  dv_S1Centered <- base::scale(dataSet[, oName_S1], center=TRUE, scale=FALSE)
  dv_S2Centered <- base::scale(dataSet[, oName_S2], center=TRUE, scale=FALSE)
  interaction <- dv_S2Centered*dataSet[, rName]
    
  lmDetails <- stats::lm(dv_S1Centered ~ 0 + dv_S2Centered + interaction) #The '0' specifies and intercept-free model.
  brief <- summary(lmDetails)
    
  coeficients <- stats::coef(brief)
  nDouble <- length(brief$residuals) 
  b1 <- coeficients["dv_S2Centered", "Estimate"]  
  b2 <- coeficients["interaction", "Estimate"]
  eSquared <- 1 - (b1+b2)
    
  details <- list(lm=lmDetails)
  aceEstimate <- NlsyLinks::CreateAceEstimate(aSquared=b2, cSquared=b1, eSquared=eSquared, caseCount=nDouble, details=details)
  return( aceEstimate )
}