Rで2つのベクトルの各要素の組み合わせから行列や表を生成

ちょっとややこしいが、こんな感じだろうか。

nums = c(1,2)

lett = letters[1:3]

matrix(apply(expand.grid(nums, lett), 1, function(istr){ paste(istr[1], istr[2])}), nrow = length(nums), dimnames = list(nums, lett))

もっとうまい方法はあるだろうか。