Visualize the output of estimate.shiftvec if the data
has only two isotope measurements.
Usage
heatmap.samples(samples)
Arguments
samples
A list of samples of discrimination factors produced by
estimate.shiftvec.
Author(s)
Robert Robere
See Also
estimate.shiftvec
Examples
## get the input from a file, where the predator measurements are simple
## raw isotopic measurements, and prey measurements are in mean/sd form
## Heatmap.samples assumes only two isotopes are measured!
## Not run: input <- build.isotopeframe("predator measurements.dat",
"prey measurements.dat")
## End(Not run)
## now get 2000 samples of discrimination factors
## Not run: out <- estimate.shiftvec(input$pred.frame, input$prey.frame,
num.samples=2000, pred.distinput=FALSE, prey.distinput=TRUE)
## End(Not run)
## now visualize
## Not run: heatmap.samples(out)