sas统计分析代码内容

/*统计描述*/

/*制作频数表*/

dataal;

inputx@@;

low=;

dis=;

z=x-mod(x-low,dis);

cards;

12

34

/

procfreq;

tablesz;

run;

/*用meaNS过程计算例数、均数、标准差、标准误和95%可信区间*/

procmeansnmeanstdstderrclm;

varx;

run;

/*用means过程对频数表计算例数、均数、标准差、最大值最小值*/

procmeans;

freqf;

varx;

run;

procunivariate;

freqf;

varx;

run;

/*正态性检验*/

procunivariatenormalplot;

varx;

run;

/*t检验*/

dataaa3;

inputx;

Y=x-比较的数；

/*配对资料d=xl-x2*/

cards;

/

procunivariatenormal;

vary;

run;

/*用MEANS作配对资料两个样本均数比较的t检险*/

dataex3_6;

inputxlx2@@;

d=xl-x2;

cards,

0.8400.580

0.5910.509

0.6740.500

0.6320.316

0.6870.337

0.9780.517

0.7500.454

0.7300.512

1.2000.997

0.8700.506

procmeanstprt;

vard;

run;

/*用UNIVARIATE过程作配对资料两样本均数比较的t检验*/

procunivariatedata=ex3_6;

vard;

run,

procttestdata=ex3_6;

vard;

run;

dataex3_7;

inputx@@;

if_n_<21thenc=l;

elsec=2;

cards;

-0.70-5.602.002.800.703.5C4.005.807.10-0.50

2.50-1.601.703.000.404.504.602.506.00-1.40

3.706.505.005.200.800.203.603.406.60-1.10

6.003.802.001.602.002.201.203.101.70-2.00

/

procttest;

varx;

classc;

run;

/*方差分析*/

DATARANDOM;

INPUTXGROUP

CARDS;

3.531

4.591

4.341

2.661

3.591

3.131

2.641

2.561

3.501

3.251

3.301

4.041

3.531

3.561

3.851

4.071

3.521

3.931

4.191

2.961

1.371

3.931

2.331

2.981

4.001

3.551

2.961

4.301

4.161

2.591

2.422

3.362

4.322

2.342

2.682

2.952

1.562

3.112

1.812

1.772

1.982

2.632

2.862

2.932

2.172

2.722

2.652

2.222

2.902

2.972

2.362

2.562

2.522

2.272

2.982

3.722

2.802

3.572

4.022

2.312

2.863

2.283

2.393

2.283

2.483

2.283

3.213

2.233

2.323

2.683

2.663

2.323

2.613

3.643

2.583

3.653

2.663

3.683

2.653

3.023

3.483

2.423

2.413

2.663

3.293

2.703

3.043

2.813

1.973

1.683

0.894

1.064

1.084

1.274

1.634

1.894

1.194

2.174

2.284

1.724

1.984

1.744

2.164

3.374

2.974

1.694

0.944

2.114

2.814

2.524

1.314

2.514

1.884

1.414

3.194

1.924

2.474

1.024

2.104

3.714

/

RUN;

/*正态性检验*/

PROCUNIVARIATENORMAL;

CLASSGROUP;

VARX;

RUN;

PROCANOVA;

CLASSGROUP;

MODELX=GROUP;

MEANSGROUP/HOVTESTSNKLSDDUNNETT;

RUN;

/*hovtest方差齐性检验*/

/*方差分析*/

PROCGLM;

CLASSGROUP；

MODELX=GROUP;

MEANSGROUP/HOVTESTSNKLSDDUNNETT;

RUN;

DATAA;

INPUTXTREATBLOCK

CARDS;

0.8211

0.7312

0.4313

0.4114

0.6815

0.6521

0.5422

0.3423

0.2124

0.4325

0.5131

0.2332

0.2833

0.3134

0.2435

RUN;

PROCPRINT;RUN;

PROCANOVA;

CLASSTREATBLOCK;

MODELX=TREATBLOCK;

MEANSTREAT/SNKALPHA=0.05;

MEANSTREAT;

RUN;

DATABLOCK;

INPUTBLOCKTREATX@@;

