stepfun.hist: Create a probability density step function from a histogram...

Description Usage Arguments Value Author(s) See Also Examples

View source: R/connectivity_estimation.distributions.R

Description

This function creates a step function from the bars in a histogram object. By default, the step function will be normalized so that it integrates to 1.

Usage

1
stepfun.hist(h, ..., normalize = TRUE)

Arguments

h

an object of type histogram

...

Additional arguments for the default stepfun function.

normalize

Boolean indicating whether or not to normalize the output stepfun so that it integrates to 1. Defaults to TRUE. If FALSE, then the function will integrate to sum(h$counts)

Value

A function of class stepfun. The height of the steps will be divided by the distance between breaks and possibly the total count.

Author(s)

David M. Kaplan dmkaplan2000@gmail.com

See Also

See also d.rel.conn.dists.func, optim.rel.conn.dists.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
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.