> a = matrix(1:25,5,5)
> a
[,1] [,2] [,3] [,4] [,5]
[1,] 1 6 11 16 21
[2,] 2 7 12 17 22
[3,] 3 8 13 18 23
[4,] 4 9 14 19 24
[5,] 5 10 15 20 25
> b = sweep(a,2,colMeans(a),"/")
> b
[,1] [,2] [,3] [,4] [,5]
[1,] 0.3333333 0.750 0.8461538 0.8888889 0.9130435
[2,] 0.6666667 0.875 0.9230769 0.9444444 0.9565217
[3,] 1.0000000 1.000 1.0000000 1.0000000 1.0000000
[4,] 1.3333333 1.125 1.0769231 1.0555556 1.0434783
[5,] 1.6666667 1.250 1.1538462 1.1111111 1.0869565
其实
t(t(a)/colMeans(a))
也不错。