libname learn "C:\sas\STAT6250\Examples"; /*Identifies library of data*/

Transcription

1 CHAPTER 7 libname learn "C:\sas\STAT6250\Examples"; /*Identifies library of data*/ /*Problem 7.2*/ proc print data=learn.hosp; where Subject eq 5 or Subject eq 100 or Subject eq 150 or Subject eq 200; /*selects data based on criteria provided in the question using "OR"*/ proc print data=learn.hosp; where Subject in (5, 100, 150,200); /*selects data based on criteria provided in the question using "IN"*/ Output: the output is generated twice. The SAS System Obs AdmitDate quarter DOB DischrDate Subject 5 08/03/ /28/ /12/ /05/ /10/ /15/ /31/ /07/ /11/ /02/ /25/ /10/ /*Problem 7.4*/ data Sales_Weight; set learn.sales; select; /*Use select statement to assign Weights based on Region*/ when (Region='North') Weight=1.5; when (Region='South') Weight=1.7; when (Region='West') Weight=2.0; when (Region='East') Weight=2.0; otherwise; end; proc print data=sales_weight; Page 1 of 11

2 Output: The SAS System Obs EmpID Name Region Customer Item Quantity UnitCost TotalSales Weight George Smith North Barco Corporation 144L George Smith South Cost Cutter's George Smith North Minimart Inc. 188S George Smith North Barco Corporation 908X George Smith South Ely Corp. 122L Glenda Johnson East Food Unlimited 188X Glenda Johnson East Shop and Drop 144L George Smith South Cost Cutter's 855W Sharon Lu West Cost Cutter's Sharon Lu West Pet's are Us 100W Jason Nguyen East Roger's Spirits 122L Jason Nguyen South Spirited Spirits 407XX Glenda Johnson North Minimart Inc Glenda Johnson East Barco Corporation George Smith North Minimart Inc. 188S /*Prob 7.6*/ proc print data=learn.sales; where (Region = 'North' and Quantity lt 60) Or Customer like "%Pet's are Us%"; /*Use combination of logical statements to select based on criteria described in problem*/ Output: The SAS System Obs EmpID Name Region Customer Item Quantity UnitCost TotalSales George Smith North Barco Corporation 144L George Smith North Minimart Inc. 188S George Smith North Barco Corporation 908X Page 2 of 11

3 Obs EmpID Name Region Customer Item Quantity UnitCost TotalSales Sharon Lu West Pet's are Us 100W Glenda Johnson North Minimart Inc George Smith North Minimart Inc. 188S CHAPTER 8 Question2 Run the program here to create a temporary SAS data set (MonthSales): data monthsales; input month /* add your line(s) here */ datalines; ; Modify this program so that a new variable, SumSales, representing Sales to date, is added to the data set. Be sure that the missing value for Sales in month 3 does not result in a missing value for SumSales. SAS Code /* Chapter 8 - Problem 2 */ data monthsales; input month /* added lines to the program */ SumSales + sales; /* adding like this ignores the missing values in sales, SumSales is retained and also initialized to 0 */ format Sales dollar8.2 SumSales dollar10.2; datalines; ; proc print data = monthsales noobs; title 'Sales to date'; Page 3 of 11

4 Sales to date month sales SumSales 1 $ $4, $ $9, $9, $ $14, $ $19, $ $25, $ $32, $ $36, $ $41, $ $47, $ $53, $ $61, Question7 Use an iterative DO loop to plot the following equation: y = 3*x2 5*x + 10 Use values of x from 0 to 10, with an increment of.10. Copy the GPLOT statements from Problem 8 or use PROC PLOT to display the resulting equation. SAS Code /* Chapter 8 - Problem 7 */ data eqnplot; do x = 1 to 10 by.1; y = 3*x**2-5*x + 10; output; end; goptions reset=all ftext='arial' htext=1.0 ftitle='arial/bo' htitle=1.5 colors=(black); symbol v=none i=sm; proc gplot data = eqnplot; quit; title 'Plot of the equation y = 3 * x^2-5*x + 10'; plot y * x; Page 4 of 11

5 Plot of the equation y = 3 * x^2-5*x + 10 y Question8 Use an iterative DO loop to plot the following equation: Logit = log(p / (1 p)) Use values of p from 0 to 1 (with a point at every.05). SAS Code /* Chapter 8 - Problem 8 */ data logitplot; do p = 0 to 1 by.05; logit = log(p / ( 1 - p)); output; end; goptions reset=all ftext='arial' htext=1.0 ftitle='arial/bo' htitle=1.5 colors=(black); symbol v=none i=sm; x Page 5 of 11

