Html Version of SAS output from exemplar 5 ANALYSIS
.
To see programme that was used for it view the SAS program file.
Links in this page
Standard errors and design effects
Regression models
back to top First some exploratory analyses
| Cannabis | ||||
| q85a | Frequency | Percent | Cumulative Frequency |
Cumulative Percent |
| never used | 210 | 58.66 | 210 | 58.66 |
| tried once or twice | 90 | 25.14 | 300 | 83.80 |
| use daily | 12 | 3.35 | 312 | 87.15 |
| use weekly | 8 | 2.23 | 320 | 89.39 |
| used in last month | 11 | 3.07 | 331 | 92.46 |
| used more than a month ago | 27 | 7.54 | 358 | 100.00 |
| amphetamines | ||||
| q85b | Frequency | Percent | Cumulative Frequency |
Cumulative Percent |
| never used | 281 | 83.14 | 281 | 83.14 |
| tried once or twice | 32 | 9.47 | 313 | 92.60 |
| use daily | 1 | 0.30 | 314 | 92.90 |
| use weekly | 1 | 0.30 | 315 | 93.20 |
| used in last month | 5 | 1.48 | 320 | 94.67 |
| used more than a month ago | 18 | 5.33 | 338 | 100.00 |
| In general, would you say your health is.... | ||||
| genhelf | Frequency | Percent | Cumulative Frequency |
Cumulative Percent |
| excellent | 72 | 19.94 | 72 | 19.94 |
| very good | 139 | 38.50 | 211 | 58.45 |
| good | 120 | 33.24 | 331 | 91.69 |
| fair | 24 | 6.65 | 355 | 98.34 |
| poor | 6 | 1.66 | 361 | 100.00 |
Now survey nmeans with right kind of weights
| The SURVEYMEANS Procedure |
| Data Summary | |
| Number of Observations | 361 |
| Sum of Weights | 435.174782 |
| Statistics | ||||||
| Variable | Label | N | Mean | Std Error of Mean | Lower 95% CL for Mean |
Upper 95% CL for Mean |
| sinc sacc genhelf |
data zone income deprivation score data zone access deprivation score In general, would you say your health is.... |
361 361 361 |
18.711681 -0.067572 2.363313 |
0.694596 0.032883 0.052741 |
17.345705 -0.132240 2.259594 |
20.077657 -0.002904 2.467031 |
| The MEANS Procedure |
Mean squared weight gives value design effect would have if variables unrelated to weight
| Variable | Label | N | Mean | Std Dev | Minimum | Maximum |
| weight wtsq |
weight |
361 361 |
1.0007474 1.0974063 |
0.3101247 0.7729044 |
0.6304000 0.3974041 |
2.5298242 6.4000104 |
Now actual design factors for these variables (See SAS code for how calculated - a bit tiresome. design effect for income deprivation is slightly larger than the calculation from the mean squared weight, as we might expect since we know response depended on this..
| Obs | dfsinc | dfsacc | dfgenhelf | desinc | desacc | degenhelf |
| 1 | 1.14547 | 1.01905 | 1.06833 | 1.31211 | 1.03846 | 1.14133 |
| The SURVEYREG Procedure |
| Regression Analysis for Dependent Variable genhelf |
| Data Summary | |
| Number of Observations | 338 |
| Sum of Weights | 338.60539 |
| Weighted Mean of genhelf | 2.34349 |
| Weighted Sum of genhelf | 793.51678 |
| Fit Statistics | |
| R-square | 0.09035 |
| Root MSE | 0.8884 |
| Denominator DF | 337 |
| Class Level Information | |||
| Class Variable | Label | Levels | Values |
| living | composition of household | 5 | couple couple w kids other sing parent single |
| q85a | cannabis | 6 | never used tried once or tw use daily use weekly used in last mon used more than a |
| q85b | amphetamines | 6 | never used tried once or tw use daily use weekly used in last mon used more than a |
| ANOVA for Dependent Variable genhelf | |||||
| Source | DF | Sum of Squares | Mean Square | F Value | Pr > F |
| Model | 11 | 25.5985 | 2.327137 | 2.94 | 0.0010 |
| Error | 326 | 257.7389 | 0.790610 | ||
| Corrected Total | 337 | 283.3374 | |||
| Tests of Model Effects | |||
| Effect | Num DF | F Value | Pr > F |
| Model | 11 | 22.44 | <.0001 |
| Intercept | 1 | 361.12 | <.0001 |
| sinc | 1 | 6.09 | 0.0141 |
| q85a | 5 | 3.11 | 0.0093 |
| q85b | 5 | 15.71 | <.0001 |
|
| Estimated Regression Coefficients | ||||
| Parameter | Estimate | Standard Error | t Value | Pr > |t| |
| Intercept | 2.3111636 | 0.25331000 | 9.12 | <.0001 |
| sinc | 0.0134122 | 0.00543692 | 2.47 | 0.0141 |
| q85a never used | -0.1489457 | 0.17464000 | -0.85 | 0.3943 |
| q85a tried once or tw | -0.1176441 | 0.17512565 | -0.67 | 0.5022 |
| q85a use daily | 0.1516325 | 0.23036159 | 0.66 | 0.5108 |
| q85a use weekly | 0.1475498 | 0.30360484 | 0.49 | 0.6273 |
| q85a used in last mon | 0.8386110 | 0.29410484 | 2.85 | 0.0046 |
| q85a used more than a | 0.0000000 | 0.00000000 | . | . |
| q85b never used | -0.1427034 | 0.20778187 | -0.69 | 0.4927 |
| q85b tried once or tw | -0.3035456 | 0.21716450 | -1.40 | 0.1631 |
| q85b use daily | 0.6240744 | 0.21000573 | 2.97 | 0.0032 |
| q85b use weekly | 1.0047980 | 0.31863538 | 3.15 | 0.0018 |
| q85b used in last mon | 0.4212229 | 0.34530883 | 1.22 | 0.2234 |
| q85b used more than a | 0.0000000 | 0.00000000 | . | . |
|
| The SURVEYREG Procedure |
| Regression Analysis for Dependent Variable genhelf |
| Data Summary | |
| Number of Observations | 358 |
| Sum of Weights | 357.73923 |
| Weighted Mean of genhelf | 2.34130 |
| Weighted Sum of genhelf | 837.57366 |
| Fit Statistics | |
| R-square | 0.02247 |
| Root MSE | 0.9033 |
| Denominator DF | 357 |
| ANOVA for Dependent Variable genhelf | |||||
| Source | DF | Sum of Squares | Mean Square | F Value | Pr > F |
| Model | 1 | 6.6722 | 6.672154 | 8.18 | 0.0045 |
| Error | 356 | 290.2655 | 0.815352 | ||
| Corrected Total | 357 | 296.9376 | |||
| Tests of Model Effects | |||
| Effect | Num DF | F Value | Pr > F |
| Model | 1 | 9.18 | 0.0026 |
| Intercept | 1 | 1619.80 | <.0001 |
| canscore | 1 | 9.18 | 0.0026 |
|
| Estimated Regression Coefficients | ||||
| Parameter | Estimate | Standard Error | t Value | Pr > |t| |
| Intercept | 2.28375766 | 0.05674393 | 40.25 | <.0001 |
| canscore | 0.45467341 | 0.15003568 | 3.03 | 0.0026 |
|
| The GLM Procedure |
| Number of observations | 361 |
|
| The GLM Procedure |
| Dependent Variable: genhelf In general, would you say your health is.... |
| Weight: gweight weight adding to population totals |
| Source | DF | Sum of Squares | Mean Square | F Value | Pr > F |
| Model | 1 | 544.02811 | 544.02811 | 8.18 | 0.0045 |
| Error | 356 | 23667.40531 | 66.48148 | ||
| Corrected Total | 357 | 24211.43342 |
| R-Square | Coeff Var | Root MSE | genhelf Mean |
| 0.022470 | 348.2522 | 8.153617 | 2.341297 |
| Source | DF | Type I SS | Mean Square | F Value | Pr > F |
| canscore | 1 | 544.0281087 | 544.0281087 | 8.18 | 0.0045 |
| Source | DF | Type III SS | Mean Square | F Value | Pr > F |
| canscore | 1 | 544.0281087 | 544.0281087 | 8.18 | 0.0045 |
| Parameter | Estimate | Standard Error | t Value | Pr > |t| |
| Intercept | 2.283757665 | 0.05180502 | 44.08 | <.0001 |
| canscore | 0.454673407 | 0.15894220 | 2.86 | 0.0045 |
| The GLM Procedure |
| Number of observations | 361 |
|
| The GLM Procedure |
| Dependent Variable: genhelf In general, would you say your health is.... |
| Source | DF | Sum of Squares | Mean Square | F Value | Pr > F |
| Model | 3 | 18.5858474 | 6.1952825 | 8.01 | <.0001 |
| Error | 334 | 258.2336793 | 0.7731547 | ||
| Corrected Total | 337 | 276.8195266 |
| R-Square | Coeff Var | Root MSE | genhelf Mean |
| 0.067141 | 38.24978 | 0.879292 | 2.298817 |
| Source | DF | Type I SS | Mean Square | F Value | Pr > F |
| sinc | 1 | 6.59991804 | 6.59991804 | 8.54 | 0.0037 |
| canscore | 1 | 9.79435975 | 9.79435975 | 12.67 | 0.0004 |
| ampscore | 1 | 2.19156958 | 2.19156958 | 2.83 | 0.0932 |
| Source | DF | Type III SS | Mean Square | F Value | Pr > F |
| sinc | 1 | 5.98174624 | 5.98174624 | 7.74 | 0.0057 |
| canscore | 1 | 6.62583762 | 6.62583762 | 8.57 | 0.0037 |
| ampscore | 1 | 2.19156958 | 2.19156958 | 2.83 | 0.0932 |
| Parameter | Estimate | Standard Error | t Value | Pr > |t| |
| Intercept | 2.021651909 | 0.08727740 | 23.16 | <.0001 |
| sinc | 0.012601011 | 0.00453028 | 2.78 | 0.0057 |
| canscore | 0.492456944 | 0.16822134 | 2.93 | 0.0037 |
| ampscore | 0.472872905 | 0.28086654 | 1.68 | 0.0932 |