Problem replicating gcta function