CARDS;

110.82120.65130.51

210.73220.54230.23

310.43320.34330.28

410.41420.21430.31

510.68520.43530.24

PROCANOVA;

CLASSBLOCKTREAT;

MODELX=BLOCKTREAT;

MEANSTREAT/HOVTESTSNKLSDDUNNETT;

RUN;

/

PROCGLM;

CLASSBLOCKTREAT;

MODELX=BLOCKTREAT;

MEANSTREAT/HOVTESTSNKLSDDUNNETT;

RUN;

/*roc曲线分析*/

dataex21;

inputgroupnum;

doi=ltonum;

inputvalue@@;

output;end;

cards;

125

6.513.512.86.213.914.79.59.06.916.8

13.310.812.214.9L3.712.85.311.812.4

7.613.311.911.212.312.7

020

8.56.44.61.79.75.34.95.73.86.5

6.35.43.34.78.66.35.94.84.55.2

;

run;

proclogistiedescending;

modelgroup=value/scaLe=noneoutroc=rocl;

run;

procprint;

run;

procgplot;

plot_SENSIT_*_lMSPEC_;

run;

DATACHISQGROUP;

INPUTCOLROWX00;

CARDS；

1146126

2118228

;

PROCFREQ;

WEIGHTX;

TABLECOL*ROW/CHISQ;

RUN;

DATACH1SQPA工RED;/*请注意配对结果中没有进行校正，所以对于需要校正的情况，需要手动计

算。但是别的一些软件可以实现校正。*/

INPUTCOLROWxee；

CARDS；

11111212

2122233

PROCFREQ;

WEIGHTX;

TABLECOL*ROW/AGREE;

RUN;

/*行*列表资料注意资料录入格式*/

/*非参数检验过程步格式*/

DATAPAIRED;

INPUTX1X2@@;

D=X1-X2;

CARDS;

6080

142152

195243

8082

242240

220220

190205

2538

212243

3844

236200

95100

PROCUNIVARIATE;

VARD;

RUN;

PROCNPARlWAY<options>;

BYvariables;

CLASSvariable;

EXACTStatistic-options</computation-options>;

FREQvariable;

OUTPUT<OUT=SAS-data-set><options>;

VARvariables;

RUN;

DATAGROUPED;

INPUTXGROUPS

CARDS;

2.7813.232

3.2313.502

4.2014.042

4.8714.152

5.1214.282

6.2114.342

7.1814.472

8.0514.642

8.5614.752

9.6014.822

4.952

5.102

PROCNPAR1WAY;

CLASSGROUP;

VARX;

RUN;

DATAmysas.A1127;

INPUTIDAGEINDEX;

CARDS;

1133.54

2113.01

393.09

462.48

582.56

6103.36

7123.18

872.65

9

PROCCORR;

VARAGEINDEX;

RUN;

PROCREG;

MODELINDEX=AGESELECTZON=STEPWISE

SLE=0.10SLS=0.10;

PLOTINDEX*AGE;

RUN;

proclogistiedescending;

modely=xSELECTION=STEPWISE

SLE=0.10SLS=0.10;

run;

/*条件逻辑回归*/

dataa;

inputiyxl-x6;

yy=2-y；

proephreg;

modelyy*y(0)=x/riskLimitsselection=stepwise

sle=0.lsls=0.lties=dLScrete;

STRATA=i;

run;

/*生存分析km法*/

PROCLIFETESTMETHOD=knPLOTS=(S)；

TIMEtime*censor(1);

STRATAgroup;

RUN;

dataa;

inputtimegroupcensor

cards;

111011510600025400

3211124101020028100

4210174111650029501

4210186101890145800

71102781020000

;

PROCLIFETESTMETHOD=knPLOTS=(S);

TIMEtime*censor(1);

STRATAgroup;

RUN;

/*生存分析寿命表法km*/

DATAsampl2_3;

INPUTtcensorf

CARDS;

0176009

1168108

2142203

3118301

419400

5105024

PROCPRINT;

PROCLIFETESTMETHOD=ltPLOTS=(s)WIDTH=1;

TIMEt*censor(0);

