Provavelmente isto resolve seu problema

    A <- matrix(1:16, ncol = 4)
    apply(A, 2, function(x) x/mean(x))

A segunda linha pode ser descrita como "Aplicar a função a segunda dimensão(colunas) da matriz"

== Edit ==

Também existe a seguinte opção que retorna o mesmo resultado

    sweep(A, 2, colMeans(A), "/", FALSE)