cancel
Turn on suggestions
Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.
Showing results for 
Search instead for 
Did you mean: 
SuzanneDorinski
‎05-10-2023
  • 175 Posts
  • 410 Likes Given
  • 60 Solutions
  • 135 Likes Received

Re: Reading in 5GB birth dataset from https://www.cdc.gov/nchs/data_ac...

by Lapis Lazuli | Level 10 in SAS Programming
‎02-18-2020 09:50 PM
1 Like
‎02-18-2020 09:50 PM
1 Like
I'll share part of the program in this post, and attach the whole thing.  I borrowed some code from the SAS Dummy blog, to figure out the name of the file within the zipped file.   I ran that, and then copied the really long file name into the INFILE statement in the data step.     * borrowing code from https://blogs.sas.com/content/sasdummy/2014/01/29/using-filename-zip/ ; %let ziploc = /folders/myfolders/NCHS Vital Statistics/Nat2018us.zip; /* Assign a fileref wth the ZIP method */ filename inzip zip "&ziploc"; /* Read the "members" (files) from the ZIP file */ data contents(keep=memname); length memname $200; fid=dopen("inzip"); if fid=0 then stop; memcount=dnum(fid); do i=1 to memcount; memname=dread(fid,i); output; end; rc=dclose(fid); run; /* create a report of the ZIP contents */ title "Files in the ZIP file"; proc print data=contents noobs N; run; title; *options obs=100 ; options obs=max; *options nocenter ; **------------------------------------------------ ; ** by Jean Roth Thu Oct 12 11:09:27 EDT 2017 ** This program reads the 2017 NCHS Natality Detail Data File ; ** Report errors to jroth@nber.org ; ** This program is distributed under the GNU GPL. ; ** See end of this file and ** http://www.gnu.org/licenses/ for details. ; ** ----------------------------------------------- ; * The following line should contain the directory where the SAS file is to be stored ; *libname library "/folders/myfolders/NCHS Vital Statistics/"; * The following line should contain the complete path and name of the raw data file. On a PC, use backslashes in paths as in C:\ ; *FILENAME datafile pipe "7z e /homes/data/natality/2017/natl2017.zip -so "; * The following line should contain the name of the SAS dataset ; %let dataset = natl2018; DATA &dataset ; INFILE inzip(Nat2018PublicUS.c20190509.r20190717.txt) zip truncover LRECL = 20000 ; attrib dob_yy length=4 label="Birth Year"; attrib dob_mm length=3 label="Birth Month 01 January"; attrib dob_tt length=4 label="Time of Birth 0000-2359 Time of Birth"; When I tested the code, the options obs=100 statement was uncommented.  Once I was sure that the code was correct, I commented out that line, then typed in the options obs=max statement.    Since the data set is 2 GB, I decided that I didn't want to have a permanent version, so I commented out the LIBNAME statement.  I commented out the FILENAME statement provided by NBER, because I don't  think it will work on my home computer.   I'm using the FILENAME statement from the SAS blog instead.     I didn't have to change anything else in the NBER program after the INFILE statement in the data step.   ... View more

Re: Reading in 5GB birth dataset from https://www.cdc.gov/nchs/data_ac...

by Lapis Lazuli | Level 10 in SAS Programming
‎02-18-2020 12:27 AM
2 Likes
‎02-18-2020 12:27 AM
2 Likes
The page where you get the birth data files points users to a page on the National Bureau of Economic Research (NBER) website. If you go to that page on the NBER website and scroll down, you'll see a table for the United States birth data and documentation.  Jean Roth at NBER has posted SAS, Stata, and SPSS code to read the ASCII file.  She has also posted the file as a Stata file, a SAS data set. and CSV.   The bad news is that she has not done that for the 2018 file.  However, I think the record layout for the 2017 file is the same as the record layout for the 2018 file.  So, http://data.nber.org/natality/2017/natl2017.sas should help you get started reading the ASCII file into SAS.    I was able to modify the 2017 program to read the zipped version of the 2018 file.  One odd note:  while the PDF for the 2018 file shows the record length as 1330, the record length is really 1345.  The NBER program for 2017 shows 15 variables in columns 1330 to 1345, but those columns are all missing in the 2018 file. ... View more

Re: Pie Chart Slice Colors

