options ls=72; data one; input chd $ ldl chol age smoke $; ldllimit = 'HighLDL'; if ldl lt 200 then ldllimit = 'LowLDL'; logchol = log(chol); label chd = Coronary Heart Disease Status ldl = Low Density Lippoprotein ldllimit = How or Low LDL chol = Cholesterol logchol = Log(Cholesterol) age = Patient Age smoke = Smoking Status; cards; Healthy 99.00000 245.00000 60.00000 Nonsmoke Healthy 123.00000 193.00000 60.00000 Nonsmoke Healthy 95.00000 178.00000 68.00000 Nonsmoke Healthy 155.00000 213.00000 61.00000 Nonsmoke Healthy 192.00000 255.00000 69.00000 Smoker Healthy 101.00000 161.00000 61.00000 Nonsmoke Healthy 181.00000 248.00000 65.00000 Nonsmoke Healthy 166.00000 255.00000 60.00000 Nonsmoke Healthy 134.00000 225.00000 60.00000 Smoker Healthy 91.00000 180.00000 66.00000 Nonsmoke Healthy 188.00000 248.00000 63.00000 Nonsmoke Healthy 129.00000 189.00000 62.00000 Nonsmoke Healthy 155.00000 197.00000 68.00000 Smoker Healthy 109.00000 166.00000 67.00000 Nonsmoke Healthy 118.00000 193.00000 67.00000 Smoker Healthy 201.00000 268.00000 60.00000 Nonsmoke Healthy 155.00000 221.00000 66.00000 Nonsmoke Healthy 139.00000 189.00000 60.00000 Nonsmoke Healthy 227.00000 290.00000 63.00000 Smoker Healthy 101.00000 145.00000 63.00000 Nonsmoke Healthy 84.00000 193.00000 60.00000 Smoker Healthy 103.00000 162.00000 60.00000 Nonsmoke Healthy 128.00000 178.00000 67.00000 Nonsmoke Healthy 186.00000 242.00000 62.00000 Smoker Healthy 156.00000 239.00000 61.00000 Smoker Healthy 127.00000 181.00000 67.00000 Nonsmoke Healthy 164.00000 248.00000 64.00000 Nonsmoke Healthy 197.00000 281.00000 60.00000 Nonsmoke Healthy 171.00000 226.00000 63.00000 Smoker Healthy 132.00000 237.00000 60.00000 Smoker Healthy 110.00000 173.00000 62.00000 Nonsmoke Healthy 192.00000 284.00000 62.00000 Nonsmoke Healthy 107.00000 155.00000 60.00000 Nonsmoke Healthy 135.00000 243.00000 61.00000 Nonsmoke Healthy 106.00000 204.00000 62.00000 Nonsmoke Healthy 120.00000 169.00000 63.00000 Nonsmoke Healthy 107.00000 169.00000 63.00000 Nonsmoke Healthy 134.00000 196.00000 68.00000 Smoker Healthy 149.00000 227.00000 63.00000 Nonsmoke Healthy 142.00000 202.00000 64.00000 Nonsmoke Healthy 136.00000 200.00000 63.00000 Nonsmoke Healthy 87.00000 180.00000 61.00000 Nonsmoke Healthy 177.00000 239.00000 65.00000 Nonsmoke Healthy 77.00000 220.00000 61.00000 Nonsmoke Healthy 140.00000 225.00000 61.00000 Nonsmoke Healthy 122.00000 175.00000 64.00000 Nonsmoke Healthy 140.00000 196.00000 69.00000 Nonsmoke Healthy 65.00000 130.00000 62.00000 Nonsmoke Healthy 137.00000 208.00000 68.00000 Smoker Healthy 96.00000 166.00000 65.00000 Nonsmoke Healthy 145.00000 205.00000 65.00000 Nonsmoke Healthy 131.00000 215.00000 60.00000 Smoker Healthy 165.00000 229.00000 62.00000 Smoker Healthy 210.00000 286.00000 60.00000 Smoker Healthy 196.00000 258.00000 66.00000 Nonsmoke Healthy 182.00000 259.00000 65.00000 Smoker Healthy 166.00000 234.00000 61.00000 Smoker Healthy 152.00000 216.00000 65.00000 Nonsmoke Healthy 194.00000 285.00000 67.00000 Nonsmoke Healthy 109.00000 217.00000 63.00000 Nonsmoke Healthy 109.00000 170.00000 63.00000 Nonsmoke Healthy 145.00000 214.00000 61.00000 Nonsmoke Healthy 189.00000 256.00000 61.00000 Smoker Healthy 145.00000 190.00000 69.00000 Nonsmoke Healthy 196.00000 260.00000 68.00000 Nonsmoke Healthy 138.00000 207.00000 65.00000 Smoker Healthy 157.00000 226.00000 65.00000 Smoker Healthy 214.00000 280.00000 61.00000 Nonsmoke Healthy 119.00000 195.00000 60.00000 Smoker Healthy 146.00000 198.00000 67.00000 Nonsmoke Healthy 145.00000 227.00000 62.00000 Nonsmoke Healthy 146.00000 210.00000 64.00000 Nonsmoke Healthy 61.00000 174.00000 63.00000 Smoker Healthy 145.00000 199.00000 60.00000 Nonsmoke Healthy 137.00000 184.00000 60.00000 Smoker Healthy 238.00000 304.00000 60.00000 Nonsmoke Healthy 138.00000 182.00000 62.00000 Nonsmoke Healthy 213.00000 288.00000 66.00000 Nonsmoke Healthy 164.00000 238.00000 65.00000 Nonsmoke Healthy 112.00000 174.00000 65.00000 Nonsmoke Healthy 121.00000 181.00000 62.00000 Nonsmoke Healthy 105.00000 185.00000 67.00000 Smoker Healthy 150.00000 239.00000 65.00000 Smoker Healthy 99.00000 169.00000 62.00000 Nonsmoke Healthy 66.00000 112.00000 61.00000 Smoker Healthy 168.00000 261.00000 60.00000 Smoker Healthy 146.00000 231.00000 67.00000 Smoker Healthy 121.00000 171.00000 61.00000 Smoker Healthy 119.00000 178.00000 63.00000 Nonsmoke Healthy 152.00000 274.00000 65.00000 Nonsmoke Healthy 132.00000 201.00000 61.00000 Nonsmoke Healthy 158.00000 228.00000 66.00000 Nonsmoke Healthy 120.00000 172.00000 62.00000 Nonsmoke Healthy 117.00000 206.00000 63.00000 Nonsmoke Healthy 113.00000 188.00000 68.00000 Nonsmoke Healthy 144.00000 213.00000 63.00000 Smoker Healthy 201.00000 246.00000 61.00000 Smoker Healthy 88.00000 199.00000 62.00000 Nonsmoke Healthy 130.00000 190.00000 66.00000 Nonsmoke Healthy 169.00000 237.00000 60.00000 Nonsmoke Healthy 197.00000 293.00000 67.00000 Nonsmoke Healthy 155.00000 238.00000 60.00000 Nonsmoke Healthy 133.00000 211.00000 63.00000 Nonsmoke Healthy 134.00000 195.00000 64.00000 Smoker Healthy 143.00000 215.00000 63.00000 Nonsmoke Healthy 160.00000 219.00000 60.00000 Nonsmoke Healthy 160.00000 230.00000 65.00000 Nonsmoke Healthy 179.00000 253.00000 64.00000 Nonsmoke Healthy 128.00000 200.00000 64.00000 Nonsmoke Healthy 159.00000 175.00000 60.00000 Nonsmoke Healthy 174.00000 245.00000 64.00000 Nonsmoke Healthy 224.00000 319.00000 60.00000 Nonsmoke Healthy 103.00000 168.00000 62.00000 Nonsmoke Healthy 162.00000 255.00000 62.00000 Nonsmoke Healthy 154.00000 226.00000 64.00000 Nonsmoke Healthy 183.00000 266.00000 62.00000 Nonsmoke Healthy 184.00000 275.00000 61.00000 Nonsmoke Healthy 241.00000 323.00000 65.00000 Nonsmoke Healthy 139.00000 194.00000 62.00000 Smoker Healthy 208.00000 273.00000 68.00000 Nonsmoke Healthy 163.00000 230.00000 60.00000 Nonsmoke Healthy 115.00000 196.00000 63.00000 Nonsmoke Healthy 164.00000 235.00000 68.00000 Nonsmoke Healthy 148.00000 208.00000 60.00000 Smoker Healthy 233.00000 324.00000 64.00000 Nonsmoke Healthy 160.00000 238.00000 62.00000 Nonsmoke Healthy 119.00000 200.00000 64.00000 Nonsmoke Healthy 138.00000 195.00000 65.00000 Nonsmoke Healthy 156.00000 262.00000 63.00000 Nonsmoke Healthy 191.00000 260.00000 61.00000 Nonsmoke Healthy 136.00000 224.00000 64.00000 Nonsmoke Healthy 140.00000 225.00000 67.00000 Nonsmoke Healthy 174.00000 232.00000 64.00000 Nonsmoke Healthy 119.00000 192.00000 66.00000 Nonsmoke Healthy 179.00000 262.00000 60.00000 Nonsmoke Healthy 138.00000 214.00000 61.00000 Nonsmoke Healthy 184.00000 268.00000 61.00000 Nonsmoke Healthy 223.00000 299.00000 63.00000 Smoker Healthy 133.00000 225.00000 67.00000 Smoker Healthy 178.00000 271.00000 61.00000 Smoker Healthy 111.00000 187.00000 67.00000 Nonsmoke Healthy 157.00000 231.00000 61.00000 Nonsmoke Healthy 157.00000 239.00000 62.00000 Nonsmoke Healthy 144.00000 230.00000 69.00000 Nonsmoke Healthy 55.00000 184.00000 63.00000 Nonsmoke Healthy 175.00000 257.00000 63.00000 Nonsmoke Healthy 138.00000 219.00000 60.00000 Smoker Healthy 203.00000 288.00000 62.00000 Nonsmoke Healthy 142.00000 220.00000 62.00000 Nonsmoke Healthy 136.00000 234.00000 64.00000 Nonsmoke Healthy 117.00000 173.00000 67.00000 Nonsmoke Healthy 129.00000 181.00000 66.00000 Nonsmoke Healthy 132.00000 190.00000 61.00000 Smoker Healthy 167.00000 250.00000 66.00000 Nonsmoke Healthy 126.00000 196.00000 64.00000 Nonsmoke Healthy 164.00000 225.00000 63.00000 Nonsmoke Healthy 82.00000 170.00000 63.00000 Nonsmoke Healthy 159.00000 240.00000 65.00000 Nonsmoke Healthy 127.00000 198.00000 62.00000 Nonsmoke Healthy 113.00000 174.00000 60.00000 Nonsmoke Healthy 158.00000 225.00000 61.00000 Nonsmoke Healthy 145.00000 213.00000 61.00000 Nonsmoke Healthy 194.00000 237.00000 69.00000 Smoker Healthy 124.00000 198.00000 69.00000 Nonsmoke Healthy 138.00000 231.00000 69.00000 Smoker Healthy 137.00000 185.00000 60.00000 Nonsmoke Healthy 138.00000 196.00000 63.00000 Smoker Healthy 203.00000 282.00000 64.00000 Nonsmoke Healthy 166.00000 227.00000 68.00000 Nonsmoke Healthy 152.00000 232.00000 61.00000 Nonsmoke Healthy 202.00000 286.00000 60.00000 Nonsmoke Healthy 106.00000 161.00000 62.00000 Nonsmoke Healthy 262.00000 333.00000 63.00000 Nonsmoke Healthy 155.00000 230.00000 61.00000 Smoker Healthy 215.00000 318.00000 67.00000 Nonsmoke Healthy 154.00000 224.00000 62.00000 Nonsmoke Healthy 170.00000 239.00000 61.00000 Nonsmoke Healthy 157.00000 250.00000 61.00000 Nonsmoke Healthy 130.00000 188.00000 66.00000 Nonsmoke Healthy 147.00000 220.00000 62.00000 Smoker Healthy 125.00000 191.00000 66.00000 Nonsmoke Healthy 130.00000 190.00000 62.00000 Nonsmoke Healthy 144.00000 233.00000 69.00000 Nonsmoke Healthy 97.00000 160.00000 62.00000 Nonsmoke Healthy 163.00000 241.00000 67.00000 Nonsmoke Healthy 158.00000 237.00000 64.00000 Nonsmoke Healthy 151.00000 234.00000 67.00000 Nonsmoke BadHeart 93.00000 195.00000 63.00000 Nonsmoke BadHeart 174.00000 249.00000 65.00000 Nonsmoke BadHeart 184.00000 275.00000 62.00000 Nonsmoke BadHeart 151.00000 231.00000 61.00000 Smoker BadHeart 226.00000 313.00000 60.00000 Smoker BadHeart 130.00000 198.00000 65.00000 Nonsmoke BadHeart 175.00000 254.00000 64.00000 Nonsmoke BadHeart 106.00000 217.00000 66.00000 Nonsmoke BadHeart 186.00000 252.00000 64.00000 Nonsmoke BadHeart 132.00000 211.00000 68.00000 Nonsmoke BadHeart 79.00000 165.00000 60.00000 Nonsmoke BadHeart 160.00000 255.00000 65.00000 Nonsmoke BadHeart 84.00000 208.00000 64.00000 Nonsmoke BadHeart 125.00000 199.00000 66.00000 Nonsmoke BadHeart 147.00000 211.00000 66.00000 Nonsmoke BadHeart 143.00000 213.00000 68.00000 Nonsmoke BadHeart 131.00000 188.00000 63.00000 Nonsmoke BadHeart 96.00000 167.00000 67.00000 Smoker BadHeart 168.00000 257.00000 68.00000 Nonsmoke BadHeart 166.00000 240.00000 66.00000 Nonsmoke BadHeart 103.00000 140.00000 69.00000 Nonsmoke BadHeart 150.00000 232.00000 66.00000 Smoker BadHeart 123.00000 240.00000 68.00000 Nonsmoke BadHeart 211.00000 254.00000 62.00000 Nonsmoke BadHeart 85.00000 174.00000 65.00000 Nonsmoke BadHeart 137.00000 224.00000 65.00000 Nonsmoke BadHeart 167.00000 250.00000 62.00000 Smoker BadHeart 187.00000 269.00000 62.00000 Nonsmoke BadHeart 222.00000 291.00000 66.00000 Smoker BadHeart 239.00000 327.00000 61.00000 Smoker BadHeart 216.00000 292.00000 63.00000 Nonsmoke BadHeart 97.00000 188.00000 66.00000 Nonsmoke BadHeart 186.00000 257.00000 68.00000 Nonsmoke BadHeart 132.00000 202.00000 66.00000 Nonsmoke BadHeart 131.00000 283.00000 61.00000 Nonsmoke BadHeart 132.00000 186.00000 63.00000 Nonsmoke BadHeart 225.00000 311.00000 68.00000 Nonsmoke BadHeart 179.00000 250.00000 67.00000 Nonsmoke BadHeart 194.00000 255.00000 60.00000 Nonsmoke BadHeart 176.00000 254.00000 61.00000 Nonsmoke BadHeart 164.00000 227.00000 61.00000 Nonsmoke BadHeart 166.00000 245.00000 64.00000 Nonsmoke BadHeart 241.00000 275.00000 64.00000 Nonsmoke BadHeart 141.00000 194.00000 65.00000 Nonsmoke BadHeart 169.00000 238.00000 62.00000 Nonsmoke BadHeart 201.00000 253.00000 60.00000 Smoker BadHeart 144.00000 202.00000 69.00000 Nonsmoke BadHeart 189.00000 276.00000 66.00000 Nonsmoke BadHeart 226.00000 292.00000 64.00000 Smoker BadHeart 225.00000 289.00000 67.00000 Nonsmoke BadHeart 180.00000 251.00000 63.00000 Nonsmoke BadHeart 207.00000 286.00000 67.00000 Smoker BadHeart 208.00000 280.00000 66.00000 Smoker BadHeart 223.00000 306.00000 63.00000 Smoker BadHeart 125.00000 173.00000 60.00000 Nonsmoke BadHeart 149.00000 208.00000 66.00000 Nonsmoke BadHeart 141.00000 198.00000 60.00000 Nonsmoke BadHeart 201.00000 291.00000 66.00000 Nonsmoke BadHeart 117.00000 186.00000 61.00000 Smoker BadHeart 121.00000 198.00000 62.00000 Nonsmoke BadHeart 168.00000 250.00000 62.00000 Nonsmoke BadHeart 176.00000 241.00000 67.00000 Nonsmoke BadHeart 157.00000 210.00000 64.00000 Nonsmoke BadHeart 176.00000 238.00000 66.00000 Nonsmoke BadHeart 139.00000 197.00000 64.00000 Nonsmoke BadHeart 152.00000 218.00000 65.00000 Smoker BadHeart 185.00000 233.00000 67.00000 Nonsmoke BadHeart 232.00000 340.00000 68.00000 Nonsmoke BadHeart 158.00000 224.00000 65.00000 Nonsmoke BadHeart 145.00000 228.00000 64.00000 Nonsmoke BadHeart 159.00000 203.00000 67.00000 Nonsmoke BadHeart 146.00000 211.00000 69.00000 Nonsmoke BadHeart 166.00000 257.00000 65.00000 Nonsmoke BadHeart 170.00000 248.00000 65.00000 Nonsmoke BadHeart 163.00000 220.00000 63.00000 Smoker BadHeart 181.00000 234.00000 63.00000 Nonsmoke BadHeart 251.00000 318.00000 67.00000 Smoker BadHeart 152.00000 237.00000 65.00000 Nonsmoke BadHeart 148.00000 258.00000 67.00000 Nonsmoke BadHeart 114.00000 175.00000 60.00000 Nonsmoke BadHeart 135.00000 192.00000 68.00000 Nonsmoke BadHeart 146.00000 213.00000 66.00000 Nonsmoke BadHeart 211.00000 274.00000 65.00000 Nonsmoke BadHeart 114.00000 162.00000 60.00000 Smoker BadHeart 194.00000 263.00000 61.00000 Nonsmoke BadHeart 142.00000 204.00000 62.00000 Nonsmoke BadHeart 205.00000 274.00000 65.00000 Nonsmoke BadHeart 141.00000 225.00000 60.00000 Nonsmoke BadHeart 86.00000 166.00000 68.00000 Nonsmoke BadHeart 215.00000 292.00000 60.00000 Nonsmoke BadHeart 119.00000 191.00000 62.00000 Nonsmoke BadHeart 97.00000 210.00000 67.00000 Smoker BadHeart 145.00000 231.00000 63.00000 Smoker BadHeart 107.00000 173.00000 60.00000 Nonsmoke BadHeart 141.00000 202.00000 63.00000 Smoker BadHeart 122.00000 199.00000 61.00000 Nonsmoke BadHeart 160.00000 221.00000 64.00000 Smoker BadHeart 122.00000 179.00000 61.00000 Nonsmoke BadHeart 178.00000 247.00000 64.00000 Nonsmoke BadHeart 117.00000 164.00000 60.00000 Nonsmoke BadHeart 173.00000 243.00000 60.00000 Nonsmoke BadHeart 122.00000 200.00000 64.00000 Smoker BadHeart 124.00000 195.00000 63.00000 Smoker BadHeart 205.00000 293.00000 61.00000 Nonsmoke BadHeart 177.00000 249.00000 63.00000 Nonsmoke BadHeart 68.00000 187.00000 64.00000 Nonsmoke BadHeart 148.00000 225.00000 63.00000 Nonsmoke BadHeart 106.00000 151.00000 62.00000 Smoker BadHeart 163.00000 261.00000 62.00000 Smoker BadHeart 171.00000 222.00000 63.00000 Smoker BadHeart 304.00000 321.00000 64.00000 Nonsmoke BadHeart 157.00000 223.00000 60.00000 Nonsmoke BadHeart 183.00000 263.00000 60.00000 Smoker BadHeart 228.00000 280.00000 69.00000 Smoker BadHeart 150.00000 220.00000 62.00000 Nonsmoke BadHeart 128.00000 204.00000 66.00000 Nonsmoke BadHeart 111.00000 273.00000 69.00000 Smoker BadHeart 233.00000 340.00000 63.00000 Nonsmoke BadHeart 235.00000 294.00000 60.00000 Smoker BadHeart 141.00000 198.00000 62.00000 Smoker BadHeart 182.00000 239.00000 60.00000 Smoker BadHeart 113.00000 188.00000 60.00000 Nonsmoke BadHeart 143.00000 211.00000 66.00000 Smoker BadHeart 175.00000 223.00000 63.00000 Smoker BadHeart 129.00000 218.00000 60.00000 Nonsmoke BadHeart 139.00000 207.00000 65.00000 Nonsmoke BadHeart 205.00000 275.00000 65.00000 Nonsmoke BadHeart 122.00000 167.00000 64.00000 Nonsmoke BadHeart 156.00000 234.00000 63.00000 Nonsmoke BadHeart 140.00000 209.00000 67.00000 Nonsmoke BadHeart 146.00000 251.00000 64.00000 Nonsmoke BadHeart 154.00000 210.00000 67.00000 Nonsmoke BadHeart 122.00000 196.00000 61.00000 Nonsmoke BadHeart 139.00000 224.00000 63.00000 Nonsmoke BadHeart 140.00000 201.00000 62.00000 Nonsmoke BadHeart 183.00000 258.00000 62.00000 Nonsmoke BadHeart 153.00000 228.00000 61.00000 Nonsmoke BadHeart 156.00000 212.00000 67.00000 Nonsmoke BadHeart 152.00000 242.00000 61.00000 Smoker BadHeart 197.00000 273.00000 62.00000 Nonsmoke BadHeart 226.00000 278.00000 62.00000 Nonsmoke BadHeart 166.00000 218.00000 67.00000 Nonsmoke BadHeart 98.00000 168.00000 64.00000 Nonsmoke BadHeart 105.00000 178.00000 65.00000 Smoker BadHeart 163.00000 229.00000 64.00000 Nonsmoke BadHeart 120.00000 198.00000 65.00000 Nonsmoke BadHeart 186.00000 234.00000 68.00000 Smoker BadHeart 134.00000 198.00000 61.00000 Nonsmoke BadHeart 130.00000 195.00000 68.00000 Nonsmoke BadHeart 220.00000 306.00000 66.00000 Nonsmoke data two; set one; smokechd = 'smokechd'; if smoke = 'Nonsmoke' and chd = 'BadHeart' then smokechd = 'NoSmchd'; if smoke = 'Smoker' and chd = 'Healthy' then smokechd = 'smokeOK'; if smoke = 'Nonsmoke' and chd = 'Healthy' then smokechd = 'NoSmOK'; proc sort data = two; by smokechd; proc univariate data=two plot; var chol ldl; by smokechd; proc means data=two; var chol ldl; by smokechd; proc glm data = two; class smokechd; model ldl chol logchol = smokechd; means smokechd / lsd tukey scheffe; contrast 'Smk-CHD v NoSmk-OK ' smokechd -1 0 0 1; contrast 'Smk-CHD v NoSmk-ill' smokechd 0 -1 0 1; contrast 'Smk-CHD v Smk-OK ' smokechd 0 0 -1 1; contrast 'Smk-OK v NoSmk-ill' smokechd 0 -1 1 0; contrast 'Smk-OK v NoSmk-OK ' smokechd -1 0 1 0; contrast 'NoSmk-ill v NoSmk-OK' smokechd -1 1 0 0; estimate 'Smk-CHD v NoSmk-OK ' smokechd -1 0 0 1; estimate 'Smk-CHD v NoSmk-ill' smokechd 0 -1 0 1; estimate 'Smk-CHD v Smk-OK ' smokechd 0 0 -1 1; estimate 'Smk-OK v NoSmk-ill' smokechd 0 -1 1 0; estimate 'Smk-OK v NoSmk-OK ' smokechd -1 0 1 0; estimate 'NoSmk-ill v NoSmk-OK' smokechd -1 1 0 0; proc freq data=two; tables chd*smoke chd*ldllimit/ chisq; proc npar1way data = two wilcoxon; class smokechd; var chol ldl; data three; set two; levene = 0; if smokechd = 'smokechd' then levene = abs(ldl - 163); if smokechd = 'NoSmchd' then levene = abs(ldl - 152); if smokechd = 'smokeOK' then levene = abs(ldl - 145); if smokechd = 'NoSmOK' then levene = abs(ldl - 145); label levene = Levene Statistic for Variance; proc means data=three; var levene; by smokechd; proc npar1way data = three wilcoxon; class smokechd; var levene; run;