## File Name: mle.reliability.R
## File Version: 0.08
#########################################################
# MLE reliability
mle.reliability <- function(meas, se.meas ){
meas[ abs(meas) > 1E10 ] <- NA
v1 <- stats::var(meas, na.rm=TRUE)
v2 <- mean( se.meas^2, na.rm=TRUE )
rel <- ( v1 - v2 ) / v1
return(rel)
}
#########################################################
#########################################################
# EAP reliability
eap.reliability <- function(meas, se.meas ){
meas[ abs(meas) > 1E10 ] <- NA
v1 <- stats::var(meas, na.rm=TRUE)
v2 <- mean( se.meas^2, na.rm=TRUE )
rel <- v1 / ( v1 + v2 )
return(rel)
}
#########################################################
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.