Generate Artificial data under the GUILDS model

Description

Using this function it is possible to generate a community dataset consisting of two separate abundance vectors for each guild, where the data generated adhere to the Guilds model.

Usage

generate.Guilds(theta, alpha_x, alpha_y, J)

Arguments

Value

Author(s)

Thijs Janzen

Examples

generate.Guilds(theta=200,alpha_x = 0.005, alpha_y = 0.001,J=20000);

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(GUILDS)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/GUILDS/generate.Guilds.Rd_%03d_medium.png", width=480, height=480)
> ### Name: generate.Guilds
> ### Title: Generate Artificial data under the GUILDS model
> ### Aliases: generate.Guilds
> 
> ### ** Examples
> 
> generate.Guilds(theta=200,alpha_x = 0.005, alpha_y = 0.001,J=20000);
$guildX
  [1] 1112 1046  812  661  630  611  603  589  543  541  512  420  416  397  396
 [16]  386  335  324  266  236  229  227  225  218  202  200  142  137  136  122
 [31]  108  107   93   92   91   91   81   80   77   77   77   76   74   73   72
 [46]   72   71   68   64   62   61   56   56   56   54   54   52   52   52   52
 [61]   50   48   42   41   40   39   38   37   33   32   32   30   30   30   28
 [76]   28   28   28   28   28   27   26   26   25   25   24   23   22   21   21
 [91]   20   18   18   17   16   15   14   13   12   11   11   11   10    9    9
[106]    9    8    8    8    7    7    7    7    7    6    6    6    5    5    5
[121]    5    4    4    4    4    4    3    3    3    3    2    2    2    2    2
[136]    2    2    2    2    2    2    2    1    1    1    1    1    1    1    1
[151]    1    1    1    1    1

$guildY
 [1] 578 421 358 284 237 211 198 152 144 137 135 125 109  69  67  63  59  51  47
[20]  46  45  42  30  29  24  22  20  17  16  14  13  10   9   9   8   7   6   5
[39]   5   5   4   4   3   3   3   3   2   2   2   2   2   2   2   1   1   1   1
[58]   1

> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>