GNU-R >apply(array(NA,dim=c(2,1)), 2 ,tcrossprod) [,1] [1,] NA [2,] NA [3,] NA [4,] NA fastR > apply(array(NA,dim=c(2,1)), 2 ,tcrossprod) Error in tcrossprod(x, y) : non-conformable arguments