@@ -47,22 +47,22 @@ TODO: add colony creation image
4747
4848# Initiating with genomes and parameters
4949
50- We initiate simulation by simulating honyebee genomes. We can quickly generate
50+ We initiate the simulation by simulating honeybee genomes. We can quickly generate
5151genomes using AlphaSimR's ` quickHaplo() ` or with SIMplyBee's ` TODO() ` with
5252appropriate honeybee genome parameters and demography. We will use
53- ` quickHaplo() ` and simulate genomes of two founding individuals. The genomes
54- will be represented by only three chromosomes and 1000 segregating sites per
55- chromosome. Honeybees have 16 chromosomes and way more segregating sites per
56- chromosome, but we want a quick simulation here.
53+ ` quickHaplo() ` and simulate genomes of two founding individuals.
54+ In this example, the genomes will be represented by only three chromosomes
55+ and 1000 segregating sites per chromosome. Honeybees have 16 chromosomes
56+ and far more segregating sites per chromosome, but we want a quick simulation here.
5757
5858``` {r founder genomes}
5959founderGenomes <- quickHaplo(nInd = 2, nChr = 3, segSites = 100)
6060```
6161
6262Now we are ready to setup global simulation parameters using ` SimParamBee ` .
63- ` SimParamBee ` enables tweaking many simulation settings, which you can read
63+ ` SimParamBee ` enables the user to tweak many simulation settings, which you can read
6464about in the help pages (` help(SimParamBee) ` and ` help(SimParam) ` ). Here we will
65- only point out that we recommend saving this object with the name ` SP ` . Namely,
65+ point out that we recommend saving this object with the name ` SP ` . Namely,
6666all SIMplyBee functions will use this object, if you don't directly specify
6767` simParamBee ` argument.
6868
@@ -100,8 +100,8 @@ baseDrones <- createDrones(x = baseQueens[1], nInd = 5)
100100baseDrones
101101```
102102
103- A keen reader will notice that drones' ploidy is two. We will address this in
104- the next sections.
103+ A keen reader will notice that drones' ploidy is two. We will address this in
104+ the next sections.
105105
106106# Creating and building up a colony
107107
@@ -151,11 +151,11 @@ colony
151151
152152Lets explore this colony. In every colony we have different groups of
153153individuals (castes). These include: queen, fathers, workers, drones, and virgin
154- queens. Queen controls the colony, workers do all the hard work, drones
155- disseminate queen's genes, and one of the virgin queens will eventually replace
156- the queen. We also store fathers, which represent drones that the queen mated
157- with. Storing fathers enables us to generate colony members on demand. Let's
158- count how many individuals we have for each caste.
154+ queens. The queen controls the colony, workers do all the hard work, drones
155+ disseminate queen's genes, and one of the virgin queens will eventually replace the
156+ queen. We also store fathers, which represent drones that the queen mated with.
157+ Storing fathers enables us to generate colony members on demand.
158+ Let's count how many individuals we have for each caste.
159159
160160``` {r colony numbers}
161161nQueens(colony)
@@ -183,14 +183,16 @@ view the complete colony data structure with `str()`.
183183str(colony)
184184```
185185
186- To explore the various accessor or modificator functions, study the SIMplyBee's
186+ To explore the various accessor or modification functions, study the SIMplyBee's
187187index of functions.
188188
189189``` {r help, eval = FALSE}
190190help(SIMplyBee)
191191```
192192
193- # Sex and caste
193+ TODO: add colony creation image- not working yet
194+ # knitr::include_graphics("~ /Desktop/GitHub/SIMplyBee/SIMplyBee/vignettes/images/founderpop_diagram2.pdf", auto_pdf = getOption("knitr.graphics.auto_pdf",FALSE))
195+
194196
195197Sex of the caste members is as you would expect.
196198
@@ -240,7 +242,7 @@ These drones are usually discarded by workers. SIMplyBee doesn't store these
240242unviable drones, but it does store their number.
241243
242244We can retrieve information about csd variation with ` getCsdAlleles() ` . For
243- details on where is the csd locus and the number of distinct alleles see
245+ details on where the csd locus is and the number of distinct alleles, see
244246` help(SimParamBee) ` . The output below shows two sequences (in two rows) of 0s
245247and 1s representing respectively ancestral and mutation alleles along the csd
246248locus for the queen. You can see that the two sequences are different, meaning
0 commit comments