cchar, v. 1.6: Count characters in sequence data.
cpg, v. 0.7: Compute the CpG content of DNA sequences.
cutSeq, v. 0.11: Cut regions from molecular sequences.
generateQuerySbjct, v. 0.4: Generate pairs of homologous DNA sequences.
gd, v. 0.12: Calculate genetic diversity (pi, S, and Tajima's D) from aligned DNA sequences with or without sliding window.
getSeq, v. 0.4: Get specific sequences from a FASTA file containing multiple entries.
ms2dna, v. 1.16: Generate samples of homologous DNA sequences evolved under defined evolutionary scenarios by converting the output of Richard Hudson's coalescent simulation program ms. As of version 1.11, it can also deal with output generated by Gary Chen's fast coalescent simulator MaCS using the pipeline macs [options] | msformatter | ms2dna -a.
randomizeSeq, v. 0.8: Randomize sequences.
sequencer, v. 1.14: Simulate shotgun sequencing with paired (as of version 1.11) or unpaired reads and a user-defined error rate.
simK, v. 0.4: Simulate pair of sequences with given number of substitutions/site (K).
The rpanel package has two aims. The first is to provide tools which make the construction of gui control panels for R applications as easy as possible. The package is aimed at those who know R but are not familiar with the technicalities of the various gui construction systems which are available. The rpanel package is built on rtcltk and manages the process of communication so that controls can be constructed directly by R simple function calls.
Selection is one of the factors that most influence the shape of genealogical trees. Here we report results of simulations of the infinite-sites version of Moran's model of population genetics aiming at quantifying how the presence of selection affects the branching pattern (topology) of binary genealogical trees. In particular, we consider a scenario of purifying or negative selection in which all mutations are deleterious and each new mutation reduces the fitness of the individual by the same fraction. Analysis of five statistical measures of tree balance or symmetry borrowed from taxonomy indicates that the genealogical trees of samples of populations in which selection is actuating are in the average more asymmetric than neutral trees and that this effect is enhanced by increasing the sample size. However, a quantitative evaluation of the power of these balance measures to detect a tree topology significantly distinct from the neutral one indicates that they are not useful as tests of neutrality of mutations.
M. Bramson, R. Durrett, и G. Swindle. Ann. Probab., 17 (2):
444--481(1989)This paper examines a version of the contact process with a large range. Particles die at rate 1, and a particle is created at an empty site $x$ at rate $łambda$ times the fraction of occupied sites in $y:||x-y||M$. This contact process is dominated by a branching random walk with death rate 1 and birth rate $łambda$, and it is shown that in many ways these two processes are very similar when $M$ is large. In particular, as $M\toınfty$, the critical value for the contact process converges to 1, which is the critical value for branching random walks. The authors obtain precise rates for this convergence, in every dimension, enabling them to describe the ``crossover'' from contact process to branching process behavior in terms of the survival probability of a process started from a single particle. The proofs of the main results use many estimates for branching random walks, further detailing the nature of this crossover behavior..
D. Mollison. Proceedings of the Sixth Berkeley Symposium on Mathematical Statistics and Probability (Univ. California, Berkeley, Calif., 1970/1971), Vol. III: Probability theory, стр. 579--614. Berkeley, Calif., Univ. California Press, (1972)
N. Berestycki. Ensaios Matematicos, 16, cite arxiv:0909.3985
Comment: Lecture notes, to appaear in the collection "Ensaios Matematicos". 17
figures.(2009)
A. Kolmogorov. Izv. Akad. Nauk SSSR, Ser. Math, (1937)Computes density of fairly general Johnson-Mehl crystals and the probability that a point is not in a crystal yet..