Pode-se fazer o que a pergunta pede com o `xtabs` depois de reformatar os dados de formato largo para formato longo.

    df2 <- reshape2::melt(df, id.vars = 'LOCATION')
    df2 <- xtabs(value ~ variable + LOCATION, df2)
    df2 <- as.data.frame.matrix(df2)
    #         EU28       JPN      USA
    #jan.1999 86.7  99.76086 90.87153
    #fev.1999 85.5  99.76086 91.62469
    #mar.1999 85.9 101.47350 91.56094