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)