pollen: Pollen data from Mosimann 1962

Description Usage Format Details References See Also Examples

Description

Data from Mosimann 1962 detailing forest pollen counts

Usage

1

Format

A matrix with four columns and 76 rows.

Details

The rows each sum to 100; the values are counts of four different types of pollen. Each row corresponds to a different level in the core; the levels are in sequence with the first row being most recent and the last row being the oldest.

References

J. E. Mosimann 1962. “On the compound multinomial distribution, the multivariate beta-distribution, and correlations among proportions”. Biometrika, volume 49, numbers 1 and 2, pp65-82.

See Also

serum

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
data(pollen)

func <- function(x,l){
ifelse(any(l<0),Inf,
 lfactorial(sum(x)) -sum(lfactorial(x))
    +lgamma(sum(l)) +sum(lgamma  (x+l))
    -sum(lgamma(l)) -lgamma(sum  (x+l))
)
}

start_vec <- c(51.6, 1, 5.3 , 2)

optim(start_vec , function(l){ -sum(apply(pollen , 1, FUN=func,l=l))})

hyperdirichlet documentation built on May 31, 2017, 5:18 a.m.

Search within the hyperdirichlet package
Search all R packages, documentation and source code