by Lapis Lazuli | Level 10 in Graphics Programming
‎02-16-2020 03:38 PM
‎02-16-2020 03:38 PM
I don't use the graph template language, but I was able to use an attribute map data set with several SASHELP data sets.   As the log notes, SGPIE is pre-production for SAS 9.4M6.  I tried to control the color for the "Other" slice when using AgeAtStart for the SASHELP.HEART data set, but couldn't get the syntax to work.   proc format; value AgeGrp 1 - 29 ='Age < 30' 30 - 59 ='Age 30 - 59' 60 - 100 ='Age 60+' ; run; data attrmap; length ID $ 8 FillColor $ 9 value $ 11; infile datalines dlm='|'; input id $ fillcolor $ value $; datalines; myagegrp|green|Age < 30 myagegrp|red|Age 30 - 59 myagegrp|blue|Age 60+ ; run; ods listing gpath='/folders/myfolders/ODS Graphics/SGPIE' style=HTMLBlue; ods graphics / noborder imagename='SGPIE_example_' reset=index(1); proc sgpie data=sashelp.bmimen dattrmap=attrmap; pie age / startangle=90 direction=clockwise attrid=myagegrp datalabelattrs=(size=12) dataskin=matte; format age AgeGrp. ; title 'SASHELP.BMIMEN data set'; run; proc sgpie data=sashelp.class dattrmap=attrmap; pie age / startangle=90 direction=clockwise attrid=myagegrp datalabelattrs=(size=12) dataskin=matte; format age AgeGrp. ; title 'SASHELP.CLASS data set'; run; proc sgpie data=sashelp.heart dattrmap=attrmap; pie AgeAtStart / startangle=90 direction=clockwise attrid=myagegrp datalabelattrs=(size=12) dataskin=matte; format AgeAtStart AgeGrp. ; title 'SASHELP.HEART data set'; title2 'AgeAtStart'; run; proc sgpie data=sashelp.heart dattrmap=attrmap; pie AgeAtDeath / startangle=90 direction=clockwise attrid=myagegrp datalabelattrs=(size=12) dataskin=matte; format AgeAtDeath AgeGrp. ; title 'SASHELP.HEART data set'; title2 'AgeAtDeath'; run; proc sgpie data=sashelp.heart dattrmap=attrmap; pie AgeCHDdiag / startangle=90 direction=clockwise attrid=myagegrp datalabelattrs=(size=12) dataskin=matte; format AgeCHDdiag AgeGrp. ; title 'SASHELP.HEART data set'; title2 'AgeCHDdiag'; run; These are the pie charts: Pie chart of age in SASHELP.BMIMEN Pie chart of age in SASHELP.CLASS Pie chart of AgeAtStart in SASHELP.HEART Pie chart of AgeAtDeath in SASHELP.HEART Pie chart of AgeCHDdiag in SASHELP.HEART ... View more

Re: SAS University Ed JupyterLab issue - Unable to import and use pyto...

by Lapis Lazuli | Level 10 in Administration and Deployment
‎12-28-2019 06:57 PM
‎12-28-2019 06:57 PM
See the section "Python limitations in SAS University Edition" towards the bottom of Coding in Python with SAS University Edition. ... View more

Re: Naming SGE files

by Lapis Lazuli | Level 10 in Graphics Programming
‎08-30-2019 12:10 PM
‎08-30-2019 12:10 PM
You can use the IMAGENAME option on the ODS GRAPHICS statement to specify the base name.  SAS then adds an underscore followed by the name of the destination.   I used the first PROC LIFETEST in the online documentation example 2 as example code below.     The code below produces kmplot_RTF.sge.     ods graphics on / imagename="kmplot" reset=index; ODS RTF file="/folders/myfolders/ODS RTF/kmplot.rtf" sge=on; ODS select SurvivalPlot; proc lifetest data=sashelp.BMT plots=survival(atrisk=0 to 2500 by 500); time T * Status(0); strata Group / test=logrank adjust=sidak; run; ODS RTF close;   ... View more

Re: Importing tab-delimited file failed

by Lapis Lazuli | Level 10 in SAS Studio
‎08-25-2019 07:31 PM
‎08-25-2019 07:31 PM
The text file that you posted has some spaces in addition to tabs.     Your code using DATALINES also has a mix of tabs and spaces.  The 75 value for Y in the first record of DATALINES starts in column 29, when the input statement is looking for the Y value in columns 18 through 21.   I used Notepad++ to see that the text file and the DATALINES version have a mixture of tabs and spaces.   ... View more

Re: Fit table on single slide Proc Report ODS Powerpoint

