Dear R experts, Suppose I have an data frame likes this: SampleID grain husk fill 1 KB001 0.076350637 0.28662970 1.45790428 2 KB001 0.118940999 0.22795690 2.94869785 3 KB001 0.150254385 0.56241033 1.09677303 4 KB002 0.188315675 0.21455241 0.96852151 5 KB002 0.180792658 0.23970933 0.97539467 6 KB002 0.180792658 0.23970933 0.97539467 Is there anyway to extract the best linear unbiased predictions (BLUPs) values? Thanks! Regards, KB