Estoy buscando una manera eficiente de calcular el producto Kronecker de dos matrices grandes. He intentado usar el métodokronecker() como sigue:

 I = diag(700)
 data = replicate(15, rnorm(120))
 test = kronecker(I,data)

Sin embargo, toma mucho tiempo para ejecutarse y luego da el siguiente error:

 Error: cannot allocate vector of size 6.8 Gb