dear @Reeza   i did not do that, i tried to insert them carefully  ... View more

Dear @Tom   thanks for your help, i removed the comments to be on the safe side. ... View more

dear @ballardw   thanks for your valuable response, exactly when i fixed this issue the code run better than before but, it did not work completely, here is when the macro fails based on the log:   ERROR: (execution) Invalid argument to function. operation : LGAMMA at line 18830 column 1 operands : _TEM1001 _TEM1001 1 row 1 col (numeric) 0 statement : ASSIGN at line 18830 column 1 traceback : module CV at line 18830 column 1 NOTE: Paused in module CV. MPRINT(GOLDEN): npar1=res1[3]; ERROR: (execution) Matrix has not been set to a value. ... View more

As I said I just copy paste it from the paper and I am not expert on coding macro.   Thank you for your time and help ... View more

I run the previous attached file first "GWNBR_Alan_SAS_Macro", then run the following code:   data gwrmodel; input centdgy centdgx Indp dep; infile cards dlm =' '; cards; 34.1358 -117.924 10.9692 12 34.1761 -118.315 10.722 8 33.8236 -118.188 10.6537 2 34.1994 -118.535 10.6964 5 34.0669 -117.751 10.952 15 33.9292 -118.21 10.6143 6 34.015 -118.06 10.7066 4 34.0672 -118.226 10.5749 2 34.0833 -118.107 10.9292 6 34.3875 -118.535 10.7737 0 34.0508 -118.454 10.7889 2 33.9236 -118.371 11.3734 3 34.6901 -118.133 11.5754 22 33.8203 -117.913 11.4777 21 33.9261 -117.951 11.6989 4 33.6275 -117.69 10.9151 3 33.6751 -117.927 10.1761 0 33.85 -116.543 10.514 2 34.0103 -117.426 10.6558 8 33.7083 -117.243 10.7527 4 33.925 -116.876 11.8476 12 33.7146 -116.234 11.6559 7 33.674 -117.321 11.5871 26 33.9756 -117.333 10.0677 2 34.1036 -117.629 10.9899 3 34.2439 -117.274 10.3513 2 34.0994 -117.504 10.0889 2 34.1072 -117.274 9.1505 3 34.0667 -117.151 11.0462 10 33.9846 -117.513 10.8081 3 34.2623 -117.189 11.0214 6 34.1439 -117.851 9.798 4 ;run; %golden(data=gwrmodel,y=dep,x=indp,lat=centdgy,long=centdgx,method=ADAPTIVE1,type=aic,gwr=poisson,out=band); %gwnbr(data=gwrmodel,y=dep,x=indp,lat=centdgy,long=centdgx,h=9.3796134,gwr=poisson,method=ADAPTIVE1,geocod=indp,out=gwr);   /***************/ here is the part of the log:   6237 %golden(data=gwrmodel,y=dep,x=indp,lat=centdgy,long=centdgx,method=ADAPTIVE1,type=aic,gwr=poi - 180 6237! sson,out=band); WARNING: Apparent invocation of macro GOLDEN not resolved. ERROR 180-322: Statement is not valid or it is used out of proper order. 6238 %gwnbr(data=gwrmodel,y=dep,x=indp,lat=centdgy,long=centdgx,h=9.3796134,gwr=poisson,method=ADA 6238! PTIVE1,geocod=indp,out=gwr); NOTE: IML Ready NOTE: Line generated by the invoked macro "GWNBR". 11 par=choose(par<1E-10,1E- 10,par); - 22 76 ERROR 22-322: Syntax error, expecting one of the following: #, ##, &, (|, *, **, +, ',', -, /, //, :, <, <=, <>, =, >, ><, >=, @, [, ^=, |, ||. ERROR 76-322: Syntax error, statement will be ignored. NOTE: Line generated by the invoked macro "GWNBR". 13 hess=choose(abs(hess)<1E -     Thank you ... View more

Dear @Reeza   thanks for your prompt reply, please find it attached ... View more

Dear @FreelanceReinh, Yas i just got it eventually where I should check the log not the list for the significant variables, that was my mistake.   I need that to see the candidate variables to be included in the model, which your code is enough for me. however, if you able to manage including all variables and get the significant variables in present of the others will ease a alot of time and become good contribution to the PROC GENMOD from my point view. It just doing similar the option (Selection = forward) in logistic regression.   thanks again for your help and cooperation ... View more

Dear @RW9 thanks for your explanation, i will try it ... View more

