%*===================== Macro Argument List Start =====================;
%macro NLMeans(
version, /* Any specified text as first argument displays the */
/* current version of the NLMeans macro in the log. */
instore=, /* Specifies the fitted model that was saved using the */
/* STORE statement in the modeling procedure. The OUT= */
/* option in the STORE statement saves the model in a */
/* file known as an item store. This is the preferred */
/* method for providing the required model information. */
/* However, if the modeling procedure does not offer the */
/* STORE statement, then you might be able to use the */
/* inest= and incovb= parameters. */
inest=, /* Specifies the data set of parameter estimates saved */
/* using an ODS OUTPUT statement in the modeling */
/* procedure. The parameter estimates of the model should */
/* be stored in a variable named ESTIMATE in this data */
/* set. An error is issued if the ESTIMATE variable is */
/* not found. */
incovb=, /* Specifies the data set containing the variance- */
/* covariance matrix of model parameters saved using an */
/* ODS OUTPUT statement in the modeling procedure. */
/* Typically, an option such as COVB is required in the */
/* modeling procedure to make this matrix available for */
/* saving. If the saved data set contains numeric */
/* variables other than those containing the covariance */
/* matrix, they should be removed before specifying the */
/* data set in this macro parameter in order to avoid a */
/* compatibility error. The incovb= data set should have */
/* the same number of observations (rows) and variables */
/* (columns) as the number of rows in the inest= data set */
/* in order to be compatible. */
coef=, /* Specifies the data set of estimate coefficients saved */
/* using an ODS OUTPUT statement in the modeling */
/* procedure. The required table must first be produced */
/* by including the E option in the LSMEANS, SLICE, or */
/* ESTIMATE statement in the modeling procedure. The */
/* table name is typically Coef and can be saved in a */
/* data set by specifying: ods output Coef=data-set-name; */
link=, /* Specifies the link function used in the modeling */
/* procedure, typically in the LINK= option in the MODEL */
/* statement. Note that for models fit by PROC PROBIT */
/* which specifies the link with the DIST= option, */
/* specify link=probit or cumprobit when DIST=NORMAL; */
/* link=logit or cumlogit when DIST=LOGISTIC; or link=cll */
/* or cumcll when DIST=EXTREME or GOMPERTZ. */
diff=all, /* To estimate and test all pairwise differences of means,*/
/* specify diff=all. Sequential differences of means */
/* (mean1-mean2, mean2-mean3, ...) are provided when */
/* diff=seq. To obtain all differences with a control */
/* level, specify diff=number, where number is the */
/* position of the control level in the ordered list of */
/* levels shown in the results from the modeling */
/* procedure. For example, if level A is the control in a */
/* variable whose levels are shown in the order A, B, C, */
/* then specify diff=1. When there are multiple sets of */
/* estimates from the modeling procedure (presented in */
/* multiple tables), such as when the SLICE / E SLICEBY= */
/* statement is used or when multiple LSMEANS /E, */
/* SLICE /E, or ESTIMATE /E statements are used, then you */
/* can either specify a single type to apply to all */
/* estimate sets, or a list of types to apply different */
/* types to the sets. If diff= is omitted, all pairwise */
/* differences are done for all sets (diff=all). */
contrasts=, /* Specifies an optional data of coefficients defining */
/* contrasts among the estimates in each set of estimates.*/
/* The specified data set must contain a variable named */
/* SET and k variables named K1, K2, ... , Kk, where k is */
/* the number of Row variables in coef= data set. */
/* Optionally, a LABEL variable can be included to label */
/* each of the specified contrasts. For a given */
/* row in the data set, the value of SET indicates */
/* the set of estimates that the contrast in that row */
/* applies to. The SET value must match one of the values */
/* of the LMatrix variable in the coef= data set. The K */
/* variables correspond to the Row variables in the coef= */
/* data set which define each estimate. For multinomial */
/* models, k=mf where m is the number of estimates */
/* requested in each response function and f is the */
/* number of response functions. For link=glogit models, */
/* k=mf+m. The additional set of m variables correspond */
/* to the last response level. See the description of */
/* options=ratio for limitations on the contrast */
/* coefficients when mean ratios are desired. If */
/* contrasts= and diff= are both specified, diff= is */
/* ignored. */
df=, /* Specifies the degrees of freedom to be used in the */
/* tests and confidence intervals computed for the */
/* estimated functions. If omitted, large-sample Wald */
/* statistics are given. The degrees of freedom for */
/* testing a linear combination of parameters in a linear */
/* model would typically be the number of observations */
/* used in fitting the model minus the number of */
/* parameters estimated in the model - essentially, the */
/* error degrees of freedom. */
alpha=0.05, /* Specifies the alpha level to be used in computing */
/* confidence limits. If omitted, alpha=0.05. */
title=, /* Specifies a title for the table(s) of results. The */
/* title must not contain quotes (" or '), ampersands (&),*/
/* commas, or parentheses. If omitted, */
/* title=Nonlinear Function Estimate. */
options=NOJOINT NONAMES NOREVERSE NORATIO DIFINFNS NOAPPEND /* */
/* Use options= to enable or disable any of several */
/* binary options: JOINT|NOJOINT, NAMES|NONAMES, */
/* REVERSE|NOREVERSE, RATIO|NORATIO, DIFINFNS|DIFALL. */
/* JOINT combines multiple sets of estimates into a */
/* single set before differencing or applying contrasts. */
/* NAMES displays the names of the model parameters used */
/* by the NLEstimate macro. REVERSE reverses the */
/* direction of differencing done by diff=. REVERSE is */
/* ignored when contrasts= is specified. RATIO requests */
/* estimation of ratios rather than differences of pairs */
/* of means. RATIO can be used with diff= or contrasts=, */
/* but with contrasts= each contrast must specify exactly */
/* one 1 coefficient to select the numerator, one -1 */
/* coefficient to select the denominator, and 0 */
/* coefficients otherwise. For multinomial models, DIFALL */
/* applies the differencing type specified in diff= */
/* across all estimates rather than within the response */
/* functions as done by DIFINFNS. This option is ignored */
/* for nonmultinomial models. APPEND merges results from */
/* all estimate sets into a single results data set named */
/* EST_ALL. NOAPPEND saves results sets in separate data */
/* sets (EST1, EST2, ...). The last results set is always */
/* saved in data set EST as well. If not specified, */
/* options=nojoint nonames noreverse noratio difinfns */
/* noappend. */
) / minoperator;
;*=========================== Macro Start =============================;
%local notesopt;
%let notesopt = %sysfunc(getoption(notes));
%let timenlm = %sysfunc(datetime());
%let _version=1.04;
%if &version ne %then %put NOTE: &sysmacroname macro Version &_version..;
%if %index(%upcase(&version),DEBUG) %then %do;
options notes mprint
%if %index(%upcase(&version),DEBUG2) %then mlogic symbolgen;
;
ods select all;
%put _user_;
%let dbgopt=debug,;
%if %index(%upcase(&version),DEBUG2) %then %let dbgopt=debug2,;
%end;
%else %do;
options nonotes nomprint nomlogic nosymbolgen;
ods exclude all;
%if &version ne %then %let dbgopt=v,;
%else %let dbgopt=;
%end;
/* Check for newer version */
%let _notfound=0;
filename _ver url 'http://ftp.sas.com/techsup/download/stat/versions.dat'
termstr=crlf;
data _null_;
infile _ver end=_eof;
input name:$15. ver;
if upcase(name)="&sysmacroname" then do;
call symput("_newver",ver); stop;
end;
if _eof then call symput("_notfound",1);
run;
options notes;
%if &syserr ne 0 or &_notfound=1 %then
%put NOTE: Unable to check for newer version of &sysmacroname macro.;
%else %if %sysevalf(&_newver > &_version) %then %do;
%put NOTE: A newer version of the &sysmacroname macro is available at;
%put NOTE- this location: http://support.sas.com/ ;
%end;
%if %index(%upcase(&version),DEBUG)=0 %then options nonotes;;
/* Check inputs */
%if (&instore= and &inest= and &incovb=) or
(&instore= and (&inest= or &incovb=)) %then %do;
%put ERROR: Either INSTORE= or both INEST= and INCOVB= must be specified.;
%goto exit;
%end;
%if &coef= or %sysfunc(exist(&coef)) ne 1 %then %do;
%put ERROR: COEF= data set not specified or not found.;
%goto exit;
%end;
%if &contrasts ne and %sysfunc(exist(&contrasts)) ne 1 %then %do;
%put ERROR: CONTRASTS= data set not found.;
%goto exit;
%end;
/* Verify COEF= data set has necessary variables */
%let status=ok;
%let dsid=%sysfunc(open(&coef));
%if &dsid %then %do;
%if %sysfunc(varnum(&dsid,LMATRIX))=0 or %sysfunc(varnum(&dsid,ROW1))=0 or
(%sysfunc(varnum(&dsid,PARAMETER))=0 and %sysfunc(varnum(&dsid,EFFECT))=0)
%then %do;
%put ERROR: Invalid COEF= data set. At least one variable Parameter (or;
%put ERROR- Effect), LMatrix, or Row1 not found.;
%let status=input_err;
%end;
%let rc=%sysfunc(close(&dsid));
%if &status=input_err %then %goto exit;
%end;
%else %do;
%put ERROR: Could not open COEF= data set.;
%goto exit;
%end;
%let validopts=JOINT NOJOINT NAMES NONAMES REVERSE NOREVERSE RATIO NORATIO
DIFINFNS DIFALL APPEND NOAPPEND;
%let joint=0; %let names=0; %let rdiff=0; %let ratio=0; %let difall=0;
%let append=0;
%let i=1;
%do %while (%scan(&options,&i) ne %str() );
%let option&i=%upcase(%scan(&options,&i));
%if &&option&i=JOINT %then %let joint=1;
%else %if &&option&i=NAMES %then %let names=1;
%else %if &&option&i=REVERSE %then %let rdiff=1;
%else %if &&option&i=RATIO %then %let ratio=1;
%else %if &&option&i=DIFALL %then %let difall=1;
%else %if &&option&i=APPEND %then %let append=1;
%else %do;
%let chk=%eval(&&option&i in &validopts);
%if not &chk %then %do;
%put ERROR: Valid values of OPTIONS= are &validopts..;
%goto exit;
%end;
%end;
%let i=%eval(&i+1);
%end;
%if %index(%quote(&title),%str(%")) or %index(%quote(&title),%str(%'))
%then %do;
%put ERROR: Do not use quotes (%str(%") or %str(%')) in TITLE=.;
%goto exit;
%end;
%if &df= %then %let dfopt=;
%else %let dfopt=, df=&df;
%if &alpha= %then %let alphaopt=;
%else %let alphaopt=, alpha=α
%if %quote(&title)= %then %let ttlopt=;
%else %let ttlopt=, title=%quote(&title);
%if &names=1 %then %let namopt=;
%else %let namopt=, listnames=no;
%if &rdiff %then %do;
%let locode=-1; %let hicode=1;
%end;
%else %do;
%let locode=1; %let hicode=-1;
%end;
/* Determine inverse link function */
%let link=%quote(%upcase(&link));
%if &link= %then %do;
%put ERROR: The LINK= option is required.;
%goto exit;
%end;
%else %if &link=LOGIT or &link=CUMLOGIT or &link=CLOGIT %then
%let ilink="logistic("||trim(left(exp))||")";
%else %if &link=PROBIT or &link=CUMPROBIT or &link=CPROBIT %then
%let ilink="probnorm("||trim(left(exp))||")";
%else %if &link=CLL or &link=CLOGLOG or &link=CUMCLL or &link=CCLL
or &link=CCLOGLOG or &link=CUMCLOGLOG %then
%let ilink="(1-exp(-exp("||trim(left(exp))||")))";
%else %if &link=LOGLOG or &link=CUMLOGLOG %then
%let ilink="exp(-exp(-("||trim(left(exp))||")))";
%else %if &link=LOG or &link=GLOGIT %then
%let ilink="exp("||trim(left(exp))||")";
%else %if %substr(&link,1,3)=POW %then %do;
%let power=%quote(%substr(&link,6));
%if &power=0 %then %let ilink="exp("||trim(left(exp))||")";
%else %let ilink="exp(log("||trim(left(exp))||")/&power)";
%end;
%else %if %upcase(&link)=IDENTITY %then %do;
options notes;
%put NOTE: This macro is not needed when LINK=IDENTITY.;
%put NOTE- Use results of LSMEANS, SLICE or ESTIMATE.;
options nonotes;
%goto exit;
%end;
%else %do;
%put ERROR: Valid LINK= values are LOG, LOGIT, CUMLOGIT, PROBIT, ;
%put ERROR- CUMPROBIT, CLL, CUMCLL, LOGLOG, CUMLOGLOG, GLOGIT, ;
%put ERROR- or POWERx where x is a number.;
%goto exit;
%end;
%if &difall and &link=GLOGIT %then %do;
%put ERROR: OPTIONS=DIFALL cannot be used with LINK=GLOGIT.;
%goto exit;
%end;
/* Determine number of response functions for multinomial models */
%let nfn=1;
%let cumlinks=CUMLOGIT CLOGIT CUMPROBIT CPROBIT CUMCLL CCLL CCLOGLOG
CUMCLOGLOG CUMLOGLOG GLOGIT;
%let cumchk=%eval(&link in &cumlinks);
%if &cumchk %then %do;
data _null_;
set &coef end=eof;
if index(Parameter,"Intercept") or index(Effect,"Intercept") then nfn+1;
else do;
call symput('nfn',cats(nfn));
stop;
end;
run;
%end;
/* Get number of sets of estimates */
%if &joint %then %do;
%let nsets=1; %let set1=1;
proc transpose data=&coef out=_coef(drop=_name_);
by LMatrix;
var row:;
run;
proc transpose data=_coef out=_coef(drop=_name_) prefix=Row;
var col:;
run;
data _coef; set _coef;
LMatrix=1;
run;
%end;
%else %do;
proc freq data=&coef nlevels;
table LMatrix / out=_setnums noprint;
ods output nlevels=_nlvl;
run;
data _null_;
set _nlvl;
call symput('nsets',cats(nlevels));
run;
data _null_;
set _setnums;
call symput(cats('set',_n_),cats(lmatrix));
run;
data _coef;
set &coef;
run;
%end;
%if &contrasts= %then %do;
%let diff=%upcase(&diff);
%let ndifvals=%sysfunc(countw(&diff));
%if &ndifvals ne 1 and &ndifvals ne &nsets %then %do;
%put ERROR: There are &nsets sets of estimates. Specify a single value;
%put ERROR- in DIFF= to apply to all sets, or one value for each set.;
%put ERROR- Valid values are ALL, SEQ, or a positive integer.;
%goto exit;
%end;
%let i=1;
%do %while (%scan(&diff,&i) ne %str() );
%let diferr=0;
%let setdif&i=%scan(&diff,&i);
%if &&setdif&i ne ALL and &&setdif&i ne SEQ %then %do;
%let diferr=1;
%if %sysevalf(&&setdif&i>0 and &&setdif&i<1e8) %then %do;
%if %sysevalf(%sysfunc(mod(&&setdif&i,1))=0) %then %let diferr=0;
%end;
%if &diferr %then %do;
%put ERROR: Invalid DIFF= value found. Valid values are ALL, SEQ or;
%put ERROR- a positive integer.;
%goto exit;
%end;
%end;
%let i=%eval(&i+1);
%end;
%end;
/* Process each set of estimates in COEF= data set */
%do s=1 %to &nsets;
/* Construct mean expressions for NLEstimate macro */
proc transpose data=_coef
%if &joint=0 %then (where=(LMatrix=&&set&s));
out=_ct(keep=col: where=(col1 ne .));
by lmatrix;
var row:;
run;
data _glognum(keep=fn1) _glogden(keep=sum);
set _ct nobs=nlsm end=last;
array lb (*) col:;
array fn (&nfn) $32767;
length exp sum sum1 $32767;
do i=1 to &nfn; *initialize &nfn fns for each row in _ct;
if lb(1) ne 0 then fn(i)=catx("*",lb(1),cats("B_p",i));
end;
do i=2 to dim(lb); *cols;
do j=1 to &nfn; *fns;
if lb(i) ne 0 then
fn(j)=catx("+",fn(j),catx("*",lb(i),cats("B_p",i+j-1)));
end;
end;
do i=1 to &nfn;
exp=fn(i); fn(i)=&ilink;
if i=1 then sum1=fn(i);
else sum1=catx("+",sum1,fn(i));
end;
sum="(1+"||trim(left(sum1))||")";
call symput("nlsm",cats(nlsm));
output _glognum;
if _n_<=nlsm/&nfn then output _glogden;
%if &link=GLOGIT %then %do;
if last and &nfn>1 then do;
fn1="1";
do i=1 to nlsm/&nfn;
output _glognum;
end;
end;
%end;
run;
data _prexp;
length exp $32767;
%if &link=GLOGIT %then %do;
set %do i=1 %to &nfn+1; _glogden %end;
;
merge _glognum;
exp=catx("/",fn1,sum);
%if &ratio %then exp="("||trim(left(exp))||")";;
%end;
%else %do;
set _glognum;
exp=fn1;
%end;
keep exp;
run;
proc transpose data=_prexp out=_lbt(keep=col:);
var exp;
run;
/* With DIFALL: don't restrict differencing to be within response functions */
%if &difall %then %let nfn=1;
/* Number of estimates per response function in multinomial models */
%let nperfn=%sysevalf(&nlsm/&nfn);
/* Comparisons: pairwise diffs, diff with control, or input contrasts */
%if &link=GLOGIT and &nfn>1 %then %do;
%let nk=%sysevalf((&nfn+1)*&nperfn);
%let badglgcoef=0;
%if %sysevalf(%sysfunc(mod(&nperfn,1)) ne 0) %then %let badglgcoef=1;
%else %do;
data
%do f=1 %to &nfn;
%let first=%eval((&f-1)*&nperfn+1);
%let last=%eval(&f*&nperfn);
_glg&f(keep=row&first - row&last)
%end;
;
set _coef
%if &joint=0 %then (where=(LMatrix=&&set&s));
;
%do i=1 %to &nfn-1;
if mod(_n_,&nfn)=&i then output _glg&i;
%end;
if mod(_n_,&nfn)=0 then output _glg&nfn;
run;
%do f=2 %to &nfn;
proc compare base=_glg1 compare=_glg&f noprint out=_glgcomp outnoequal;
var row1-row&nperfn;
%let first=%eval((&f-1)*&nperfn+1);
%let last=%eval(&f*&nperfn);
with row&first-row&last;
run;
%let uneq=0;
data _null_;
dsid=open("_glgcomp");
nobs=attrn(dsid, "nobs");
if nobs>0 then call symput ("uneq",1);
rc=close(dsid);
run;
%if &uneq %then %let badglgcoef=1;
%end;
%end;
%if &badglgcoef %then %do;
%put ERROR: For LINK=GLOGIT, the requested estimates must be estimated in;
%put ERROR- each of the &nfn logits. If the ESTIMATE statement was used,;
%put ERROR- the CATEGORY=JOINT option is required.;
%put ERROR: The above message is for estimate set &&set&s...;
%goto endset;
%end;
%end;
%else %let nk=&nlsm;
%if &nperfn=1 %then %do;
%put ERROR: Only one estimate - no differencing possible;
%put ERROR: The above message is for estimate set &&set&s...;
%goto endset;
%end;
%let label=0;
%if &contrasts= %then %do;
%if &ndifvals=1 %then %let setdif&s=&setdif1;
%if %sysevalf(&&setdif&s>0 and &&setdif&s<1e8) %then %do;
%if &&setdif&s>&nperfn %then %do;
%put ERROR: The DIFF= value, &&setdif&s, for estimate set &&set&s must be an;
%put ERROR- integer between 1 to the the number of estimates, &nperfn..;
%goto exit;
%end;
%end;
data _cntrs(keep=k:);
array k (&nk) (&nk*0);
do f=1 to &nfn
%if &link=GLOGIT and &nfn>1 %then +1;
;
first=(f-1)*&nperfn+1;
last=f*&nperfn;
%if &&setdif&s=ALL %then %do;
do i=first to last-1;
do j=i+1 to last;
do h=1 to last;
if h=i then k(h)=&locode;
else if h=j then k(h)=&hicode;
else k(h)=0;
end;
output;
end;
end;
%end;
%else %if &&setdif&s=SEQ %then %do;
do i=first to last-1;
do h=1 to last;
if h=i then k(h)=&locode;
else if h=i+1 then k(h)=&hicode;
else k(h)=0;
end;
output;
end;
%end;
%else %do;
do j=first to last;
do h=1 to last;
if h=(first-1)+&&setdif&s then k(h)=&locode;
else if h=j then k(h)=&hicode;
else k(h)=0;
end;
if j ne (first-1)+&&setdif&s then output;
end;
%end;
end;
run;
%let noobs=0;
data _null_;
dsid=open("_cntrs");
nobs=attrn(dsid, "nobs");
if nobs=0 then call symput ("noobs",1);
rc=close(dsid);
run;
%if &noobs %then %do;
%put ERROR: Could not produce differencing coefficients. If a multinomial;
%put ERROR- model was used you might need to specify OPTIONS=DIFALL.;
%put ERROR: The above message is for estimate set &&set&s...;
%goto endset;
%end;
%end;
%else %do;
%let misscoef=0; %let miss_set=0;
data _cntrs;
set &contrasts;
array k (*) k:; nk=dim(k);
if set=&&set&s;
if (nk-nmiss(of k:))<&nk then do;
call symput("misscoef",1);
stop;
end;
run;
data _null_;
dsid=open("_cntrs");
nobs=attrn(dsid, "nobs");
if nobs=0 then call symput ("miss_set",1);
if varnum(dsid,"LABEL") then call symput("label",1);
rc=close(dsid);
run;
%if &misscoef %then %do;
%put ERROR: Too few contrast coefficients found in CONTRASTS= data set;
%put ERROR- for estimate set &&set&s (LMatrix=&&set&s) in COEF= data set.;
%put ERROR- No analysis done for estimate set &&set&s...;
%goto endset;
%end;
%if &miss_set %then %do;
%put ERROR: No contrast coefficients found in CONTRASTS= data found for;
%put ERROR- estimate set &&set&s (LMatrix=&&set&s) in COEF= data set.;
%goto endset;
%end;
%if &ratio %then %do;
%let badcoef=0;
data _null_;
set _cntrs;
badcoef=0; cnt1=0; cntm1=0;
array k (*) k:;
do i=1 to dim(k);
if k(i) not in(-1,0,1,.) then badcoef=1;
if k(i)=1 then cnt1+1;
if k(i)=-1 then cntm1+1;
end;
if cnt1 ne 1 or cntm1 ne 1 then badcoef=1;
if badcoef then call symput("badcoef","1");
run;
%if &badcoef %then %do;
%put ERROR: OPTIONS=RATIO requires all contrasts to contain a single 1;
%put ERROR- for the numerator, a single -1 for the denominator, and;
%put ERROR- zeros otherwise.;
%put ERROR: The above message is for estimate set &&set&s...;
%goto endset;
%end;
%end;
%end;
/* Construct fdata= data set for NLEstimate macro */
data _fd;
if _n_=1 then set _lbt;
set _cntrs;
array lb (*) col:;
array k (*) k:;
length f $32767;
%if &label=0 %then length label $32767;;
do i=1 to dim(lb);
%if &label=0 %then label=catx(" ",label,k(i));;
%if &ratio %then %do;
if k(i)=1 then num=lb(i);
if k(i)=-1 then den=lb(i);
f=catx("/",num,den);
%end;
%else
if k(i) ne 0 then f=catx("+",f,catx("*",k(i),lb(i)));
;
end;
%if &ratio and &label=0 %then %do;
negpos=index(label,"-");
substr(label,negpos,1)="/";
%end;
keep label f;
run;
/* Call NLEstimate macro to estimate/test the contrasts */
%if &instore ne %then
%nlestimate(&dbgopt instore=&instore, fdata=_fd &dfopt &alphaopt
&namopt &ttlopt);
%else %if &inest ne and &incovb ne %then
%nlestimate(&dbgopt inest=&inest, incovb=&incovb, fdata=_fd &dfopt &alphaopt
&namopt &ttlopt);
/* Append results sets if requested or save as separate files */
%if &append %then %do;
%if %sysfunc(exist(est_all)) and &s=1 %then %do;
options notes;
%put NOTE: Data set EST_ALL already exists. Appending to it.;
%if %index(%upcase(&version),DEBUG)=0 %then options nonotes;;
%end;
proc append base=est_all data=est;
run;
%end;
%else %if &nsets > 1 %then %do;
data est&s;
set est;
run;
%end;
/* End of processing set of estimates */
%endset:
%end;
%exit:
%if %index(%upcase(&version),DEBUG) %then %do;
options nomprint nomlogic nosymbolgen;
%put _user_;
%end;
%else %do;
proc datasets nolist nowarn;
delete _coef _ct _setnums _nlvl _glognum _glogden _prexp _cntrs _fd _lbt
_rows _nrows _glg:;
run; quit;
%end;
ods select all;
options ¬esopt;
%let timenlm=%sysfunc(round(%sysevalf(%sysfunc(datetime())-&timenlm),0.01));
%put NOTE: The &sysmacroname macro used &timenlm seconds.;
%mend;
