SAS课程上机练习一

SAS课程上机练习一
SAS课程上机练习一

SAS课程上机练习一

下面的数据是一次对20岁以上人群进行心血管病随访研究的结果,数据中各变量的含义如下:

V1—随访号V2—年龄V3—开始随访时的求诊机构V4/V5—开始随访时的收缩压/舒张压V6/V7—体重/身高V8—开始随访时的胆固醇V9—社会经济地位(1-高2-中3-底)V10—原有心血管病(0-其他心脏病1-冠心病2-冠心与高心3-高心4-高心及风心5-风心6-可疑心脏病7-高血压8-正常)V11—随访结束时的求诊机构V12/V13—随访结束时的收缩压/舒张压V14—随访结束时的胆固醇V15—随访结束时的体重V16—随访结束时的最后诊断(0-未诊断1-3心肌梗塞4-7心绞痛8-9其他)V17—死亡年份(63-68)0-未死。

一、数据集练习操作

1、用data步建立SAS永久数据集(下次练习可以调用此数据集)

2、将年龄(V2)按20—30—40—50—及60岁以上分组

3、计算体重指数={体重/身高2}*100

4、将随访结束时最后诊断的结果(V16)分为四组:(1)心肌梗塞(2) 心绞痛(3)其他心脏

病(4)未诊断

5、按体重指数的大小( 0.30)分为二类:1-超重0-未超重

6、将死亡年份(V17)按是否死亡分为二类(1为死亡,0为未死亡)。

7、将原有心血管病(V10)分为患有心血管病与不患有心血管病二类(1有0无)

8、将数据集按是否死亡分为二个数据集。

9、将社会经济地位(V9)设置为哑变量。

10、计算开始随访时的收缩压与/舒张压(V4/V5)之差。

二.。计算分析

1.年龄分布的特征。

2.开始时收缩压频数分布的特征。

3·超重病人开始时胆固醇(V8)的分布的特征。

3·比较不同体重的人开始时的收缩压。

4·比较不同社会经济地位与患心血管病的差别。

5·比较心肌梗塞的病人开始时的收缩压与随访结束时的收缩压的差别。

6·比较不同年龄、不同体重(指数)间开始时的胆固醇的差别。

7·试分析年龄与血压(V4)之间的关系。

9、比较不同社会经济地位开始随访是收缩压的差别。

三、结果存盘

1.要求将你编辑好的SAS程序保存。

2.要求将所有的计算结果(OUTPUT)保存.。

3.文件名的格式为:你的姓名+你的学号

例如:王永,学号为20918003,则你存盘的文件名为:王永208918003

v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12 v13 v14 v15 v16 v17

1 4

2 1 110 65 64 147 291 2 8 4 120 78 271 146 2 68

2 5

3 1 130 72 69 167 278 1 6 2 122 68 250 165 9 67

3 53 2 120 90 70 222 342

4 8 1 132 90 304 223 2 64

4 48 4 120 80 72 229 239 4 8 2 118 68 209 227 3 65

5 53 3 118 74 6

6 134 243 3 8 5 118 56 261 138 2 66

6 58 2 122 72 69 135 210 3 8 4 130 72 245 136 2 64

7 48 4 130 90 67 165 219 3 8 4 138 86 275 166 2 63

8 60 1 124 80 74 235 203 3 8 1 160 90 271 226 3 65

9 59 4 160 100 72 206 269 5 8 3 150 100 291 198 3 67

10 40 3 120 80 69 148 185 3 8 3 110 64 241 152 2 66

11 56 3 115 80 64 147 260 3 8 4 140 80 326 152 2 68

12 58 3 140 90 63 121 312 5 8 1 120 75 234 114 2 63

13 64 2 135 85 64 189 185 1 8 4 140 78 153 168 3 66

14 57 2 110 78 70 173 282 3 8 2 144 74 236 171 2 66

15 32 1 112 70 69 171 254 2 8 4 142 96 249 179 3 64

16 59 1 140 90 65 150 303 2 8 1 205 85 302 153 2 65

17 48 1 130 80 64 147 271 4 8 3 165 85 251 163 3 64

18 47 2 115 84 67 211 304 1 8 1 155 80 278 149 9 68

19 47 2 130 80 67 147 334 1 8 3 138 85 303 147 2 0

20 28 1 120 86 70 189 328 3 8 2 128 88 300 194 6 0

21 37 3 95 55 69 190 226 3 8 3 155 105 311 191 2 0

22 54 1 141 100 65 171 363 3 7 2 180 100 276 154 2 65

23 38 1 130 90 67 170 399 2 8 2 132 86 353 167 2 0

24 52 2 125 90 65 141 199 2 8 2 152 100 234 135 6 0

25 46 1 110 70 67 159 271 3 8 3 152 88 299 164 6 65

26 51 4 120 80 70 139 261 3 8 1 130 95 285 173 6 0

27 49 1 120 80 68 194 263 3 8 1 178 76 230 196 9 0

28 46 4 110 70 66 160 242 3 8 3 130 90 254 175 6 0

29 26 1 110 80 70 206 260 3 8 3 130 76 325 231 2 0

30 35 2 120 80 72 191 321 3 8 1 130 80 334 169 6 0

31 45 2 108 80 70 155 258 5 8 4 138 88 259 182 2 0

32 57 1 130 80 69 184 167 3 8 3 155 90 237 173 2 0

33 24 3 104 75 70 157 185 5 8 2 120 80 236 166 7 0

34 64 1 144 95 66 191 244 1 8 2 198 110 227 187 7 0

35 34 3 142 102 71 176 314 1 7 1 145 100 233 176 9 0

36 30 3 110 80 71 198 234 2 8 1 100 65 227 187 2 0

37 52 4 145 90 66 183 289 3 7 3 150 88 299 143 6 0

38 56 2 125 75 65 122 329 5 8 1 140 80 253 130 2 66

39 44 1 125 90 65 156 439 1 8 1 130 80 342 152 2 0

40 45 3 130 90 73 143 243 4 8 4 158 78 249 146 6 0

41 29 2 140 95 64 148 419 5 7 1 130 85 254 139 5 0

43 46 3 134 90 73 198 271 3 8 1 125 80 219 205 4 0

