function [data, auxData, metaData, txtData, weights] = mydata_Sistrurus_catenatus
%% set metaData
metaData.phylum     = 'Chordata'; 
metaData.class      = 'Reptilia'; 
metaData.order      = 'Squamata'; 
metaData.family     = 'Viperidae';
metaData.species    = 'Sistrurus_catenatus'; 
metaData.species_en = 'Massasauga'; 

metaData.ecoCode.climate = {'BSk','Cfa','Dfa'};
metaData.ecoCode.ecozone = {'THn'};
metaData.ecoCode.habitat = {'0iTg','0iFm'};
metaData.ecoCode.embryo  = {'Tv'};
metaData.ecoCode.migrate = {};
metaData.ecoCode.food    = {'biCv'};
metaData.ecoCode.gender  = {'Dg'};
metaData.ecoCode.reprod  = {'O'};

metaData.T_typical  = C2K(26); % K, body temp
metaData.data_0     = {'ab'; 'ap'; 'am'; 'Lb'; 'Li'; 'Wwi'; 'Ri'}; 
metaData.data_1     = {'t-L'}; 

metaData.COMPLETE = 2.5; % using criteria of LikaKear2011

metaData.author   = {'Bas Kooijman'};    
metaData.date_subm = [2022 21 31];              
metaData.email    = {'bas.kooijman@vu.nl'};            
metaData.address  = {'VU University, Amsterdam'};   

metaData.curator     = {'Starrlight Augustine'};
metaData.email_cur   = {'starrlight@tecnico.ulisboa.pt'}; 
metaData.date_acc    = [2022 12 31];

%% set data
% zero-variate data

data.ab = 107;     units.ab = 'd';    label.ab = 'age at birth';     bibkey.ab = 'AnAge';
  temp.ab = C2K(30); units.temp.ab = 'K'; label.temp.ab = 'temperature';
data.tp = 1100;     units.tp = 'd';    label.tp = 'time since birth at puberty';  bibkey.tp = 'AnAge';
  temp.tp = C2K(30); units.temp.tp = 'K'; label.temp.tp = 'temperature';
data.am = 20*365;  units.am = 'd';    label.am = 'life span';                    bibkey.am = 'AnAge';   
  temp.am = C2K(25); units.temp.am = 'K'; label.temp.am = 'temperature'; 

data.Lb  = 17;    units.Lb  = 'cm';  label.Lb  = 'SVL atbirth';      bibkey.Lb  = 'HileKing2017';
data.Li  = 70;    units.Li  = 'cm';  label.Li  = 'ultimate SVL';  bibkey.Li  = 'HileKing2017';
  comment.Li = 'inferred from tL data';

data.Wwi = 384;   units.Wwi = 'g';  label.Wwi = 'ultimate wet weight';             bibkey.Wwi = 'AnAge';
  
data.Ri  = 8/365;   units.Ri  = '#/d'; label.Ri  = 'maximum reprod rate';          bibkey.Ri  = 'ADW';   
  temp.Ri = C2K(25);   units.temp.Ri = 'K'; label.temp.Ri = 'temperature';
  comment.Ri = 'ADW: 5-20 young per clutch, 1 clutch per yr assumed';

% univariate data

