About MikeZdeb

MikeZdeb · ‎10-05-2015

Hi, doesn't give you a matrix but gives you the values in a data set. data have; input id $10. latitude longitude; datalines; 1000000999 -75.0000 50.0000 1000000888 -90.0000 30.0000 1000000777 -75.0000 50.0500 1000000666 -75.0000 50.0400 1000000555 -80.0000 35.0000 1000000444 -75.0300 50.0000 ; data want; format latitude longitude la2 lo2 10.5 distance 15.5; set have nobs=howmany; do j=_n_+1 to howmany; set have (rename=(id=id2 latitude=la2 longitude=lo2)) point=j; distance = geodist(latitude,longitude,la2,lo2,'m'); output; end; run; data set WANT ... ID ID2 and distance in miles 1 1000000999 1000000888 11391.48746 2 1000000999 1000000777 0.89814 3 1000000999 1000000666 0.71851 4 1000000999 1000000555 410.93468 5 1000000999 1000000444 2.08114 6 1000000888 1000000777 11391.48746 7 1000000888 1000000666 11391.48746 8 1000000888 1000000555 11738.41243 9 1000000888 1000000444 11393.56860 10 1000000777 1000000666 0.17963 11 1000000777 1000000555 411.32648 12 1000000777 1000000444 2.26632 13 1000000666 1000000555 411.24805 14 1000000666 1000000444 2.20145 15 1000000555 1000000444 409.06250 If you really want a matrix with those zeroes in the diagonal, first change the loop in the data step as follows ... from ... do j=_n_+1 to howmany; to ... do j=_n_ to howmany; then do the following ... proc transpose data=want out=want (drop=_name_) prefix=d; var distance; id id2; by id notsorted; run; and you get ... Obs id d1000000999 d1000000888 d1000000777 d1000000666 d1000000555 d1000000444 1 1000000999 0.00000 11391.48746 0.89814 0.71851 410.93468 2.08114 2 1000000888 . 0.00000 11391.48746 11391.48746 11738.41243 11393.56860 3 1000000777 . . 0.00000 0.17963 411.32648 2.26632 4 1000000666 . . . 0.00000 411.24805 2.20145 5 1000000555 . . . . 0.00000 409.06250 6 1000000444 . . . . . 0.00000

MikeZdeb · ‎10-05-2015

Hi. It's good to know the variable type before using functions. But if you don't care, you could use a CAT function with that INPUT function and SAS would give the same result whether the variable is CHAR or NUM ... data _null_; format num char 10.2; retain x '12345' y 12345; char = input(cat(x),5.); num = input(cat(y),5.); put x= char= / y= num=; run; the LOG .(no "nasty" messages) ... 191 data _null_; 192 format num char 10.2; 193 retain x '12345' y 12345; 194 char = input(cat(x),5.); 195 num = input(cat(y),5.); 196 put x= char= / y= num=; 197 run; x=12345 char=12345.00 y=12345 num=12345.00

MikeZdeb · ‎10-05-2015

Hi, from your posting, it looked as if the extra stuff was using an INDEX function. I was trying to point out that you can use any SAS character function with numeric variables if you just use the CAT function with the numeric variables.

MikeZdeb · ‎10-05-2015

Hi. Be careful with character variable lengths. The results of using character functions all have default lengths. If you run PROC CONTENTS on data set TEMP, you'll see that the variable W created with the CATS has a length of 200 and that none of the other functions (TRIM or STRIP) change the length of the resulting variables V, Y, and Z. The lengths will still be the length of any variables involved in the concatenation plus the lengths of any character constants. In this example that's 8, not because that's the default lenght of a character variabe, but because it's 6 (length of X) plus 2 (those two $);

MikeZdeb · ‎10-05-2015

Hi, another idea ... data x; input student :$1. grade @@; datalines; A 50 B 60 C 70 D 33 ; data y; set x; grade + 10*(student in ('A' 'B' 'C')); run; data set Y ... Obs student grade 1 A 60 2 B 70 3 C 80 4 D 33