FREQf;

RUN;

/*cox回归*/

dataa;

inputnumsexagestagebloodxraychmthrpcensorday;

cards;

114522011578

2036220111549

315722101938

4045201004717

5042201114111

6039210111245

7138211114435

8145221013750

9130201013958

10045210112581

11045310113572

12157211012938

13057220111932

14149221113205

15133210113451

16051221012363

r

PROCPHREG;

Modelday*censor(0)=sexagestagebloodxraychmthrp/

risklimits;

RUN;

/*判别分析*/

datad2;

doi=ltol6;

dogrp=lto2；

inputidxlx200;output;end;end;

cards;

19.865.18110.662.07

213.333.73212.534.45

1513.335.961511.203.42

16•*169.333.63

run;

procdiscrimdata=d2

outstat=outanovamanovapool=testlisterrposterr;

classgrp;

varxlx2;

priorsprop;

run;

/*判别分析*/

dataexl84;

inputxl-x4g;

cards;

6.0-11.5二9901

-11.0-18.525-363

90.2-17.01732

-4.0-15.013541

0.0-14.020352

0.5-11.5二9373

-10.0-19.021-423

0.0-23.05-351

20.0-22.08-203

-100.0-21.47-151

-100.0-21.515-402

13.0-17.21822

-5.0-18.515181

10.0-18.014501

-8.0-14.016561

0.6-13.026213

-40.0-20.022-503

procdiscrim;

classg;

varxl-x4;

priorsprop;run;

/*样品聚类*/

dataexl9_3;

inputxl-xlO;

cards;

4625521381.680.358.114

41

35122035102.761.436.843

31

52252027842.190.544.113

31

3272024511.930.4711.459

61

3822032472.560.8011.685

51

51313037102.920.3711.602

21

4091031942.510.4011.405

51

34172046583.670.4611.353

31

5029050193.950.4713.4510

81

42202074825.890.1213.110

02

57301538002.990.1910.762

2■

36152024781.950.2510.000

01

3712038273.010.8210.504

41

5232029842.350.1611.153

31

52321037492.950.7211.4511

101

42273049413.890.7313.807

61

44272039483.110.3313.6516

141

4021533602.640.3711.400

01

3821529362.310.6911.401

11

44272068515.390.9912.287

62

4327039263.090.4711.950

0■

2610343813.450.5211.807

51

37182071425.620.8511.815

52

2892026122.060.3711.651

11

2593026382.080.7812.251

11

34142043223.400.4115.005

51

50322028622.250.698.802

21

f

proccluster

method=average;

varxl-xlO;

run;

proctree;

run;

/*指标聚类*/

procvarclusdata=exl9_3;

varxl-xlO;

run;

/*主成分分析*/

dataex20_l;

inputxl-x4;

cards;

402.0520

101.5530

1203.01350

2504.5180

1203.5950

101.51250

401.01940

2704.01360

2803.51160

1703.0960

1803.51440

1302.03050

2201.51720

1601.53560

2202.51430

1402.02020

2202.01410

401.0100

201.01260

1202.0200

procprincomp;

varxl-x4;

run;

/*因子分析*/

dataex20_2;

inputxl-x9;

cards;

4.3438999.061.2325.4693.153.56

97.5161.66

3.4527188.280.8523.5594.312.44

97.9473.33

4.38385103.971.2126.5492.534.02

98.4876.79

4.1837799.481.1926.8993.862.92

99.4163.16

4.32378102.011.1927.6393.181.99

99.7180.00

4.1334997.551.1027.3490.634.38

99.0363.16

4.5736191.661.1424.8990.602.73

99.6973.53

4.3120962.180.5231.7491.673.65

99.4861.11

4.0642583.270.9326.5693.813.09

99.4870.73

4.4345892.390.9524.2691.124.21

99.7679.07

4.1349695.431.0328.7593.433.50

99.1080.49

4.1051492.991.0726.3193.244.22

100.0078.95

4.1149080.900.9726.9093.684.97

99.7780.53

3.5334479.660.6831.8794.773.59

100.0081.97

4.1650890.981.0129.4395.752.77

