dual.mark.transmission: Fraction of eggs marked for male and female mark transmission

In dmkaplan2000/ConnMatTools: Tools for Working with Connectivity Data

Description

Estimates the fraction of eggs produced at the source site that are the result of crossing parents, one or both of which have been genotyped. Based on the assumption that probability of breeding between pairs of individuals is completely independent of whether or not one or more of those individuals was genotyped.

Usage

dual.mark.transmission(p.female, p.male = p.female)

Arguments

p.female	Fraction of all adult females genotyped in the source population
p.male	Fraction of all adult males genotyped in the source population. Defaults to be equal to p.female

Value

A list with the following elements:

prob.matrix

2x2 matrix with probabilities for producing offspring with male or female known or unknown parents

p

fraction of all eggs produced at source site that will come from at least one genotyped parent

p.female.known

Fraction of eggs with a single known female parent among all eggs that have one or more known parents

p.male.known

Fraction of eggs with a single known male parent among all eggs that have one or more known parents

p.two.known.parents

Fraction of eggs with two known parents among all eggs that have one or more known parents

Author(s)

David M. Kaplan [e-mail]

References

Kaplan DM, Cuif M, Fauvelot C, Vigliola L, Nguyen-Huu T, Tiavouane J and Lett C (in press) Uncertainty in empirical estimates of marine larval connectivity. ICES Journal of Marine Science. doi:10.1093/icesjms/fsw182.

See Also

Other connectivity estimation: d.rel.conn.beta.prior(), d.rel.conn.dists.func(), d.rel.conn.finite.settlement(), d.rel.conn.multinomial.unnorm(), d.rel.conn.multiple(), d.rel.conn.unif.prior(), optim.rel.conn.dists(), r.marked.egg.fraction()

Examples

library(ConnMatTools)
data(damselfish.lods)

# Histograms of simulated LODs
l <- seq(-1,30,0.5)
h.in <- hist(damselfish.lods$in.group,breaks=l)
h.out <- hist(damselfish.lods$out.group,breaks=l)

# PDFs for marked and unmarked individuals based on simulations
d.marked <- stepfun.hist(h.in)
d.unmarked <- stepfun.hist(h.out)

# Fraction of adults genotyped at source site
p.adults <- 0.25

# prior.shape1=1 # Uniform prior
prior.shape1=0.5 # Jeffreys prior

# Fraction of eggs from one or more genotyped parents
p <- dual.mark.transmission(p.adults)$p

# PDF for relative connectivity
D <- d.rel.conn.dists.func(damselfish.lods$real.children,
                           d.unmarked,d.marked,p,
                           prior.shape1=prior.shape1)

# Estimate most probable value for relative connectivity
phi.mx <- optim.rel.conn.dists(damselfish.lods$real.children,
                                    d.unmarked,d.marked,p)$phi

# Estimate 95% confidence interval for relative connectivity
Q <- q.rel.conn.dists.func(damselfish.lods$real.children,
                           d.unmarked,d.marked,p,
                           prior.shape1=prior.shape1)

# Plot it up
phi <- seq(0,1,0.001)
plot(phi,D(phi),type="l",
     xlim=c(0,0.1),
     main="PDF for relative connectivity",
     xlab=expression(phi),
     ylab="Probability density")

abline(v=phi.mx,col="green",lty="dashed")
abline(v=Q(c(0.025,0.975)),col="red",lty="dashed")

dmkaplan2000/ConnMatTools documentation built on Feb. 8, 2020, 4:41 p.m.