Atualização do script com bfactor() analysis