[R] (no subject)

[unknown]

There is a function 'simMD()' in 'popgen' library which 
simulates a sample of genotype data as follows:
> library(popgen)
> x <- simMD(20, 2, 2, p = NULL, c(0.09, 0.05), ac = 2, beta = 1)
> x
  , , 1
 
      [,1] [,2]
 [1,]    1    1
 [2,]    1    1
      ... 
[37,]    1    2
[38,]    2    2
[39,]    2    2
[40,]    2    2
 
, , 2
 
      [,1] [,2]
 [1,]    2    2
 [2,]    1    2
 [3,]    1    2
      ...
[38,]    2    1
[39,]    1    2
[40,]    1    2
>
How can I repeat this function, for example, 1000 times to generate 1000 samples and assign each output to distinct 'vector' Xi, where i=1,2,...,1000
The goal is to generate a large number of samples using this function and then use them in further analysis.
Any suggestions would be appreciated
Sitdikov Mansor