6 proc gplot data = logitplot; quit; title 'Logit Plot'; plot logit * p; Logit Plot logit p Page 6 of 11

7 Question12 You place money in a fund that returns a compound interest of 4.25% annually. You deposit $1,000 every year. How many years will it take to reach $30,000? Use Do While or Do Until to calculate the number of years. SAS Code /* Chapter 8 - Problem 12 */ data years; retain year 0; /* initializing the year variable as 0 */ retain total 0; /* initializing the total variable as 0 */ deposit = 1000; interest =.0425; end; do until (total gt 30000); total = (deposit + total) * (1 + interest); /* calculating the compound interest */ year = year + 1; output; format total dollar10.2 deposit dollar8.2; proc print data = years noobs; title 'DO Until Statement to calculate the Compound Interest'; var interest deposit total year; DO Until Statement to calculate the Compound Interest interest deposit total year $ $1, $ $2, $ $3, $ $4, $ $5, $ $6, $ $8, $ $9, $ $11, $ $12, $ $14, $ $15, $ $17, $ $19, Page 7 of 11

8 $ $21, $ $23, $ $25, $ $27, $ $29, $ $31, Therefore, 20 years will be required for the initial deposit to reach $30,000 when the compound interest is 4.25% annually and the deposit of $1,000 is made every year. CHAPTER 9 *9-3; options yearcutoff=1910; data year1910_2006; Date mmddyy8.; format Date date9.; datalines; 01/01/11 02/23/05 03/15/15 05/09/06 ; options yearcutoff=1920; /* Good idea to set yearcutoff back to the default after you change it */ title "Listing of YEAR1910_2006"; proc print data=year1910_2006 noobs; *9-4; data ages; set learn.hosp; AgeJan1 = int(yrdif(dob,'01jan2006'd,'actual')); AgeToday = int(yrdif(dob,today(),'actual')); title "First 10 Observations from Ages"; proc print data=ages(obs=10) noobs; *9-5; data freq; set learn.hosp(keep=admitdate); Day = weekday(admitdate); Month = month(admitdate); Year = year(admitdate); proc format; value days 1='Sun' 2='Mon' 3='Tue' 4='Wed' 5='Thu' 6='Fri' Page 8 of 11

9 7='Sat'; value months 1='Jan' 2='Feb' 3='Mar' 4='Apr' 5='May' 6='Jun' 7='Jul' 8='Aug' 9='Sep' 10='Oct' 11='Nov' 12='Dec'; title "Frequencies for Hospital Admissions"; proc freq data=freq; tables Day Month Year / nocum nopercent; format Day days. Month months.; *9-6; data days; set learn.medical(keep=visitdate); Day = weekday(visitdate); proc format; value days 1='Sun' 2='Mon' 3='Tue' 4='Wed' 5='Thu' 6='Fri' 7='Sat'; title "Frequencies for Visit Dates"; proc freq data=days; tables Day / nocum nopercent; format Day days.; /* Alternative solution not using date functions */ proc freq data=learn.medical; tables VisitDate / nocum nopercent; format VisitDate monname.; ***Note: the MONNAME format prints a month name given a SAS data value; *9-7; title "Admissions before July 15, 2002"; proc print data=learn.hosp; where AdmitDate le '01Jul2002'd and AdmitDate is not missing; Question8 /* Chapter 9 - Problem 8 */ data date_08; input day month year; date = mdy(month, day, year); /* function mdy to create the desired date */ Page 9 of 11

10 datalines; ; quit; format date mmddyy10.; proc print data = date_08 noobs; title 'Date created from the values for Problem 8'; Date created from the values for Problem 8 day month year date Question /25/ /01/ /21/1946 Repeat Problem 8, except use the following data. If there is a missing value for the day, substitute the 15th of the month SAS Code /* Chapter 9 - Problem 9 */ data date_09; input day month year; if missing(day) then date = mdy(month, 15, year); /* substituting 15 for the missing date */ else date = mdy(month, day, year); format date mmddyy10.; datalines; ; proc print data = date_09 noobs; title 'Date created from the values for Problem 9'; Page 10 of 11

11 Date created from the values for Problem 9 day month year date /25/ /15/ /12/2006 Page 11 of 11