44 45 4 150 104 71 187 278 1 7 1 210 110 368 189 2 0

45 51 3 120 90 66 163 226 3 8 1 130 80 271 167 1 0

46 44 2 120 90 72 211 188 3 8 2 138 88 240 196 2 68

47 32 2 108 78 66 151 235 4 8 2 120 70 226 155 6 0

48 34 2 130 90 61 120 317 2 5 2 122 68 248 118 6 0

49 30 2 120 84 68 170 258 3 8 1 165 110 359 181 4 0

50 40 1 112 80 69 167 334 3 8 2 100 70 306 154 6 0

51 34 3 124 88 66 195 345 3 8 1 130 85 309 183 4 0

52 43 3 118 72 71 149 224 3 8 2 120 70 209 147 2 0

53 39 2 164 110 66 245 220 4 7 1 145 95 258 194 2 0

54 45 2 110 80 69 170 347 3 8 3 150 90 296 186 2 0

55 41 2 115 80 68 145 339 5 8 1 160 80 254 149 1 0

56 57 1 130 90 68 188 353 3 8 1 160 80 230 172 7 0

57 53 3 110 80 67 150 235 3 8 4 120 76 221 170 9 0

58 55 2 125 90 65 163 235 2 8 2 170 94 255 139 6 0

59 57 2 210 110 67 165 220 1 4 2 150 94 178 174 6 0

60 38 3 115 90 70 187 385 3 8 2 142 108 334 202 9 0

61 45 4 110 80 67 209 240 3 8 2 156 108 251 205 1 0

62 33 1 130 90 68 200 188 3 8 4 125 70 210 191 7 0

63 61 2 160 100 68 160 241 2 7 3 170 110 235 163 6 0

64 36 3 100 75 72 164 241 3 8 4 125 72 271 175 2 0

65 37 3 130 88 67 178 295 4 8 3 170 90 367 194 1 0

66 51 3 125 85 73 198 283 4 8 4 175 85 309 211 9 0

67 46 1 110 80 69 178 277 3 8 4 150 88 311 186 2 0

68 51 3 138 100 72 208 296 3 7 1 135 95 286 199 8 0

69 60 3 130 84 69 122 243 2 8 4 166 102 291 113 0 65

70 57 2 110 80 71 224 158 1 8 3 135 80 202 224 0 68

71 63 4 130 80 67 143 243 3 8 2 166 90 273 121 0 64

72 63 3 115 65 69 196 278 2 8 1 150 65 248 192 0 66

73 68 3 120 80 63 109 215 1 8 2 136 76 251 108 0 68

74 57 3 145 85 66 140 308 1 8 3 150 80 247 134 0 63

75 64 1 150 90 70 147 226 2 8 4 140 85 192 145 0 63

76 63 2 115 75 67 180 303 3 8 2 122 80 289 126 0 65

77 62 1 120 80 68 174 535 2 8 5 146 76 268 156 0 66

78 55 1 140 82 69 145 199 3 8 2 176 106 218 155 0 65

79 50 3 150 90 71 170 326 1 7 1 195 98 204 178 0 64

80 39 1 114 72 65 156 187 5 8 3 168 110 199 170 0 66

81 50 3 150 115 71 220 283 3 7 1 164 120 275 239 0 67

82 42 3 105 78 67 166 195 3 8 2 112 80 218 168 0 67

83 53 4 100 80 71 199 209 2 8 1 110 65 220 185 0 63

84 56 1 150 90 72 233 284 1 8 5 134 74 170 225 0 68

85 62 4 166 90 66 130 258 3 7 3 150 70 185 126 0 64

87 57 2 110 68 71 166 300 3 8 1 105 60 262 143 0 68

88 43 1 120 85 70 134 220 3 8 4 140 85 354 141 0 63

89 65 1 170 105 67 183 214 5 3 4 190 100 148 161 0 65

90 49 2 120 90 69 139 273 3 8 3 160 104 308 163 0 64

91 47 1 110 70 70 130 203 3 8 1 125 75 198 127 0 66

92 59 1 110 80 70 167 220 3 8 1 125 75 275 152 0 67

93 53 4 120 90 62 166 253 3 8 4 134 78 316 168 0 67

94 43 1 120 80 72 171 198 3 8 3 152 86 190 178 0 65

95 52 3 140 86 67 128 300 1 8 1 185 95 260 144 0 64

96 65 4 125 85 67 164 228 5 8 1 130 70 240 152 0 64

97 50 4 110 70 67 162 239 3 8 3 70 70 233 155 0 68

98 33 4 106 80 67 151 191 2 8 4 128 70 207 142 0 66

99 49 4 120 80 63 142 283 3 8 5 116 88 311 160 0 67 100 62 2 100 65 69 141 224 3 8 2 170 78 242 125 0 65 101 45 4 130 90 70 200 220 3 8 2 130 90 203 168 0 67 102 53 3 125 88 67 167 226 3 8 3 160 100 212 165 0 67 103 30 2 160 85 71 155 187 3 8 1 130 95 184 156 0 63 104 25 1 110 80 74 190 235 3 8 2 116 90 280 210 0 63 105 26 3 118 80 65 120 328 2 8 3 132 90 354 137 0 63 106 69 1 160 90 67 185 314 4 8 3 170 85 191 159 0 66 107 60 2 140 80 68 170 356 1 6 5 192 74 250 150 0 0 108 45 3 130 100 72 168 252 5 3 1 150 100 234 176 0 0 109 55 1 110 80 75 198 358 3 8 3 110 80 264 177 0 0 110 63 3 190 100 66 187 207 5 3 2 148 70 210 188 0 0 111 52 1 170 100 65 164 218 5 7 1 110 70 261 130 0 0 112 53 1 200 140 67 197 210 3 3 1 215 100 139 164 0 64 113 55 3 118 82 69 124 265 3 8 5 132 90 284 124 0 66 114 48 3 120 85 68 161 267 2 8 1 110 75 258 162 0 0 115 50 3 105 70 65 161 325 4 8 3 125 75 186 114 0 0 116 44 1 130 80 69 202 246 3 8 1 130 60 176 165 0 64 117 49 4 120 80 69 189 295 3 8 1 140 75 305 203 0 0 118 61 1 150 90 69 142 247 3 6 2 150 90 218 174 0 67 119 42 4 120 85 67 192 250 5 8 1 162 98 268 207 0 0 120 47 1 110 80 71 228 250 5 8 2 128 88 249 207 0 0 121 40 2 100 70 68 169 260 2 6 1 112 70 269 151 0 0 122 33 3 125 88 67 149 220 3 8 1 135 80 225 166 0 0 123 29 3 130 90 70 173 280 3 8 1 145 90 308 180 0 0 124 43 3 120 80 71 164 260 3 8 1 125 85 271 176 0 0 125 51 4 130 90 69 193 290 2 8 3 120 80 254 178 0 0 126 57 3 160 90 63 144 280 4 7 3 170 80 321 152 0 0 127 30 3 115 86 70 172 210 2 8 1 108 80 201 187 0 0 128 44 2 120 90 67 178 260 2 8 3 130 96 242 173 0 0 129 44 2 120 90 72 196 240 3 8 4 150 84 240 193 0 0