% time-length data
data.tL = [ ... % time since birth (d), SVL (cm)
203.035	16.719
203.035	17.544
203.035	19.430
203.376	16.248
206.791	16.837
207.132	18.723
207.132	19.548
207.132	20.255
207.473	20.845
208.156	18.016
209.181	18.723
209.522	15.540
210.205	19.430
210.888	16.601
211.229	14.361
211.229	17.073
211.912	13.183
211.912	19.312
211.912	19.666
212.253	18.134
213.961	15.187
213.961	16.955
214.302	14.479
214.302	20.845
214.643	19.194
214.985	18.251
215.326	20.609
216.351	18.369
217.033	19.430
217.033	21.788
217.716	16.601
219.082	19.666
222.155	12.829
222.155	17.426
222.155	17.898
222.155	18.487
222.155	18.487
222.496	14.244
222.496	19.666
223.862	19.194
224.203	15.187
224.545	20.373
224.886	15.658
225.228	18.369
225.228	19.430
225.569	20.845
226.593	19.194
226.935	17.898
227.959	14.361
227.959	15.305
227.959	17.191
227.959	21.198
227.959	18.369
228.300	20.138
228.983	16.365
230.008	19.194
230.008	18.487
230.690	12.122
231.032	20.373
231.715	17.544
231.715	19.312
231.715	20.845
232.056	16.130
232.739	17.426
233.422	16.130
233.422	21.552
234.105	18.251
234.788	22.024
234.788	23.202
235.129	18.959
235.470	20.255
235.470	21.316
237.178	20.963
237.519	23.202
238.885	20.609
239.226	22.141
239.226	20.845
239.568	19.548
239.568	20.373
240.250	18.251
240.250	17.544
240.933	22.377
241.275	23.910
241.275	18.487
241.616	23.202
241.616	20.020
241.958	21.316
242.982	21.081
242.982	17.073
243.323	19.194
243.665	18.841
244.006	22.731
244.347	23.202
244.347	19.902
244.689	21.788
245.372	17.426
245.713	20.845
245.713	18.723
246.737	19.430
247.079	20.255
247.420	17.898
249.127	19.784
249.127	21.788
249.810	20.963
250.835	19.902
251.176	20.609
251.176	23.556
251.517	21.906
251.517	22.731
252.883	23.084
253.566	22.259
253.566	21.316
254.932	22.259
255.273	24.263
255.615	20.373
255.615	18.251
256.639	23.792
256.639	21.434
258.346	23.556
258.346	22.613
258.346	20.373
259.370	23.792
259.370	18.841
259.712	19.548
260.395	22.377
260.395	18.959
260.736	21.788
260.736	19.666
261.077	20.491
261.077	16.837
262.443	23.202
263.809	23.674
263.809	22.731
264.492	24.381
265.516	18.959
266.199	20.727
266.882	20.255
267.223	22.731
267.223	21.316
267.564	21.081
267.906	19.077
269.272	21.906
270.637	23.792
272.686	23.320
273.369	23.792
273.710	22.259
274.052	19.548
274.393	20.491
274.734	22.731
275.759	24.028
275.759	24.263
275.759	21.552
276.442	22.849
276.783	25.324
276.783	20.845
277.466	20.020
277.466	19.902
277.807	24.028
278.490	20.845
279.173	23.674
279.173	22.377
279.514	24.028
279.856	20.845
279.856	19.312
280.539	19.077
281.222	21.788
281.222	23.084
281.904	20.373
282.246	22.731
282.587	24.028
283.270	24.853
286.002	24.145
286.002	23.320
287.709	18.487

365+93.096	21.906
365+93.778	23.910
365+103.338	21.552
365+104.363	25.560
365+111.874	24.145
365+115.971	22.495
365+116.654	21.788
365+117.337	23.438
365+117.337	21.552
365+117.678	22.613
365+119.385	24.145
365+119.727	22.495
365+120.410	23.084
365+120.751	21.081
365+121.434	28.271
365+121.434	24.617
365+123.141	25.560
365+123.141	23.202
365+124.165	29.921
365+124.165	22.259
365+125.190	22.024
365+125.531	24.617
365+126.214	25.914
365+126.897	23.556
365+126.897	21.434
365+133.725	26.149
365+133.725	26.974
365+136.457	24.971
365+138.505	23.910
365+139.188	24.499
365+140.554	24.145
365+141.237	25.324
365+141.237	26.385
365+146.017	27.917
365+146.017	28.625
365+146.358	22.259
365+147.041	25.324
365+147.382	23.556
365+148.065	25.796
365+148.065	26.857
365+148.065	27.682
365+148.065	28.153
365+148.407	20.845
365+149.431	24.617
365+151.138	26.739
365+151.480	24.735
365+153.187	23.910
365+153.528	28.389
365+153.528	26.503
365+156.601	25.560
365+160.357	25.442
365+160.357	32.868
365+161.039	28.035
365+161.722	25.088
365+163.429	29.804
365+164.454	28.861
365+165.137	26.739
365+165.137	30.511
365+165.819	29.804
365+166.161	28.861
365+166.502	28.035
365+167.527	26.267
365+168.551	25.324
365+168.551	27.210
365+168.892	29.096
365+169.234	32.986
365+171.624	32.043
365+174.697	24.145
365+175.038	26.149
365+175.038	29.214
365+178.794	25.560
365+180.842	25.560
365+181.525	30.747
365+183.574	31.218
365+185.281	32.161
365+185.964	29.450
365+188.012	26.149
365+188.012	34.754
365+188.354	25.088
365+188.695	29.096
365+189.378	28.035
365+193.134	33.576
365+194.158	32.279
365+194.841	28.625
365+194.841	31.218
365+195.865	33.458
365+197.914	28.389
365+202.352	36.287
365+203.035	36.994
365+206.791	32.868
365+207.815	31.454
365+207.815	30.393
365+209.863	32.397
365+210.888	30.275
365+211.912	39.705
365+212.936	27.800
365+213.619	36.876
365+216.009	33.694
365+216.692	36.523
365+217.033	38.998
365+217.375	31.218
365+218.399	32.279
365+219.765	35.108
365+220.789	36.758
365+222.155	37.819
365+222.496	36.640
365+223.179	30.275
365+224.886	43.124
365+226.593	34.754
365+226.935	38.762
365+232.056	38.527
365+232.739	39.352
365+234.105	29.568
365+234.105	45.835
365+234.788	34.283
365+234.788	34.283
365+234.788	34.990
365+234.788	39.352
365+238.885	30.275
365+238.885	32.279
365+241.275	32.397
365+241.616	30.747
365+242.640	37.112
365+242.982	36.287
365+245.372	41.120
365+247.079	37.701
365+247.420	39.234
365+247.762	40.648
365+248.445	28.035
365+248.445	37.937
365+248.445	35.462
365+248.445	34.872
365+248.445	33.929
365+248.445	33.576
365+248.445	32.279
365+249.810	36.169
365+250.152	39.823
365+251.859	42.063
365+251.859	36.876
365+252.542	33.222
365+252.883	38.762
365+253.225	32.279
365+254.590	40.295
365+256.297	37.701
365+256.639	35.933
365+256.639	36.523
365+256.980	42.652
365+258.005	37.583
365+258.346	39.941
365+258.687	38.409
365+259.370	36.405
365+260.102	39.587
365+264.833	37.348
365+268.247	29.804
365+268.589	34.637
365+269.954	39.470
365+276.442	35.933

730+111.191	42.299
730+114.605	35.462
730+114.947	36.405
730+115.288	42.063
730+115.971	40.884
730+116.654	43.949
730+117.337	38.055
730+120.410	41.709
730+120.751	36.640
730+121.434	34.990
730+121.775	35.344
730+122.117	39.234
730+122.800	31.218
730+122.800	43.124
730+122.800	41.356
730+122.800	38.762
730+123.141	42.534
730+123.483	40.648
730+124.848	43.006
730+124.848	38.291
730+124.848	39.705
730+125.190	36.876
730+125.531	37.466
730+125.873	45.363
730+126.214	34.754
730+126.555	45.835
730+126.897	43.477
730+126.897	40.413
730+127.238	41.827
730+128.263	35.462
730+129.970	41.002
730+132.701	40.413
730+133.725	35.815
730+135.091	37.819
730+136.115	36.994
730+136.457	34.872
730+137.140	40.177
730+138.847	43.949
730+139.188	36.876
730+140.554	40.177
730+141.237	43.949
730+142.944	39.587
730+147.041	43.595
730+147.382	39.823
730+147.382	37.819
730+147.724	36.758
730+148.748	44.774
730+148.748	44.303
730+150.114	39.234
730+150.114	36.758
730+150.455	42.063
730+151.138	44.774
730+151.821	43.006
730+156.601	45.835
730+156.942	37.348
730+158.649	41.120
730+163.771	43.006
730+165.137	41.591
730+166.502	44.303
730+167.868	39.823
730+169.917	43.242
730+171.965	43.006
730+173.331	43.242
730+174.697	47.367
730+178.794	40.648
730+178.794	39.352
730+180.159	42.770
730+181.525	41.945
730+181.525	40.766
730+188.354	43.949
730+206.108	46.778
730+208.498	50.432
730+214.643	48.310
730+214.985	42.417
730+217.033	44.303
730+218.399	46.896
730+219.765	44.420
730+231.032	45.363
730+234.788	43.360
730+242.640	42.652
730+242.982	45.835
730+245.372	45.128
730+249.810	51.493
730+255.956	49.018
730+258.687	46.424
730+258.687	50.668
730+262.443	47.132
730+262.785	50.079
730+263.126	45.363
730+263.467	46.071
730+275.076	52.554
730+282.587	45.010

1095+87.633	46.306
1095+120.068	51.022
1095+120.068	50.432
1095+120.751	48.546
1095+121.434	52.908
1095+121.434	47.367
1095+121.434	46.778
1095+121.434	46.778
1095+121.434	44.656
1095+122.117	47.721
1095+122.800	50.668
1095+122.800	49.607
1095+122.800	45.599
1095+122.800	43.006
1095+123.141	44.185
1095+125.190	44.067
1095+127.580	51.493
1095+131.677	47.014
1095+137.481	45.835
1095+139.188	49.018
1095+139.871	51.493
1095+141.920	52.083
1095+144.651	45.835
1095+144.651	43.006
1095+147.724	56.326
1095+150.114	51.965
1095+150.114	46.306
1095+155.235	47.367
1095+158.649	46.660
1095+158.991	52.672
1095+160.015	54.912
1095+161.039	49.607
1095+161.381	50.550
1095+165.819	50.079
1095+169.234	53.379
1095+169.575	52.554
1095+170.941	51.847
1095+174.355	48.900
1095+177.769	46.306
1095+178.111	48.664
1095+180.159	56.916
1095+181.184	48.782
1095+181.525	49.725
1095+182.208	48.075
1095+183.232	53.733
1095+183.574	52.083
1095+184.256	49.607
1095+184.256	47.250
1095+184.256	51.375
1095+184.598	48.664
1095+185.964	52.790
1095+189.378	46.542
1095+189.378	50.550
1095+200.645	52.790
1095+200.986	48.664
1095+202.693	47.839
1095+204.742	56.208
1095+208.498	45.835
1095+209.522	52.083
1095+212.936	52.436
1095+212.936	51.493
1095+213.278	54.794
1095+229.325	53.379
1095+230.690	50.079
1095+231.715	56.090
1095+235.812	54.322
1095+242.299	47.957
1095+248.103	52.554
1095+256.980	57.269
1095+265.857	53.851
1095+267.564	50.668
1095+271.320	48.193

1460+94.461	53.026
1460+120.068	49.607
1460+120.068	48.664
1460+124.165	54.558
1460+124.165	53.733
1460+126.555	48.900
1460+133.042	49.607
1460+141.920	49.253
1460+144.992	50.079
1460+151.480	49.371
1460+154.211	56.208
1460+158.991	52.200
1460+164.795	50.904
1460+165.137	51.965
1460+167.868	51.611
1460+169.234	52.554
1460+172.648	50.550
1460+187.329	54.322
1460+193.134	52.672
1460+208.839	51.022
1460+214.302	64.224
1460+235.470	50.904
1460+238.543	51.375
1460+250.835	57.269];
data.tL(:,1) = data.tL(:,1) - 203; % set origing at birth
n=size(data.tL,1); for i=2:n; if data.tL(i,1)<=data.tL(i-1,1);data.tL(i,1)=data.tL(i-1,1)+1e-3;end;end
units.tL  = {'d', 'cm'};  label.tL = {'time since birth', 'SVL'};  
temp.tL   = C2K(25);  units.temp.tL = 'K'; label.temp.tL = 'temperature';
bibkey.tL = 'HileKing2017';

