1 The SAS System 13:15 Sunday, May 19, 1996 NOTE: Copyright (c) 1989-1996 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Release 6.11 TS040 Licensed to UNIVERSITY OF WISCONSIN, Site 0002176032. This message is contained in the SAS news file, and is presented upon initialization. Edit the files "news" in the "misc/base" directory to display site-specific news and information in the program log. The command line option "-nonews" will prevent this display. NOTE: AUTOEXEC processing beginning; file is /afs/.cs.wisc.edu/s/sas-6.11/alpha_osf20/sas611/autoexec.sas. NOTE: SAS initialization used: real time 0.81 seconds cpu time 0.36 seconds NOTE: AUTOEXEC processing completed. 1 options ls=79 ps=200; 2 3 data species; 4 infile 'wheat.dat' firstobs=2; 5 input ploidy species $ cross $ access plant tr1 tr2 tr3; 6 NOTE: The infile 'wheat.dat' is: File Name=/afs/cs.wisc.edu/p/stat/src/statistics/pda/data/Wheat/wheat.dat, Owner Name=yandell,Group Name=yandell, Access Permission=r--r--r--, File Size (bytes)=127873 NOTE: 2685 records were read from the infile 'wheat.dat'. The minimum record length was 24. The maximum record length was 53. NOTE: The data set WORK.SPECIES has 2685 observations and 8 variables. NOTE: DATA statement used: real time 2.88 seconds cpu time 0.70 seconds 7 proc mixed data = species; 8 class ploidy species cross access plant; 9 model tr1 = ploidy species(ploidy) / solution; 10 random cross(ploidy*species) access(cross*ploidy*species); 11 estimate 'grand mean' intercept 1 ploidy .25 .25 .25 .25 12 species(ploidy) .25 .05 .05 .05 .05 .05 .25 .0625 .0625 .0625 .0625; 13 lsmeans ploidy species(ploidy); 14 make 'lsmeans' out = lsmmix; 15 16 /* look at random effects only for species sativa */ 17 NOTE: 88 observations are not included because of missing values. NOTE: The data set WORK.LSMMIX has 15 observations and 6 variables. NOTE: At least one W.D format was too small for the number to be printed. The decimal may be shifted by the "BEST" format. NOTE: The PROCEDURE MIXED printed page 1. NOTE: PROCEDURE MIXED used: real time 9.23 seconds cpu time 7.33 seconds 18 data sativa; set species; 19 if ploidy = 6; 20 NOTE: The data set WORK.SATIVA has 555 observations and 8 variables. NOTE: DATA statement used: real time 0.70 seconds cpu time 0.20 seconds 21 proc mixed data = sativa; 22 class cross access plant; 23 model tr1 = ; 24 random cross access(cross) / solution; 25 estimate 'sativa' intercept 1; 26 estimate 'cross ci2525' intercept 1 | cross 1 27 access(cross) .2 .2 .2 .2 .2; 28 estimate 'cross ci3816' intercept 1 | cross 0 1 29 access(cross) 0 0 0 0 0 .2 .2 .2 .2 .2; 30 estimate 'cross pi291990' intercept 1 | cross 0 0 1 31 access(cross) 0 0 0 0 0 0 0 0 0 0 .2 .2 .2 .2 .2; 32 estimate 'cross pi291991' intercept 1 | cross 0 0 0 1 33 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 .2 .2 .2 .2 .2; 34 estimate 'access 1 ci2525' intercept 1 | cross 1 access(cross) 1; 35 estimate 'access 2 ci2525' intercept 1 | cross 1 access(cross) 0 1; 36 estimate 'access 3 ci2525' intercept 1 | cross 1 access(cross) 0 0 1; 37 estimate 'access 4 ci2525' intercept 1 | cross 1 access(cross) 0 0 0 1; 38 estimate 'access 5 ci2525' intercept 1 | cross 1 access(cross) 0 0 0 0 1; 39 estimate 'access 1 ci3816' intercept 1 | cross 0 1 40 access(cross) 0 0 0 0 0 1; 41 estimate 'access 2 ci3816' intercept 1 | cross 0 1 42 access(cross) 0 0 0 0 0 0 1; 43 estimate 'access 3 ci3816' intercept 1 | cross 0 1 44 access(cross) 0 0 0 0 0 0 0 1; 45 estimate 'access 4 ci3816' intercept 1 | cross 0 1 46 access(cross) 0 0 0 0 0 0 0 0 1; 47 estimate 'access 5 ci3816' intercept 1 | cross 0 1 48 access(cross) 0 0 0 0 0 0 0 0 0 1; 49 estimate 'access 1 pi291990' intercept 1 | cross 0 0 1 50 access(cross) 0 0 0 0 0 0 0 0 0 0 1; 51 estimate 'access 2 pi291990' intercept 1 | cross 0 0 1 52 access(cross) 0 0 0 0 0 0 0 0 0 0 0 1; 53 estimate 'access 3 pi291990' intercept 1 | cross 0 0 1 54 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 1; 55 estimate 'access 4 pi291990' intercept 1 | cross 0 0 1 56 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 1; 57 estimate 'access 5 pi291990' intercept 1 | cross 0 0 1 58 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1; 59 estimate 'access 1 pi291991' intercept 1 | cross 0 0 0 1 60 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1; 61 estimate 'access 2 pi291991' intercept 1 | cross 0 0 0 1 62 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1; 63 estimate 'access 3 pi291991' intercept 1 | cross 0 0 0 1 64 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1; 65 estimate 'access 4 pi291991' intercept 1 | cross 0 0 0 1 66 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1; 67 estimate 'access 5 pi291991' intercept 1 | cross 0 0 0 1 68 access(cross) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1; NOTE: 4 observations are not included because of missing values. NOTE: At least one W.D format was too small for the number to be printed. The decimal may be shifted by the "BEST" format. NOTE: The PROCEDURE MIXED printed page 2. NOTE: PROCEDURE MIXED used: real time 1.43 seconds cpu time 0.63 seconds NOTE: The SAS System used: real time 15.68 seconds cpu time 9.45 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414