130 38 2 106 80 67 181 210 5 8 3 110 80 208 179 0 0 131 35 1 124 90 68 189 320 3 8 4 138 88 338 201 0 0 132 34 3 126 85 64 165 310 2 8 1 130 90 296 167 0 0 133 30 1 104 70 69 161 300 3 8 1 135 85 290 178 0 0 134 51 3 140 92 69 170 310 1 8 3 172 100 282 163 0 0 135 51 4 120 80 70 200 260 2 8 2 144 80 283 169 0 0 136 57 3 108 76 66 161 200 2 8 3 106 75 150 161 0 0 137 46 2 115 90 69 189 300 3 8 2 110 60 169 150 0 0 138 50 2 132 88 69 220 220 3 8 3 182 100 204 246 0 0 139 25 4 120 84 72 180 220 3 8 3 130 75 222 179 0 0 140 55 2 120 92 65 154 310 2 8 2 126 72 265 138 0 0 141 47 2 130 75 67 145 260 3 8 4 138 76 289 163 0 0 142 49 1 130 80 64 162 230 5 8 1 195 105 220 194 0 0 143 47 3 130 95 64 163 280 4 8 1 130 70 190 152 0 0 144 34 3 105 65 65 137 220 5 8 2 118 82 306 152 0 0 145 47 2 120 80 68 152 220 3 8 1 145 76 270 160 0 0 146 46 1 120 80 68 152 270 2 8 3 150 100 326 166 0 0 147 45 2 120 82 71 171 240 1 8 4 140 84 238 169 0 0 148 37 1 150 105 69 205 220 3 7 3 180 110 214 187 0 0 149 36 3 120 90 67 188 220 5 8 4 125 82 189 183 0 0 150 42 1 110 70 71 162 190 3 8 3 155 90 174 187 0 0 151 57 3 130 90 71 181 260 2 8 3 190 104 288 185 0 0 152 43 2 120 80 69 201 300 1 8 4 138 82 204 198 0 0 153 58 3 110 78 69 175 300 3 8 4 148 80 229 198 0 0 154 34 1 110 70 72 157 220 3 8 1 135 80 243 158 0 0 155 29 3 148 98 69 203 260 5 8 2 150 110 245 194 0 0 156 43 2 110 75 67 143 190 3 8 3 124 90 204 139 0 0 157 48 3 122 94 73 198 250 1 8 3 170 120 276 211 0 0 158 41 3 120 80 64 147 220 4 8 1 170 110 333 160 0 0 159 54 3 120 78 69 137 210 3 8 4 160 80 329 143 0 0 160 47 1 110 80 66 187 210 3 8 1 120 85 279 188 0 0 161 35 1 110 80 70 154 250 3 8 1 110 75 256 140 0 0 162 48 2 142 90 67 163 220 4 8 1 180 80 223 167 0 0 163 49 1 112 78 64 149 240 3 8 4 125 72 253 153 0 0 164 25 3 115 78 73 180 160 3 8 1 120 80 216 200 0 0 165 48 4 150 110 66 189 220 1 3 3 180 105 224 178 0 0 166 42 1 120 90 75 207 210 3 8 3 140 94 219 224 0 0 167 34 3 152 102 68 185 300 3 7 3 195 125 280 180 0 0 168 39 3 130 92 67 144 210 3 8 2 128 88 204 151 0 0 169 46 3 108 75 71 140 260 3 8 3 120 75 228 148 0 0 170 30 2 106 80 70 173 260 4 8 5 118 80 254 186 0 0 171 64 4 104 74 63 146 300 3 8 1 130 80 249 143 0 0 172 47 3 104 80 73 177 390 3 8 1 110 70 339 189 0 0 173 36 3 135 80 69 155 290 3 8 2 120 84 238 157 0 0

174 30 4 120 80 70 130 260 3 8 4 135 80 260 157 0 0 175 48 3 130 92 70 198 230 3 8 3 140 100 235 204 0 0 176 34 1 110 80 66 155 250 3 8 1 180 100 278 121 0 0 177 26 3 110 75 66 136 230 3 8 3 118 76 183 134 0 0 178 22 4 120 90 69 192 240 5 8 1 150 105 269 201 0 0 179 26 2 88 60 70 178 300 2 8 2 112 76 308 206 0 0 180 54 2 100 64 71 173 300 4 8 1 125 75 259 202 0 0 181 53 3 138 90 66 215 310 4 8 2 148 88 336 202 0 0 182 48 2 120 75 70 203 220 3 8 1 130 70 273 213 0 0 183 42 3 115 75 71 204 220 2 8 3 140 80 225 224 0 0 184 42 3 122 78 66 128 190 2 8 2 112 66 196 127 0 0 185 60 4 128 80 70 141 230 3 8 3 140 80 224 122 0 0 186 35 2 100 78 68 141 220 3 8 2 144 90 239 138 0 0 187 52 2 120 70 68 135 250 3 8 1 135 80 300 137 0 0 188 44 2 150 95 61 127 260 3 3 1 155 95 312 130 0 0 189 41 2 142 96 66 145 300 2 8 4 148 85 315 143 0 0 190 28 2 115 65 66 150 220 2 8 3 124 86 276 149 0 0 191 36 3 118 84 63 153 300 2 8 2 124 80 252 149 0 0 192 54 3 148 90 71 140 220 1 8 1 230 105 265 143 0 0 193 49 2 110 75 70 141 150 3 8 4 170 98 221 147 0 0 194 49 4 105 75 69 144 250 3 8 1 110 70 254 144 0 0 195 61 3 122 78 64 133 180 3 8 2 170 72 229 142 0 0 196 50 1 115 80 66 148 300 2 8 1 115 65 273 152 0 0 197 23 1 110 70 69 137 120 3 8 2 112 76 198 153 0 0 198 20 3 130 80 66 150 210 5 0 1 130 85 274 158 0 0 199 46 3 140 84 66 138 130 4 6 2 148 88 160 157 0 0 200 36 1 100 70 70 157 260 3 8 3 120 86 251 152 0 0