%% set weights for all real data
weights = setweights(data, []);

%% set pseudodata and respective weights
[data, units, label, weights] = addpseudodata(data, units, label, weights);

%% pack auxData and txtData for output
auxData.temp = temp;
txtData.units = units;
txtData.label = label;
txtData.bibkey = bibkey;
txtData.comment = comment;

%% Discussion points
D1 = 'Males are assumed not to differ from females';
D2 = 'Temperatures are guessed and for tL assumed to vary as 15+25*sin(2*pi*(t+100)/365)';
metaData.discussion = struct('D1',D1, 'D2',D2);

%% Facts
F1 = 'Sex is determined by incubation temperature';
metaData.bibkey.F1 = 'retileDB'; 
metaData.facts = struct('F1',F1);

%% Links
metaData.links.id_CoL = '4XRRX'; % Cat of Life
metaData.links.id_ITIS = '174304'; % ITIS
metaData.links.id_EoL = '790627'; % Ency of Life
metaData.links.id_Wiki = 'Sistrurus_catenatus'; % Wikipedia
metaData.links.id_ADW = 'Sistrurus_catenatus'; % ADW
metaData.links.id_Taxo = '93130'; % Taxonomicon
metaData.links.id_WoRMS = ''; % WoRMS
metaData.links.id_ReptileDB = 'genus=Sistrurus&species=miliarius'; % ReptileDB
metaData.links.id_AnAge = 'Sistrurus_catenatus'; % AnAge


