Dynamic simulation periods are specified in Time's definition. This is usually a list of numbers or labels, typically in some unit of time (days, weeks, months, etc.). Use the ÒDynamic()Ó function in your variables to perform dynamic simulation. 1000 0 0 1 4 13 0 1 2 0 0 0 Risks from farmed and wild salmon v4 28.6.2004 Jouni Tuomisto This is the version 4 of the model calculating risks and benefits of farmed salmon. (c) Copyright KTL (National Public Health Institute, Finland). <ref>[http://ytoswww/yhteiset/Huippuyksikko/Tutkimus/Viljelylohi/Materiaali/Viljelylohi.rmd Reference Manager database</ref> <ref>[http://ytoswww/yhteiset/Huippuyksikko/Tutkimus/Viljelylohi/ Directory for data and models</ref> Jouni Tuomisto 9. tamta 2004 20:14 jtue 8. syyta 2008 17:08 48,24 1,19,28,905,585,17 2,102,90,553,461 Trebuchet MS, 13 0,Model Risks_from_farmed_an,2,2,0,1,E:\Tee\Farmed_salmon-3c.ANA 97,1,1,0,2,1,2794,4312,0 2,25,65,696,600 Pollutant health risk avoided cases/a Pollutant health risk is calculated assuming additivity between the pollutants. However, dioxin risks are not considered because they were not considered in Hites. After unit conversion, numbers are calculated for Western Europe as cases per year. Note that negative numbers mean increased risk unlike in previous versions of the model. var a:= -H1906*H1905/1000*H1910; var b:= (if pollutant1='Dioxin' then 0 else a); var c:= sum(b,pollutant1); c 256,336,1 48,24 1,1,1,1,1,1,0,,1, 2,291,123,476,224 2,357,128,628,450,0,MEAN Graphtool:0 Distresol:10 Diststeps:1 Cdfresol:5 Cdfsteps:1 Symbolsize:6 Baroverlap:0 Linestyle:1 Frame:1 Grid:1 Ticks:1 Mesh:1 Scales:1 Rotation:45 Tilt:0 Depth:70 Frameauto:1 Showkey:1 Xminimum:0 Xmaximum:1 Yminimum:0 Ymaximum:1 Zminimum:0 Zmaximum:1 Xintervals:0 Yintervals:0 Includexzero:0 Includeyzero:0 Includezzero:0 Statsselect:[1,1,1,1,1,0,0,0] Probindex:[0.05,0.25,0.5,0.75,0.95] [H1899,H1898] [1,0,0,0] Health effect of fish avoided cases/a Numbers are calculated for Western Europe as avoided deaths per year. Note that positive numbers mean increased benefit unlike in previous versions of the model. -H1909*min([H1908,benefit_limit])*H1911 384,336,1 48,24 1,1,1,1,1,1,0,,1, 2,80,222,476,224 2,95,7,589,375,0,MEAN [H1898,Salmon1] [1,0,0,0] [Reg_poll,1,Year3,1,Salmon1,1,Recommendation1,1] Net health effect avoided cases/a Net health effect of pollutant cancer and omega-3 cardiac benefit. H1912+H1900 320,408,1 48,24 1,1,1,1,1,1,0,,1, 2,102,90,476,420 2,537,56,726,259,0,MIDM [H1898,Salmon1,1,2] [0,0,0,0] [Salmon,1,Reg_poll,1,Recommendation,1,Sys_localindex('PROBABILITY'),1] Fish advisories ktluser 11. tamta 2004 9:20 48,24 56,208,1 48,24 1,40,0,610,544,17 100,1,1,0,2,9,4744,6798,7 Based on linear cancer risk extrapolation Epa_model 288,40,1 56,32 The model by EPA is well-respected and sound However, defending his study in an interview with IntraFish yesterday afternoon, David Carpenter of the State University of New York at Albany said that the differences between wild and farmed salmon PCBs levels are not insignificant. ÒI donÕt agree with [Gallo] and I donÕt think others would agree with him either.Ó Carpenter said that the EPA risk assessment model that he and his colleagues used to determine that salmon posed a health risk is a well-respected and sound one. ÒIt is not just something that we made up. It is a time-tested measureÉitÕs a yard stick that we have.Ó <ref>http://www.intrafish.com/articlea.php?articleID=41070&s=1</ref> Epa_model 288,128,1 52,36 2,102,90,630,401 Applies only to non-commercial fishing A set of four volumes that provides guidance for assessing health risks associated with the consumption of chemically contaminated non-commercial fish and wildlife. EPA developed the series of documents to help state, local, regional, and tribal environmental health officials who are responsible for developing and managing fish consumption advisories. <ref>http://www.epa.gov/waterscience/fish/guidance.html</ref> Developed_for_spceci 528,192,1 56,28 Point of view is that of a local authority: how to give advice to a fisherman about the consumption of his prey. Developed_for_spceci 360,432,1 68,56 This is not a public health problem but a special case where the authority has a restricted responsibility Applies_only_to_non_+Developed_for_spceci 528,296,1 64,56 Precautionary principle is relevant in this case 200,296,1 48,38 [Constant Developed_for_spceci] The use of the EPA model is problematic for farmed fish. Hites et al, 2004 do not discuss this issue. Point_of_view_is_tha+This_is_not_a_public 528,432,1 64,56 2,102,90,476,373 Based on 1/100000 additional lifetime cancer risk assuming additivity and using linearised multistage model Epa_model 528,80,1 72,64 Developed for spcecial high-exposure subgroups such as tribes and non-commercial fishermen, who eat a lot of fish anyway Epa_model 360,296,1 84,56 [Constant Precautionary_princi] Should we use EPA screening values, FDA action levels or something else? Epa_model+Fda_model 64,288,1 52,52 FDA action level model 5 64,448,1 48,24 2,136,146,416,303,0,MIDM [] EPA fish advisory model <ref>U.S.EPA. Guidance for Assessing Chemical Contaminant Data for Use in Fish Advisory. Volume 2: Risk Assessment and Fish Consumption Limits, 3rd Edition. 2000. Table 3-1. [http://www.epa.gov/waterscience/fish/guidance.html Open access Internet file] [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/EPAFishAdvisory/ Intranet file]</ref> jtue 28. Junta 2004 18:03 48,24 64,80,1 48,29 1,40,0,517,300,17 Advised fish consumption 2^(Floor(logten(Epa_model)/logten(2))) 56,176,1 48,24 2,48,219,743,303,1,MIDM [Location1,Undefined] EPA fish advisory model meals/month CRmm variable in the U.S.EPA advisory model. <ref>U.S.EPA. Guidance for Assessing Chemical Contaminant Data for Use in Fish Advisory. Volume 2: Risk Assessment and Fish Consumption Limits, 3rd Edition. 2000. Table 3-1. [http://www.epa.gov/waterscience/fish/guidance.html Open access Internet file] [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/EPAFishAdvisory/ Intranet file]</ref> index Effect:=['Cancer','Non-cancer']; var MS:= 0.227; var Tap:= 365.25/12; var CSF:=Potency[potency='Ca (CSF)']; var RfD:= Potency[Potency='Non-Ca (RfD)']; var CRlimCa:= ARL*BW/sum(CSF*in1*(Poll_salmon_hites/1000),Pollutant1); var a:= In1*RfD/(Poll_salmon_hites/1000); var b:= if isnan(a) then 1 else a; var c:= if b>0 then b else 1; var CRlimNonCa:= min(c,Pollutant1)*BW; var CRlim:= array(Effect,[CRlimCa,CRlimNonCa]); var CRmm:= CRlim*Tap/MS; CRmm 56,104,1 48,29 2,43,74,632,412 2,120,130,416,303,0,MIDM [] Pollutants in salmon 1 56,32,1 48,24 Poll_salmon_hites Include pollutants Table(Pollutant1)( 1,1,0,1) 216,160,1 48,24 2,469,131,476,224 2,242,231,416,303,0,MIDM Potency 1 216,104,1 48,24 Potency ARL probability Acceptable risk level <ref>U.S.EPA. Guidance for Assessing Chemical Contaminant Data for Use in Fish Advisory. Volume 2: Risk Assessment and Fish Consumption Limits, 3rd Edition. 2000. Table 3-1. [http://www.epa.gov/waterscience/fish/guidance.html Open access Internet file] [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/EPAFishAdvisory/ Intranet file]</ref> 10u 216,32,1 48,24 2,471,122,476,382 Other parts ktluser 11. tamta 2004 9:20 48,24 760,96,1 48,24 1,0,0,1,1,1,0,,0, 1,193,72,636,526,17 Pollutant ['Dieldrin','Toxaphene','Dioxin','PCB'] 504,64,1 48,12 ['Dieldrin','Toxaphene','Dioxin','PCB'] kg 70 504,392,1 48,24 1,1,1,1,1,1,0,0,0,0 Location ['Scotland','Faroe Islands','Frankfurt','Edinburgh','Norway','Paris','London','Oslo','East Canada','Boston','Maine','San Francisco','West Canada','Toronto','Los Angeles','Vancouver','Washington DC','Seattle','Chicago','New York','Washington St','Chile','SE AK Chinook','Denver','BC Chinook','BC Sockeye','Oregon Chinook','SE AK Sockeye','New Orleans','BC Coho','Kodiak AK Sockeye','SE AK Coho','Kodiak AK Coho','BC Pink','Kodiak AK Pink','SE AK Pink','SE AK Chum','BC Chum','Kodiak AK Chum'] 504,32,1 48,12 ['Scotland','Faroe Islands','Frankfurt','Edinburgh','Norway','Paris','London','Oslo','East Canada','Boston','Maine','San Francisco','West Canada','Toronto','Los Angeles','Vancouver','Washington DC','Seattle','Chicago','New York','Washington St','Chile','SE AK Chinook','Denver','BC Chinook','BC Sockeye','Oregon Chinook','SE AK Sockeye','New Orleans','BC Coho','Kodiak AK Sockeye','SE AK Coho','Kodiak AK Coho','BC Pink','Kodiak AK Pink','SE AK Pink','SE AK Chum','BC Chum','Kodiak AK Chum'] ['Farmed salmon','Wild salmon','Market salmon'] 504,96,1 48,13 ['Farmed salmon','Wild salmon','Market salmon'] Loki v 2 20.1.2004 Jouni Tuomisto En ole ehtinyt aiemmin lokia kirjoittaa, joten nyt yleiskuvaus mallista. Viljelylohi on tehty arvioimaan, onko Hites et al (Science 9.1.2004) riskinarviointi hyvin tehty. Ensikommenttina KTL:ssŠ oli se, ettŠ kalan terveyshyšdyt on unohdettu. NiinpŠ rakensimme mallin, joka 1) kŠyttŠŠ samaa EPAn riskimallia saasteiden terveyshaittojen laskemiseen kuin Hites (PCB:n, dieldriinin ja toksafeenin (mutta ei dioksiinien) aiheuttama yhdistetty syšŠriski olettaen additiivisuuden ja linearised multistage-mallin eli suoraan kŠyttŠmŠllŠ EPAn CSF-arvoja) ja 2) laskee myšs omega-3-rasvahappojen tuoman hyšdyn sydŠnkuolemariskiin. Vertailu tehtiin 1) olettamalla lohensyšntiŠ 0.25 - 32 amerikkalaista annosta kuukaudessa ja laskemalla syšpŠriski ja/tai sydŠnhyšty sekŠ 2) olettamalla jokin lohensyšnti (esim. 20 g/d) ja lisŠksi jotain oletuksia muista omega-3-lŠhteistŠ sekŠ niiden muutoksista jos lohensyšnti muuttuisi. Vaikutuksen lisŠksi tehtiin argumenttianalyysi (oma moduli) jossa katsottiin importance analysis eli rank-korrelaatio lŠhtšmuuttujien ja lopputuleman vŠlille. TŠssŠ oli mukana erilaisia pŠŠtšksiŠ, mm. pitŠisikš katsoa saasteiden syšpŠhaittaa vai nettovaikutusta?, PitŠisikš katsoa terveysvastetta lainkaan vai pelkkŠŠ altistusta? ja MillŠ viljelty lohi pitŠisi korvata? PŠŠtškset otettiin mukaan analyysiin siten, ettŠ kullekin pŠŠtšsvaihtoehdolle oletettiin sama todennŠkšisyys, ja ne otettiin mukaan satunnaismuuttujina (ikŠŠn kuin me yrittŠisimme arvioida, mikŠ on ŠŠnestyksen tulos kun tŠstŠ ŠŠnestetŠŠn). TŠhŠn liittyen jŠin pohtimaan sitŠ, pitŠisikš meidŠn olettaa pienempi todennŠkšisyys huonoille vaihtoehdoille (kuten epŠtodennŠkšisille presidenttiehdokkaille annetaan vŠhemmŠn aikaa televisiossa) mutta en pŠŠtynyt tŠssŠ mihinkŠŠn lopputulokseen, ja niin tasajako jŠi malliin. LisŠksi on tehty VOI-analyysi (oma moduli). TŠssŠ yritin rakentaa VOI-funktiota, joka olisi suoraan laskenut mielenkiinnon kohteena olevan tuloksen (helpottaisi mallinrakennusta jatkossa ja tekisi erilaisten VOIn laskemisen kŠtevŠksi), mutta ongelmaksi muodostui se, ettŠ mean-funktio toimi oikein vain, kun se laskettiin variablesta. Jos yritti laskea sen tilapŠisestŠ, solmun sisŠllŠ olevasta muuttujasta, tuloksena oli yleensŠ mid. NiinpŠ tyydyttiin laskemaan homma kŠsipelillŠ kuten Particle VOI -mallissa. Conclusions from Hites 2004 sisŠltŠŠ sitaatteja ja argumentteja keskustelusta, joka on Hitesin myštŠ kŠynnistynyt. What should be the scope of the assessment oli aikeissa olla moduli, josta eri pŠŠtšsvaihtoehdot olisivat sinne kirjatun argumentoinnin seurauksena nousseet, mutta sitŠ ei ollut aikaa tyšstŠŠ kovin pitkŠlle. Confounder analysis -moduli sisŠltŠŠ pohdintaa siitŠ, millaiset tekijŠt voivat vaikuttaa Hitesin lopputulokseen ja kvalitatiivista argumentointia niiden mahdollisista vaikutuksista tulosten tulkintaan. Help on yksinkertaisesti kopio Help v4 mallista. 20.1. alkoi ongelmaksi tulla se, ettŠ malli ršnsysi liikaa, ja oli hankalaa saada indeksit tŠsmŠŠmŠŠn lŠhtšarvojen ja lopputuloksen kesken. NiinpŠ pŠŠtin tehdŠ uuden version 3, josta kaikki ršnsyt on poistettu ja jonka tarkoituksena on toimia mallina Science-artikkelia varten. Kaikki laajemmat tarkastelut siis sŠŠstetŠŠn mallin seuraaviin versioihin. NiinpŠ versio 2:een jŠtetŠŠn kaikki ršnsyt, josta niitŠ sitten voi tarpeen mukaan kopioida takaisin kŠytšssŠ olevaan malliversioon. NŠin ehkŠ pysyy selvŠnŠ se, mitŠ Science-vastineessa on ja mitŠ ei ole. Fishing and farming, Arguments on fish pollutants, Total pollutant exposure, What should be the scope of risk assessment?, Conclusions from Hites 2004, ja Confounder analysis ovat semmoiset modulit jotka nyt poistetaan versiosta 3. Samoin poistetaan pŠŠmodulista solmut Risk or net health effect?, Acceptable risk ja Health effects or exposures? sekŠ nŠiden input nodet. 0 504,128,1 48,12 2,463,67,476,399 65535,54067,19661 Loki v3 20.1.2003 Jouni Tuomisto Versiosta 3 siis on tehty riisuttu versio Science-juttua varten. Lue tarkemmin Loki v 2:sta. Nyt versiosta 3 poistetaan aiemmin kuvatun lisŠksi Other parts -modulista indeksejŠ, joita ei kŠytetŠ missŠŠn. NŠmŠ ovat Viljelyalue, Kalastusalue, Ostokaupunki, Lohilaji ja Saaste. Argument analysis -modulista poistetaan solmut We should not consider concentrations..., Acceptable exposure, Va2, Health or exposure?, What is salmon replacement?, Va5, Va3, Va3 inputs, Va3 importance eli kaikki solmut. TŠrkeyssolmu luodaan uudelleen, mutta nyt se voidaan tehdŠ suoraan Outcome-solmulle ilman indeksimuunnoksia. NiinpŠ koko Argument analysis -moduli poistetaan ja asia siirretŠŠn VOIs-moduliin, joka nimetŠŠn uudelleen VOI and importance analysis. Food intake -modulista poistetaan Salmon consumption, Salmon replacement, Food change, Salmon amount, Oil increase, Food intake, Wild salmon compensation, Va1, Food_rec, Source1. Eli kaikki solmut keskittyvŠt nyt vain loheen, eikŠ muita omega-3-lŠhteitŠ huomioida. Ne tulevat mukaan malliin raja-arvossa, joka kuvaa hyšdyllisen lisŠsaannin rajaa ja siis sisŠltŠŠ absoluuttisen fysiologisen rajan, josta on vŠhennetty muusta ravinnosta tuleva mŠŠrŠ. TŠmŠ solmu tehdŠŠn Annosvastemoduliin. VOIs-modulista poistetaan Va16, Va12, VOI, VOI1 ja VOI-laskenta tehdŠŠn suoraan Outcome-solmusta. 21.1.2004 Jouni Tuomisto Malli muuttui eilen siten, ettŠ nyt lasketaan VOI kahdelle eri kysymykselle: pitŠisikš suositella viljellyn lohen enimmŠissaanniksi 1 annos/kk ja pitŠisikš rajoittaa enemmŠn kalanrehun saastepitoisuuksia. NŠmŠ kaksi nostetaan esiin, koska edellinen on suora vastine Hitesin ym. argumenttiin, ja jŠlkimmŠinen on korostamassa sitŠ, ettŠ asetettu kysymys mŠŠrŠŠ sen, mikŠ tieto on tŠrkeŠŠ ja mikŠ ei. Other parts -modulista poistetaan solmut Acceptable exposure increase ja Amount or replacement, ja ARL siirretŠŠn Fish advirories -moduliin sekŠ Potency Exposure-response function for pollutant risk -moduliin. NŠistŠ moduleista poistetaan vastaavat aliakset. Unit- ja Description-kentŠt pŠivitetŠŠn koko mallissa, ja viitteitŠ listŠtŠŠn sikŠli kuin ne ovat helposti kŠsillŠ. Kuitenkin viitteet on vielŠ pistettŠvŠ kuntoon, nyt muotoilut eivŠt ole kunnossa. 0 504,152,1 48,12 2,212,144,476,344 65535,54067,19661 Compensating fish amount g/d H1900/H1909/H1910 504,224,1 48,24 [H1898,Salmon1] Probability of decision var a:= sum(H1900,Salmon1); Probability(a[H1898='BAU']+0.0001>a[H1898='Restrict farmed salmon use']) 504,336,1 48,24 2,115,372,476,291 Probability of decision var a:= sum(H1901,Salmon1); Probability(a[H1898='BAU']>a[H1898='Change farmed to wild salmon']) 504,280,1 48,24 Outcomes index a:= ['Net effect of salmon recommendation','Net effect of feed regulation','Cancer effect of recommendation']; var b:= array(a,[H1901,0,H1900]); var c:= b[H1898='Restrict farmed salmon use']-b[H1898='BAU']; var d:= (if regulate_pollutants_=1 then c[H1899='More actions'] else c[H1899='BAU']); var e:= H1901[H1899='More actions'] - H1901[H1899='BAU']; var f:= (if recommend= 1 then e[H1898='Restrict farmed salmon use'] else e[H1898='BAU']); var g:= (if a= 'Net effect of feed regulation' then f else d); var h:= sum(g,Salmon1); h 232,264,1 48,24 2,452,264,476,469 2,767,202,367,474,0,MEAN Pollutant or net health effect? probability A chance node that collapses the decision about whether the proper endpoint metric is pollutant risk or net health effect. Bernoulli( .5 ) 336,56,1 48,29 1,1,1,1,1,1,0,0,0,0 2,585,196,416,303,0,MIDM 2,168,178,416,303,0,SAMP Mortality by recommendation cases/a Net health effect indexed by only consumption recommendation. if Regulate_pollutants_=1 then outcome3[H1899='More actions'] else outcome3[H1899='BAU'] 168,152,1 48,32 2,102,90,476,293 2,499,269,416,303,0,MIDM 2,415,126,518,378,0,MEAN [Salmon1,H1898] Regulate pollutants? probability A chance node that collapses the decision about regulating fish feed. bernoulli(0.5) 56,152,1 48,24 Mortality by feed regulation cases/a Net health effect indexed by only fish feed regulation. if recommend = 1 then outcome3[H1898='Restrict farmed salmon use'] else outcome3[H1898='BAU'] 280,152,1 48,32 2,408,141,646,274 2,499,269,416,303,0,MIDM 2,77,202,518,378,0,MEAN [Salmon1,H1898] Recommend? probability A chance node that collapses the decision about consumption recommendations for farmed salmon. bernoulli(0.5) 384,152,1 48,24 1,1,1,1,1,1,0,0,0,0 2,298,233,476,224 2,80,145,416,303,0,SAMP Lifetime cancer+CHD mortality prevented by salmon cases/a Net health effect indexed by the two decisions of concern: a) whether to recommend salmon consumption restrictions and b) whether to apply stricter regulations for fish feed. The definition contained also this row: if a>ARL then 1 else 0 But it was removed when we cut the acceptable concentration concept out of the model. var a:= (if Poll_or_net = 1 then H1900 else H1901); sum(a,Salmon1) 224,56,1 48,46 2,102,90,476,277 2,499,269,416,303,0,MIDM 2,723,592,518,176,0,MEAN [H1899,H1898] [1,0,0,0] Log v4 28.6.2004 Jouni Tuomisto This version is an update of the model that was used for calculating the results for Tuomisto paper submitted to Science on January 28, 2004. Only argumentation and comments have been clarified and added. No substantive changes have been made to definitions. The descriptions of the variables described in the Table 1 in 'Supportive online material' have not been changed. Therefore, this version produces the results that were presented in the manuscript. The main conclusions of this study were added as arguments on the top level of the model. The argumentation about the pollutant model selection was clarified (see module Fish advisories). 11.1.2005 Jouni Tuomisto The model was given a unified resource name (URN) and metadata. The only addition since 28.6.2004 is the node Urn and this note. 0 504,184,1 48,12 2,306,93,476,540 [Alias Log_v4] 65535,54067,19661 Benefit-risk diagram for farmed salmon index benefit_risk: ['Benefits','Risks']; var a:= array(benefit_risk,[H1912,H1900]); a[Salmon1='Farmed salmon'H1899='BAU'] 72,264,1 48,42 2,20,7,551,791,1,MEAN [Sys_localindex('BENEFIT_RISK'),H1898,Undefined,Undefined,1] [0,0,0,0] [Recommendation,1,Sys_localindex('BR'),1,Sys_localindex('STEP'),1] VOI analysis for farmed salmon 72,352,1 48,32 Rows for sql 1 var a:= slice(Result_reporting,Result_reporting.decision,4); '"'&3000+run&'";"'&1&'";"'&a&'";"'&run&'"' 400,288,1 48,24 2,102,90,490,397 2,568,87,416,303,0,MIDM 2,I,4,2,0,0,4,0,$,0,"ABBREV",0 1..4000 56,32,1 48,24 [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108,109,110,111,112,113,114,115,116,117,118,119,120,121,122,123,124,125,126,127,128,129,130,131,132,133,134,135,136,137,138,139,140,141,142,143,144,145,146,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,201,202,203,204,205,206,207,208,209,210,211,212,213,214,215,216,217,218,219,220,221,222,223,224,225,226,227,228,229,230,231,232,233,234,235,236,237,238,239,240,241,242,243,244,245,246,247,248,249,250,251,252,253,254,255,256,257,258,259,260,261,262,263,264,265,266,267,268,269,270,271,272,273,274,275,276,277,278,279,280,281,282,283,284,285,286,287,288,289,290,291,292,293,294,295,296,297,298,299,300,301,302,303,304,305,306,307,308,309,310,311,312,313,314,315,316,317,318,319,320,321,322,323,324,325,326,327,328,329,330,331,332,333,334,335,336,337,338,339,340,341,342,343,344,345,346,347,348,349,350,351,352,353,354,355,356,357,358,359,360,361,362,363,364,365,366,367,368,369,370,371,372,373,374,375,376,377,378,379,380,381,382,383,384,385,386,387,388,389,390,391,392,393,394,395,396,397,398,399,400,401,402,403,404,405,406,407,408,409,410,411,412,413,414,415,416,417,418,419,420,421,422,423,424,425,426,427,428,429,430,431,432,433,434,435,436,437,438,439,440,441,442,443,444,445,446,447,448,449,450,451,452,453,454,455,456,457,458,459,460,461,462,463,464,465,466,467,468,469,470,471,472,473,474,475,476,477,478,479,480,481,482,483,484,485,486,487,488,489,490,491,492,493,494,495,496,497,498,499,500,501,502,503,504,505,506,507,508,509,510,511,512,513,514,515,516,517,518,519,520,521,522,523,524,525,526,527,528,529,530,531,532,533,534,535,536,537,538,539,540,541,542,543,544,545,546,547,548,549,550,551,552,553,554,555,556,557,558,559,560,561,562,563,564,565,566,567,568,569,570,571,572,573,574,575,576,577,578,579,580,581,582,583,584,585,586,587,588,589,590,591,592,593,594,595,596,597,598,599,600,601,602,603,604,605,606,607,608,609,610,611,612,613,614,615,616,617,618,619,620,621,622,623,624,625,626,627,628,629,630,631,632,633,634,635,636,637,638,639,640,641,642,643,644,645,646,647,648,649,650,651,652,653,654,655,656,657,658,659,660,661,662,663,664,665,666,667,668,669,670,671,672,673,674,675,676,677,678,679,680,681,682,683,684,685,686,687,688,689,690,691,692,693,694,695,696,697,698,699,700,701,702,703,704,705,706,707,708,709,710,711,712,713,714,715,716,717,718,719,720,721,722,723,724,725,726,727,728,729,730,731,732,733,734,735,736,737,738,739,740,741,742,743,744,745,746,747,748,749,750,751,752,753,754,755,756,757,758,759,760,761,762,763,764,765,766,767,768,769,770,771,772,773,774,775,776,777,778,779,780,781,782,783,784,785,786,787,788,789,790,791,792,793,794,795,796,797,798,799,800,801,802,803,804,805,806,807,808,809,810,811,812,813,814,815,816,817,818,819,820,821,822,823,824,825,826,827,828,829,830,831,832,833,834,835,836,837,838,839,840,841,842,843,844,845,846,847,848,849,850,851,852,853,854,855,856,857,858,859,860,861,862,863,864,865,866,867,868,869,870,871,872,873,874,875,876,877,878,879,880,881,882,883,884,885,886,887,888,889,890,891,892,893,894,895,896,897,898,899,900,901,902,903,904,905,906,907,908,909,910,911,912,913,914,915,916,917,918,919,920,921,922,923,924,925,926,927,928,929,930,931,932,933,934,935,936,937,938,939,940,941,942,943,944,945,946,947,948,949,950,951,952,953,954,955,956,957,958,959,960,961,962,963,964,965,966,967,968,969,970,971,972,973,974,975,976,977,978,979,980,981,982,983,984,985,986,987,988,989,990,991,992,993,994,995,996,997,998,999,1000,1001,1002,1003,1004,1005,1006,1007,1008,1009,1010,1011,1012,1013,1014,1015,1016,1017,1018,1019,1020,1021,1022,1023,1024,1025,1026,1027,1028,1029,1030,1031,1032,1033,1034,1035,1036,1037,1038,1039,1040,1041,1042,1043,1044,1045,1046,1047,1048,1049,1050,1051,1052,1053,1054,1055,1056,1057,1058,1059,1060,1061,1062,1063,1064,1065,1066,1067,1068,1069,1070,1071,1072,1073,1074,1075,1076,1077,1078,1079,1080,1081,1082,1083,1084,1085,1086,1087,1088,1089,1090,1091,1092,1093,1094,1095,1096,1097,1098,1099,1100,1101,1102,1103,1104,1105,1106,1107,1108,1109,1110,1111,1112,1113,1114,1115,1116,1117,1118,1119,1120,1121,1122,1123,1124,1125,1126,1127,1128,1129,1130,1131,1132,1133,1134,1135,1136,1137,1138,1139,1140,1141,1142,1143,1144,1145,1146,1147,1148,1149,1150,1151,1152,1153,1154,1155,1156,1157,1158,1159,1160,1161,1162,1163,1164,1165,1166,1167,1168,1169,1170,1171,1172,1173,1174,1175,1176,1177,1178,1179,1180,1181,1182,1183,1184,1185,1186,1187,1188,1189,1190,1191,1192,1193,1194,1195,1196,1197,1198,1199,1200,1201,1202,1203,1204,1205,1206,1207,1208,1209,1210,1211,1212,1213,1214,1215,1216,1217,1218,1219,1220,1221,1222,1223,1224,1225,1226,1227,1228,1229,1230,1231,1232,1233,1234,1235,1236,1237,1238,1239,1240,1241,1242,1243,1244,1245,1246,1247,1248,1249,1250,1251,1252,1253,1254,1255,1256,1257,1258,1259,1260,1261,1262,1263,1264,1265,1266,1267,1268,1269,1270,1271,1272,1273,1274,1275,1276,1277,1278,1279,1280,1281,1282,1283,1284,1285,1286,1287,1288,1289,1290,1291,1292,1293,1294,1295,1296,1297,1298,1299,1300,1301,1302,1303,1304,1305,1306,1307,1308,1309,1310,1311,1312,1313,1314,1315,1316,1317,1318,1319,1320,1321,1322,1323,1324,1325,1326,1327,1328,1329,1330,1331,1332,1333,1334,1335,1336,1337,1338,1339,1340,1341,1342,1343,1344,1345,1346,1347,1348,1349,1350,1351,1352,1353,1354,1355,1356,1357,1358,1359,1360,1361,1362,1363,1364,1365,1366,1367,1368,1369,1370,1371,1372,1373,1374,1375,1376,1377,1378,1379,1380,1381,1382,1383,1384,1385,1386,1387,1388,1389,1390,1391,1392,1393,1394,1395,1396,1397,1398,1399,1400,1401,1402,1403,1404,1405,1406,1407,1408,1409,1410,1411,1412,1413,1414,1415,1416,1417,1418,1419,1420,1421,1422,1423,1424,1425,1426,1427,1428,1429,1430,1431,1432,1433,1434,1435,1436,1437,1438,1439,1440,1441,1442,1443,1444,1445,1446,1447,1448,1449,1450,1451,1452,1453,1454,1455,1456,1457,1458,1459,1460,1461,1462,1463,1464,1465,1466,1467,1468,1469,1470,1471,1472,1473,1474,1475,1476,1477,1478,1479,1480,1481,1482,1483,1484,1485,1486,1487,1488,1489,1490,1491,1492,1493,1494,1495,1496,1497,1498,1499,1500,1501,1502,1503,1504,1505,1506,1507,1508,1509,1510,1511,1512,1513,1514,1515,1516,1517,1518,1519,1520,1521,1522,1523,1524,1525,1526,1527,1528,1529,1530,1531,1532,1533,1534,1535,1536,1537,1538,1539,1540,1541,1542,1543,1544,1545,1546,1547,1548,1549,1550,1551,1552,1553,1554,1555,1556,1557,1558,1559,1560,1561,1562,1563,1564,1565,1566,1567,1568,1569,1570,1571,1572,1573,1574,1575,1576,1577,1578,1579,1580,1581,1582,1583,1584,1585,1586,1587,1588,1589,1590,1591,1592,1593,1594,1595,1596,1597,1598,1599,1600,1601,1602,1603,1604,1605,1606,1607,1608,1609,1610,1611,1612,1613,1614,1615,1616,1617,1618,1619,1620,1621,1622,1623,1624,1625,1626,1627,1628,1629,1630,1631,1632,1633,1634,1635,1636,1637,1638,1639,1640,1641,1642,1643,1644,1645,1646,1647,1648,1649,1650,1651,1652,1653,1654,1655,1656,1657,1658,1659,1660,1661,1662,1663,1664,1665,1666,1667,1668,1669,1670,1671,1672,1673,1674,1675,1676,1677,1678,1679,1680,1681,1682,1683,1684,1685,1686,1687,1688,1689,1690,1691,1692,1693,1694,1695,1696,1697,1698,1699,1700,1701,1702,1703,1704,1705,1706,1707,1708,1709,1710,1711,1712,1713,1714,1715,1716,1717,1718,1719,1720,1721,1722,1723,1724,1725,1726,1727,1728,1729,1730,1731,1732,1733,1734,1735,1736,1737,1738,1739,1740,1741,1742,1743,1744,1745,1746,1747,1748,1749,1750,1751,1752,1753,1754,1755,1756,1757,1758,1759,1760,1761,1762,1763,1764,1765,1766,1767,1768,1769,1770,1771,1772,1773,1774,1775,1776,1777,1778,1779,1780,1781,1782,1783,1784,1785,1786,1787,1788,1789,1790,1791,1792,1793,1794,1795,1796,1797,1798,1799,1800,1801,1802,1803,1804,1805,1806,1807,1808,1809,1810,1811,1812,1813,1814,1815,1816,1817,1818,1819,1820,1821,1822,1823,1824,1825,1826,1827,1828,1829,1830,1831,1832,1833,1834,1835,1836,1837,1838,1839,1840,1841,1842,1843,1844,1845,1846,1847,1848,1849,1850,1851,1852,1853,1854,1855,1856,1857,1858,1859,1860,1861,1862,1863,1864,1865,1866,1867,1868,1869,1870,1871,1872,1873,1874,1875,1876,1877,1878,1879,1880,1881,1882,1883,1884,1885,1886,1887,1888,1889,1890,1891,1892,1893,1894,1895,1896,1897,1898,1899,1900,1901,1902,1903,1904,1905,1906,1907,1908,1909,1910,1911,1912,1913,1914,1915,1916,1917,1918,1919,1920,1921,1922,1923,1924,1925,1926,1927,1928,1929,1930,1931,1932,1933,1934,1935,1936,1937,1938,1939,1940,1941,1942,1943,1944,1945,1946,1947,1948,1949,1950,1951,1952,1953,1954,1955,1956,1957,1958,1959,1960,1961,1962,1963,1964,1965,1966,1967,1968,1969,1970,1971,1972,1973,1974,1975,1976,1977,1978,1979,1980,1981,1982,1983,1984,1985,1986,1987,1988,1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001,2002,2003,2004,2005,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016,2017,2018,2019,2020,2021,2022,2023,2024,2025,2026,2027,2028,2029,2030,2031,2032,2033,2034,2035,2036,2037,2038,2039,2040,2041,2042,2043,2044,2045,2046,2047,2048,2049,2050,2051,2052,2053,2054,2055,2056,2057,2058,2059,2060,2061,2062,2063,2064,2065,2066,2067,2068,2069,2070,2071,2072,2073,2074,2075,2076,2077,2078,2079,2080,2081,2082,2083,2084,2085,2086,2087,2088,2089,2090,2091,2092,2093,2094,2095,2096,2097,2098,2099,2100,2101,2102,2103,2104,2105,2106,2107,2108,2109,2110,2111,2112,2113,2114,2115,2116,2117,2118,2119,2120,2121,2122,2123,2124,2125,2126,2127,2128,2129,2130,2131,2132,2133,2134,2135,2136,2137,2138,2139,2140,2141,2142,2143,2144,2145,2146,2147,2148,2149,2150,2151,2152,2153,2154,2155,2156,2157,2158,2159,2160,2161,2162,2163,2164,2165,2166,2167,2168,2169,2170,2171,2172,2173,2174,2175,2176,2177,2178,2179,2180,2181,2182,2183,2184,2185,2186,2187,2188,2189,2190,2191,2192,2193,2194,2195,2196,2197,2198,2199,2200,2201,2202,2203,2204,2205,2206,2207,2208,2209,2210,2211,2212,2213,2214,2215,2216,2217,2218,2219,2220,2221,2222,2223,2224,2225,2226,2227,2228,2229,2230,2231,2232,2233,2234,2235,2236,2237,2238,2239,2240,2241,2242,2243,2244,2245,2246,2247,2248,2249,2250,2251,2252,2253,2254,2255,2256,2257,2258,2259,2260,2261,2262,2263,2264,2265,2266,2267,2268,2269,2270,2271,2272,2273,2274,2275,2276,2277,2278,2279,2280,2281,2282,2283,2284,2285,2286,2287,2288,2289,2290,2291,2292,2293,2294,2295,2296,2297,2298,2299,2300,2301,2302,2303,2304,2305,2306,2307,2308,2309,2310,2311,2312,2313,2314,2315,2316,2317,2318,2319,2320,2321,2322,2323,2324,2325,2326,2327,2328,2329,2330,2331,2332,2333,2334,2335,2336,2337,2338,2339,2340,2341,2342,2343,2344,2345,2346,2347,2348,2349,2350,2351,2352,2353,2354,2355,2356,2357,2358,2359,2360,2361,2362,2363,2364,2365,2366,2367,2368,2369,2370,2371,2372,2373,2374,2375,2376,2377,2378,2379,2380,2381,2382,2383,2384,2385,2386,2387,2388,2389,2390,2391,2392,2393,2394,2395,2396,2397,2398,2399,2400,2401,2402,2403,2404,2405,2406,2407,2408,2409,2410,2411,2412,2413,2414,2415,2416,2417,2418,2419,2420,2421,2422,2423,2424,2425,2426,2427,2428,2429,2430,2431,2432,2433,2434,2435,2436,2437,2438,2439,2440,2441,2442,2443,2444,2445,2446,2447,2448,2449,2450,2451,2452,2453,2454,2455,2456,2457,2458,2459,2460,2461,2462,2463,2464,2465,2466,2467,2468,2469,2470,2471,2472,2473,2474,2475,2476,2477,2478,2479,2480,2481,2482,2483,2484,2485,2486,2487,2488,2489,2490,2491,2492,2493,2494,2495,2496,2497,2498,2499,2500,2501,2502,2503,2504,2505,2506,2507,2508,2509,2510,2511,2512,2513,2514,2515,2516,2517,2518,2519,2520,2521,2522,2523,2524,2525,2526,2527,2528,2529,2530,2531,2532,2533,2534,2535,2536,2537,2538,2539,2540,2541,2542,2543,2544,2545,2546,2547,2548,2549,2550,2551,2552,2553,2554,2555,2556,2557,2558,2559,2560,2561,2562,2563,2564,2565,2566,2567,2568,2569,2570,2571,2572,2573,2574,2575,2576,2577,2578,2579,2580,2581,2582,2583,2584,2585,2586,2587,2588,2589,2590,2591,2592,2593,2594,2595,2596,2597,2598,2599,2600,2601,2602,2603,2604,2605,2606,2607,2608,2609,2610,2611,2612,2613,2614,2615,2616,2617,2618,2619,2620,2621,2622,2623,2624,2625,2626,2627,2628,2629,2630,2631,2632,2633,2634,2635,2636,2637,2638,2639,2640,2641,2642,2643,2644,2645,2646,2647,2648,2649,2650,2651,2652,2653,2654,2655,2656,2657,2658,2659,2660,2661,2662,2663,2664,2665,2666,2667,2668,2669,2670,2671,2672,2673,2674,2675,2676,2677,2678,2679,2680,2681,2682,2683,2684,2685,2686,2687,2688,2689,2690,2691,2692,2693,2694,2695,2696,2697,2698,2699,2700,2701,2702,2703,2704,2705,2706,2707,2708,2709,2710,2711,2712,2713,2714,2715,2716,2717,2718,2719,2720,2721,2722,2723,2724,2725,2726,2727,2728,2729,2730,2731,2732,2733,2734,2735,2736,2737,2738,2739,2740,2741,2742,2743,2744,2745,2746,2747,2748,2749,2750,2751,2752,2753,2754,2755,2756,2757,2758,2759,2760,2761,2762,2763,2764,2765,2766,2767,2768,2769,2770,2771,2772,2773,2774,2775,2776,2777,2778,2779,2780,2781,2782,2783,2784,2785,2786,2787,2788,2789,2790,2791,2792,2793,2794,2795,2796,2797,2798,2799,2800,2801,2802,2803,2804,2805,2806,2807,2808,2809,2810,2811,2812,2813,2814,2815,2816,2817,2818,2819,2820,2821,2822,2823,2824,2825,2826,2827,2828,2829,2830,2831,2832,2833,2834,2835,2836,2837,2838,2839,2840,2841,2842,2843,2844,2845,2846,2847,2848,2849,2850,2851,2852,2853,2854,2855,2856,2857,2858,2859,2860,2861,2862,2863,2864,2865,2866,2867,2868,2869,2870,2871,2872,2873,2874,2875,2876,2877,2878,2879,2880,2881,2882,2883,2884,2885,2886,2887,2888,2889,2890,2891,2892,2893,2894,2895,2896,2897,2898,2899,2900,2901,2902,2903,2904,2905,2906,2907,2908,2909,2910,2911,2912,2913,2914,2915,2916,2917,2918,2919,2920,2921,2922,2923,2924,2925,2926,2927,2928,2929,2930,2931,2932,2933,2934,2935,2936,2937,2938,2939,2940,2941,2942,2943,2944,2945,2946,2947,2948,2949,2950,2951,2952,2953,2954,2955,2956,2957,2958,2959,2960,2961,2962,2963,2964,2965,2966,2967,2968,2969,2970,2971,2972,2973,2974,2975,2976,2977,2978,2979,2980,2981,2982,2983,2984,2985,2986,2987,2988,2989,2990,2991,2992,2993,2994,2995,2996,2997,2998,2999,3000,3001,3002,3003,3004,3005,3006,3007,3008,3009,3010,3011,3012,3013,3014,3015,3016,3017,3018,3019,3020,3021,3022,3023,3024,3025,3026,3027,3028,3029,3030,3031,3032,3033,3034,3035,3036,3037,3038,3039,3040,3041,3042,3043,3044,3045,3046,3047,3048,3049,3050,3051,3052,3053,3054,3055,3056,3057,3058,3059,3060,3061,3062,3063,3064,3065,3066,3067,3068,3069,3070,3071,3072,3073,3074,3075,3076,3077,3078,3079,3080,3081,3082,3083,3084,3085,3086,3087,3088,3089,3090,3091,3092,3093,3094,3095,3096,3097,3098,3099,3100,3101,3102,3103,3104,3105,3106,3107,3108,3109,3110,3111,3112,3113,3114,3115,3116,3117,3118,3119,3120,3121,3122,3123,3124,3125,3126,3127,3128,3129,3130,3131,3132,3133,3134,3135,3136,3137,3138,3139,3140,3141,3142,3143,3144,3145,3146,3147,3148,3149,3150,3151,3152,3153,3154,3155,3156,3157,3158,3159,3160,3161,3162,3163,3164,3165,3166,3167,3168,3169,3170,3171,3172,3173,3174,3175,3176,3177,3178,3179,3180,3181,3182,3183,3184,3185,3186,3187,3188,3189,3190,3191,3192,3193,3194,3195,3196,3197,3198,3199,3200,3201,3202,3203,3204,3205,3206,3207,3208,3209,3210,3211,3212,3213,3214,3215,3216,3217,3218,3219,3220,3221,3222,3223,3224,3225,3226,3227,3228,3229,3230,3231,3232,3233,3234,3235,3236,3237,3238,3239,3240,3241,3242,3243,3244,3245,3246,3247,3248,3249,3250,3251,3252,3253,3254,3255,3256,3257,3258,3259,3260,3261,3262,3263,3264,3265,3266,3267,3268,3269,3270,3271,3272,3273,3274,3275,3276,3277,3278,3279,3280,3281,3282,3283,3284,3285,3286,3287,3288,3289,3290,3291,3292,3293,3294,3295,3296,3297,3298,3299,3300,3301,3302,3303,3304,3305,3306,3307,3308,3309,3310,3311,3312,3313,3314,3315,3316,3317,3318,3319,3320,3321,3322,3323,3324,3325,3326,3327,3328,3329,3330,3331,3332,3333,3334,3335,3336,3337,3338,3339,3340,3341,3342,3343,3344,3345,3346,3347,3348,3349,3350,3351,3352,3353,3354,3355,3356,3357,3358,3359,3360,3361,3362,3363,3364,3365,3366,3367,3368,3369,3370,3371,3372,3373,3374,3375,3376,3377,3378,3379,3380,3381,3382,3383,3384,3385,3386,3387,3388,3389,3390,3391,3392,3393,3394,3395,3396,3397,3398,3399,3400,3401,3402,3403,3404,3405,3406,3407,3408,3409,3410,3411,3412,3413,3414,3415,3416,3417,3418,3419,3420,3421,3422,3423,3424,3425,3426,3427,3428,3429,3430,3431,3432,3433,3434,3435,3436,3437,3438,3439,3440,3441,3442,3443,3444,3445,3446,3447,3448,3449,3450,3451,3452,3453,3454,3455,3456,3457,3458,3459,3460,3461,3462,3463,3464,3465,3466,3467,3468,3469,3470,3471,3472,3473,3474,3475,3476,3477,3478,3479,3480,3481,3482,3483,3484,3485,3486,3487,3488,3489,3490,3491,3492,3493,3494,3495,3496,3497,3498,3499,3500,3501,3502,3503,3504,3505,3506,3507,3508,3509,3510,3511,3512,3513,3514,3515,3516,3517,3518,3519,3520,3521,3522,3523,3524,3525,3526,3527,3528,3529,3530,3531,3532,3533,3534,3535,3536,3537,3538,3539,3540,3541,3542,3543,3544,3545,3546,3547,3548,3549,3550,3551,3552,3553,3554,3555,3556,3557,3558,3559,3560,3561,3562,3563,3564,3565,3566,3567,3568,3569,3570,3571,3572,3573,3574,3575,3576,3577,3578,3579,3580,3581,3582,3583,3584,3585,3586,3587,3588,3589,3590,3591,3592,3593,3594,3595,3596,3597,3598,3599,3600,3601,3602,3603,3604,3605,3606,3607,3608,3609,3610,3611,3612,3613,3614,3615,3616,3617,3618,3619,3620,3621,3622,3623,3624,3625,3626,3627,3628,3629,3630,3631,3632,3633,3634,3635,3636,3637,3638,3639,3640,3641,3642,3643,3644,3645,3646,3647,3648,3649,3650,3651,3652,3653,3654,3655,3656,3657,3658,3659,3660,3661,3662,3663,3664,3665,3666,3667,3668,3669,3670,3671,3672,3673,3674,3675,3676,3677,3678,3679,3680,3681,3682,3683,3684,3685,3686,3687,3688,3689,3690,3691,3692,3693,3694,3695,3696,3697,3698,3699,3700,3701,3702,3703,3704,3705,3706,3707,3708,3709,3710,3711,3712,3713,3714,3715,3716,3717,3718,3719,3720,3721,3722,3723,3724,3725,3726,3727,3728,3729,3730,3731,3732,3733,3734,3735,3736,3737,3738,3739,3740,3741,3742,3743,3744,3745,3746,3747,3748,3749,3750,3751,3752,3753,3754,3755,3756,3757,3758,3759,3760,3761,3762,3763,3764,3765,3766,3767,3768,3769,3770,3771,3772,3773,3774,3775,3776,3777,3778,3779,3780,3781,3782,3783,3784,3785,3786,3787,3788,3789,3790,3791,3792,3793,3794,3795,3796,3797,3798,3799,3800,3801,3802,3803,3804,3805,3806,3807,3808,3809,3810,3811,3812,3813,3814,3815,3816,3817,3818,3819,3820,3821,3822,3823,3824,3825,3826,3827,3828,3829,3830,3831,3832,3833,3834,3835,3836,3837,3838,3839,3840,3841,3842,3843,3844,3845,3846,3847,3848,3849,3850,3851,3852,3853,3854,3855,3856,3857,3858,3859,3860,3861,3862,3863,3864,3865,3866,3867,3868,3869,3870,3871,3872,3873,3874,3875,3876,3877,3878,3879,3880,3881,3882,3883,3884,3885,3886,3887,3888,3889,3890,3891,3892,3893,3894,3895,3896,3897,3898,3899,3900,3901,3902,3903,3904,3905,3906,3907,3908,3909,3910,3911,3912,3913,3914,3915,3916,3917,3918,3919,3920,3921,3922,3923,3924,3925,3926,3927,3928,3929,3930,3931,3932,3933,3934,3935,3936,3937,3938,3939,3940,3941,3942,3943,3944,3945,3946,3947,3948,3949,3950,3951,3952,3953,3954,3955,3956,3957,3958,3959,3960,3961,3962,3963,3964,3965,3966,3967,3968,3969,3970,3971,3972,3973,3974,3975,3976,3977,3978,3979,3980,3981,3982,3983,3984,3985,3986,3987,3988,3989,3990,3991,3992,3993,3994,3995,3996,3997,3998,3999,4000] Rows for sql Calculate rows for the SQL result database. It assumes 1000 iterations. var a:= ['Business as usual','Recommend restrictions to salmon consumption','Stricter limits for fish feed pollutants','Restrictions to salmon consumption AND stricter fish feed limits']; var b:= floor((in2-1)/1000); '"'&in2&'";"1";"'&slice(a,a,b+1)&'"' 400,344,1 48,24 (param1) Doresult var a:= slice(Net_health_effects_i,net_health_effects_i.decision,4); '"'&3000+run&'";"'&1&'";"'&a&'";"'&run&'"' 56,88,1 48,24 2,102,90,476,224 param1 Result reporting This node is used to report results from this model. Just replace the variable name on the first row with the variable you want to report and calculate. This node calculates the result only for farmed salmon, and it looks all four possible deicisions along one dimension (not 2*2 table). var a:= sample(H1909); a:= a[Salmon1='Farmed salmon']; index decision:= ['Business as usual','Recommend restrictions','Stricter rules for feed','Both']; a:= array(decision,[ slice(slice(a,H1898,1),H1899,1), slice(slice(a,H1898,2),H1899,1), slice(slice(a,H1898,1),H1899,2), slice(slice(a,H1898,2),H1899,2)]); index statistics:= ['Mean','SD','0.01','0.025','0.05','0.25','0.5 (Median)','0.75','0.95','0.975','0.99']; array(statistics,[mean(a), sdeviation(a), getfract(a,0.01), getfract(a,0.025), getfract(a,0.05), getfract(a,0.25), getfract(a,0.5), getfract(a,0.75), getfract(a,0.95), getfract(a,0.975), getfract(a,0.99)]) 232,328,1 48,24 2,102,90,476,485 2,362,39,568,303,0,MIDM [Sys_localindex('DECISION'),Sys_localindex('STATISTICS'),Undefined,Undefined,Undefined,1] 2,D,4,2,0,0,4,0,$,0,"ABBREV",0 [Pollutant1,1,Sys_localindex('STATISTICS'),1,Sys_localindex('DECISION'),1] Pollutant exposure µg/kg/d Pollutant exposure per body weight per day. H1903*H1904/1000/BW 256,272,1 48,24 1,1,1,1,1,1,0,,1, 2,287,149,476,224 2,72,47,653,399,0,MIDM [Salmon1,Pollutant1] [1,0,0,0] Fish feed ktluser 11. Janta 2004 12:08 48,24 168,144,1 48,24 1,472,152,516,377,17 2,40,50,576,600 The concentrations of pollutants in fish feed have been reducing Rideout said that major feed companies have been able reduce toxin levels in their fish meal over the past several years by using substitute ingredients and less contaminated fish. ÒIt has been an issue that the industry is responding to. Feed companies have been working overtime to use high quality meals that are very low in contaminants and have twice the amount of omega-3sÉThe line is trending downward. We donÕt like having this in our product.Ó <ref>http://www.intrafish.com/articlea.php?articleID=41061&s=1</ref> Feed_backgr 328,48,1 64,36 2,102,90,720,332 Pollutant concentration in fish feed var a:= (if H1899='More actions' then impr_in_feed else 0); Feed_backgr*(1-a) 192,128,1 48,29 1,1,1,1,1,1,0,,1, 2,291,326,636,303,0,SAMP [Undefined,H1899,Undefined,Undefined,Undefined,1] [1,0,0,0] What has been done and what should be done to reduce pollutants in fish feed? H1902 192,240,1 60,44 Should we change fish feed instead of giving fish consumption advisories? Impr_in_feed+H1898+H1899 64,248,1 48,55 [Alias Should_we_change_fi1] Pollutant levels in fish feed after lower limits - Pollutant concentrations in fish feed can be further reduced from current levels. The estimate of 0-100 % with a gradual decrease in probability density is based on author judgement. It reflects rather a theoretical range of improvement than a realistic estimate. Triangular( 0, 0, 1 ) 64,128,1 48,38 2,102,90,476,409 2,515,277,416,303,0,MIDM 2,216,226,416,303,1,CDFP 52425,39321,65535 Fish feed background - This is a dummy variable only, because the actual concentrations in fish feed are not needed in the current model. 1 192,48,1 48,24 52425,39321,65535 Lower limits for pollutants in fish feed? Two options are assumed for fish feed regulations: 1) business as usual (BAU) with current legislation, and 2) More restrictive regulations for fish feed, resulting in reduction of pollutant levels in feed and consequently in salmon. (This is irrespective of any trends unrelated to the decision). ['BAU','More actions'] 280,144,1 48,32 1,1,1,1,1,1,0,,1, ['BAU','More actions'] Exposure- response function for omega3 jtue 12. Janta 2004 8:51 48,24 504,336,1 48,42 1,76,122,586,421,17 Exposure- response function for health benefit probability/(g/d) Exposure-response function where also the uncertainty about the population that benefits from omega-3 is taken into account. var a:= Erf_hcrude*(if All_or_chd=1 then 1 else F_chd_pati); array(cause_of_death3,[a]) 200,248,1 52,44 1,1,1,1,1,1,0,,1, 2,291,175,476,224 2,136,146,489,288,0,MEAN Graphtool:0 Distresol:10 Diststeps:1 Cdfresol:5 Cdfsteps:1 Symbolsize:6 Baroverlap:0 Linestyle:1 Frame:1 Grid:1 Ticks:1 Mesh:1 Scales:1 Rotation:45 Tilt:0 Depth:70 Frameauto:1 Showkey:1 Xminimum:0 Xmaximum:1 Yminimum:0 Ymaximum:1 Zminimum:0 Zmaximum:1 Xintervals:0 Yintervals:0 Includexzero:0 Includeyzero:0 Includezzero:0 Statsselect:[1,1,1,1,1,0,0,0] Probindex:[0,0.01,0.05,0.25,0.5,0.75,0.95,0.99,1] [0,0,0,0] Benefits: effects of omega-3 fatty acids on cardiovascular mortality Erf_hcrude 336,152,1 52,56 Does omega 3 help other people than CHD patients? All_or_chd 64,56,1 48,46 Does omega-3 help CHD patients or everyone? probability A large part of omega-3 benefit literature is based on studies on cardiac patients. This node reflects the uncertainty whether there is cardiac health benefit for everyone or only coronary heart disease (CHD) patients. The estimate is not based on data but the aim is to maximise uncertainty. Bernoulli( 0.5 ) 64,160,1 48,38 2,102,90,476,333 52425,39321,65535 Fraction of CHD patients among deaths fraction Fraction of coronary heart disease patients among the deaths. Current estimate is based on the fraction of cardiac deaths from total deaths in EEA countries, although there are cardiac deaths among non-CHD patients, and there are CHD patients with other causes of death. <ref>[http://www.who.int WHO data]</ref> 1.5717M/3.8664M 64,248,1 48,38 2,102,90,476,434 Dose-response of health benefit probability/(g/d) Dose-response function comes from secondary prevention trials reviewed by Din 2004 Table 1. The relative risk reductions are divided by the omega3 exposure in each study. A continuous distribution is used, and each study result is used as a quintile point for the distribution. Another review is Marckmann and Gronbaek 1999 that concluded that 0.6-0.9 g/d of omega-3 results in 40-60 % decrease in coronary heart disease mortality. The low estimate from this result was used (40% per 0.9 g/d). <ref>Din JN, Newby DE, Flapan AD. Science, medicine, and the future - Omega 3 fatty acids and cardiovascular disease - fishing for a natural treatment. British Medical Journal 2004; 328(7430):30-35. [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/Din_Omega3andCVD_BMJ2004.pdf Intranet file]</ref> <ref>Marckmann P, Gronbaek M. Fish consumption and coronary heart disease mortality. A systematic review of prospective cohort studies. European Journal of Clinical Nutrition 1999; 53(8):585-590.</ref> -Fractiles( [0/3.5,.325/1.5,.482/1.8,.297/0.85, 0.4/0.9 ]) 200,152,1 52,32 2,102,90,512,527 2,72,82,416,303,1,PDFP 52425,39321,65535 Highest omega3 dose with health benefit g/d Describes the amount of fish that is still beneficial when added to diet. After this limit, no extra benefit is assumed from omega-3 fatty acids. The value reflects both the physiological need of omega-3 and the current intake of omega-3 from other sources than salmon. Both parts of the estimate are complicated, and the latter varies from country to country. This might have implications to the decision if we could give country-wise recommendations of feed regulations. The estimate is based on author judgement. A rough idea about the magnitude comes from Din 2004, where the trials had 0.85 - 1.8 g/d of omega-3 with benefit but 3.5 g/d showed no benefit in a small trial where the population used reasonable amount of fish anyway. If the physiological limit is lower, the slope of the exposure-response function should be steeper. Another data comes from Albert 1998: the benefit may be limited to omega-3 doses <4.9 g/mo = 0.16 g/d. Markmann and Gronbaek concluded that 0.6-0.9 g/d is beneficial. <ref>Din JN, Newby DE, Flapan AD. Science, medicine, and the future - Omega 3 fatty acids and cardiovascular disease - fishing for a natural treatment. British Medical Journal 2004; 328(7430):30-35. [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/Din_Omega3andCVD_BMJ2004.pdf Intranet file]</ref> <ref>Albert CM, Hennekens CH, O'Donnell CJ, Ajani UA, Carey VJ, Willett WC et al. Fish consumption and risk of sudden cardiac death. Jama-Journal of the American Medical Association 1998; 279(1):23-28.</ref> <ref>Marckmann P, Gronbaek M. Fish consumption and coronary heart disease mortality. A systematic review of prospective cohort studies. European Journal of Clinical Nutrition 1999; 53(8):585-590.</ref> Triangular( .2, .5, 1 ) 496,64,1 48,38 2,135,16,476,598 52425,39321,65535 Cause of death3 ['Cardiovascular'] 200,304,1 56,12 Exposure- response function for pollutant risk Pieta 16. tamta 2004 1:32 48,24 136,336,1 48,51 1,141,229,418,300,17 Potency of pollutants (mg/kg/d)^Æ1 Potency of pollutants. Cancer slope factors (CSF) are used for cancer and Reference Dose (RfD) values are used for non-cancer endpoints. The data comes from the U.S.EPA. <ref>U.S.EPA. Guidance for Assessing Chemical Contaminant Data for Use in Fish Advisory. Volume 2: Risk Assessment and Fish Consumption Limits, 3rd Edition. 2000. Table 3-1. [http://www.epa.gov/waterscience/fish/guidance.html Open access Internet file]</ref> <ref>[http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/EPAFishAdvisory/ Intranet file]</ref> Table(Pollutant1,Self)( 16,50u, 1.1,250u, 156K,0, 2,20u ) ['Ca (CSF)','Non-Ca (RfD)'] 64,40,1 48,24 2,249,11,476,457 2,480,276,416,303,0,MIDM 2,103,144,416,303,0,MIDM [Alias Potency1] 65535,52427,65534 [Self,Pollutant1] [Self,Pollutant1] [1,0,0,0] Exposure-response function for pollutant risk (mg/kg/d)-1 The response assessment is restricted to cancer endpoints, because it is the more sensitive endpoint. potency[Potency='Ca (CSF)'] 64,128,1 48,38 1,1,1,1,1,1,0,,1, 2,102,90,476,399 2,499,251,655,303,0,MIDM [0,0,0,0] Is the exposure-response function affected by the target population and its background cancer risk? Should this be taken into account in the model? H1906 232,128,1 72,72 2,436,19,476,224 Salmon intake jtue 16. Janta 2004 12:54 48,24 320,208,1 48,24 1,81,109,605,297,17 Alternatives for farmed salmon: wild salmon, other fish, canola oil, flaxseed oil. ÒWeÕre telling people that if they want to reduce their risk of cancer, they should not eat more than one meal of farmed salmon a month,Ó Carpenter said. He added that the cancer risk from the toxins effectively cancels out the benefits of omega-3 fatty acids found in farmed salmon, which have not been proven to prevent or reduce the risk of cancer. ÒThere are other places to get omega-3s - wild salmon, other fish, canola oil, flaxseed oil,Ó Carpenter said. <ref>[http://www.intrafish.com/articlea.php?articleID=41061&s=1 Intrafish.com press release 9 Jan 2004]</ref> F 120,56,1 60,44 2,102,90,609,347 Current average consumption of salmon g/d Data comes from EPIC study looking at fish consumption in 10 European countries by gender. We take the minimum, the unweighed average and the maximum of these values in the distribution to represent uncertainty in population average fatty fish intake. All fatty fish is assumed to be salmon. <ref>Welch AA, Lund E, Amiano P, Dorronsoro M. Variability in fish consumption in 10 European countries. In: Riboli E, Lambert R, editors. Nutrition and lifestyle: opportunities for cancer prevention. Lyon: International Agency for Research on Cancer, 2002: 221-222. [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/RiboliNutritionLifestyle_IARC156_2002.pdf PDF of article] [http://ytoswww/yhteiset/Huippuyksikko/Tutkimus/Viljelylohi/Materiaali/ConsumptionOfFish.xls Data in Excel]</ref> Triangular( 7.5, 15.3, 31 ) 352,48,1 48,38 2,376,70,476,496 2,0,0,793,492,0,MEAN [Chance Welch_et_al_2002] [1,0,0,0] Fraction of farmed from total salmon use fraction Fraction of farmed salmon of total salmon consumption in Western Europe. The current estimate is based on author judgement after discussions with people from the Finnish Game and Fisheries Research Institute. Uniform( .8, 1 ) 248,48,1 52,44 2,102,90,476,367 Salmon intake g/d Intake of farmed and wild salmon after the two decisions (regulate fish feed pollutants / recommend restrictions for farmed salmon use) has been made. Although market salmon exists in the index, it is not used in this version of the model. Wild salmon use after restricting farmed salmon use has a triangular probability distribution. Min assumes the same relative decrease as in farmed salmon; mode assumes no change; max assumes that wild salmon intake increases so much that it totally compensates the decrease in farmed salmon use. Estimates are based on author judgement. The wild salmon production capacity is probably much than the max used for the variable. This overestimation causes bias towards smaller costs due to salmon use restrictions. Table(Salmon1,H1899,H1898)( F,(A*F), (F+Pollutant_scare),((A*F)+Pollutant_scare), W,Triangular((A*W),W,(W+((1-A)*F))), W,Triangular((A*W),W,(W+((1-A)*F))), 0,0, 0,0 ) 248,216,1 48,24 1,1,1,1,1,1,0,,1, 2,516,77,476,224 2,58,276,500,232,0,MIDM 2,248,258,443,303,1,PDFP [H1898,Salmon1] [Undefined,H1898,Undefined,Undefined,1] [Index Salmon1] [1,0,0,0] Farmed salmon baseline g/d Average farmed salmon consumption in Western Europe in the base case. Fraction_farmed*crude_salmon 248,144,1 48,29 2,469,178,476,300 Wild salmon baseline g/d Average wild salmon consumption in Western Europe in the base case. (1-Fraction_farmed)*crude_salmon 352,144,1 48,24 Farmed salmon use after recommendation fraction Farmed salmon use per baseline after a recommendation to restrict farmed salmon use. A uniform distribution between 1 American meal/month (227 g) and no change to baseline. If baseline is lower than the recommendation, no change occurs. var a:= 1*227/(365.25/12); var b:= min([a/f,1]); uniform(b,1) 128,144,1 56,36 2,264,94,476,252 2,552,65,424,320,0,MIDM Salmon consumption after feed limits g/d Change in farmed salmon use when fish feed is more strictly regulated. Consumer may consume more salmon, when pollutant problems are handled. However, there is a possibility of bad reputation ('There is a big problem, because authorities have to intervene'). The range overlaps zero to reflect this uncertainty. The expectation is slightly positive. The estimate is based on author judgement. Triangular( -1, 0.5, 1 ) 128,216,1 52,32 2,151,377,416,303,1,PDFP Welch et al 2002 <ref>Welch AA, Lund E, Amiano P, Dorronsoro M. Variability in fish consumption in 10 European countries. In: Riboli E, Lambert R, editors. Nutrition and lifestyle: opportunities for cancer prevention. Lyon: International Agency for Research on Cancer, 2002: 221-222. [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/RiboliNutritionLifestyle_IARC156_2002.pdf PDF of article] [http://ytoswww/yhteiset/Huippuyksikko/Tutkimus/Viljelylohi/Materiaali/ConsumptionOfFish.xls Data in Excel]</ref> 0 472,48,1 48,24 2,102,90,476,379 2,40,50,416,303,0,MIDM 65535,52427,65534 [Chance Crude_salmon] Recommend restricted farmed salmon consumption? A decision about whether a general recommendation should be given to restrict the consumption of European farmed salmon to one meal (227 g) per month or not (business as usual, BAU). ['BAU','Restrict farmed salmon use'] 424,144,1 76,32 1,1,1,1,1,1,0,,1, 2,102,90,476,354 ['BAU','Restrict farmed salmon use'] Omega3 content in salmon g/g Omega-3 fatty acid content in salmon. Estimate is from Din 2004 Table 2. In a previous version, there were other fish types as well: min is the lowest value, max is the highest value and mode is the average of the two. <ref>Din JN, Newby DE, Flapan AD. Science, medicine, and the future - Omega 3 fatty acids and cardiovascular disease - fishing for a natural treatment. British Medical Journal 2004; 328(7430):30-35. [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/Din_Omega3andCVD_BMJ2004.pdf Intranet file]</ref> array(Year3,[Uniform(0.0128,0.0215)]) 464,216,1 48,32 1,1,1,1,1,1,0,0,1,0 2,102,90,476,445 2,106,70,416,303,0,MIDM 2,472,482,416,303,0,MIDM 52425,39321,65535 [1,0,0,0] Omega3 exposure g/d Omega-3 fatty acid intake from salmon. H1904*H1907 384,272,1 48,24 1,1,1,1,1,1,0,,1, 2,102,90,476,224 2,57,94,723,303,0,MIDM [H1898,Salmon1] [1,0,0,0] [H1899,1,Year3,1,Salmon1,1,H1898,1] Pollutants in salmon jtok 16. tamta 2004 22:14 48,24 168,208,1 48,24 1,439,232,556,308,17 Pollutants in salmon Hites 2004 µg/kg Pollutant concentration data from Hites et al, Fig 2. <ref>Hites RA, Foran JA, Carpenter DO, Hamilton MC, Knuth BA, Schwager SJ. Global assessment of organic contaminants in farmed salmon. Science 2004; 303(5655):226-229. [http://ytoswww/yhteiset/Huippuyksikko/Kirjallisuus/Fish_and_health/HitesRA%26al_Science2004.pdf Intranet file]</ref> Table(Location1,Pollutant1)( 5.607,160.008,2.94m,50.904, 5.607,190.0095,2.5725m,48.177, 6.8085,136.6735,2.52m,47.268, 3.738,126.673,2.31m,49.995, 4.539,113.339,2.31m,41.814, 4.9395,150.0075,1.995m,37.269, 3.6045,93.338,2.52m,46.359, 4.4055,153.341,2.1m,34.542, 3.204,66.67,1.7325m,39.087, 2.5365,43.3355,1.89m,34.1784, 3.3375,70.0035,1.4175m,29.997, 2.8035,43.3355,1.449m,34.542, 2.5365,36.6685,1.5225m,34.542, 2.403,36.6685,1.3545m,32.724, 2.0025,30.0015,1.554m,26.361, 1.7355,32.0016,1.1025m,21.816, 2.0025,58.6696,840u,12.726, 1.2015,23.3345,1.344m,22.725, 1.2015,23.3345,945u,18.18, 1.2015,26.668,787.5u,15.453, 0.9345,14.6674,1.344m,21.816, 1.068,18.6676,892.5u,21.816, 1.2015,53.336,1.575m,7.272, 0.801,18.6676,1.155m,14.544, 0.534,20.6677,430.5u,15.453, 0.6675,23.3345,315u,7.272, 0.534,26.668,283.5u,8.181, 0.4806,20.001,315u,6.363, 0.4806,13.334,735u,9.09, 0.4005,16.6675,178.5u,4.545, 0.267,13.334,210u,3.636, 0.4005,14.0007,105u,3.636, 0.4005,15.3341,105u,3.636, 0.4005,12.0006,126u,3.636, 0.4005,10.0005,105u,2.727, 0.4005,10.0005,105u,1.818, 0.267,6.667,105u,1.818, 0.267,6.667,105u,1.818, 0.267,6.667,105u,1.818 ) 56,48,1 48,29 2,354,132,476,475 2,17,12,416,638,0,MIDM [Alias Pollutants_in_salmo2] 65535,52427,65534 [Pollutant1,Location1] [Pollutant1,Location1] , , , [1,0,0,0] Pollutants in salmon µg/kg Pollutant concentrations in salmon Poll_i_types*H1902 312,112,1 48,24 1,1,1,1,1,1,0,,1, [Salmon1,Pollutant1] [1,0,0,0] [Reg_poll,2,Pollutant1,1,Salmon1,1] Salmon type Data from Hites classified based on salmon type (farmed, wild, market) Table(Location1,Self)( 'Farmed salmon','Europe', 'Farmed salmon','Europe', 'Market salmon','Europe', 'Market salmon','Europe', 'Farmed salmon','Europe', 'Market salmon','Europe', 'Market salmon','Europe', 'Market salmon','Europe', 'Farmed salmon','North America', 'Market salmon','North America', 'Farmed salmon','North America', 'Market salmon','North America', 'Farmed salmon','North America', 'Market salmon','North America', 'Market salmon','North America', 'Market salmon','North America', 'Market salmon','North America', 'Market salmon','North America', 'Market salmon','North America', 'Market salmon','North America', 'Farmed salmon','North America', 'Farmed salmon','South America', 'Wild salmon','North America', 'Market salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America', 'Wild salmon','North America' ) ['Type','Region'] 56,112,1 48,24 2,102,90,476,395 2,72,82,416,614,0,MIDM 52425,39321,65535 [Self,Location1] [Self,Location1] Pollutant concentration in f/w/m salmon µg/kg Dieldrin, toxaphene, PCB, and dioxin concentrations in farmed, wild, and market salmon. Triangular probability distribution is used for each pollutant and salmon type. Estimates are based on data from Hites 2004. Parameters min, mode, and max are the minimum, average, and maximum of Hites' data, respectively. var typ:= salmon_type[salmon_type='Type']; var a:= (if Salmon1=typ then poll_salmon_hites else 0); var b:= (if Salmon1=typ then 1 else 0); var c:= max(a,location1); var d:= sum(a,location1)/sum(b,location1); var e:= (if Salmon1=typ then poll_salmon_hites else 1M); var f:= min(e,location1); triangular(f,d,c) 192,113,1 48,38 2,36,56,476,309 2,489,195,416,303,1,PDFP Other parts jtue 28. Junta 2004 18:03 48,24 488,40,1 48,24 1,0,1,1,1,1,0,,0, 1,40,0,517,300,17 Pollutants per types and region µg/kg Triangular probability distribution for concentrations indexed by salmon type AND THE THREE REGIONS based on data from Hites 2004. Min is the minimum of Hites, Max is the maximum of Hites, and Mode is the average of Hites. Currently not used, but probably we should look at European values, because the consumption and population data comes from Europe. var typ:= salmon_type[salmon_type='Type']; var reg:= salmon_type[salmon_type='Region']; var a:= (if typ=Salmon1 and reg=Region then poll_salmon_hites else 0); var b:= sum((if typ=Salmon1 and reg=Region then 1 else 0),location1); var c:= max(a,location1); var d:= if b>0 then sum(a,location1)/b else 0; var e:= (if Salmon1=salmon_type then poll_salmon_hites else 1M); var f:= min(e,location1); d 56,32,1 48,29 2,102,90,476,402 2,72,82,451,390,0,MIDM [Salmon1,Pollutant1] Region The three regions considered in Hites et al 2004. ['Europe','North America','South America'] 56,72,1 48,12 Concentrationparameters for model description µg/kg Triangular probability distribution for concentrations indexed by salmon type based on data from Hites et al, 2004. Min is the minimum, Mode is the average, and Max is the maximum calculated for each salmon type (farmed, wild, and market) separately. var typ:= salmon_type[salmon_type='Type']; var a:= (if Salmon1=typ then poll_salmon_hites else 0); var b:= (if Salmon1=typ then 1 else 0); var c:= max(a,location1); var d:= sum(a,location1)/sum(b,location1); var e:= (if Salmon1=typ then poll_salmon_hites else 1M); var f:= min(e,location1); index x:=['Min','Mode','Max']; array(x,[f,d,c]) 56,136,1 48,38 2,102,90,476,365 2,561,214,416,303,0,MIDM 1,D,4,2,0,0 Total mortality W Europe cases/a Total mortality in European Economic Area countries (386.63 million inhabitants) <ref>[http://www.who.int WHO data]</ref> array(Year3,[3.8664M]) 136,424,1 48,32 1,1,1,1,1,1,0,,1, 2,102,90,476,330 65535,52427,65534 [1,1,0,1] CHD mortality W Europe cases/a Coronary heart disease mortality in European Economic Area countries (386.63 million inhabitants). The estimate consists of acute myocardial infarction and other ischaemic heart diseases (ICD 10: 270, 279). <ref>[http://www.who.int WHO data]</ref> array(Year3,[615.3k]) 504,416,1 48,32 1,1,1,1,1,1,0,,1, 2,102,90,476,224 65535,52427,65534 [1,1,0,1] Should we change fish feed instead of giving fish consumption advisories? 1 568,145,1 48,61 Should_we_change_fis Log v4 1 760,144,1 52,12 65535,54067,19661 Loki_v4 Pollutant risk is much smaller than the net health benefit of farmed salmon H1901+H1900 256,488,1 60,51 65535,65532,19661 Scientific uncertainties related to recommendations are unimportant H1898+V3 424,51,1 64,51 65535,65532,19661 Some scientific and political uncertainties related to feed limits are important H1899+V1 280,48,1 72,42 65535,65532,19661 URN:NBN:fi-fe20042774 DC-attribute with refinement Scheme (if any) Value Title Risk benefit analysis of eating farmed salmon Creator.personalName Tuomisto, Jouni T Creator.personalName Tuomisto, Jouko Creator.personalName Tainio, Marko Creator.personalName Niittynen, Marjo Creator.personalName Verkasalo, Pia Creator.personalName Vartiainen, Terttu Creator.personalName Kiviranta, Hannu Creator.personalName Pekkanen, Juha Subject risk benefit analysis Subject persistent organic pollutants Subject omega-3 fatty acids Subject MeSH polychlorinated biphenyls Subject MeSH salmon Subject MeSH risk assessment Subject MeSH fatty acids, omega-3 Subject UDC 614 Public health and hygiene. Description.abstract In their Report ÒGlobal assessment of organic contaminants in farmed salmon,Ó R. A. Hites and co-workers analyzed wild and farmed salmon samples from North and South America and Europe for organic pollutants (9 Jan. 2004, p. 226). The authors conclude that, because of chemical contaminants, farmed salmon should not be eaten more often than 0.25 to 1 times per month. However, the model used does not take into account any beneficial effects of eating fish. We analyzed both risks and benefits. We also performed a value-of-information analysis to see which uncertainties were relevant for decision-making. This is the version 4 of the model calculating risks and benefits of farmed salmon. (c) Copyright Kansanterveyslaitos (KTL; National Public Health Institute, Finland). Publisher Kansanterveyslaitos Date.issued W3C-DTF 2004-07-23 Type DCMIType Software Format IMT text/xml Format.medium computerFile Format 115 kB Identifier http://www.ktl.fi/risk Identifier URN URN:NBN:fi-fe20042774 Language ISO639-2 en Relation.hasPart URL http://www.sciencemag.org/cgi/content/full/305/5683/476 Rights Copyright Kansanterveyslaitos, 2004 0 744,176,1 80,12 2,102,90,523,439 65535,54067,19661 Year3 [2000] 504,456,1 48,12 Inputs for RDB The node Variables_to_be_saved in the module RDB-connection.ANA from Heande should contain the following: Var_name Probabilistic? 'H1898' 0 'H1899' 0 'H1900' 1 'H1901' 1 'H1902' 1 'H1903' 1 'H1904' 1 'H1905' 1 'H1906' 0 'H1907' 1 'H1908' 1 'H1909' 1 'H1910' 0 'H1911' 0 'H1912' 1 760,32,1 48,24 65535,54067,19661