by Lapis Lazuli | Level 10 in ODS and Base Reporting
‎08-25-2019 07:02 PM
‎08-25-2019 07:02 PM
https://communities.sas.com/t5/ODS-and-Base-Reporting/Changing-ODS-Powerpoint-Slide-Title-Size/td-p/239776 shows how to change the size of the title for PowerPoint.   You can use PROC PRINT to print the first 20 rows, then another PROC PRINT to print the remaining 21 rows.     proc sql outobs=41; create table car_test as select make, model from sashelp.cars; quit; run; ods _all_ close; run; * use style template to change size of title font. see * https://communities.sas.com/t5/ODS-and-Base-Reporting/Changing-ODS-Powerpoint-Slide-Title-Size/td-p/239776 * for more details. ; ods path (prepend) work.templat(update); proc template; define style styles.mystyle; parent = styles.powerpointlight; class systemtitle / fontsize=13pt fontweight=bold textalign=left; end; run; options orientation=landscape /*papersize=letter*/ nodate nonumber topmargin=.025in bottommargin=.001in leftmargin=.5in rightmargin=.5in; ods noproctitle; ods escapechar="^"; title; footnote; ods powerpoint file="/folders/myfolders/ODS POWERPOINT examples/test.pptx" style=mystyle; ***Create first slide; ods powerpoint layout=twocontent; title "Test title"; proc print data=car_test(obs=20) noobs style(header)=[fontweight=bold backgroundcolor=cxd3add9 color=black just=c fontsize=13pt]; var make model / style(column)=[just=c fontsize=13pt verticalalign=middle height=9pt]; run; proc print data=car_test(firstobs=21 obs=41) noobs style(header)=[fontweight=bold backgroundcolor=cxd3add9 color=black just=c fontsize=13pt]; var make model / style(column)=[just=c fontsize=13pt verticalalign=middle height=9pt]; run; ods powerpoint close; ods listing; run; Split 41 row table into 2 pieces to display in PowerPoint ... View more

Re: writing to separate excel sheets in an excel file using separate p...

by Lapis Lazuli | Level 10 in ODS and Base Reporting
‎08-06-2019 09:40 PM
1 Like
‎08-06-2019 09:40 PM
1 Like
Quick example that produces 2 files from the SASHELP.BASEBALL dataset:   %let path=/folders/myfolders/ODS Excel examples/; ods excel file="&path.one_sheet_per_city.xlsx" options(sheet_name="#byval(team)" index='on' embedded_titles='on' frozen_headers='on'); proc sort data=sashelp.baseball out=sorted_by_team; by team; run; options nobyline; proc print data=sorted_by_team noobs; title "#byval(team) players in 1986"; by team; pageby team; run; ods excel close; proc sort data=sorted_by_team out=sorted_by_league_division; by league division name; run; ods excel file="&path.one_sheet_per_league_and_division.xlsx" options(sheet_name="#byval(league) #byval(division)" index='off' embedded_titles='on' frozen_headers='on' start_at='2,2' autofilter='all'); proc print data=sorted_by_league_division noobs; title "Players listed in #byval(league) #byval(division) in 1986"; by league division; pageby division; run; ods excel close; I'd suggest searching on "ods excel Chevell Parker" to get a bunch of papers with examples, and a few videos of his talks.   ... View more

Re: Colon modifier seems not work

by Lapis Lazuli | Level 10 in SAS Studio
‎08-03-2019 11:18 PM
‎08-03-2019 11:18 PM
I copied your program to Notepad++, and can see that the data lines are tab delimited.   This code works for me:   data activity; length age $7 gender $1 race $7 poverty $3 limitation $5 ; infile datalines dlm='09'x dsd missover; input Age Gender Race Poverty Limitation Counts ; datalines; 10–14 M whites low Major 5 10–14 M whites med Major 21 10–14 M whites hig Major 21 10–14 M other low Major 53 10–14 M other med Major 102 10–14 M other hig Major 52 10–14 F whites low Major 1 10–14 F whites med Major 11 10–14 F whites hig Major 15 10–14 F other low Major 21 10–14 F other med Major 52 10–14 F other hig Major 28 15–19 M whites low Major 1 15–19 M whites med Major 13 15–19 M whites hig Major 21 15–19 M other low Major 36 15–19 M other med Major 86 15–19 M other hig Major 40 15–19 F whites low Major 6 15–19 F whites med Major 8 15–19 F whites hig Major 17 15–19 F other low Major 42 15–19 F other med Major 73 15–19 F other hig Major 26 10–14 M whites low Other 2 10–14 M whites med Other 14 10–14 M whites hig Other 6 10–14 M other low Other 16 10–14 M other med Other 36 10–14 M other hig Other 21 10–14 F whites low Other 2 10–14 F whites med Other 2 10–14 F whites hig Other 3 10–14 F other low Other 17 10–14 F other med Other 33 10–14 F other hig Other 12 15–19 M whites low Other 2 15–19 M whites med Other 8 15–19 M whites hig Other 7 15–19 M other low Other 28 15–19 M other med Other 42 15–19 M other hig Other 23 15–19 F whites low Other 3 15–19 F whites med Other 7 15–19 F whites hig Other 2 15–19 F other low Other 31 15–19 F other med Other 36 15–19 F other hig Other 29 10–14 M whites low None 102 10–14 M whites med None 336 10–14 M whites hig None 285 10–14 M other low None 827 10–14 M other med None 1518 10–14 M other hig None 667 10–14 F whites low None 87 10–14 F whites med None 304 10–14 F whites hig None 302 10–14 F other low None 846 10–14 F other med None 1452 10–14 F other hig None 687 15–19 M whites low None 91 15–19 M whites med None 326 15–19 M whites hig None 313 15–19 M other low None 888 15–19 M other med None 1509 15–19 M other hig None 725 15–19 F whites low None 115 15–19 F whites med None 344 15–19 F whites hig None 286 15–19 F other low None 817 15–19 F other med None 1500 15–19 F other hig None 716 ; run; proc print data=activity; run; ... View more