MikeZdeb · ‎10-05-2015

Hi, just curious as to why you would still want to check for variable type when you can use the INDEX function for either NUMERIC or CHARACTER variables if you use the CAT function within the INDEX function ... INDEX(CAT(var), "search string"). Seems as if checking for variable type with VTYPE is extra effort.

MikeZdeb · ‎10-04-2015

Hi. You can search a numeric array with the WHICHN function (no loop needed ... there's a WHICHC for character data). To look for one value ... * note NUMERIC emplid; data have1; input emplid @@; datalines; 1111 2222 3333 ; data have2; input mgrid_1-mgrid_5; datalines; 101 1111 1212 2124 14322 201 2211 2222 3434 5656 301 3322 2124 5434 3333 401 3355 1199 55 77 ; data want (keep=emplid mgrid); retain emplid 1111; array grid(*) mgrid_1-mgrid_5; set have2; found = whichn(emplid, of grid(*)); if found then do; mgrid=grid(found-1); output; end; run; If you want to search for more than one EMPLID ... data want (keep=emplid mgrid); * lookup array (HAVE1) ... places all values in an array; array lookup(3) _temporary_ do j=1 to howmany; set have1 nobs=howmany; lookup(j) = emplid; end; * check grid (HAVE2) for values in lookup (HAVE1); array grid(5) mgrid_1-mgrid_5; do until (last); set have2 end=last; do j=1 to howmany; found = whichn(lookup(j), of grid(*)); if found then do; emplid=lookup(j); mgrid=grid(found-1); output; end; end; end; run; data set WANT ... Obs emplid mgrid 1 1111 101 2 2222 2211 3 3333 5434

MikeZdeb · ‎10-03-2015

Very nice. I think that there are a couple of extra spaces (leading and within) to get rid of so maybe add ... new_unit=prxchange('s/(^|\s)(?=\.)/$1 0/',-1,prxchange('s/([\/])/ $1 /i',-1,prxchange('s/(d?.\d+)(?=\w)/$1 /',-1,unit))); new_unit=compbl(left(new_unit));

MikeZdeb · ‎10-02-2015

Hi, another idea ... data x; input id1 $ id2 $; cards; ds246 jk900 kl78 gh456 jk900 ds246 ko000 rt234 gh456 kl789 nm879 th45 ; data y; set x; y1 = id1; y2 = id2; array y(2) ; call sortc(of y(*)); run; proc sort data=y out=y (drop=y: ) nodupkey; by y: ; run; data set y ... Obs id1 id2 1 ds246 jk900 2 kl78 gh456 3 gh456 kl789 4 ko000 rt234 5 nm879 th45 If changing the order of ID1 and ID2 doesn't matter, simpler code ... data y; set x; array id(2) ; call sortc(of id(*)); run; proc sort data=y nodupkey; by id: ; run;

MikeZdeb · ‎10-02-2015

Hi, another suggestion ... data x; input type course :$5. cm meetday:$2. meettime; datalines; 2 MAT20 215 MW 09 4 MAT20 154 MW 08 4 ACG15 214 T 12 4 ACG15 214 MW 12 ; data y; do new_column=1 by 1 until (last.cm); set x; by course cm notsorted; if first.cm and last.cm then call missing (new_column); output; end; run; data set y ... new_ Obs column type course cm meetday meettime 1 . 2 MAT20 215 MW 9 2 . 4 MAT20 154 MW 8 3 1 4 ACG15 214 T 12 4 2 4 ACG15 214 MW 12

MikeZdeb · ‎10-02-2015

Hi, one idea ... Always use a CAT function (both for character and numeric) and just search for characters. Here's a data step example, should also work in SQL, don't have to check for type ... data x; input x1 x2 :$5.; datalines; 55255 A-AAA 66662 BBB-B 99999 CCCCC ; data y; set x; if index(cat(x1),'2') then found_num = 1; else found_num = 0; if index(cat(x2),'-') then found_char = 1; else found_char = 0; run; No problems in the LOG ... 38 data y; 39 set x; 40 if index(cat(x1),'2') then found_num = 1; else found_num = 0; 41 if index(cat(x2),'-') then found_char = 1; else found_char = 0; 42 run; NOTE: There were 3 observations read from the data set WORK.X. NOTE: The data set WORK.Y has 3 observations and 4 variables. and it works ... found_ found_ Obs x1 x2 num char 1 55255 A-AAA 1 1 2 66662 BBB-B 1 1 3 99999 CCCCC 0 0

MikeZdeb · ‎10-01-2015

Hi, made some assumptions as to how your data looks. Solution also creates the variable as a real date. If you don't care about that and just want the variable name, use the second data step. data x; input id d200101-d200103; datalines; 1 1 . 1 2 . . 1 3 . 1 1 ; data y; format first1 yymon7.; set x; array d(*) d: ; first1 = input(substr(vname(d(whichn( 1,of d: ))),2,6),yymmn6.); run; Obs date id d200101 d200102 d200103 1 2001JAN 1 1 . 1 2 2001MAR 2 . . 1 3 2001FEB 3 . 1 1 Alternative, FIRST1 as a character variable (variable name with the first occurrence of a 1) ... data y; set x; array d(*) d: ; first1 = vname(d(whichn( 1,of d: ))); run;

MikeZdeb · ‎10-01-2015

Hi. There is a way to add a legend that contains just the text you would like to display (map with legend attached) ... NOTE: revised version eliminates data step that separates the projected data set into the MAP and ANNOTATE data sets and uses WHERE statements in subsequent steps to use the correct portion of the projected data set Data Position; Input @1 City $20. @30 y @42 x; x=atan(1)/45 * x; y=atan(1)/45 * y; Nr=_N_; Datalines; New York 40.714352 74.0059731 Ann Arbor 42.2808256 83.7430378 Santiago de Compostela 42.8782132 8.5448445 Dalian 38.914003 -121.614682 ; * create a macro variable with text to be used in a LEGEND statement; proc sql noprint; select catt('j=l "',nr,': ',city,'"') into :cities separated by ' ' from position; quit; * combine positions and world map (no Antarctica); Data combo; Set position maps.world (Drop=x y Rename=(long=x lat=y) where=(cont ne 97)); run; Proc GProject Data=Combo Out=Proj Project=Gall; ID Cont ID; Run; Data Anno; Retain xsys ysys '2' hsys '3' function 'label' style 'calibri' size 2.5 color 'steel' when 'a'; Set proj; where ^missing(city); text=cat(nr); run; goptions reset=all ftext='calibri' htext=2.5 gunit=pct; * use the macro variable in the LEGEND statement; legend1 origin=(3,20)pct mode=share value=none shape=bar(.00001,.00001)pct label=(&cities); Pattern1 v=ms r=77 c=PALG; Title1 'MAP WITH MODIFIED LEGEND' ls=2; Proc GMap Data=proj Map=proj; where missing(city); ID Cont ID; Choro ID / Anno=Anno COutLine=white legend=legend1 statistic=first; Run; Quit;

MikeZdeb · ‎10-01-2015

Hi. Here's an example using GPLOT (admittedly more stuff than SGPLOT). It produced the attached chart (default color scheme used) ... * your data; data x; input class overall pct_over v_222 pct_222 v_233 pct_233 v_244 pct_244; datalines; 4 3 1 0 0 0 0 3 2 5 22 15 60 45 31 15 40 30 7 50 35 58 38 44 33 33 40 8 130 27 77 32 50 40 20 15 9 70 18 18 13 15 11 18 9 10 33 8 13 8 17 8 9 7 13 9 4 6 4 25 20 7 4 14 11 2 2 1 7 3 1 1 16 18 6 0 0 9 8 4 5 17 12 3 0 0 4 2 5 5 19 8 2 0 0 6 3 1 1 20 6 1 0 0 3 1 2 3 ; * rearrange your data for use in PROC GCHART; * three variables in new data set; * CLASS; * _NAME_ (the name of the PCT variable from data set x); * COL1 (values of percentages in data set x); proc transpose data=x out=y; by class; var pct: ; run; * rescale the percentages from 0-100 to 0-1 (so you can use the PERCENT format in GCHART); data y; set y; col1 = col1/100; run; * specify a destination for the PNG file produced by GCHART; filename gout 'z:\chart.png'; * specify details for the chart; goptions reset=all gunit=pct ftext='calibri' htext=2 gsfname=gout dev=png xpixels=1500 ypixels=1000; * override the default axis look in GCHART; axis1 label=none value=none; axis2 label=none offset=(2,2)pct; axis3 label=(h=3 a=90 'ACCOUNT(%)') order=0 to .50 by .05 minor=none; * move the legend inside the chart frame and tweak defaults; legend1 position=(top inside right) across=1 label=none shape=bar(3,4) value=('V_222' 'V_233' 'V_244' 'OVERALL') mode=protect offset=(-5,-5)pct; * adds whitespace at top; title1 ls=2; * adds a centered label for the horizontal access; footnote1 h=3 j=c 'CLASS'; * draw the chart (produces a PNG file); proc gchart data=y; vbar _name_ / subgroup=_name_ sumvar=col1 group=class maxis=axis1 gaxis=axis2 raxis=axis3 space=0 legend=legend1; format col1 percent.; run; quit;

MikeZdeb · ‎10-01-2015

Hi. Not sure what you mean by a Manhattan plot since your exaple does not look like what I understand such a plot to be. If you want a standard Manhattan plot, take a look at ... Creating Manhattan Plots with SAS/Graph http://www.sascommunity.org/wiki/Creating_Manhattan_Plots_with_SAS/Graph First paragraph in this link ... A Manhattan plot is a type of bar chart that is sometimes used in the visualization of genome-wide association data. Some examples can be found by using Google to search for "Manhattan plots" then clicking on "Images" on the upper left-hand side of the search results. You can see Manhattan plots used in publications by clicking on either of the following links: example 1 or example 2. has couple links to Manhattan plots in publications.

Online Status	Offline
Date Last Visited	‎06-28-2018 10:46 AM

Re: Create variables as per separators in string

Re: How to create sub-totals ?

Re: Running a function for each unique value

Re: Aggregating dataset by date, including missing dates

Re: Aggregating dataset by date, including missing dates

Re: Aggregating dataset by date, including missing dates

Re: Aggregating dataset by date, including missing dates

Re: How to write generic libname to multiple databases?

Re: delete/retein observations by frequency (n) of variables

Re: delete/retein observations by frequency (n) of variables

Re: is there simpler way to convert char flag to numeric number

Re: INPUT function requires a character argument

Re: Renaming all variables in a table

Re: Report on missing values of each coloumn

Re: Remove the last character of my string

Re: Proc GMap: Caption and custom legend

Re: Saving PROC GPLOT graph to folder

Re: re: Catx Function Error

How to calculate geodesic distance in SAS

Re: Using Proc Distance to Calculate Geodesic Distance

Re: Converting character format to numeric format

Re: Applying SAS functions on variables based on data types

Re: Trim and left together

Re: Can you help me please immediately

Re: Applying SAS functions on variables based on data types

Re: how to find the value I am looking for from my dataset

Re: re: Insert Spaces

Re: Identifying duplicated of pair of id

Re: if first. and last observations

Re: Applying SAS functions on variables based on data types

Re: Filling variable with names of other variables

Re: Proc GMap: Caption and custom legend

Re: Manhattan-Like chart generation

Re: Manhattan-Like chart generation