data( Blanciforti86 )
estResult <- aidsEst( c( "pFood1", "pFood2", "pFood3", "pFood4" ),
c( "wFood1", "wFood2", "wFood3", "wFood4" ), "xFood",
data = Blanciforti86, method = "LA:L" )
wMeans <- colMeans( Blanciforti86[ , c( "wFood1", "wFood2",
"wFood3", "wFood4" ) ] )
aidsEla( estResult$coef, wMeans, formula = "Ch" )
Run the code above in your browser using DataLab