Use of LSMEANS PROC CORR contrast iml covariance auto correlation (http://nue.okstate.edu/Research_Methods/covariance.html) 1. compute glm for gsndvi versus yield (by loc yr) 2. run rep trt anova by loc year for yield 3. significance of trt? 4. identify an outlier in Perkins, 1998 (treatment for yield) data one;, input rep 1 trt 2-3 loc year gsndvi gsinsey yield; cards; 101 Perkins_N*P 1998 0.54804946 0.005708849 1.26600 102 Perkins_N*P 1998 0.532290864 0.005544697 1.19200 103 Perkins_N*P 1998 0.521282328 0.005430024 1.21600 104 Perkins_N*P 1998 0.569138276 0.005928524 1.56400 105 Perkins_N*P 1998 0.72662996 0.007569062 1.36500 106 Perkins_N*P 1998 0.691646312 0.007204649 2.01100 107 Perkins_N*P 1998 0.774145064 0.008064011 2.73100 108 Perkins_N*P 1998 0.757762796 0.007893362 1.88700 109 Perkins_N*P 1998 0.825133876 0.008595145 2.60700 110 Perkins_N*P 1998 0.795067084 0.008281949 2.65600 111 Perkins_N*P 1998 0.78000468 0.008125049 2.18500 112 Perkins_N*P 1998 0.789461288 0.008223555 2.18500 201 Perkins_N*P 1998 0.548470076 0.00571323 0.99300 202 Perkins_N*P 1998 0.513399404 0.00534791 1.21600 203 Perkins_N*P 1998 0.742809172 0.007737596 1.51400 204 Perkins_N*P 1998 0.545634544 0.005683693 1.14200 205 Perkins_N*P 1998 0.64729308 0.006742636 1.66300 206 Perkins_N*P 1998 0.781694396 0.00814265 2.20900 207 Perkins_N*P 1998 0.628350856 0.006545321 1.53900 208 Perkins_N*P 1998 0.848775396 0.00884141 2.43300 209 Perkins_N*P 1998 0.7003052 0.007294846 1.93600 210 Perkins_N*P 1998 0.745318364 0.007763733 2.08500 211 Perkins_N*P 1998 0.819578844 0.00853728 2.75600 212 Perkins_N*P 1998 0.746079824 0.007771665 2.28400 301 Perkins_N*P 1998 0.539477596 0.005619558 0.94300 302 Perkins_N*P 1998 0.53071718 0.005528304 0.72000 303 Perkins_N*P 1998 0.650389684 0.006774893 1.26600 304 Perkins_N*P 1998 0.534248904 0.005565093 1.16700 305 Perkins_N*P 1998 0.65773596 0.006851416 1.66300 306 Perkins_N*P 1998 0.745688216 0.007767586 1.63800 307 Perkins_N*P 1998 0.671362468 0.006993359 1.34100 308 Perkins_N*P 1998 0.692712356 0.007215754 1.88700 309 Perkins_N*P 1998 0.803689712 0.008371768 2.20900 310 Perkins_N*P 1998 0.602954352 0.006280775 0.99300 311 Perkins_N*P 1998 0.8234079 0.008577166 1.96100 312 Perkins_N*P 1998 0.844852064 0.008800542 2.65600 101 Perkins_S*N 1998 0.568507352 0.005921952 1.04119 102 Perkins_S*N 1998 0.621468708 0.006473632 1.89645 103 Perkins_S*N 1998 0.672036904 0.007000384 1.15275 104 Perkins_S*N 1998 0.585730852 0.006101363 1.63616 105 Perkins_S*N 1998 0.596289764 0.006211352 1.63616 106 Perkins_S*N 1998 0.570740968 0.005945218 0.96682 107 Perkins_S*N 1998 0.561545432 0.005849432 1.22712 108 Perkins_S*N 1998 0.504697004 0.00525726 0.66934 109 Perkins_S*N 1998 0.542494428 0.005650984 1.71053 110 Perkins_S*N 1998 0.469539308 0.004891034 1.11556 111 Perkins_S*N 1998 0.642499508 0.006692703 1.00400 112 Perkins_S*N 1998 0.695569644 0.007245517 1.22712 113 Perkins_S*N 1998 0.548818172 0.005716856 1.22712 114 Perkins_S*N 1998 0.542059308 0.005646451 1.00400 115 Perkins_S*N 1998 0.642064388 0.006688171 1.37586 116 Perkins_S*N 1998 0.610264368 0.006356921 0.96682 201 Perkins_S*N 1998 0.536003888 0.005583374 1.52460 202 Perkins_S*N 1998 0.765181592 0.007970642 1.71053 203 Perkins_S*N 1998 0.65668442 0.006840463 1.52460 204 Perkins_S*N 1998 0.774979044 0.008072698 2.30549 205 Perkins_S*N 1998 0.583062116 0.006073564 1.00400 206 Perkins_S*N 1998 0.591728256 0.006163836 1.71053 207 Perkins_S*N 1998 0.780396288 0.008129128 2.15675 208 Perkins_S*N 1998 0.537947424 0.005603619 1.63616 209 Perkins_S*N 1998 0.45859604 0.004777042 1.37586 210 Perkins_S*N 1998 0.747022584 0.007781485 1.59897 211 Perkins_S*N 1998 0.500715656 0.005215788 1.00400 212 Perkins_S*N 1998 0.667475396 0.006952869 1.15275 213 Perkins_S*N 1998 0.583845332 0.006081722 1.37586 214 Perkins_S*N 1998 0.494087328 0.005146743 1.41304 215 Perkins_S*N 1998 0.51091922 0.005322075 1.26430 216 Perkins_S*N 1998 0.56657832 0.005901858 1.07838 301 Perkins_S*N 1998 0.5355035 0.005578161 1.07838 302 Perkins_S*N 1998 0.688905056 0.007176094 1.89645 303 Perkins_S*N 1998 0.623383236 0.006493575 1.82208 304 Perkins_S*N 1998 0.642412484 0.006691797 1.82208 305 Perkins_S*N 1998 0.53651878 0.005588737 1.30149 306 Perkins_S*N 1998 0.673987692 0.007020705 1.22712 307 Perkins_S*N 1998 0.729545264 0.00759943 1.63616 308 Perkins_S*N 1998 0.76042428 0.007921086 309 Perkins_S*N 1998 0.601308148 0.006263627 1.45023 310 Perkins_S*N 1998 0.486828076 0.005071126 1.15275 311 Perkins_S*N 1998 0.682458028 0.007108938 1.52460 312 Perkins_S*N 1998 0.651513744 0.006786602 1.67334 313 Perkins_S*N 1998 0.498162952 0.005189197 1.04119 314 Perkins_S*N 1998 0.497633556 0.005183683 1.00400 315 Perkins_S*N 1998 0.6422892 0.006690513 1.71053 316 Perkins_S*N 1998 0.6422892 0.006690513 2.19394 101 Tipton_S*N 1998 0.74711686 0.007115399 3.19794 102 Tipton_S*N 1998 0.819107464 0.007801023 4.12758 103 Tipton_S*N 1998 0.844467708 0.00804255 4.20195 104 Tipton_S*N 1998 0.852430404 0.008118385 5.09440 105 Tipton_S*N 1998 0.70839118 0.006746583 2.34268 106 Tipton_S*N 1998 0.812740208 0.007740383 3.27231 107 Tipton_S*N 1998 0.855294944 0.008145666 3.42105 108 Tipton_S*N 1998 0.871430644 0.008299339 5.28032 109 Tipton_S*N 1998 0.777916104 0.007408725 2.82609 110 Tipton_S*N 1998 0.81637346 0.007774985 3.75572 111 Tipton_S*N 1998 0.843778768 0.008035988 3.86728 112 Tipton_S*N 1998 0.868435568 0.008270815 5.50343 113 Tipton_S*N 1998 0.708441944 0.006747066 3.08639 114 Tipton_S*N 1998 0.766958332 0.007304365 3.64417 115 Tipton_S*N 1998 0.794827768 0.007569788 4.42506 116 Tipton_S*N 1998 0.851096036 0.008105677 4.87128 201 Tipton_S*N 1998 0.798772856 0.007607361 3.08639 202 Tipton_S*N 1998 0.806307684 0.007679121 4.01602 203 Tipton_S*N 1998 0.81263868 0.007739416 3.75572 204 Tipton_S*N 1998 0.875665812 0.008339674 5.65218 205 Tipton_S*N 1998 0.68514852 0.006525224 2.71453 206 Tipton_S*N 1998 0.840116508 0.00800111 3.30950 207 Tipton_S*N 1998 0.813088304 0.007743698 3.12357 208 Tipton_S*N 1998 0.881438404 0.008394651 4.61099 209 Tipton_S*N 1998 0.705113276 0.006715365 3.16076 210 Tipton_S*N 1998 0.838230988 0.007983152 3.34668 211 Tipton_S*N 1998 0.79976638 0.007616823 3.83009 212 Tipton_S*N 1998 0.86920428 0.008278136 5.50343 213 Tipton_S*N 1998 0.665908964 0.00634199 3.53261 214 Tipton_S*N 1998 0.786429952 0.007489809 3.38387 215 Tipton_S*N 1998 0.802181296 0.007639822 4.61099 216 Tipton_S*N 1998 0.839195504 0.007992338 5.13158 301 Tipton_S*N 1998 0.75774104 0.007216581 3.56980 302 Tipton_S*N 1998 0.769046908 0.007324256 3.75572 303 Tipton_S*N 1998 0.84124782 0.008011884 4.53661 304 Tipton_S*N 1998 0.887254508 0.008450043 5.65218 305 Tipton_S*N 1998 0.755543684 0.007195654 2.67735 306 Tipton_S*N 1998 0.757327676 0.007212645 3.27231 307 Tipton_S*N 1998 0.86452674 0.008233588 3.83009 308 Tipton_S*N 1998 0.88602892 0.008438371 5.09440 309 Tipton_S*N 1998 0.740140436 0.007048957 3.08639 310 Tipton_S*N 1998 0.790273512 0.007526414 3.53261 311 Tipton_S*N 1998 0.791448336 0.007537603 4.53661 312 Tipton_S*N 1998 0.882569716 0.008405426 5.46625 313 Tipton_S*N 1998 0.764311352 0.007279156 3.08639 314 Tipton_S*N 1998 0.739349968 0.007041428 4.64817 315 Tipton_S*N 1998 0.828730868 0.007892675 4.01602 316 Tipton_S*N 1998 0.851552912 0.008110028 5.50343 103 Perkins_N*P 1999 0.4997874 0.0043841 0.80068 106 Perkins_N*P 1999 0.723743664 0.006348629 2.27878 109 Perkins_N*P 1999 0.777234416 0.006817846 2.95018 112 Perkins_N*P 1999 0.67630108 0.005932466 2.18225 203 Perkins_N*P 1999 0.58546978 0.0051357 1.48759 206 Perkins_N*P 1999 0.63474712 0.005567957 1.70827 209 Perkins_N*P 1999 0.575048656 0.005044286 2.35217 212 Perkins_N*P 1999 0.666663172 0.005847923 2.11363 303 Perkins_N*P 1999 0.535605028 0.00469829 0.94332 306 Perkins_N*P 1999 0.59399088 0.005210446 1.38210 309 Perkins_N*P 1999 0.699007092 0.006131641 2.40270 312 Perkins_N*P 1999 0.65004884 0.005702183 2.58665 101 Stillwater_222 1999 0.689710028 0.006696214 1.33614 102 Stillwater_222 1999 0.7581399 0.007360582 1.52734 103 Stillwater_222 1999 0.698586476 0.006782393 1.57693 104 Stillwater_222 1999 0.319712988 0.00310401 110 Stillwater_222 1999 0.550573156 0.00534537 0.82960 201 Stillwater_222 1999 0.570704708 0.005540822 0.76743 202 Stillwater_222 1999 0.75172188 0.007298271 1.20995 203 Stillwater_222 1999 0.664712384 0.006453518 1.53420 204 Stillwater_222 1999 0.670107872 0.006505902 2.21685 210 Stillwater_222 1999 0.543357416 0.005275315 0.84265 301 Stillwater_222 1999 0.616160244 0.005982138 0.79309 302 Stillwater_222 1999 0.674118228 0.006544837 1.44035 303 Stillwater_222 1999 0.69435856 0.006741345 1.96333 304 Stillwater_222 1999 0.594549284 0.005772323 0.82923 310 Stillwater_222 1999 0.494464432 0.004800626 0.70623 401 Stillwater_222 1999 0.610242612 0.005924686 0.80548 402 Stillwater_222 1999 0.650933584 0.006319744 1.20187 403 Stillwater_222 1999 0.76615336 0.007438382 1.82096 404 Stillwater_222 1999 0.609488404 0.005917363 2.04443 410 Stillwater_222 1999 0.48952582 0.004752678 0.76830 102 Efaw_301 1999 0.462055244 0.003756547 0.88776 108 Efaw_301 1999 0.635457816 0.005166324 1.09017 109 Efaw_301 1999 0.766204124 0.006229302 1.35902 110 Efaw_301 1999 0.83468476 0.006786055 2.62536 111 Efaw_301 1999 0.83530118 0.006791067 2.32269 114 Efaw_301 1999 0.828868656 0.00673877 2.64657 202 Efaw_301 1999 0.563111864 0.004578145 1.11416 208 Efaw_301 1999 0.687374884 0.005588414 1.23548 209 Efaw_301 1999 0.712329016 0.005791293 1.55491 210 Efaw_301 1999 0.803733224 0.006534416 2.94043 211 Efaw_301 1999 0.807627548 0.006566078 2.94440 214 Efaw_301 1999 0.826772828 0.00672173 3.96662 302 Efaw_301 1999 0.393422316 0.003198555 0.81442 308 Efaw_301 1999 0.644443044 0.005239374 0.88954 309 Efaw_301 1999 0.728341432 0.005921475 1.97478 310 Efaw_301 1999 0.736246112 0.005985741 2.41972 311 Efaw_301 1999 0.56588938 0.004600727 0.98890 314 Efaw_301 1999 0.682951164 0.005552448 2.59189 101 Efaw_AA 1999 0.656851216 0.006702563 2.18774 105 Efaw_AA 1999 0.795647244 0.008118849 3.67204 106 Efaw_AA 1999 0.770562576 0.007862883 3.82269 107 Efaw_AA 1999 0.790273512 0.008064015 3.82473 108 Efaw_AA 1999 0.780940188 0.007968777 3.71793 109 Efaw_AA 1999 0.80444392 0.008208611 3.86486 110 Efaw_AA 1999 0.71270612 0.007272511 3.32468 201 Efaw_AA 1999 0.747160372 0.007624085 1.98488 205 Efaw_AA 1999 0.728899836 0.007437753 2.97733 206 Efaw_AA 1999 0.770932428 0.007866657 3.41553 207 Efaw_AA 1999 0.772926728 0.007887007 3.47529 208 Efaw_AA 1999 0.765943052 0.007815745 3.31966 209 Efaw_AA 1999 0.807366476 0.008238433 3.42392 210 Efaw_AA 1999 0.787445232 0.008035155 3.97166 301 Efaw_AA 1999 0.658475664 0.006719139 2.33457 305 Efaw_AA 1999 0.735832748 0.007508497 2.73304 306 Efaw_AA 1999 0.747907328 0.007631707 3.22713 307 Efaw_AA 1999 0.700965132 0.007152705 2.67572 308 Efaw_AA 1999 0.736144584 0.007511679 3.07488 309 Efaw_AA 1999 0.761163984 0.007766979 3.32304 310 Efaw_AA 1999 0.780019184 0.007959379 3.82681 101 Haskell_801 1999 0.795139604 0.006571402 1.52413 102 Haskell_801 1999 0.76843774 0.006350725 1.43461 108 Haskell_801 1999 0.866129432 0.007158094 2.47284 109 Haskell_801 1999 0.756965076 0.00625591 0.55233 110 Haskell_801 1999 0.82435066 0.006812815 2.23725 111 Haskell_801 1999 0.807511516 0.006673649 3.14677 112 Haskell_801 1999 0.859515608 0.007103435 1.97058 201 Haskell_801 1999 0.715585164 0.005913927 1.35846 202 Haskell_801 1999 0.78816318 0.006513745 1.71865 208 Haskell_801 1999 0.853474692 0.00705351 3.09049 209 Haskell_801 1999 0.676939256 0.005594539 0.82189 210 Haskell_801 1999 0.837418764 0.006920816 2.41896 211 Haskell_801 1999 0.87511466 0.007232353 2.92302 212 Haskell_801 1999 0.8625687 0.007128667 1.68540 301 Haskell_801 1999 0.728254408 0.006018631 1.70223 302 Haskell_801 1999 0.80788862 0.006676765 2.22456 308 Haskell_801 1999 0.85611442 0.007075326 2.31451 309 Haskell_801 1999 0.484362396 0.004002995 0.41723 310 Haskell_801 1999 0.846766592 0.006998071 2.98003 311 Haskell_801 1999 0.865991644 0.007156956 2.86654 312 Haskell_801 1999 0.859573624 0.007103914 1.83821 401 Haskell_801 1999 0.741808396 0.006130648 1.59411 402 Haskell_801 1999 0.642426988 0.005309314 2.58256 408 Haskell_801 1999 0.861357616 0.007118658 2.52094 409 Haskell_801 1999 0.655785172 0.005419712 0.46283 410 Haskell_801 1999 0.854932344 0.007065557 3.00649 411 Haskell_801 1999 0.861074788 0.007116321 2.88147 412 Haskell_801 1999 0.855323952 0.007068793 0.83485 options linesize = 90; options pagesize = 80; data one; infile 'd:\oficial\osu\lt\502\lah50288.sas'; input rep trt depth ph nh4n no3n minn totn orgn no3 ke actn orgc melp cec; data two; set one; /* ------------------------------------------------------------- depth = 1 = (0-6"), depth = 2 = (6-12"), depth = 3 = (12-18") depth = 4 = (18-24"), depth = 5 = (24-36"), depth = 6 = (36-48") depth = 7 = (48-60"), depth = 8 = (60-72"), depth = 9 = (72-84") depth = 10 = (84-96"), depth = 11 = (96-108"), depth = 12 = (108-120") g/cm3 * 62.43 = lb/ft3 lb/ft3 * 43560ft3/ac = lb soil/ac (0-12") (21780ft3/ac 0-6") lb soil/ac * 453.592g/lb = g soil/ac g soil/ac * ug N/g = ug N/ac ug N/ac * 0.000001g/ug = g N/ac g N/ac * 0.002204623 lb/g = lb N/ac Pb*(62.43*21780*453.592*0.000001*0.002204623)* no3n = lb no3n/ac (cancel 453.592 g/lb with 0.002204623 lb/g) Pb*(62.43*21780*0.000001) * no3n = lb no3n/ac (6" total depth) Pb*(1.3597254)*no3n = lb no3n/ac (6" total depth) Pb*(2.7194508)*no3n = lb no3n/ac (12" total depth) --------------------------------------------------------------- */ if depth = 1 then no3lb = no3n * 1.4 * 1.3597254; if depth = 2 then no3lb = no3n * 1.4 * 1.3597254; if depth = 3 then no3lb = no3n * 1.55 * 1.3597254; if depth = 4 then no3lb = no3n * 1.55 * 1.3597254; if depth = 5 then no3lb = no3n * 1.55 * 2.7194508; if depth = 6 then no3lb = no3n * 1.55 * 2.7194508; if depth = 7 then no3lb = no3n * 1.525 * 2.7194508; if depth = 8 then no3lb = no3n * 1.525 * 2.7194508; if depth = 9 then no3lb = no3n * 1.525 * 2.7194508; if depth = 10 then no3lb = no3n * 1.525 * 2.7194508; if depth = 1 then nh4lb = nh4n * 1.4 * 1.3597254; if depth = 2 then nh4lb = nh4n * 1.4 * 1.3597254; if depth = 3 then nh4lb = nh4n * 1.55 * 1.3597254; if depth = 4 then nh4lb = nh4n * 1.55 * 1.3597254; if depth = 5 then nh4lb = nh4n * 1.55 * 2.7194508; if depth = 6 then nh4lb = nh4n * 1.55 * 2.7194508; if depth = 7 then nh4lb = nh4n * 1.525 * 2.7194508; if depth = 8 then nh4lb = nh4n * 1.525 * 2.7194508; if depth = 9 then nh4lb = nh4n * 1.525 * 2.7194508; if depth = 10 then nh4lb = nh4n * 1.525 * 2.7194508; proc print; data two; set two; if trt = 2 then delete; if trt = 8 then delete; proc glm; classes rep trt depth; model no3n nh4n = rep trt rep*trt depth trt*depth; test h = rep trt e = rep*trt; means depth trt trt*depth; run; data three; set two; proc sort; by rep trt depth; proc transpose data = three out = outn prefix = ndepth; by rep trt; var no3lb; run; data outn; set outn; no3lbt = ndepth1 + ndepth2 + ndepth3 + ndepth4 + ndepth5 + ndepth6 + ndepth7 + ndepth8 + ndepth9 + ndepth10; proc print data = outn (obs=32); proc sort data = outn; by trt; proc means; by trt; proc transpose data = three out = outa prefix = adepth; by rep trt; var nh4lb; run; data outa; set outa; nh4lbt = adepth1 + adepth2 + adepth3 + adepth4 + adepth5 + adepth6 + adepth7 + adepth8 + adepth9 + adepth10; proc print data = outa (obs=32); proc sort data = outa; by trt; proc means; by trt; data new1; set outn; if trt = 2 then delete; if trt = 8 then delete; proc glm data = new1; classes rep trt; model no3lbt = rep trt; means trt; run; data new2; set outa; if trt = 2 then delete; if trt = 8 then delete; proc glm data = new2; classes rep trt; model nh4lbt = rep trt; means trt; run; /* Regression + line (individual estimates, and mean) and confidence limits data one; input red nir NDVI yield height; cards; .5 .6 .7 4000 32 .4 .7 .8 5000 39 .5 .8 .9 5500 41 .4 .5 .7 3500 30 .3 .5 .6 3400 29 .2 .4 .5 2200 25 .24 .45 .55 2400 28 .29 .46 .58 2600 29 proc reg data = one; model yield=height/p clm; output out = dog p = pyd l95m=l95m u95m=u95m; proc gplot; plot pyd*height yield*height='x' l95m*height u95m*height/overlay; run; Program for Determining the Significance Between Slope and Intercept Components from 2 Independent Regressions data one; input exp rep x y; if rep = 1 then intc_dif=0; if rep = 2 then intc_dif=1; slop_dif=intc_dif*x; cards; 2 1 3.3 39 2 1 2.4 48 2 1 3.1 42 2 1 3.35 40 2 2 3 30 2 2 3.1 26 2 2 3.2 20 2 2 3.15 22 data two; set one; proc sort; by rep; proc reg; model y = x intc_dif slop_dif; run; proc reg; by rep; model y=x; run; Surface Response Model linear and quadratic relationships of x and y with z and a linear interaction term. Z = x x2 y y2 xy data one; red nir NDVI yield height; cards; .5 .6 .7 4000 32 .4 .7 .8 5000 39 .5 .8 .9 5500 41 .4 .5 .7 3500 30 .3 .5 .6 3400 29 .2 .4 .5 2200 25 .24 .45 .55 2400 28 .29 .46 .58 2600 29 proc rsreg data = one out = two; model NDVI = red nir /predict; proc g3grid data = two out = three; grid red*nir=NDVI/spline; proc g3d data = three gout=new; plot red*nir=NDVI; run; /* Linear-Plateau Program */ data one; input x y; cards; 76.14909677 63.55625138 75.96631579 65.4660614 77.238 62.26415791 79.48367647 71.42115047 62.68826667 58.17183011 65.84746667 64.03178591 61.45 61.71603593 59.4864 41.77450512 61.43911765 60.84082216 64.93566667 78.62898724 61.66886111 62.61466573 69.6395 53.31725572 67.4186 48.80491021 58.4233 68.06579264 83.23 58.37199282 69.35263158 65.88514318 77.76564706 54.55210173 52.274425 33.41021734 72.732 60 78.36666667 78.60202271 71.68488889 67.10545084 65.1 57.79182613 68.76423529 64.05385342 71.25847826 50.73669537 67.353 50.12004621 70.2999375 57.03644636 64.77823529 38.82080477 83.12588462 . 73.58125 58.13735681 76.48976471 78.65050657 88.662 68.4277427 87.41186667 . 77.7387 54.07236738 76.35694444 51.45835383 59.85 59.95821758 76.21136842 78.32628929 77.62115789 67.82357647 76.415625 70.37309863 62.9145 51.64352563 78.27915789 63.36663219 69.146 51.77236091 75.60282353 69.15255398 80.54402941 72.78673663 58.10171053 43.85799351 76.35269444 81.92505801 77.19766667 54.12075644 83.1105 61.28358395 78.4921875 71.0019481 80.42631579 77.25710421 74.589 48.91214741 proc nlin data = one best = 3; parms b0=10 to 50 by 10 b1=0.5 to 1.1 by .1 njoint=20 to 100 by 10; if x 1.e-10 then scale=pr[j]; if abs(p[j,i]) < 1.e-10 then p[j,i]=0; end; p[,i]=p[,i]/scale; end; print p; run; http://nue.okstate.edu/Research_Methods/RCBD_CRD.html (has interaction contrasts) means trt/duncan lsd alpha level