NOTE: Copyright (c) 1999-2001 by SAS Institute Inc., Cary, NC, USA.
NOTE: SAS (r) Proprietary Software Release 8.2 (TS2M0)
Licensed to US GEOLOGICAL SURVEY, Site 0034757169.
NOTE: This session is executing on the WIN_98 platform.
NOTE: SAS initialization used:
real time 3.68 seconds
1 /**********************************************************************************/
2 /* Name: logexpos.sas */
3 /* Author: Terry Shaffer (terry_shaffer@usgs.gov) */
4 /* SAS Version: 8 */
5 /* Supporting files: aic_mac.sas */
6 /* Date: 4 Dec 2002 */
7 /* */
8 /* This program illustrates 1) fitting of logistic-exposure nest-survival */
9 /* models with PROC GENMOD, 2) computation of AIC model-selection criteria, */
10 /* and 3) computation of model-averaged regression coefficients and unconditional */
11 /* standard errors. The example involves evaluation of eight candidate models. */
12 /* */
13 /**********************************************************************************/
14
15 /* Include the file containing the macros for computing model-selection criteria */
16 /* and model-averaged regression coefficients. The user will need to modify the */
17 /* %Include statement to point to the location of the aic_mac.sas file on their */
18 /* computer. Detailed descriptions of the macros are given in the macro file. */
19 %Include "i:\klett\macros\aic_mac.sas";
339
340 /* Read in one observation for each interval of exposure on each nest. */
341 /* expos = interval length; survive=1 (if the nest survives the interval, */
342 /* survive=0 otherwise; nest_ht is a continuous explanatory variable; */
343 /* parastat and patsize are categorical explanatory variables. */
344 data chat;
345 length patsize $ 5;
346 input expos nest_ht parastat patsize survive;
347 if nest_ht=. then delete;
348 trials=1;
349 cards;
NOTE: The data set WORK.CHAT has 292 observations and 6 variables.
NOTE: DATA statement used:
real time 0.38 seconds
661 ;
662
663 /* following code computes the effective sample size for computing AICc */
664 Data N_eff(keep=n_eff);
665 Set Chat end=lastobs;
666 if survive=0 then n_eff+1;
667 else if survive=1 then n_eff+expos;
668 if lastobs then do;
669 put "The effective sample size for computing AICc is " n_eff ;
670 output;
671 end;
672 run;
The effective sample size for computing AICc is 485
NOTE: There were 292 observations read from the data set WORK.CHAT.
NOTE: The data set WORK.N_EFF has 1 observations and 1 variables.
NOTE: DATA statement used:
real time 0.00 seconds
673
674 /* following code fits a logistic-exposure constant-survival model */
675 /* and creates three data sets containing information about the model */
676 /* and the results. */
677 proc genmod data=chat;
678 a=1/expos;
679 fwdlink link = log((_mean_**a)/(1-_mean_**a));
680 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
681 model survive/trials = / dist=bin;
682 ods output modelfit=modelfit;
683 ods output modelinfo=modelinfo;
684 ods output ParameterEstimates=ParameterEstimates;
685 title 'logistic-exposure, constant-survival model';
686 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 2 observations and 8 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.76 seconds
687
688 /* following code invokes the aicc macro to read the data sets created by GENMOD */
689 /* and compute aic values. AIC results will be stored in the data set named in */
690 /* the dsn= statement. Parameter estimates will be stored in the data set named */
691 /* in the estdsn= statement. model_dimension is set to 2 because the highest */
692 /* order interaction in the suite of candidate models is 2 (parastat*patsize). */
693 %aicc(dsn=chat_aic,estdsn=estimate,model=constant survival,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.10 seconds
NOTE: There were 2 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 1 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 1 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 1 observations and 6 variables.
NOTE: DATA statement used:
real time 0.11 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 1 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 1 observations and 9 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: BASE data set does not exist. DATA file is being copied to BASE file.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: The data set WORK.CHAT_AIC has 1 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.11 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: BASE data set does not exist. DATA file is being copied to BASE file.
NOTE: There were 1 observations read from the data set WORK.TEMP6.
NOTE: The data set WORK.ESTIMATE has 1 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.00 seconds
694
695 /* following code fits a logistic-exposure model with a parastat main effect */
696 proc genmod data=chat;
697 class parastat;
698 a=1/expos;
699 fwdlink link = log((_mean_**a)/(1-_mean_**a));
700 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
701 model survive/trials = parastat / dist=bin;
702 ods output modelfit=modelfit;
703 ods output modelinfo=modelinfo;
704 ods output ParameterEstimates=ParameterEstimates;
705 title 'logistic-exposure, parasitism status main effect';
706 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 4 observations and 9 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.17 seconds
707
708 /* Summarize and store the results from the above model */
709 %aicc(dsn=chat_aic,estdsn=estimate,model=parastat main effect,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 4 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 3 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 3 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 3 observations and 6 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 3 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 3 observations and 9 variables.
NOTE: DATA statement used:
real time 0.06 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: 1 observations added.
NOTE: The data set WORK.CHAT_AIC has 2 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.10 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: There were 3 observations read from the data set WORK.TEMP6.
NOTE: 3 observations added.
NOTE: The data set WORK.ESTIMATE has 4 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.05 seconds
710
711 /* following code fits a logistic-exposure model with a patsize main effect */
712 proc genmod data=chat;
713 class patsize;
714 a=1/expos;
715 fwdlink link = log((_mean_**a)/(1-_mean_**a));
716 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
717 model survive/trials = patsize / dist=bin;
718 ods output modelfit=modelfit;
719 ods output modelinfo=modelinfo;
720 ods output ParameterEstimates=ParameterEstimates;
721 title 'logistic-exposure, patch size main effect';
722 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 4 observations and 9 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.16 seconds
723
724 /* Summarize and store the results from the above model */
725 %aicc(dsn=chat_aic,estdsn=estimate,model=patsize main effect,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 4 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 3 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 3 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 3 observations and 6 variables.
NOTE: DATA statement used:
real time 0.11 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 3 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 3 observations and 9 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: 1 observations added.
NOTE: The data set WORK.CHAT_AIC has 3 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.04 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: There were 3 observations read from the data set WORK.TEMP6.
NOTE: 3 observations added.
NOTE: The data set WORK.ESTIMATE has 7 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.00 seconds
726
727 /* following code fits a logistic-exposure model with a parastat and patsize main effects
727! */
728 proc genmod data=chat;
729 class parastat patsize;
730 a=1/expos;
731 fwdlink link = log((_mean_**a)/(1-_mean_**a));
732 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
733 model survive/trials = parastat patsize / dist=bin;
734 ods output modelfit=modelfit;
735 ods output modelinfo=modelinfo;
736 ods output ParameterEstimates=ParameterEstimates;
737 title 'logistic-exposure, parastat and patsize main effect';
738 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 6 observations and 9 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.17 seconds
739
740 /* Summarize and store the results from the above model */
741 %aicc(dsn=chat_aic,estdsn=estimate,model=parastat and patsize main
741! effect,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 6 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 5 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 5 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 5 observations and 6 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 5 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 5 observations and 9 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: 1 observations added.
NOTE: The data set WORK.CHAT_AIC has 4 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.00 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: There were 5 observations read from the data set WORK.TEMP6.
NOTE: 5 observations added.
NOTE: The data set WORK.ESTIMATE has 12 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.05 seconds
742
743 /* following code fits a logistic-exposure model with a parastat and patsize main effects
743! */
744 /* and their interaction
744! */
745 proc genmod data=chat;
746 class parastat patsize;
747 a=1/expos;
748 fwdlink link = log((_mean_**a)/(1-_mean_**a));
749 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
750 model survive/trials = parastat patsize parastat*patsize / dist=bin;
751 ods output modelfit=modelfit;
752 ods output modelinfo=modelinfo;
753 ods output ParameterEstimates=ParameterEstimates;
754 title 'logistic-exposure, main effects and interaction';
755 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 10 observations and 10 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.16 seconds
756
757 /* Summarize and store the results from the above model */
758 %aicc(dsn=chat_aic,estdsn=estimate,model=main effects and interaction,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 10 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 9 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 9 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 9 observations and 6 variables.
NOTE: DATA statement used:
real time 0.11 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 9 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 9 observations and 9 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: 1 observations added.
NOTE: The data set WORK.CHAT_AIC has 5 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.05 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: There were 9 observations read from the data set WORK.TEMP6.
NOTE: 9 observations added.
NOTE: The data set WORK.ESTIMATE has 21 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.00 seconds
759
760 /* following code fits a logistic-exposure model for nest_ht, a continuous covariate.
760! */
761 proc genmod data=chat;
762 a=1/expos;
763 fwdlink link = log((_mean_**a)/(1-_mean_**a));
764 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
765 model survive/trials = nest_ht / dist=bin;
766 ods output modelfit=modelfit;
767 ods output modelinfo=modelinfo;
768 ods output ParameterEstimates=ParameterEstimates;
769 title 'logistic-exposure, nest_ht';
770 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 3 observations and 8 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.16 seconds
771
772 /* Summarize and store the results from the above model */
773 %aicc(dsn=chat_aic,estdsn=estimate,model=nest height only,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 3 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 2 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 2 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 2 observations and 6 variables.
NOTE: DATA statement used:
real time 0.10 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 2 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 2 observations and 9 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: 1 observations added.
NOTE: The data set WORK.CHAT_AIC has 6 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.05 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: There were 2 observations read from the data set WORK.TEMP6.
NOTE: 2 observations added.
NOTE: The data set WORK.ESTIMATE has 23 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.00 seconds
774
775 /* following code fits a logistic-exposure model with effects of parastat, patsize */
776 /* and nest_ht. */
777 proc genmod data=chat;
778 class parastat patsize;
779 a=1/expos;
780 fwdlink link = log((_mean_**a)/(1-_mean_**a));
781 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
782 model survive/trials = parastat patsize nest_ht / dist=bin;
783 ods output modelfit=modelfit;
784 ods output modelinfo=modelinfo;
785 ods output ParameterEstimates=ParameterEstimates;
786 title 'logistic-exposure, main effects and nest-height';
787 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 7 observations and 9 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.15 seconds
788
789 /* Summarize and store the results from the above model */
790 %aicc(dsn=chat_aic,estdsn=estimate,model=main effects and nest-height,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 7 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 6 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 6 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 6 observations and 6 variables.
NOTE: DATA statement used:
real time 0.10 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 6 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 6 observations and 9 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: 1 observations added.
NOTE: The data set WORK.CHAT_AIC has 7 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.06 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: There were 6 observations read from the data set WORK.TEMP6.
NOTE: 6 observations added.
NOTE: The data set WORK.ESTIMATE has 29 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.00 seconds
791
792 /* following code fits a logistic-exposure model with effects of parastat, patsize */
793 /* parastat*patsize, and nest_ht.
793! */
794 proc genmod data=chat;
795 class parastat patsize;
796 a=1/expos;
797 fwdlink link = log((_mean_**a)/(1-_mean_**a));
798 invlink ilink = (exp(_xbeta_)/(1+exp(_xbeta_)))**expos;
799 model survive/trials = parastat patsize parastat*patsize nest_ht / dist=bin;
800 ods output modelfit=modelfit;
801 ods output modelinfo=modelinfo;
802 ods output ParameterEstimates=ParameterEstimates;
803 title 'logistic-exposure, main effects, interaction, and nest-height';
804 run;
NOTE: Algorithm converged.
NOTE: The scale parameter was held fixed.
NOTE: The data set WORK.PARAMETERESTIMATES has 11 observations and 10 variables.
NOTE: The data set WORK.MODELINFO has 8 observations and 3 variables.
NOTE: The data set WORK.MODELFIT has 5 observations and 4 variables.
NOTE: PROCEDURE GENMOD used:
real time 0.15 seconds
805
806 /* Summarize and store the results from the above model */
807 %aicc(dsn=chat_aic,estdsn=estimate,model=main effects + interaction +
807! nest-height,model_dimension=2);
NOTE: There were 5 observations read from the data set WORK.MODELFIT.
NOTE: The data set WORK.TEMP has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 1 observations read from the data set WORK.N_EFF.
NOTE: There were 8 observations read from the data set WORK.MODELINFO.
NOTE: The data set WORK.TEMP2 has 1 observations and 3 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 11 observations read from the data set WORK.PARAMETERESTIMATES.
NOTE: The data set WORK.TEMP4 has 10 observations and 2 variables.
NOTE: The data set WORK.TEMP4A has 10 observations and 2 variables.
NOTE: The data set WORK.TEMP4B has 10 observations and 6 variables.
NOTE: DATA statement used:
real time 0.11 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: There were 10 observations read from the data set WORK.TEMP4B.
NOTE: The data set WORK.TEMP6 has 10 observations and 9 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: There were 1 observations read from the data set WORK.TEMP.
NOTE: There were 1 observations read from the data set WORK.TEMP2.
NOTE: The data set WORK.TEMP3 has 1 observations and 10 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: Appending WORK.TEMP3 to WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.TEMP3.
NOTE: 1 observations added.
NOTE: The data set WORK.CHAT_AIC has 8 observations and 10 variables.
NOTE: PROCEDURE APPEND used:
real time 0.00 seconds
NOTE: Appending WORK.TEMP6 to WORK.ESTIMATE.
NOTE: There were 10 observations read from the data set WORK.TEMP6.
NOTE: 10 observations added.
NOTE: The data set WORK.ESTIMATE has 39 observations and 9 variables.
NOTE: PROCEDURE APPEND used:
real time 0.05 seconds
808
809 /* Compute AIC model selection criteria and Akaike weights*/
810 %DeltaAic(DataIn=chat_aic,DataOut=chat_delta_aicc,Vari=aicc);
NOTE: The data set WORK.MIN has 1 observations and 1 variables.
NOTE: PROCEDURE UNIVARIATE used:
real time 0.00 seconds
NOTE: There were 8 observations read from the data set WORK.CHAT_AIC.
NOTE: There were 1 observations read from the data set WORK.MIN.
NOTE: The data set WORK.AIC2 has 8 observations and 13 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: There were 8 observations read from the data set WORK.AIC2.
NOTE: The data set WORK.AIC2 has 8 observations and 13 variables.
NOTE: PROCEDURE SORT used:
real time 0.11 seconds
NOTE: The data set WORK.WDENOMINATOR has 1 observations and 1 variables.
NOTE: PROCEDURE UNIVARIATE used:
real time 0.00 seconds
NOTE: There were 1 observations read from the data set WORK.WDENOMINATOR.
NOTE: There were 8 observations read from the data set WORK.AIC2.
NOTE: The data set WORK.CHAT_DELTA_AICC has 8 observations and 15 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: There were 8 observations read from the data set WORK.CHAT_DELTA_AICC.
NOTE: PROCEDURE PRINT used:
real time 0.04 seconds
811
812 /* Compute model-averaged estimates and unconditional standard errors */
813 %Modelavg(estimate=estimate, Akaike=chat_delta_aicc, DataOut=Model_Averaged_Estimates,
813! model_dimension=2);
NOTE: 29 observations with duplicate key values were deleted.
NOTE: There were 39 observations read from the data set WORK.ESTIMATE.
NOTE: The data set WORK.ALLCOMBOS has 10 observations and 3 variables.
NOTE: PROCEDURE SORT used:
real time 0.06 seconds
NOTE: There were 39 observations read from the data set WORK.ESTIMATE.
NOTE: The data set WORK.ESTIMATE has 39 observations and 9 variables.
NOTE: PROCEDURE SORT used:
real time 0.04 seconds
NOTE: There were 39 observations read from the data set WORK.ESTIMATE.
NOTE: The data set WORK.ESTIMATOR has 39 observations and 10 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 39 observations read from the data set WORK.ESTIMATOR.
NOTE: The data set WORK.ESTIMATOR has 39 observations and 10 variables.
NOTE: PROCEDURE SORT used:
real time 0.00 seconds
NOTE: The data set WORK.NMODELS has 1 observations and 1 variables.
NOTE: PROCEDURE UNIVARIATE used:
real time 0.04 seconds
NOTE: There were 1 observations read from the data set WORK.NMODELS.
NOTE: There were 10 observations read from the data set WORK.ALLCOMBOS.
NOTE: The data set WORK.ALLCOMBOS2 has 80 observations and 5 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 80 observations read from the data set WORK.ALLCOMBOS2.
NOTE: The data set WORK.ALLCOMBOS2 has 80 observations and 5 variables.
NOTE: PROCEDURE SORT used:
real time 0.05 seconds
NOTE: There were 80 observations read from the data set WORK.ALLCOMBOS2.
NOTE: There were 39 observations read from the data set WORK.ESTIMATOR.
NOTE: The data set WORK.ESTIMATE2 has 80 observations and 9 variables.
NOTE: DATA statement used:
real time 0.00 seconds
NOTE: There were 8 observations read from the data set WORK.CHAT_DELTA_AICC.
NOTE: The data set WORK.W1 has 8 observations and 2 variables.
NOTE: PROCEDURE SORT used:
real time 0.04 seconds
NOTE: There were 8 observations read from the data set WORK.W1.
NOTE: The data set WORK.W2 has 8 observations and 3 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 8 observations read from the data set WORK.W2.
NOTE: There were 80 observations read from the data set WORK.ESTIMATE2.
NOTE: The data set WORK.ESTIMATE3 has 80 observations and 11 variables.
NOTE: DATA statement used:
real time 0.04 seconds
NOTE: There were 80 observations read from the data set WORK.ESTIMATE3.
NOTE: The data set WORK.MODEL_AVG has 10 observations and 4 variables.
NOTE: PROCEDURE SUMMARY used:
real time 0.05 seconds
NOTE: There were 80 observations read from the data set WORK.ESTIMATE3.
NOTE: The data set WORK.ESTIMATE3 has 80 observations and 11 variables.
NOTE: PROCEDURE SORT used:
real time 0.05 seconds
NOTE: There were 10 observations read from the data set WORK.MODEL_AVG.
NOTE: There were 80 observations read from the data set WORK.ESTIMATE3.
NOTE: The data set WORK.ESTIMATE4 has 80 observations and 13 variables.
NOTE: DATA statement used:
real time 0.10 seconds
NOTE: There were 80 observations read from the data set WORK.ESTIMATE4.
NOTE: The data set WORK.OUTDATA has 10 observations and 6 variables.
NOTE: PROCEDURE SUMMARY used:
real time 0.05 seconds
NOTE: There were 10 observations read from the data set WORK.OUTDATA.
NOTE: The data set WORK.MODEL_AVERAGED_ESTIMATES has 10 observations and 7 variables.
NOTE: DATA statement used:
real time 0.05 seconds
NOTE: There were 10 observations read from the data set WORK.MODEL_AVERAGED_ESTIMATES.
NOTE: PROCEDURE PRINT used:
real time 0.00 seconds