Changes made in revision 1 of fakedata_mixsir.m

--- /tmp/cvx_h2mH3P	2012-05-19 21:27:29.000000000 -0700
+++ /tmp/cvx_R5kJu2	2012-05-19 21:27:29.000000000 -0700
@@ -1,73 +0,0 @@
-%This program generates fake data for testing the MixSIR program.
-%Written by BX Semmens on 6/30/07
-
-%HOW IT WORKS...
-%set isotopes to 2 (this can be easily changed)
-%randomly draw the number of sources (2 - 10)
-%generate iso sig values for each source (randomize, but constrain to high and low vals for C and N)
-%generate SD in sources (randomize, but constrain to high and low values C and N)
-
-%do a random draw to get % contribs for each source
-
-
-%from that, define the pred iso signature with variance
-
-%DATA:
-%mean and variance of sources
-%mean and variance of pred
-
-%CREDIBLE INTERVALS FOR TWO ISOTOPES
-%CARBON : -35 to -15, SD between 1 and 2
-%NITROGEN: 0 to 15, SD between 1 and 2
-
-
-%CAVIATES
-% 1) these simulations are based on the mixture exclusively derived from
-% this sources... this assumption is always made in mixing models, but rarely true in
-% reality.
-%--------------------------------------------------------------------------
-clear all;
-Cmax=-15;Cmin=-30;CmaxSD=2;CminSD=1;
-Nmax=15;Nmin=0;NmaxSD=2;NminSD=1;
-
-%set max number of sources
-mxsrc = 5;
-%specify the amount of data you want
-data_size=50;
-
-%number of sources
-src_num = round(rand * (mxsrc-2) + 2); %sets min to two sources, max to mxsrc
-
-%now randomly (uniform over interval) draw iso signatures for each source
-for i = 1:src_num
-    src_C(i) = rand * (Cmax-Cmin) + Cmin;
-    src_Csd(i) = rand * (CmaxSD-CminSD) + CminSD;
-    src_N(i) = rand * (Nmax-Nmin) + Nmin;
-    src_Nsd(i) = rand * (NmaxSD-NminSD) + NminSD;
-    
-end
-
-%randomly generate the contribs from each source
-contribs = rand_props(src_num);
-
-%now iterate to generate fake mix data
-mix = zeros(data_size,2); 
-random = randn(data_size,2);
-for i=1:data_size
-    for j=1:src_num
-        mix(i,1)= mix(i,1) + (contribs(j) * (src_C(j)+src_Csd(j).*random(i,1)));
-        mix(i,2)= mix(i,2)+ (contribs(j) * (src_N(j)+src_Nsd(j).*random(i,2)));
-    end
-end
-
-%now write the data files for MixSIR
-frac = zeros(src_num,2); %zeros out fractionation
-dlmwrite('mean_frac.txt',frac,'\t')
-dlmwrite('SD_frac.txt',frac,'\t')
-dlmwrite('mix_data.txt',mix,'\t')
-%dlmwrite('mean_pred.txt',mean(mix),'\t')
-%dlmwrite('SD_pred.txt',std(mix),'\t')
-dlmwrite('mean_source.txt',cat(2,src_C',src_N'),'\t')
-dlmwrite('SD_source.txt',cat(2,src_Csd',src_Nsd'),'\t')
-
-