Re: Adding a box around a part of a graphic title

by Lapis Lazuli | Level 10 in Graphics Programming
‎08-03-2019 01:29 PM
‎08-03-2019 01:29 PM
I think you'll need to use an annotation data set to put a box around just a part of a title.   ods graphics / noborder; proc sgplot data=sashelp.class; scatter y=weight x=height / group=sex; title j=l box=1 "A Height and Weight for SASHELP.CLASS"; footnote j=l 'This version uses TITLE statement with BOX option'; run; title; %sganno; data anno; %sgtext(label=" A", border="true", drawspace="layoutpercent", xc1="2", yc1="105", anchor="center", justify="center"); %sgtext(label=" Height and weight for SASHELP.CLASS", border="false", xc1="24", yc1="105", widthunit="percent", width=50, justify="left"); run; * pad= option is adding space at the top of the graphic for * the title from the annotation data set ; proc sgplot data=sashelp.class sganno=anno pad=(top=10pct); scatter y=weight x=height / group=sex; footnote j=l 'This version uses annotation data set for the title.'; run; ods graphics / reset=all; The screen shot below shows the graphic when using an annotation data set for the title. Use an annotation data set to put a box around part of a title in a graphic ... View more

Re: Custom PROC Template with PROC Report

by Lapis Lazuli | Level 10 in SAS Programming
‎08-02-2019 08:53 PM
1 Like
‎08-02-2019 08:53 PM
1 Like
Use an ODS PATH statement when creating custom templates.   libname process1 '/folders/myfolders/templates'; ods path template.process1(update) sashelp.tmplmst(read); proc template; define style MyStyle /store=process1.tmplts; Replace Table from Output / frame=BOX Rules=groups cellpadding=2pt borderspacing=0.75pt borderwidth=0.75pt borderstyle=solid; Replace Body from Document / topmargin=1in bottommargin=1in leftmargin=1in rightmargin=1in; parent=styles.sasdocPrinter; replace fonts / 'TitleFont2'=("Times New Roman", 10.1pt) 'TitleFont'=("Times New Roman", 10.1pt) 'StrongFont'=("Times New Roman", 10.1pt) 'EmphasisFont'=("Times New Roman", 10.1pt) 'FixedEmphasisFont'=("Courier New", 10.1pt) 'FixedStrongFont'=("Courier New", 10.1pt) 'FixedHeadingFont'=("Courier New", 10.1pt) 'BatchFixedFont'=("Courier New", 10.1pt) 'FixedFont'=("Courier New", 10.1pt) 'headingEmphasisFont'=("Times New Roman", 10.1pt) 'headingFont'=("Times New Roman", 10.1pt) 'docFont'=("Times New Roman", 10.1pt); end; run; Then use an ODS PATH statement in your other programs to let SAS know that you have your own style templates.   options nodate; libname process1 '/folders/myfolders/templates'; * MyStyle template has been saved to process1 library in a separate program ; * use ods path statement to make custom templates available ; ods path process1.tmplts(update) sashelp.tmplmst(read); ods pdf style=MyStyle file='/folders/myfolders/templates/Example using custom template.pdf' ; proc print data=sashelp.class; run; ods pdf close; ods path reset; The screen shot below shows the PROC PRINT using the custom style.   PROC PRINT using custom style ... View more

Re: sgmap bubble plot using group option