SAS上机练习题(二)参考答案

6$6? ?д乬??? ?1乬 data a; input x@@; cards; 142.3 148.8 142.7 144.4 144.7 145.1 143.3 154.2 152.3 142.7 156.6 137.9 143.9 141.2 139.3 145.8 142.2 137.9 141.2 150.6 142.7 151.3 142.4 141.5 141.9 147.9 125.8 139.9 148.9 154.9 145.7 140.8 139.6 148.8 147.8 146.7 132.7 149.7 154.0 158.2 138.2 149.8 151.1 140.1 140.5 143.4 152.9 147.5 147.7 162.6 141.6 143.6 144.0 150.6 138.9 150.8 147.9 136.9 146.5 130.6 142.5 149.0 145.4 139.5 148.9 144.5 141.8 148.1 145.4 134.6 130.5 145.2 146.2 146.4 142.4 137.1 141.4 144.0 129.4 142.8 132.1 141.8 143.3 143.8 134.7 147.1 140.9 137.4 142.5 146.6 135.5 146.8 156.3 150.0 147.3 142.9 141.4 134.7 138.5 146.6 134.5 135.1 141.9 142.1 138.1 134.9 146.7 138.5 139.6 139.2 148.8 150.3 140.7 143.5 140.2 143.6 138.7 138.9 143.5 139.9 134.4 133.1 145.9 139.2 137.4 142.3 160.9 137.7 142.9 126.8 ; proc means data=a n mean median clm qrange std cv stderr maxdec=2; var x; proc univariate data=a normal; histogram x; var x; run; ??? The MEANS Procedure Analysis Variable : x N Mean Median Lower 95% CL for Mean Upper 95% CL for Mean Quartile Range Std Dev Coeff of Variation Std Error ? Ё? 95% ??? 95% ??? ? ? ? ? ? 130 143.22142.75142.10144.337.80 6.43 4.490.56

SAS上机练习题(全部,含参考答案)

重庆医科大学--卫生统计学统计软件包 SAS上机练习题(一) 1、SAS常用的窗口有哪三个?请在三个基本窗口之间切换并记住这些命令或功能键。 2、请在PGM窗口中输入如下几行程序,提交系统执行,并查看OUTPUT窗和LOG窗中内容,注意不同颜色的含义;并根据日志窗中的信息修改完善程序。 3、将第2题的程序、结果及日志保存到磁盘。 4、试根据如下例1的程序完成后面的问题: 表1 某班16名学生3门功课成绩表如下 问题: 1)建立数据集; 2)打印至少有1门功课不及格同学的信息;(提示,使用if语句) 参考程序: data a; input id sh wl bl; cards; 083 68 71 65 084 74 61 68

085 73 75 46 086 79 80 79 087 75 71 68 084 85 85 87 085 78 79 75 086 80 76 79 087 85 80 82 088 77 71 75 089 67 73 71 080 75 81 70 118 70 54 75 083 70 66 84 084 62 73 65 099 82 70 79 ; run; data b; set a; if sh<60 or wl<60 or bl<60then output; run; proc print data=b; var id sh wl bl;

run; 5、根据下列数据建立数据集 表2 销售数据 开始时间终止时间费用 2005/04/28 25MAY2009 $123,345,000 2005 09 18 05OCT2009 $33,234,500 2007/08/12 22SEP2009 $345,600 20040508 30JUN2009 $432,334,500 提示:(格式化输入;数据之间以空格分隔,数据对齐;注意格式后面的长度应以前一个位置结束开始计算,如果读入错误,可试着调整格式的宽度;显示日期需要使用输出格式) 开始时间,输入格式yymmdd10. 终止时间,输入格式date10. 费用,输入格式dollar12. 参考程序: data a; input x1 yymmdd10. x2 date10. x3 dollar13.; cards; 2005/04/28 25MAY2009 $123,345,000 2005 09 18 05OCT2009 $33,234,500 2007/08/12 22SEP2009 $345,600

SAS上机实习-8页文档资料

