Hi Ksharp, Sorry for not being clear. "Mutual Information" is a statistics from Information Theory. Its is explained in Wikipedia. Regards, Amit ... View more

Hi Peter, I appreciate all your help. Thanks for replying to my question. There are so many options in the proc glmselect for selection=lasso, I am lost. In your experience, what are some of the options that have worked for binary response variable. If you could give me some sample code that would be great. Regards, Amit ... View more

It makes sense. I will take it into consideration. Thanks for all your help. Regards, Amit ... View more

Hi , Thanks for the reply. I will look into the reference manual. From your comment regarding 1.3M observations, I gather that logistics regression will fail to provide good predictive models. ( I am using train and validate datasets for model building) 1> Is logistics regression not suited for large datasets? 2> What methods would provide better predictive power for large datasets? I appreciate all your help. Regards, Amit ... View more

Hi Rob, I just deleted those two zipcodes and it worked. But I get some warning in proc gmap. Can you tell me what the issue is. Thanks, Amit 023 proc gmap data=map map=map; 1024 id state; 1025 choro state / anno=dot nolegend; 1026 run; WARNING: Some observations were discarded when charting STATE. Only first matching observation was used. Use STATISTIC= option for summary statistics. NOTE: PROBLEM IN OBSERVATION 1 - DATA SYSTEM REQUESTED, BUT VALUE IS NOT ON GRAPH 'Y' ... View more

Hi Rob, There are more than a million data points. The error happens at the proc gproject. How can I isolate the problem zip codes. Does proc gproject has a way to output error records. Thanks, Amit ... View more

Hi Rob, This is the entire code: Thanks, Amit * Temp contains customer zip; proc sort data=temp nodupkey; by zip; run; data longlat; merge temp(in=mine) sashelp.zipcode(rename=(x=long y=lat)keep=x y zip); by zip; if mine; x=atan(1)/45*long; y=atan(1)/45*lat; x=-x; keep zip x y; if x eq . then delete; run; /* Create an annotate data set to place a symbol at the ZIP code locations. */ data anno; set longlat; retain xsys ysys '2' function 'label' size .75 flag 1 when 'a'; style='special'; text='M'; color='red'; output; run; /* Combine the map data set with the annotate data set. */ data all; set maps.states(where=(state not in(2 15 72))) anno; run; /* Project the combined data set. */ proc gproject data=all out=allp; id state; run; quit; /* Separate the projected data set into a map and an annotate data set. */ data map dot; set allp; if flag=1 then output dot; else output map; run; /* Define the pattern for the map. */ pattern1 v=me c=black r=50; /* Define the title for the map. */ title 'Customer ZIP Code locations'; /* Generate the map and place the symbols at ZIP code locations. */ proc gmap data=map map=map; id state; choro state / anno=dot nolegend; run; quit; ... View more