... 0 SISUS: Stable Isotope Sourcing using Sampling http://StatAcumen.com/sisus, by Erik Barry Erhardt Starting Wed May 23 15:46:37 2007 ... 0.006 Reading workbook: SISUS_v01-07_Phillips2001_Mink_7-2.xls reading up to 14 sheets, as requested: 1r 2r 3r 4 5 6 7 8 9 10 11 12 13 14 ... 13.656 Process input data and assign variables ### SISUS Model Settings ### analysis.name = Phillips in Oecologia (2001), p.169 Table 2 Mink (n=7, k=2) filename.prefix = P2001Mink_7-2 mixtures.group.name = No Group m = 1 Mixtures n = 7 Sources k = 2 Isotopes Include SIMPLEX (sum of proportions is 1) constraint -- for Mixing Model NOT correcting source isotopic signatures for discrimination NOT using SOURCE CONCENTRATIONS NOT using ASSIMILATION EFFICIENCIES NOT including additional LINEAR CONSTRAINTS NOT using PRIOR, only uniform samples from solution polytope NOT using ISOTOPE STDDEV NOT using CONCENTRATION STDDEV NOT using EFFICIENCY STDDEV Mixtures are not a time series ### SISUS Model Assumptions (under development) ### Mixing Model using Simplex No discrimination or isotope values already corrected. Elemental concentrations are the same for all sources and elements. Assimilation efficiencies are the same for all sources and elements. Solutions within solution polytope are equally likely. No variance in isotope ratios in sources. No variance in elemental concetrations in sources. No variance in assimilation efficiencies. ... 13.742 Set random seed = 76543 ... 13.743 Draw 0 samples from MVN isotope distribution ... 13.838 Plot Convex Hull ==> process Mixture: Mink ... 20.796 Draw uniform samples from convex solution polytope Drawing 1 sample(s) of size 10000 for a total of 10000 samples: 1 (v) Solution Polytope with 5 vertices ... 64.132 1 of 1 = 100 % of mvn samples lead to feasible solutions ... 64.133 Back-solve from biomass to isotopes contributions ... 74.347 Display Mariginal histograms, Scatterplot, Numerical summeries, and Write out samples ... 74.349 Uniform Sample ... 74.351 Biomass Marginal Histogram ... 76.248 Biomass Scatterplot Matrix ... 234.122 Biomass Numerical Summaries ... 234.403 Write out samples ... 234.806 SISUS Complete