<html><table id="201179" bgcolor="#ffffff" border="0" height="400" width="100%"><tbody><tr><td valign="top" width="100%"><font face="Verdana" size="1"><font face="Verdana" size="2"><font size="4">Dear all,<br><br>I have used following script to see the the linear trend in data. </font></font><font face="Verdana" size="4">I am attaching one fig. for the reference.<br>Yellow line is representing mean <br>Green line represents trend<br><br>
</font><font face="Verdana" size="2"><font size="4"><font color="#000080">1. This script is working fine but in the result kinks are coming in the line representing trend, How to remove<br>it and make the trend line smooth ?<br>2. How to check the statsistical significance level in the trend (5% amd 1% significance level) ?<br><br>Script :<br><font color="#003366">*assuming this is the time limit<br>'reinit'<br>'enable print print.meta'<br>'open test1.ctl'<br>'set vpage 0 8 7.5 11'<br>'set grid off'<br>'set grads off'<br>'set t 1 57'<br>'set lat 18.96'<br>'set lon 72.83'<br>'set z 1'<br>'define yvar = ave(rf,t=1,t=57)'<br>'set vrange 0 1000'<br>'set ylabs 0| 100| 200| 300| 400| 500| 600| 700| 800| 900| 1000|'<br>'d rf'<br>* make a variable 1, 2, ..., 24 to calculate m and c of y = mx + c<br>'xvar = rf*0 + 1'<br>'xvar = const(xvar,1,-rf)'<br>'xvar = 58 - sum(xvar,t+0,t=57)'<br>'set t 1'<br>'define p = sum(xvar,t=1,t=57)'<br>'define q = sum(rf,t=1,t=57)'<br>'define r = sum(xvar*rf,t=1,t=57)'<br>'define s = sum(xvar*xvar,t=1,t=57)'<br>'define rn = ave(xvar*rf,t=1,t=57)'<br>'define n = r/rn'<br>'define amp = (n*r - p*q)/(n*s - p*p)'<br>'define dd = (q - amp*p)/n'<br>'set t 1 57'<br>'define yvar = ave(rf,t=1,t=57)'<br>'define aa= (amp*xvar + dd)'<br>'set cmark 0'<br>'d aa'<br>'set cmark 0'<br>'d yvar'<br>'draw ylab mm'<br>'draw xlab Year' <br>'draw title 1SD_Ahm blk_Act_TS ylw_mean grn_trend'<br>'close 1' <br><br>**__________________________________________________________________________________<br><br>'open test2.ctl'<br>'set vpage 0 8 4 7.5'<br>'set grid off'<br>'set grads off'<br>'set t 1 57'<br>'set lat 18.96'<br>'set lon 72.83'<br>'set z 1'<br>'define yvar = ave(rf,t=1,t=57)'<br>*'set vrange 0 1000'<br>'set ylabs 0| 100| 200| 300| 400| 500| 600| 700| 800| 900| 1000|'<br>'d rf'<br>* make a variable 1, 2, ..., 24 to calculate m and c of y = mx + c<br>'xvar = rf*0 + 1'<br>'xvar = const(xvar,1,-rf)'<br>'xvar = 58 - sum(xvar,t+0,t=57)'<br>'set t 1'<br>'define p = sum(xvar,t=1,t=57)'<br>'define q = sum(rf,t=1,t=57)'<br>'define r = sum(xvar*rf,t=1,t=57)'<br>'define s = sum(xvar*xvar,t=1,t=57)'<br>'define rn = ave(xvar*rf,t=1,t=57)'<br>'define n = r/rn'<br>'define amp = (n*r - p*q)/(n*s - p*p)'<br>'define dd = (q - amp*p)/n'<br>'set t 1 57'<br>'define yvar = ave(rf,t=1,t=57)'<br>'define aa= (amp*xvar + dd)'<br>'set cmark 0'<br>'d aa'<br>'set cmark 0'<br>'d yvar'<br>'draw ylab mm'<br>'draw xlab Year' <br>'draw title 2SD_Ahm blk_Act_TS ylw_mean grn_trend'<br>'close 1' <br><br>**_________________________________________________________________________________<br><br>'open test3.ctl'<br>'set vpage 0 8 0.5 4'<br>'set grid off'<br>'set grads off'<br>'set t 1 57'<br>'set lat 18.96'<br>'set lon 72.83'<br>'set z 1'<br>'define yvar = ave(rf,t=1,t=57)'<br>'set vrange 0 1000'<br>'set ylabs 0| 100| 200| 300| 400| 500| 600| 700| 800| 900| 1000|'<br>'d rf'<br>* make a variable 1, 2, ..., 24 to calculate m and c of y = mx + c<br>'xvar = rf*0 + 1'<br>'xvar = const(xvar,1,-rf)'<br>'xvar = 58 - sum(xvar,t+0,t=57)'<br>'set t 1'<br>'define p = sum(xvar,t=1,t=57)'<br>'define q = sum(rf,t=1,t=57)'<br>'define r = sum(xvar*rf,t=1,t=57)'<br>'define s = sum(xvar*xvar,t=1,t=57)'<br>'define rn = ave(xvar*rf,t=1,t=57)'<br>'define n = r/rn'<br>'define amp = (n*r - p*q)/(n*s - p*p)'<br>'define dd = (q - amp*p)/n'<br>'set t 1 57'<br>'define yvar = ave(rf,t=1,t=57)'<br>'define aa= (amp*xvar + dd)'<br>'set cmark 0'<br>'d aa'<br>'set cmark 0'<br>'d yvar'<br>'draw ylab mm'<br>'draw xlab Year' <br>'draw title 3SD_Ahm blk_Act_TS ylw_mean grn_trend'<br>'close 1' <br>'print'<br>'disable print'<br>'!gxps -c -i print.meta -o TS_Ahm_June_123SD_magnitude_mm.ps'<br>'!rm print.meta'</font><br><br></font></font><br><font size="4">Thanks,<br></font></font></font> <b><font size="3"><font color="#103712">Dhaval
</font><br><font color="#1c5d20">
Indian Centre for Climate and Societal Impacts Research (ICCSIR)
<br>
Ahmedabad Education Society compound,
<br>
Opp. H.L. Commerce College,
<br>
Navrangpura,
<br>
Ahmedabad 380 009,
<br>
Gujarat,
<br>
India
<br>
Contact No. : +91(079)40045462
</font><br><br></font></b>
<br>
<br>
</td></tr></tbody></table></html>