... 0 SISUS: Stable Isotope Sourcing using Sampling 
      http://StatAcumen.com/sisus, by Erik Barry Erhardt 
      Starting  Wed May 23 00:07:24 2007 
... 0.006 Reading workbook:  SISUS_v01-07_PhillipsGregg2001_Wolf_3-2.xls 
   reading up to  14 sheets, as requested:   1r  2r  3r  4  5  6  7  8  9r  10r  11  12  13  14
... 20.73 Process input data and assign variables 

### SISUS Model Settings ### 
  analysis.name =  Phillips and Gregg in Oecologia (2001) p.175 Table 3 Wolf (n=3, k=2) 
  filename.prefix =  PG2001Wolf_3-2 
  mixtures.group.name =  No group 

  m =  1  Mixtures 
  n =  3  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 
  Drawing MVN samples from Isotope distribution, given ISOTOPE STDDEV, if specified 
  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. 
  ** Experimental ** Variance in isotope ratios, following independent normal distributions with specified center and standard deviation. 
  No variance in elemental concetrations in sources. 
  No variance in assimilation efficiencies. 

... 20.878 Set random seed = 16894 
... 20.88 Draw 0 samples from MVN isotope distribution 
... 20.977 Plot Convex Hull 

==> process Mixture:  Wolf 
... 24.162 Draw uniform samples from convex solution polytope 
           Drawing 1 sample(s) of size 10000 for a total of  10000 samples: 
           1  
           Unique Solution 
... 24.35 1  of  1  =  100 % of mvn samples lead to feasible solutions 
... 24.352 Back-solve from biomass to isotopes contributions 
... 24.354 Display Mariginal histograms, Scatterplot, Numerical summeries, and Write out samples 
... 24.356    Uniform Sample 
... 24.358       Biomass Numerical Summaries 
... 24.426    Write out samples 

... 24.433 SISUS Complete