中科院研究生院统计分析与SAS实现第1次上机实习题 一、定量资料上机实习题 要求: (1)先判断定量资料所对应的实验设计类型; (2)假定资料满足参数检验的前提条件,请选用相应设计的定量资料的方差分析,并用SAS软件实现统计计算; (3)摘录主要计算结果并合理解释,给出统计学结论和专业结论。 【练习1】取4窝不同种系未成年的大白鼠,每窝3只,随机分配到三个实验组中,分别注射不同剂量雌激素,经过一定时间后处死大白鼠测子宫重量,资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义?若剂量间差别有统计学意义,请作两两比较。 表1 未成年大白鼠注射不同剂量雌激素后的子宫重量 子宫重量(g) 窝别 剂量(μg/100g):0.2 0.4 0.8 合计 1 106 116 145 367 2 42 68 115 225 3 70 111 133 314 4 42 63 87 192 合计260 358 480 1098 定量资料的随机区组设计(区组因素:窝别;实验因素:剂量)【SAS程序】: DATA PGM15G; DO A=1TO4; /*A为窝别*/ DO B=1TO3; /*B为雌激素剂量*/ INPUT X @@; /*X为子宫重量*/ OUTPUT; END; END; CARDS; 106 116 145 42 68 115 70 111 133 42 63 87 RUN; ods html; /*将结果输出成网页格式,SAS9.0以后版本可用*/ PROC GLM DATA=PGM15G; CLASS A B; MODEL X=A B / SS3; MEANS A B; /*给出因素A、B各水平下的均值和标准差*/ MEANS B / SNK; /*对因素B(即剂量)各水平下的均值进行两两比较*/ RUN; ODS HTML CLOSE; 【练习2】一位工程师研究由钻头压力产生的冲力。考察了A(钻孔速度)和B(进料

统计实验与SAS上机简易过程步

数据统计分析一般可遵循以下思路: (1)先确定研究目的,根据研究目的选择方法。不同研究目的采用的统计方法不同,常见的研究目的主要有三类:①差异性研究,即比较组间均数、率等的差异,可用的方法有t检验、方差分析、χ2检验、非参数检验等。②相关性分析,即分析两个或多个变量之间的关系,可用的方法有相关分析。③影响性分析,即分析某一结局发生的影响因素,可用的方法有线性回归、logistic回归、Cox 回归等。 (2)明确数据类型,根据数据类型进一步确定方法:①定量资料可用的方法有t检验、方差分析、非参数检验、线性相关、线性回归等。②分类资料可用的方法有χ2检验、对数线性模型、logistic回归等。下图简要列出了不同研究目的、不同数据类型常用的统计分析方法。 (3)选定统计方法后,需要利用统计软件具体实现统计分析过程。SAS中,不同的统计方法对应不同的命令,只要方法选定,便可通过对应的命令辅之以相应的选项实现统计结果的输出。 (4)统计结果的输出并非数据分析的完成。一般统计软件都会输出很多结果,需要从中选择自己需要的部分,并做出统计学结论。但统计学结论不同于专业结论,最终还需要结合实际做出合理的专业结论。 第一部分:统计描述

1.定量资料的统计描述指标及SAS实现; (1)数据分布检验:PROC UNIVARIATE ①基本格式: ②语句格式示例: 1.PROC UNIVARIATE normal;/*normal选项表示进行正态性检验*/ 2.CLASS group;/*指定group为分组变量*/ 3.VAR weight;/*指定分析变量为weight*/ 4.RUN; ③结果:正态性检验(tests for normality)结果,常用的是Shapiro-Wilk 检验和Kolmogorov-Smirnov检验。当例数小于2000时,采用Shapiro-Wilk检验W值为标准;当例数大于2000时,SAS中不显示Shapiro-Wilk检验结果,采用Kolmogorov-Smirnov检验D值为判断标准。正态性检验的P≤0.05提示不服从正态分布,P>0.05提示服从正态分布。 注:若服从正态分布,进行PROC MEANS过程步;若不服从则计算百分位数,转(3) (2)数据描述(符合正态分布的数据):PROC MEANS ①基本格式: 关键字(可以无视):不写任何关键字时默认输出n,mean,std,max,min; n:有效数据记录数(有效样本量) median:中位数 mean:均数 qrange:四分位数间距 std:标准差 var:方差 clm:95%可信区间 max、min:最大、最小值 ②语句格式示例: 1.PROC MEANS n mean std median qrange clm;/*关调用proc means过程, 要求输出的指标有例数、均值、标准差、中位数、四分位数间距、95% 可信区间*/ 2.CLASS group;/*指定group为分组变量*/ 3.VAR weight;/*指定分析变量为weight*/ 4.Run; ③结果以“均数±标准差”表示 (3)偏正态分布的统计描述: ①基本思想:计算中位数和百分位数,并且用“中位数(Q1~Q3)”表示 ②语句格式示例: 1.proc univariate data=aa; 2.var x; 3.output out=c pctlpre=P pctlpts=0 to 100 by 2.5;/*计算0到100

SAS上机练习题(全部-含参考答案).docx

重庆医科大学“卫生统计学统计软件包 SAS上机练习题(一) 1、SAS常用的窗口有哪三个?请在三个基本窗口之间切换并记住这些命令或功能键。 2、请在PGM窗口中输入如下几行程序,提交系统执行,并查看OUTPUT窗和LOG窗中内容,注意不同颜色的含义;并根据口志窗屮的信息修改完善程序。 DATS EXO 101; INPUTT NAME $ AGE SEX; CARDS; XIAOMIN 19 1 LIDONG 20 1 NANA 18 2 PROD PRONT DATS=EX1; RUN; PROC PRINT DATA=EX 1; VAR NAME AGE; RUN; 3、将第2题的程序、结果及tl志保存到磁盘。 4、试根据如下例1的程序完成后面的问题: 表1某班16名学生3门功课成绩表如下 问题: 1)建立数据集; 2)打印至少有1门功课不及格同学的信息;(提示,使用if语句) 参考程序: data a; input id sh wl bl; cards; 083 68 71 65 084 74 61 68 085 73 75 46 086 79 80 79 087 75 71 68 084 85 85 87 085 78 79 75 086 80 76 79 087 85 80 82 088 77 71 75 089 67 73 71 080 75 81 70 118 70 54 75 083 70 66 84 084 62 73 65 099 82 70 79 run; data b; set a;

if sh<60 or wl<60 or bl<60 then output; run; proc print data=b; var id sh wl bl; run; 5、根据下列数据建立数据集 表2销售数据 提示:(格式化输入;数据之间以空格分隔,数据对齐;注意格式后面的长度应以前一个位置结朿开始计算,如果读入错误,可试着调整格式的宽度;显示日期需要使用输出格式) 开始时间,输入格式yymmddlO. 终止时间,输入格式血tel0. 费用,输入格式dollar 12.

SAS上机练习题(全部_含参考问题详解)

医科大学--卫生统计学统计软件包 SAS上机练习题(一) 1、SAS常用的窗口有哪三个?请在三个基本窗口之间切换并记住这些命令或功能键。 2、请在PGM窗口中输入如下几行程序,提交系统执行,并查看OUTPUT窗和LOG窗中容,注意不同颜色的含义;并根据日志窗中的信息修改完善程序。 3、将第2题的程序、结果及日志保存到磁盘。 4、试根据如下例1的程序完成后面的问题: 表1 某班16名学生3门功课成绩表如下 问题: 1)建立数据集; 2)打印至少有1门功课不及格同学的信息;(提示,使用if语句) 参考程序: data a; input id sh wl bl; cards; 083 68 71 65 084 74 61 68 085 73 75 46 086 79 80 79 087 75 71 68 084 85 85 87 085 78 79 75 086 80 76 79

087 85 80 82 088 77 71 75 089 67 73 71 080 75 81 70 118 70 54 75 083 70 66 84 084 62 73 65 099 82 70 79 ; run; data b; set a; if sh<60 or wl<60 or bl<60then output; run; proc print data=b; var id sh wl bl; run; 5、根据下列数据建立数据集 表2 销售数据 开始时间终止时间费用 2005/04/28 25MAY2009 $123,345,000 2005 09 18 05OCT2009 $33,234,500 2007/08/12 22SEP2009 $345,600 20040508 30JUN2009 $432,334,500 提示:(格式化输入;数据之间以空格分隔,数据对齐;注意格式后面的长度应以前一个位置结束开始计算,如果读入错误,可试着调整格式的宽度;显示日期需要使用输出格式) 开始时间,输入格式yymmdd10. 终止时间,输入格式date10. 费用,输入格式dollar12. 参考程序:

SAS上机实习考试

1、为比较两种方法对乳酸饮料中脂肪含量测定结果是否不同,某人随机抽取了10份乳酸饮料制品,分别用脂肪酸水解法和哥特里-罗紫法测定其结果如表3-3第(1)~(3)栏。问两法测定结果是否不同? 表两种方法对乳酸饮料中脂肪含量的测定结果(%) 编号(1) 哥特里-罗紫法 (2) 脂肪酸水解法 (3) 1 0.840 0.580 2 0.591 0.509 3 0.67 4 0.500 4 0.632 0.316 5 0.687 0.337 6 0.978 0.517 7 0.750 0.454 8 0.730 0.512 9 1.200 0.997 10 0.870 0.506 2、为研究国产四类新药阿卡波糖胶囊的降血糖效果,某医院用20名II型糖尿病病人进行同期随机对照试验。试验者将这些病人随机等分到试验组(用阿卡波糖胶囊)和对照组(用拜唐苹胶囊),分别测得试验开始前和8周后的空腹血糖,算得空腹血糖下降值见表3-4,能否认为该国产四类新药阿卡波糖胶囊与拜唐苹胶囊对空腹血糖的降糖效果不同? 表3-4 试验组和对照组空腹血糖下降值(mmol/L) 试验组 (n1=10) -0.70 -5.60 2.00 2.80 0.70 3.50 4.00 5.80 7.10 -0.50 对照组 (n2=10) 3.70 6.50 5.00 5.20 0.80 0.20 0.60 3.40 6.60 -1.10 3. 两组肿瘤患者,单纯放疗组(A)13 例,口服平消胶囊+放疗组(B)12 例,接受放疗前后,血清Sil-2R水平(U/ml)如表所示。试评价平消胶囊对接受放疗患者血清Sil-2R水平的影响? 两组肿瘤患者的血清Sil-2R水平

SAS课程上机练习一

SAS课程上机练习一 下面的数据是一次对20岁以上人群进行心血管病随访研究的结果,数据中各变量的含义如下: V1—随访号V2—年龄V3—开始随访时的求诊机构V4/V5—开始随访时的收缩压/舒张压V6/V7—体重/身高V8—开始随访时的胆固醇V9—社会经济地位(1-高2-中3-底)V10—原有心血管病(0-其他心脏病1-冠心病2-冠心与高心3-高心4-高心及风心5-风心6-可疑心脏病7-高血压8-正常)V11—随访结束时的求诊机构V12/V13—随访结束时的收缩压/舒张压V14—随访结束时的胆固醇V15—随访结束时的体重V16—随访结束时的最后诊断(0-未诊断1-3心肌梗塞4-7心绞痛8-9其他)V17—死亡年份(63-68)0-未死。 一、数据集练习操作 1、用data步建立SAS永久数据集(下次练习可以调用此数据集) 2、将年龄(V2)按20—30—40—50—及60岁以上分组 3、计算体重指数={体重/身高2}*100 4、将随访结束时最后诊断的结果(V16)分为四组:(1)心肌梗塞(2) 心绞痛(3)其他心脏 病(4)未诊断 5、按体重指数的大小( 0.30)分为二类:1-超重0-未超重 6、将死亡年份(V17)按是否死亡分为二类(1为死亡,0为未死亡)。 7、将原有心血管病(V10)分为患有心血管病与不患有心血管病二类(1有0无) 8、将数据集按是否死亡分为二个数据集。 9、将社会经济地位(V9)设置为哑变量。 10、计算开始随访时的收缩压与/舒张压(V4/V5)之差。 二.。计算分析 1.年龄分布的特征。 2.开始时收缩压频数分布的特征。 3·超重病人开始时胆固醇(V8)的分布的特征。 3·比较不同体重的人开始时的收缩压。 4·比较不同社会经济地位与患心血管病的差别。 5·比较心肌梗塞的病人开始时的收缩压与随访结束时的收缩压的差别。 6·比较不同年龄、不同体重(指数)间开始时的胆固醇的差别。 7·试分析年龄与血压(V4)之间的关系。 9、比较不同社会经济地位开始随访是收缩压的差别。 三、结果存盘 1.要求将你编辑好的SAS程序保存。 2.要求将所有的计算结果(OUTPUT)保存.。 3.文件名的格式为:你的姓名+你的学号 例如:王永,学号为20918003,则你存盘的文件名为:王永208918003

SAS上机练习题(全部-含参考标准答案)

SAS上机练习题(全部-含参考答案)

————————————————————————————————作者:————————————————————————————————日期:

重庆医科大学--卫生统计学统计软件包 SAS上机练习题(一) 1、SAS常用的窗口有哪三个?请在三个基本窗口之间切换并记住这些命令或功能键。 2、请在PGM窗口中输入如下几行程序,提交系统执行,并查看OUTPUT窗和LOG窗中内容,注意不同颜色的含义;并根据日志窗中的信息修改完善程序。 DATS EX0101; INPUTT NAME $ AGE SEX; CARDS; XIAOMIN 19 1 LIDONG 20 1 NANA 18 2 ; PROD PRONT DATS=EX1; RUN; PROC PRINT DATA=EX1; V AR NAME AGE; RUN; 3、将第2题的程序、结果及日志保存到磁盘。 4、试根据如下例1的程序完成后面的问题: 表1 某班16名学生3门功课成绩表如下 问题: 1)建立数据集; 2)打印至少有1门功课不及格同学的信息;(提示,使用if语句) 参考程序: data a; input id sh wl bl; cards; 083 68 71 65 084 74 61 68 085 73 75 46 086 79 80 79 087 75 71 68 084 85 85 87

085 78 79 75 086 80 76 79 087 85 80 82 088 77 71 75 089 67 73 71 080 75 81 70 118 70 54 75 083 70 66 84 084 62 73 65 099 82 70 79 ; run; data b; set a; if sh<60 or wl<60 or bl<60then output; run; proc print data=b; var id sh wl bl; run; 5、根据下列数据建立数据集 表2销售数据 开始时间终止时间费用 2005/04/28 25MAY2009 $123,345,000 2005 09 18 05OCT2009 $33,234,500 2007/08/12 22SEP2009 $345,600 20040508 30JUN2009 $432,334,500 提示:(格式化输入;数据之间以空格分隔,数据对齐;注意格式后面的长度应以前一个位置结束开始计算,如果读入错误,可试着调整格式的宽度;显示日期需要使用输出格式) 开始时间,输入格式yymmdd10. 终止时间,输入格式date10. 费用,输入格式dollar12.

SAS上机实验

实验内容 1试用产生标准正态分布的随机数normal(seed)产生参数为10的卡方分布随机数100个。 2根据数据集: (1)创建一个仅包含地区、销售的产品类型、销售数量和销售额的数据集。 (2)分别创建一个仅包含产品类型a100和产品类型a200和SAS数据集。 (3)选择一个人口在50000以上的部分子集。 3 以下数据来自7位同学的高考语文、数学和英语成绩,试用编程的方法计算出平均成绩在75以上的男同学的人数。 实验步骤: 实验1: data a (drop=i) ; do i=1to100by1; z=normal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2+no rmal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2+normal(0)**2; output; end; 结果如下: 实验2: (1) data biao1;

input region$ product$ quantity price; cards; es a100 150 3750 so a100 410 10250 es a100 350 8750 so a100 710 17750 es a100 750 18750 so a100 760 19000 es a100 150 3000 so a100 410 8200 es a100 350 7000 so a100 710 14200 es a100 750 15000 so a100 760 152000 es a200 165 4125 so a200 425 10425 es a200 365 9125 ne a100 200 5000 we a100 180 4500 ne a100 600 15000 we a100 780 19500 ne a100 800 20000 we a100 880 22000 ne a100 200 4000 we a100 180 3600 ne a100 600 12000 we a100 780 15600 ne a100 800 16000 we a100 880 17600 ne a200 215 5375 we a200 195 4875 ne a200 615 15375 ; 实验结果:

0九研究生SAS软件上机试卷(1)

09硕/博士研究生课程《医学统计软件应用》上机试卷 (春学期) 一·内容 有一个名为:Hypertension的SAS数据文件,文件中的数据来自于杭州市某社区居民家庭健康档案一份2000例的随机样本,数据中各变量的含义如下: V1 地区编号:①城镇②农村 V2性别①男②女 V3 年龄(岁) V4文化程度 1.文盲或半文盲2小学 3.初中 4.中专或高中 5.大专6大学或以上 V5职业 1.专业技术人员(医生、教师等) 2.行政管理、文职人员 3.轻体力劳动工人(服务行业、钳工、检验工等) 4.重体力劳动工人(搬运工、建筑工等) 5.在校学生 6. 家务 V6你的家人是否有高血压①父母双方均有②父一方有③母一方有④父母双方均没有⑤不详 V7你的家人是否有糖尿病①父母双方均有②父一方有③母一方有④父母双方均没有⑤不详 V8你的家人是否有其他疾病①肿瘤②冠心病③脑卒中④高血脂 V9. 吸烟①是②否③已戒 V10.饮酒①是②不饮③已戒酒 V11..饮茶①是②偶尔③不饮 V12.锻炼你是否经常参加体育运动(包括散步.太极拳.跑步游泳跳舞等运动持续半小时以上为1次) ①每天参加②每周3次以上③每周1-2次 ④每月1-2次⑤不参加 V13饮食习惯①以荤食为主②以素食为主③荤素各半 V14吃.腌、霉食品①经常吃(每周3次以上) ②偶尔吃③不吃 V15食物口味咸淡程度①偏咸②偏淡③一般 V16甜食①.经常吃(每周3次以上) ②偶尔吃③不吃 V17工作压力①大②一般③没压力 V18 5年内家庭重大事件(亲属病故.伤亡.离异)①有②无 V19 身高(厘米)V20—体重(公斤)V21—腰围(厘米)V22—臀围(厘米) V23是否服降压药①是②否 V24—收缩压(mmHg) V25—舒张压(mmHg)

最新SAS上机实习

S A S上机实习

中科院研究生院统计分析与SAS实现第1次上机实习题 一、定量资料上机实习题 要求: (1)先判断定量资料所对应的实验设计类型; (2)假定资料满足参数检验的前提条件,请选用相应设计的定量资料的方差分析,并用SAS软件实现统计计算; (3)摘录主要计算结果并合理解释,给出统计学结论和专业结论。 【练习1】取4窝不同种系未成年的大白鼠,每窝3只,随机分配到三个实验组中,分别注射不同剂量雌激素,经过一定时间后处死大白鼠测子宫重量,资料见表1。问剂量和窝别的各自水平下子宫重量之间的差别有无统计学意义?若剂量间差别有统计学意义,请作两两比较。 表1 未成年大白鼠注射不同剂量雌激素后的子宫重量 窝别 子宫重量(g) 剂量(μg/100g):0.2 0.4 0.8 合计 1 106 116 145 367 2 42 68 115 225 3 70 111 133 314 4 42 63 87 192 合计260 358 480 1098 定量资料的随机区组设计(区组因素:窝别;实验因素:剂量)【SAS程序】: DATA PGM15G; DO A=1TO4; /*A为窝别*/ DO B=1TO3; /*B为雌激素剂量*/ INPUT X @@; /*X为子宫重量*/ OUTPUT; END; END; CARDS; 106 116 145 42 68 115 70 111 133 42 63 87 ; RUN; ods html; /*将结果输出成网页格式,SAS9.0以后版本可用*/ PROC GLM DATA=PGM15G; CLASS A B; MODEL X=A B / SS3;

上机练习题

上机练习题 1.括弧匹配检验work1.pas 假设表达式中允许包含两种括号:圆括号和方括号,其嵌套的顺序随意,如([]())或[([][])]等为正确的匹配,[(])或([]()或( ( ) ) )均为错误的匹配。 现在的问题是,要求检验一个给定表达式中的括弧是否正确匹配? 输入一个只包含圆括号和方括号的字符串,判断字符串中的括号是否匹配,匹配就输出“OK”,不匹配就输出“Wrong”。 输入一个字符串: [([][])] 输出: OK 【输入】 输入仅一行字符(字符个数小于255) 【输出】 匹配就输出“OK”,不匹配就输出“Wrong”。 【样例】 输入(work1.in) [(]) 输出(work1.out) Wrong 2.自然数的拆分work2.pas 任何一个大于1的自然数总可以拆分成若干个自然数之和。 1:4=1+1+1+1 2:4=1+1+2 3:4=1+3 4:4=2+2 5:4=4 分析:设拆分出的数s1≤s2≤…≤sk。定义数组s为一个栈,用来存放因子。从1开始搜索因子,求和,若sum ≤n就将因子压栈;若sum =n,则输出解,出栈;若sum >n,则修改栈顶元素的值,即回溯。 【输入】 输入仅一行,包含一个自然数n。 【输出】 这个自然数拆分结果。拆分出的因子要求满足s1+s2+…+sk=n 且s1≤s2≤…≤sk。具体要求见样例。 【样例】 输入(work2.in) 4 输出(work2.out)

1:4=1+1+1+1 2:4=1+1+2 3:4=1+3 4:4=2+2 5:4=4 3.糟糕的头发 work3.pas 农夫约翰有N (1 <= N <= 80,000)头奶牛,自从奶牛们独自照看自己的头发开始,它们的发型变得越来越糟糕。约翰决定统计奶牛们所能看到其他奶牛的头顶的数目总和。 第i头奶牛加上她头发的高度为h[i] (1 <= h[i] <= 1,000,000,000),奶牛们站成一排面向东方(见图1),所以,从第i头奶牛的发尖看过去,她能看见前面所有比她矮的奶牛头顶(这些可能的奶牛依次为i+1,i+2,……)。 例如 = = = = - = 奶牛面向右 --> = = = = - = = = = = = = = = 1 2 3 4 5 6 Cow#1 可以看到的奶牛编号为#2, 3, 4 Cow#2 什么也看不到 Cow#3 可以看到的奶牛编号为#4 Cow#4 什么也看不到 Cow#5 可以看到的奶牛编号为#6 Cow#6 什么也看不到 约翰用c[i] 表示每头奶牛能看到的奶牛数量,请你帮约翰统计所有c[i](1 <= i <= N)的总和。上面例子的统计结果就是:3 + 0 + 1 + 0 + 1 + 0 = 5 时间限制:0.5 seconds 【输入】 输入的第一行为一个单独的整数n表示奶牛的数目,以下i+1行中n个整数分别代表第i头羊的高度。 【输出】只需一行,c[1]~c[n]的和。 【样例】 输入(work3.in) 6 10 3 7 4 12 2 输出(work3.out) 5

new_SAS课程上机练习四

SAS课程上机练习四 1、用两种化疗方法治疗淋巴系统肿瘤,其疗效如下表。检验两种化疗方法的缓解率 有否差别? 两种化疗方法的缓解率比较 治疗方法缓解未缓解合计缓解率(%) 单纯化疗组 2 10 12 16.7 复合化疗组14 13 27 51.9 合计16 23 39 2.用两种方法检查110例肺结核患者,结果如下表,问用两种方法检查的结果有无差别? 110例肺结核患者分别用两种方法检查结果 甲法乙法合计 + — + 50 28 78 —12 20 32 合计62 48 110 3、某市中心血站采集血清标本673例,检出乙肝核心抗体55例,其血型分布如下 表。问不同血型的乙肝核心抗体阳性率有否差别? 乙肝核心抗体的阳性率与血型关系 血型检测数乙肝核心抗体阳性数阳性率(%) A 221 13 5.88 B 169 16 9.47 O 235 23 9.79 AB 48 3 6.25 合计673 55 8.17 4、为了解年龄与总胆固醇之间的数量关系,随机抽查了20名高血脂病人治疗前的总胆固醇水平(mmol/L)与年龄,测定值见表9.1。试求胆固醇对年龄的直线回归方程和相关系数。

5、某医院用中药结合化疗(中药组)和化疗(对照组)两种疗法治疗18例白血病患者,最初的缓解期(天)如下,试绘制两组的缓解率曲线并作比较。 中药组:10,13+,26,32,89,169+,198+,205+,220,254+ 对照组:4,12,20+,32,45,69,103,254+ (注:有+号的数据为截尾值) 6、某疾控中心对该地2496名55岁以下男性做了一次吸烟和喝咖啡对心肌梗死影响的病例对照研究,调查结果见下表。 吸烟及喝咖啡与心肌梗死关系调查结果 喝咖啡(杯/天) 吸烟(支/天) 0 1-20 21-35 35以上 病例对照病例对照病例对照病例对照 0 66 123 30 52 15 12 36 13 1-2 141 179 59 45 53 22 69 25 3-4 113 106 63 65 55 16 119 30 5以上129 80 102 58 118 44 373 85 请拟合下列模型: (1)模型1将喝咖啡(cofe)和吸烟(smoke)两个变量作为有序分类变量(0,1,2,3,),拟合模型。 (2)模型2将喝咖啡(cofe)和吸烟(smoke)两个变量以0为基准设为哑变量(0,1),拟合模型 (3)模型3将喝咖啡和吸烟两个变量作为交互效应变量,拟合模型。 7、根据练习一的数据,分析影响血压的主要因素有哪些? 8、某单位研究三种不同穴位的镇痛效果,资料如表4-16。试检验针刺不同穴位的镇痛效果有无差别?

相关主题
相关文档
最新文档