[R] (no subject)

From: unknown <yunusbb_at_inbox.ru>
Date: Fri 24 Jun 2005 - 22:42:28 EST


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

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Fri Jun 24 22:45:56 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:33:02 EST