by Lapis Lazuli | Level 10 in SAS Programming
‎08-02-2019 07:40 PM
1 Like
‎08-02-2019 07:40 PM
1 Like
The GROUP= option for the BUBBLE statement is new in SAS 9.4M6.  Below is a screen shot from the "What's New" documentation.   What's New in Base SAS shows when options are added to ODS Graphics and Mapping procedures It looks like the version tag might be in the wrong spot in the online documentation.  I circled the version tag in red in the screen shot below.  I think it should be before GROUP=variable, not after the first sentence for GROUP=variable.   Version tag might be in the wrong spot in SAS online documentation I'm using SAS 9.4M6 in SAS University Edition.  I downloaded the STATES dataset from http://support.sas.com/rnd/datavisualization/mapsonline/html/archivedmaps.html and tried your code.  The SIZE=variable has to be numeric.  The areacodes variable is character.  I included the areacode variable in the annotations dataset, and then used SIZE=areacode on the BUBBLE statement, with DATALABEL=areacodes as an option.   The image below shows the map.   Map showing area codes for spots in Nevada and New York ... View more

Re: Proc Report ods excel content issue

by Lapis Lazuli | Level 10 in ODS and Base Reporting
‎07-31-2019 10:39 PM
‎07-31-2019 10:39 PM
You might be using an older version of SAS.  I remember seeing that issue with ODS Excel several years ago, perhaps in SAS 9.4 maintenance release 1?   Using the test code you provided, the last column of data in the spreadsheet is the value of x27 in column AB.   The screen shot below shows that there aren't any blank columns between x24 and x25.   SAS 9.4 Maintenance Release 6 produces correct output for ODS Excel destination ... View more

Re: PROC REPORT:Column HEADSKIP

by Lapis Lazuli | Level 10 in SAS Programming
‎07-31-2019 10:12 PM
1 Like
‎07-31-2019 10:12 PM
1 Like
The headline and headskip options only affect LISTING output.   Instead of using rules=none for the style of the report, have you tried rules=groups?   ods pdf file='/folders/myfolders/ODS PDF/Journal_style.pdf'; %let n1=17; %let n2=49; %let n3=99; * headline and headskip options only affect LISTING output! ; proc report data=sashelp.class nowindows split="*" style(report)={rules=groups frame=void} headline headskip; columns (name sex age height weight); define name /order style={cellwidth=31% asis=on just=l} "Investigator Name, and Investigator Number"; define sex /display style={cellwidth=22% just=c} spacing=3 order order=data flow "Placebo*(n=&n2)" left; define age /display style={cellwidth=22% just=c} spacing=3 order order=data flow "CMP-135 * (n=&n1)" left; define height /display style={cellwidth=11% just=c} spacing=3 order order=data flow "All Patients * (n=&n3)" left; define weight /display style={cellwidth=11% just=c} spacing=3 order order=data flow "All Patients * (n=&n3)" left; /* break after invnid/skip; */ compute after _page_/style={just=l font_face='Courier New' font_size=10pt bordertopcolor=black bordertopwidth=0.25}; line "Note:"; line "Investigator name & ID number used."; endcomp; title1 j=l 'Company/PRJ5457C'; title2 j=c 'Table 14.1/1'; title3 j=c 'Enrollment by Investigator'; title4 j=c 'Randomized Patients'; run; ods pdf close; Below is a screen shot of the output.  I used SASHELP.CLASS since you didn't provide example data. PROC REPORT with rules=groups has a line beneath the header and at the bottom of the table   You may find this PROC REPORT Tutorial by @Cynthia_sas helpful. ... View more

Re: Problem with formatting with data forest plot

by Lapis Lazuli | Level 10 in Graphics Programming
‎07-20-2019 02:58 PM
2 Likes
‎07-20-2019 02:58 PM
2 Likes
Combining what @ballardw, @Reeza, and @twildone have already suggested, use this for the data step:   data forest; input Id Subgroup $3-27 Mean Low High PValue $ mean_CI $18.; infile datalines missover; zero=0; one=1; pval_lbl='P value'; meanCI='Mean change'; ObsId=_n_; datalines; 1 GroupA (n=50)............. 2 ..Univariate.............5.3 3.1 7.1 <0.01 5.3 (3.1-7.1) 1 GroupB (n=90)............. 2 ..Univariate.............3.1 1.6 5.2 0.01 3.1 (1.6-5.2) 1 GroupC (n=35)............. 2 ..Univariate.............0.2 -1.1 2.0 0.98 0.2 (-1.1-2.0) ; run; And then you don't have to use Photoshop on the resulting graphic.   Forest plot showing both mean and confidence interval ... View more
Contact Me
Online Status
Offline
Date Last Visited
‎05-10-2023 04:54 PM
Likes given to
Likes from
Quick links