<p>Dear suman</p>
<p>You might find interestingbthe following page<br>
<a href="http://www.cawcr.gov.au/projects/verification/">http://www.cawcr.gov.au/projects/verification/</a></p>
<p>Best regards<br>
Luis</p>
<div class="gmail_quote">El 14/10/2014 07:00, "Suman Maity" <<a href="mailto:suman.buie@gmail.com">suman.buie@gmail.com</a>> escribió:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Dear Herbert<br></div>Thank you very much for your reply. But I couldn't understand from the grads code how to make the contigency table. Will you please elaborate by giving a small example how to create it and why? so that I can calculate it by hand for a small datasets.<br></div>Thanking you once again,<br><br></div>Regards<br></div>'S Maity' <br><div class="gmail_extra"><br clear="all"><div><div dir="ltr"><div><i><font style="background-color:rgb(255,153,0)"><span style="background-color:rgb(255,255,255)"><font face="comic sans ms,sans-serif"><b><font size="4">*******************************************************<br>Suman Maity</font></b></font></span></font></i></div><div><font size="1"><font style="background-color:rgb(255,153,0)"><font size="4"><span style="background-color:rgb(255,255,255)"><font face="comic sans ms,sans-serif">Research Scholar<br></font></span></font></font></font></div><div><font size="1"><font style="background-color:rgb(255,153,0)"><font size="4"><span style="background-color:rgb(255,255,255)"><font face="comic sans ms,sans-serif">Center for Ocean, Rivers, Atmosphere and Land Sciences(CORAL)<br></font></span></font></font></font></div><div><font size="1"><font style="background-color:rgb(255,153,0)"><font size="4"><span style="background-color:rgb(255,255,255)"><font face="comic sans ms,sans-serif">Indian Institute Of Technology Kharagpur.<br></font></span></font></font></font></div><div><i><b><font style="font-family:comic sans ms,sans-serif" size="4">*******************************************************</font></b><b style="font-family:comic sans ms,sans-serif"><font size="4"> </font></b></i><br></div><div style="text-align:left"><i><font color="#3366ff" size="1"><span style="font-size:12pt;font-family:'Times New Roman','serif'">
</span> </font></i></div><div style="display:inline"></div></div></div>
<br><div class="gmail_quote">On Tue, Oct 14, 2014 at 12:38 PM, Maat, Herbert ter <span dir="ltr"><<a href="mailto:herbert.termaat@wur.nl" target="_blank">herbert.termaat@wur.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div link="blue" vlink="purple" lang="EN-GB">
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d">Dear Suman,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d">My suggestion would be build a contingency table and then import the contingency table into a spreadsheet program. I have created a script to make a contingency
table. You can find the most important part of this script below.<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d">Kind regards,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d">Herbert ter Maat<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">* modeldata<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">'use d:/UserData/Herbert/Projects/KvK/KKF/FutureWeather/GrADS/fw38-LS-2006-07-15-000000-g3.ctl'<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">* Observation data<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">filen1='d:\UserData\Herbert\Projects\KvK\KKF\FutureWeather\DATA_WUR\Rain2Grads\PrecipObs_200608_day.dat'<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">t0=9<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">rok=0<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">* Threshold value of detection<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">thres=0.5<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">result=read(filen1)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">'set lat 51.5'<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">‘set lon 5’<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">'set gxout print'<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">while (rok=0)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> 'set t 't0<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> result=read(filen1)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> rok=sublin(result,1)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> rok=subwrd(rok,1)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> if (rok!=0); say 'end of file'; break; endif<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> lin=sublin(result,2)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> stat2=subwrd(lin,3)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> say stat2<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> t1=t0<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> 'set t 't1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> a=0<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> b=0<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> c=0<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> d=0<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> j=1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> stat1=stat2<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> lat=subwrd(lin,2)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> lon=subwrd(lin,1)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> while (j < 30)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> value1=subwrd(lin,j+4)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> t1=t1+24<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> 'set t 't1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> 'd gr2stn(totpcp-totpcp(t-24),'lon','lat')'<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> reslin2=sublin(result,3)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> value2=subwrd(reslin2,1)
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> j=j+1<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> say j<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> if(value1>=thres & value2>=thres); a=a+1; endif<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> if(value1>=thres & value2<thres); c=c+1; endif<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> if(value1<thres & value2>=thres); b=b+1; endif<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> if(value1<thres & value2<thres); d=d+1; endif<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> endwhile<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> record = stat1','a','b','c','d<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"> res = write(contingency.txt,record)<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Courier New";color:#1f497d">endwhile<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Verdana","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"" lang="EN-US">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"" lang="EN-US"> <a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a> [mailto:<a href="mailto:gradsusr-bounces@gradsusr.org" target="_blank">gradsusr-bounces@gradsusr.org</a>]
<b>On Behalf Of </b>Suman Maity<br>
<b>Sent:</b> vrijdag 10 oktober 2014 19:08<br>
<b>To:</b> GrADS Users Forum; regcnet<br>
<b>Subject:</b> [gradsusr] How to calculate Equitable Threat Score (ETS)<u></u><u></u></span></p><div><div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Dear all<u></u><u></u></p>
</div>
<p class="MsoNormal">This is Suman Maity from Indian Institute of Technology Kharagpur, India. Frequently we are validating our model output with some observation datasets. For the validation purpose people are using different skill scores. Equitable Threat
Score (ETS) is such a score for verification. My query is how to calculate that? for example I have 2007 (June-July-August-September) model simulated daily rainfall data and GPCP rainfall data. If I want to calculate the ETS then how to do that?<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">Please help in this regard. Any suggestion will be highly appreciable.
<u></u><u></u></p>
</div>
<p class="MsoNormal">Regards<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">'S Maity'<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><b><i><span style="font-size:13.5pt;font-family:"Comic Sans MS";background:white">*******************************************************<br>
Suman Maity</span></i></b><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:13.5pt;font-family:"Comic Sans MS";background:white">Research Scholar</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:13.5pt;font-family:"Comic Sans MS";background:white">Center for Ocean, Rivers, Atmosphere and Land Sciences(CORAL)</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:13.5pt;font-family:"Comic Sans MS";background:white">Indian Institute Of Technology Kharagpur.</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><b><i><br>
</i></b><b><i><span style="font-size:13.5pt;font-family:"Comic Sans MS"">*******************************************************
</span></i></b><u></u><u></u></p>
</div>
<p class="MsoNormal"><i><span style="font-size:7.5pt;color:#3366ff"> </span></i><u></u><u></u></p>
</div>
</div>
</div></div></div>
</div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div></div>
<br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div>