data( Blanciforti86 )
# Data on food consumption are available only for the first 32 years
Blanciforti86 <- Blanciforti86[ 1:32, ]
pNames <- c( "pFood1", "pFood2", "pFood3", "pFood4" )
wNames <- c( "wFood1", "wFood2", "wFood3", "wFood4" )
estResult <- aidsEst( pNames, wNames, "xFood",
data = Blanciforti86, method = "IL:L" )
tc <- aidsTestConsist( pNames, wNames, "xFood", Blanciforti86,
coef = estResult$coef )
tc$mPercent # great!
tc$cPercent # Oh, that is bad!
Run the code above in your browser using DataLab