%% References
bibkey = 'Wiki'; type = 'Misc'; bib = ...
'howpublished = {\url{http://en.wikipedia.org/wiki/Sistrurus_catenatus}}';
metaData.biblist.(bibkey) = ['''@', type, '{', bibkey, ', ' bib, '}'';'];
%
bibkey = 'Kooy2010'; type = 'Book'; bib = [ ...  % used in setting of chemical parameters and pseudodata
'author = {Kooijman, S.A.L.M.}, ' ...
'year = {2010}, ' ...
'title  = {Dynamic Energy Budget theory for metabolic organisation}, ' ...
'publisher = {Cambridge Univ. Press, Cambridge}, ' ...
'pages = {Table 4.2 (page 150), 8.1 (page 300)}, ' ...
'howpublished = {\url{../../../bib/Kooy2010.html}}'];
metaData.biblist.(bibkey) = ['''@', type, '{', bibkey, ', ' bib, '}'';'];
%
bibkey = 'HileKing2017'; type = 'Article'; bib = [ ... 
'doi = {10.1371/journal.pone.0172011}, ' ...
'author = {Hileman, E.T. and King, R.B. and Adamski, J.M. and Anton, T.G. and Bailey, R.L. and Baker S.J., et al.}, ' ... 
'year = {2017}, ' ...
'title = {Climatic and geographic predictors of life history variation in Eastern Massasauga (\emph{Sistrurus catenatus}): {A} range-wide synthesis}, ' ...
'journal = {PLoS ONE}, ' ...
'volume = {12(2)}, ' ...
'pages = {e0172011}'];
metaData.biblist.(bibkey) = ['''@', type, '{', bibkey, ', ' bib, '}'';'];
%
bibkey = 'AnAge'; type = 'Misc'; bib = ...
'howpublished = {\url{http://genomics.senescence.info/species/entry.php?species=Sistrurus_catenatus}}';
metaData.biblist.(bibkey) = ['''@', type, '{', bibkey, ', ' bib, '}'';'];
%
bibkey = 'ADW'; type = 'Misc'; bib = ...
'howpublished = {\url{http://animaldiversity.org/accounts/Sistrurus_catenatus/}}';
metaData.biblist.(bibkey) = ['''@', type, '{', bibkey, ', ' bib, '}'';'];
%
bibkey = 'retileDB'; type = 'Misc'; bib = ...
'howpublished = {\url{http://reptile-database.reptarium.cz/species?genus=Sistrurus&species=catenatus}}';
metaData.biblist.(bibkey) = ['''@', type, '{', bibkey, ', ' bib, '}'';'];