Dear @FreelanceReinh i have tried many times but could not got it; Please find my sample below:     data have; infile cards dlm=" "; input res trans commn commr busn totacc facc fsacc exp; cards; 0.5106 0.00592 0.00204 0.04674 0.01549 1561 8 45 14.6401 0.54689 0.00368 0.00104 0.01984 0.00878 1293 6 50 14.5009 0.55426 0.00755 0.00368 0.01085 0.0197 1283 6 27 13.7799 0.78301 0.00229 0.00246 0.03242 0.02246 1399 3 25 13.9239 0.52035 0.00555 0.00813 0.00985 0.0148 835 3 24 13.8535 0.52406 0.00045 0.00179 0.01541 0.00757 784 5 15 14.4113 0.54079 0.01462 0.05013 0.01628 0.02279 936 4 13 14.0026 0.0878 0.00312 0.00011 0.00942 0.0025 1034 5 16 14.1336 0.04896 0.00157 0.04674 0.11843 0.0066 1179 8 25 13.9206 0.41158 0.00399 0.00734 0.25723 0.02834 1893 2 31 13.7067 0.62409 0.00807 0.00253 0.03514 0.00633 1699 7 49 13.821 0.34779 0.0289 0.05937 0.04541 0.02666 1900 9 39 14.9494 0.57421 0.00915 0.00181 0.03558 0.00675 1402 3 42 13.9342 0.68064 0.00395 0.0007 0.08047 0.02416 376 0 5 12.6189 0.38445 0.01718 0.00134 0.05237 0.0705 1623 11 42 14.7249 0.71162 0.00409 0.00075 0.02756 0.0264 1433 2 14 14.0926 0.51897 0.00997 0.00074 0.04778 0.03209 5088 11 71 14.894 0.63653 0.00411 0.00374 0.0487 0.0202 2044 7 51 14.6865 0.31564 0.00522 0.00018 0.01399 0.04705 1005 3 20 13.9002 0.35732 0.00522 0.00208 0.04063 0.02346 5741 34 157 15.3313 0.75071 0.01419 0.00551 0.02774 0.01832 980 4 25 13.9193 0.30134 0 0.00008 0.00945 0.00451 671 7 27 14.004 0.24283 0.00123 0.01685 0.00348 0.00336 564 2 10 14.074 0.46881 0.00231 0.00829 0.01293 0.03569 700 1 5 13.7692 0.03949 0.01613 0 0.00055 0.00081 407 1 6 13.2234 0.04287 0.00561 0.04142 0.01772 0.01949 4436 34 99 15.0589 0.32452 0.01828 0.00561 0.08055 0.01413 9017 40 162 15.1009 0.01545 0.00422 0.00205 0.00982 0.03275 2866 39 114 15.7243 0.05259 0.00144 0.08136 0.00561 0.01251 227 1 5 13.9364 0.0216 0.00012 0.00009 0.02055 0.012 2757 11 61 14.2239 0.01964 0.00341 0.00691 0.00493 0.00128 41 0 0 14.2912 0.02977 0.00424 0.01542 0.00407 0.00895 1440 10 35 14.6687 0.49944 0.01248 0.01842 0.01607 0.00833 332 1 4 11.7664 0.57083 0.0105 0.00156 0.12008 0.0259 759 3 13 13.136 0.00176 0.0002 0.01181 0.00386 0.00292 1784 43 126 15.3832 0.01671 0.00037 0.00009 0.00096 0.00116 234 1 6 11.4916 0.02278 0.00003 0.00231 0.0088 0.00024 3272 48 142 14.6535 0.07297 0 0.01172 0 0 62 0 1 9.2823 0.44657 0.06945 0.02416 0.05744 0.12574 1536 3 22 14.0376 0.00054 0 0 0 0.00008 655 14 39 12.7758 0.69536 0.01753 0.00429 0.03601 0.01435 1934 5 37 14.3894 0.49895 0.00396 0 0.26858 0.01808 354 4 14 12.1618 0.02764 0.00178 0.04852 0.01417 0.02677 5776 30 134 15.0845 0.00177 0.0002 0.00037 0.00006 0.00003 126 2 2 11.8863 0.45239 0.01005 0.00076 0.00014 0.01328 1246 4 31 14.4243 0.13396 0.00224 0.01808 0.0079 0.00452 5140 45 240 15.2573 0 0.0708 0 0 0 1272 9 45 13.0868 0.55166 0.01672 0.00734 0.04716 0.00758 2015 12 61 14.5673 0.49025 0.01888 0.01895 0.04154 0.07337 1777 5 18 14.0023 0.10641 0.00109 0.00371 0.03058 0.00735 3680 17 111 15.4036 ;run; data _null_; array dep totacc facc fsacc; array indp res trans commn commr busn; do i=1 to dim(dep); call execute('ods output ParameterEstimates(persist=proc)=est_' || vname(dep[i]) ||';'); do j=1 to dim(indp); call execute( 'proc genmod data=have;' || 'model ' || vname(dep[i]) || ' = exp ' || vname(indp[j]) || ' /d=nb; run;'); end; call execute('ods output close;'); call execute('proc sql noprint; select parameter into :indplist_' || vname(dep[i]) || ' separated by " " from est_' || vname(dep[i]) || ' where upcase(parameter) not in ("INTERCEPT", "EXP", "DISPERSION") & .<ProbChiSq<0.05; quit;'); end; run; %put &=indplist_totacc ; %put &=indplist_facc ; %put &=indplist_fsacc; ... View more

Ithink that because of the follwoing: first let explain here what i got in the list in terms of the (parameter):   Intercept exp indp1 Dispersion   usualy the (Intercept) or the (exp) P-value turns to be significant (p-value < 0.05), so it will be included to the list even if the (indp1) insignificant.     ... View more

what i meant is that the (%put &=indplist_indp1) contains all independant variables (indp1 .. indp16) regardless of their P-value. I thought the list contains only significant variables (if the P-value < 0.05). ... View more

Dear @FreelanceReinh you are doing very good coding. now i got it and i like it although the list contains all veriables and thier values either significant or insignificant. ... View more

i am not sure about last part (%put &=indplist_dep1; /* replace dep1 by the name of the first dep. variable */ ). i rename my dependant variable to (dep, dep1 ... dep7) and the indepdendant variables to (indp, indp1 .... indp16) just to apply the code. when i run it this time it showed the same results of old one. ... View more