data one; input junk pace age sex $; cards; 1 5.38333 28 M 2 5.46667 39 M 3 5.51667 41 M 4 5.61667 42 M 5 5.93333 40 M 6 6.00000 23 M 7 6.01667 40 M 8 6.06667 42 M 9 6.08333 28 M 10 6.18333 18 M 11 6.43333 32 M 12 6.58333 32 M 13 6.61667 52 M 14 6.76667 41 M 15 6.76667 43 M 16 6.78333 50 M 17 6.86667 17 M 18 6.90000 31 M 19 6.91667 17 M 20 6.91667 51 M 21 6.96667 33 M 22 6.98333 41 M 23 7.00000 32 M 24 7.01667 55 M 25 7.10000 23 F 26 7.13333 41 M 27 7.15000 34 M 28 7.18333 38 M 29 7.23333 56 M 30 7.25000 47 M 31 7.25000 16 M 32 7.26667 26 F 33 7.28333 14 M 34 7.28333 50 M 35 7.35000 36 M 36 7.35000 37 M 37 7.41667 50 M 38 7.43333 24 M 39 7.46667 40 M 40 7.48333 42 M 41 7.53333 36 M 42 7.56667 10 F 43 7.58333 9 M 44 7.60000 44 M 45 7.60000 32 M 46 7.70000 30 M 47 7.70000 39 M 48 7.73333 39 M 49 7.76667 54 M 50 7.76667 27 F 51 7.78333 15 M 52 7.80000 15 F 53 7.80000 38 M 54 7.85000 48 M 55 7.86667 37 M 56 7.86667 10 M 57 7.90000 48 M 58 7.98333 18 F 59 8.06667 29 M 60 8.08333 26 M 61 8.08333 25 F 62 8.16667 9 M 63 8.16667 32 M 64 8.20000 14 F 65 8.28333 50 M 66 8.28333 32 M 67 8.36667 12 M 68 8.36667 42 M 69 8.36667 28 F 70 8.41667 15 F 71 8.43333 51 M 72 8.48333 43 F 73 8.48333 11 F 74 8.51667 24 F 75 8.51667 39 M 76 8.58333 50 M 77 8.65000 60 M 78 8.73333 30 F 79 8.80000 49 M 80 8.80000 46 M 81 8.85000 47 M 82 8.85000 40 F 83 8.85000 42 F 84 8.90000 23 F 85 8.91667 26 F 86 9.15000 39 M 87 9.18333 44 F 88 9.20000 40 M 89 9.21667 42 M 90 9.21667 26 F 91 9.31667 43 F 92 9.31667 15 F 93 9.33333 15 F 94 9.38333 13 M 95 9.38333 41 F 96 9.40000 42 F 97 9.40000 50 F 98 9.41667 47 M 99 9.48333 49 M 100 9.50000 61 M 101 9.51667 54 M 102 9.55000 50 M 103 9.55000 50 F 104 9.56667 14 F 105 9.56667 43 F 106 9.58333 41 F 107 9.60000 41 F 108 9.65000 61 M 109 9.70000 44 M 110 9.70000 27 F 111 9.70000 36 F 112 9.71667 35 M 113 9.8333 24 F 114 9.9000 49 M 115 9.9833 34 M 116 10.1667 52 M 117 10.1833 15 M 118 10.2333 35 M 119 10.2667 15 F 120 10.2833 38 M 121 10.3667 38 M 122 10.4000 49 F 123 10.6167 31 F 124 10.6667 34 F 125 10.7000 9 F 126 10.7333 9 F 127 10.8000 8 M 128 10.8000 41 M 129 10.8000 40 F 130 10.8667 38 F 131 10.9167 11 M 132 10.9167 49 F 133 10.9333 10 F 134 10.9333 44 M 135 11.0000 50 F 136 11.1500 61 M 137 11.1667 13 M 138 11.5167 35 F 139 11.6000 10 F 140 11.6333 30 F 141 11.6500 28 F 142 11.6833 37 F 143 11.6833 38 F 144 11.7167 49 F 145 11.8667 51 F 146 11.9000 37 F 147 11.9333 41 F 148 11.9333 35 F 149 11.9333 76 F 150 12.0333 34 F 151 12.2000 52 M 152 13.0000 51 F 153 13.0000 46 F 154 13.0000 17 F 155 14.4667 39 M 156 14.4667 6 F 157 15.1000 52 F 158 17.2667 10 F 159 17.2667 10 F 160 17.5000 81 M ; run; proc nlin data=one; parms alpha=1 theta=30 gamma=5.5; model pace = alpha * (age-theta)**2 + gamma; output out=two p=p; run; proc gplot; plot pace*age p*age/overlay; run;