98.7262.86

4.1754592.981.0826.9294.893.14

99.4182.35

4.1650795.101.0125.8294.412.80

99.3560.61

4.8654093.171.0727.5993.472.77

99.8070.21

5.0655284.381.1027.5695.153.10

98.6369.23

4.0345372.690.9026.0391.944.50

99.0560.42

4.1552986.531.0522.4091.523.84

98.5868.42

3.9451591.011.0225.4494.882.56

99.3673.91

4.1255289.141.1025.7092.653.87

95.5266.67

4.4259790.181.1826.9493.033.76

99.2873.81

3.0543778.810.8723.0594.464.03

96.2287.10

3.9447787.340.9526.7891.784.57

94.2887.34

4.1463888.571.2726.5395.161.67

94.5091.67

3.8758389.821.1622.6693.433.55

94.4989.07

4.0855290.191.1022.5390.363.47

97.8887.14

4.1455190.811.0923.0691.652.47

97.7287.13

4.0457481.361.1426.6593.741.61

98.2093.02

3.9351576.871.0223.8893.823.09

95.4688.37

3.9055580.581.1023.0894.382.06

96.8291.79

3.6255487.211.1022.5092.433.22

97.1687.77

3.7558690.311.1223.7392.472.07

97.7493.89

3.7762786.471.2423.2291.173.40

98.9889.80

*

procfactor

n=4;

run;

procfactor

n=4

rotate=quartimax;

run;

dataex21_l;

inputxl-x4yl-y4;

cards;

1210120.123.861.010.266.32.012.73

1210120.723.459.811.367.61.922.71

1040121.222.959.010.166.51.922.60

1620121.524.659.59.567.81.952.64

1690122.524.460.711.069.22.082.64

1150122.727.264.510.569.12.192.84

1150123.220.056.110.459.31.832.61

1460123.324.958.410.569.02.012.72

1190123.421.859.010.667.41.902.71

1840123.923.560.29.667.12.002.84

1250124.525.263.011.267.82.052.78

1480124.822.358.110.767.92.052.73

1310124.922.058.010.567.81.982.68

1660125.324.760.010.869.31.952.80

1580125.622.859.09.469.12.002.65

1460125.825.761.010.269.61.952.70

1240126.030.268.09.267.12.142.88

1100126.225.260.59.868.41.982.72

1250126.823.658.510.267.51.942.74

1270127.123.057.710.869.81.902.78

1300127.624.359.010.367.91.932.84

1350127.724..60.011.069.72.032.77

1250128.321.655.510.468.51.832.70

1720128.527.162.011.471.22.032.75

1480128.522.657.410.067.32.042.83

1380129.424.960.511.569.82.042.76

1170129.026.763.79.667.42.132.98

1640129.826，62.09.871.02.002.84

1640131.628.762.89.770.71.892.89

1150130.225.058.610.571.81.962.78

1430130.526.L60.710.868.62.052.77

1150130.623.454.411.869.21.962.78

1150131.425.563.210.270.42.052.84

1320131.625.658.910.970.22.062.86

1360131.727.462.010.973.51.992.70

1460132.026.361.511.171.22.172.13

1380132.225.761.410.170.11.962.83

1300132.524.557.010.871.82.022.84

1220132.727.061.310.172.22.082.80

1320132.925.260.511.273.12.012.73

1910133.130.167.09.087.12.152.97

1800133.526.562.59.871.72.072.82

1560133.624.858.510.372.21.932.79

1840134.026.060.510.473.01.982.74

1470134.328.262.011.387.22.664.03

1590134.425.560.79.669.91.992.81

1430134.126.663.011.272.22.062.90

1760134.632.566.09.987.42.612.98

1470135.327.961.810.173.32.202.78

1580135.628.165.89.873.12.052.89

1580136.528.262.011.872.92.172.92

1840137.127.662.89.572.42.112.91

1810137.428.362.59.474.22.063.00

1850138.129.562.49.772.32.124.02

2120140.034.968.89.587.92.744.15

1760140.732.064.410.274.02.174.05

1800141.032.563.89.588.22.654.08

1260141.729.165.0