char16.docx

char16.docx

《char16.docx》由会员分享，可在线阅读，更多相关《char16.docx（18页珍藏版）》请在冰豆网上搜索。

char16

例16.1计算Pearson相关系数及其它关联测度。

proccorrdata=fdata.fitnesspearsonspearmanhoeffding;

varweightoxygenruntime;

title'MeasuresofAssociationfor';

title2'aPhysicalFitnessStudy';

run;

例16.2计算并输出相关指标。

proccorrdata=fdata.a1a0001pearsonspearmankendallhoeffding;

varopprhiprloprclpr;

title‘Spearman的rho,Kendall的tau-b,Pearson和Hoeffding相关’;

run;

proccorrdata=fdata.a1a0001csscpcov;

varopprhiprlopr;

partialclpr;

title‘偏相关阵’;

run;

proccorrdata=fdata.a1a0001covalphaoutp=corrout;

varopprhiprlopr;

title‘协方差阵和相关阵’;

run;

procprintdata=corrout;

title2‘PROCCORR产生的输出数据集’;

run;

例16.3计算两个数据集中相同变量之间的相关系数。

dataa;/*数据集准备*/

mergefdata.a1a0001（keep=dateopprclpr）

fdata.szcz（keep=dateopprclprrename=（oppr=oppr_szclpr=clpr_sz））;

bydate;

run;

proccorrdata=anomisscov;

varoppr_szclpr_sz;

withopprclpr;

title2‘长方形的COV和CORR阵’;

run;

proccorrdata=acovcsscpoutp=oup;

title2‘从含有缺失值的数据集中计算CSSCP和COV’;

run;

例16.4计算Cronbach's系数Alpha。

optionsnodatepageno=1linesize=80pagesize=60;

proccorrdata=fdata.psychdatalphanocorrnomiss;

run;

例16.5输出偏相关系数到一个数据集中。

optionsnodatepageno=1linesize=120pagesize=60;

proccorrdata=fdata.fitnessspearmankendallcovnosimpleoutp=fitcorr;

varweightoxygenruntime;

partialage;

labelage='Ageofsubject'

weight='Wtinkg'

runtime='1.5miinminutes'

oxygen='O2use';

title1'PartialCorrelationsforaFitnessandExerciseStudy';

run;

例16.6按格式化值的顺序排列。

procformat;

value$sfmt'M'='male''F'='female';

procfreqdata=fdata.classorder=formatted;

tablesex;

formatsex$sfmt.;

run;

例中，计算数据集FDATA.CLASS中变量SEX的分布，并以格式化值的顺序排列。

输出结果为：

CumulativeCumulative

Sex频数PercentFrequencyPercent

-------------------------------------------------------------------------------

F947.37947.37

M1052.6319100.00

例16.7TABLE语句举例。

tablesa*（bc）;等价于Tablesa*ba*c;

tables（ab）*（cd）;等价于Tablesa*cb*ca*db*d;

tables（abc）*d;等价于Tablesa*db*dc*d;

tablesa-c;等价于Tablesabc;

tables（a-c）*d;等价于Tablesa*db*dc*d;

例16.8随机数频数分析。

dataa;

doI=1to1000;

X=int（uniform（8888）*3）+1;

Y=int（uniform（8888）*4）+1;

output;

end;

procfreqdata=a（drop=i）;

title’没有TABLES语句’;

run;

title;

procfreq;

tablesxx*y/chisq;

run;

procfreq;

tablesx*y/list;

run;

例16.9利用TABLE语句创建输出频率数据集。

optionsnodatepageno=1linesize=80pagesize=60;

procfreqdata=fdata.color;

weightcount;

tableseyeshaireyes*hair/out=freqcntoutexpectsparse;

title'EyeandHairColorofEuropeanChildren';

run;

procprintdata=freqcntnoobs;

title2'OutputDataSetfromPROCFREQ';

run;

例16.10对One-Way频率表作卡方检验。

procsortdata=fdata.color;

byregion;

run;

procfreqdata=fdata.colororder=data;

weightcount;

tableshair/nocumtestp=（301230253）;

byregion;

title'HairColorofEuropeanChildren';

run;

例16.11计算One-Way频率表的二项比例。

procfreqdata=fdata.colororder=freq;

weightcount;

tableseyes/binomialalpha=.1;

tableshair/binomial（p=.28）;

title'HairandEyeColorofEuropeanChildren';

run;

例16.122×2列联表分析。

optionsnodatepageno=1linesize=84pagesize=64;

procformat;

valueexpfmt1='HighCholesterolDiet'

0='LowCholesterolDiet';

valuerspfmt1='Yes'

0='No';

run;

procsortdata=fdata.fatcomp;

bydescendingexposuredescendingresponse;run;

procfreqdata=fdata.fatcomporder=data;

weightcount;

tablesexposure*response/chisqrelrisk;

exactpchior;

formatexposureexpfmt.responserspfmt.;

title'Case-ControlStudyofHighFat/CholesterolDiet';

run;

例16.13创建包含卡方统计量的数据集。

optionsnodatepageno=1pagesize=60;

procfreqdata=fdata.colororder=data;

weightcount;

tableseyes*hair/chisqexpectedcellchi2norownocol;

outputout=chisqdatpchilrchinnmiss;

title'Chi-SquareTestsfor3by5TableofEyeandHairColor';

run;

procprintdata=chisqdatnoobs;

title'Chi-SquareStatisticsforEyeandHairColor';

title2'OutputDataSetfromtheFREQProcedure';

run;

例16.14计算分层表的Cochran-Mantel-Haenszel统计量。

optionsnodatepageno=1linesize=80pagesize=60;

procfreqdata=fdata.migraine;

weightfrequency;

tablesgender*treatment*improve/cmhnoprint;

title1'ClinicalTrialforTreatmentofMigraineHeadaches';

run;

例16.15进行Cochran-Armitage趋势检验。

optionsnodatepageno=1linesize=80pagesize=72;

procfreqdata=fdata.pain;

weightcount;

tablesdose*adverse/trendmeasurescl;

testsmdcr;

exacttrend/maxtime=60;

title1'ClinicalTrialforTreatmentofPain';

run;

例16.16信用等级转移矩阵计算。

procfreqdata=fdata.rating;

table_col3*_col4_col4*_col5/nocolnorownopercent;

run;

例16.17分组求均值。

procsortdata=fdata.dd_c;

bydate;

procmeansdata=fdata.dd_cnoprint;

bydate;

varsum;

outputout=sumsum=sum_c;

run;

例16.18输出常用统计量到数据集中。

procmeansdata=fdata.r_month;

varrmrfr600600;

outputout=statsum=s_rms_rfs_r600600n=n_rmn_rfn_r600600

mean=M_rmM_rfM_r600600std=std_rmstd_rfstd_r600600;

run;

procmeansdata=fdata.r_monthnmeanmaxminrangestdfw=8;

varrmrf;

run;

例16.19计算特定的描述统计量。

optionsnodatepageno=1linesize=80pagesize=60;

procmeansdata=fdata.cakenmeanmaxminrangestdfw=8;

varPresentScoreTasteScore;

title'SummaryofPresentationandTasteScores';

run;

例16.20计算有CLASS变量时的描述统计量。

procmeansdata=fdata.grademaxdec=3;

classStatusYear;

types（）status*year;

varScore;

title'FinalExamGradesforStudentStatusandYearofGraduation';

run;

例16.21BY语句和CLASS变量一起使用。

optionsnodatepageno=1linesize=80pagesize=60;

procsortdata=fdata.Gradeout=GradeBySection;

bysection;

run;

procmeansdata=GradeBySectionminmaxmedian;

bysection;

classStatusYear;

varScore;

title1'FinalExamScoresforStudentStatusandYearofGraduation';

title2'WithinEachSection';

run;

例16.22CLASSDATA=选项和CLASS变量一起使用。

optionsnodatepageno=1linesize=80pagesize=60;

procmeansdata=fdata.cakerangemedianminmax

fw=7maxdec=0classdata=fdata.caketypeexclusiveprintalltypes;

classflavorlayers;

varTasteScore;

title'TasteScoreForNumberofLayersandCakeFlavor';

run;

例16.23多标签值和CLASS变量一起使用。

optionsnodatepageno=1linesize=80pagesize=64;

procformat;

value$flvrfmt

'Chocolate'='Chocolate'

'Vanilla'='Vanilla'

'Rum',

'Spice'='OtherFlavor';

valueagefmt（multilabel）

15-29='below30years'

30-50='between30and50'

51-high='over50years'

15-19='15to19'

20-25='20to25'

25-39='25to39'

40-55='40to55'

56-high='56andabove';

run;

procmeansdata=fdata.cakefw=6nminmaxmediannonobs;

classflavor/order=freq;

classage/mlforder=fmt;

typesflavorflavor*age;

varTasteScore;

formatageagefmt.flavor$flvrfmt.;

title'TasteScoreforCakeFlavorsandParticipant''sAge';

run;

例16.24预先规定的输出格式与CLASS变量一起应用。

optionsnodatepageno=1linesize=80pagesize=64;

procformat;

valuelayerfmt1='singlelayer'

2-3='multi-layer'

.='unknown';

value$flvrfmt（notsorted）

'Vanilla'='Vanilla'

'Orange','Lemon'='Citrus'

'Spice'='Spice'

'Rum',

'Mint',

'Almond'='OtherFlavor';

run;

procmeansdata=fdata.cakefw=7completetypesmissingnonobs;

classflavorlayers/preloadfmtexclusiveorder=data;

ways12;

varTasteScore;

formatlayerslayerfmt.flavor$flvrfmt.;

title'TasteScoreForNumberofLayersandCakeFlavors';

run;

例16.25计算均值的置信限。

procmeansdata=fdata.charityfw=8maxdec=2alpha=.1clmmeanstd;

classYear;

varMoneyRaisedHoursVolunteered;

title'ConfidenceLimitsforFundRaisingStatistics';

title2'1992-94';

run;

例16.26计算输出统计量。

optionsnodatepageno=1linesize=80pagesize=60;

procmeansdata=fdata.Gradenoprint;

classStatusYear;

varfinalgrade;

outputout=sumstatmean=AverageGrade

idgroup（max（score）obsout（name）=BestScore）/wayslevels;

run;

procprintdata=sumstatnoobs;

title1'AverageUndergraduateandGraduateCourseGrades';

title2'ForTwoYears';

run;

例16.27计算并存贮多个变量的统计量。

optionsnodatepageno=1linesize=80pagesize=60;

procmeansdata=fdata.Gradenoprintdescend;

classStatusYear;

varScoreFinalGrade;

outputout=Sumdata（where=（status='1'or_type_=0））

mean=median（finalgrade）=MedianGrade;

run;

procprintdata=Sumdata;

title'ExamandCourseGradesforUndergraduatesOnly';

title2'andforAllStudents';

run;

例16.28计算带有缺失值CLASS变量时的输出统计量。

optionsnodatepageno=1linesize=80pagesize=60;

procmeansdata=fdata.cakechartypenwaynoprint;

classflavor/order=freqascending;

classlayers/missing;

varTasteScore;

outputout=cakestatmax=HighScore;

run;

procprintdata=cakestat;

title'MaximumTasteScoreforFlavorandCakeLayers';

run;

例16.29利用输出统计量判别极端值。

optionsnodatepageno=1linesize=80pagesize=60;

procmeansdata=fdata.Charitynmeanrange;

classSchoolYear;

varMoneyRaisedHoursVolunteered;

outputout=Prizemaxid（MoneyRaised（name）

hoursVolunteered（name））=MostCashMostTime

max=;

title'SummaryofVolunteerWorkbySchoolandYear';

run;

procprintdata=Prize;

title'BestResults:

MostMoneyRaisedandMostHoursWorked';

run;

例16.30应用单变量过程进行多变量分析。

optionsnodatepageno=1linesize=80pagesize=72;

procunivariatedata=fdata.statepop;

varcitypop_90citypop_80;

title'UnitedStatesCensusofPopulationandHousing';

run;

例16.31规定四舍五入单位并识别极端值。

optionsnodatepageno=1linesize=80pagesize=68;

procunivariatedata=fdata.statepopfreqround=1nextrobs=2nextrval=4;

varcitypop_90;

idregionstate;

title'UnitedStatesCensusofPopulationandHousing';

run;

例16.32计算稳健估计量。

optionsnodatepageno=1linesize=80pagesize=72;

procunivariatedata=fdata.statepoprobustscaletrimmed=6.25winsorized=.1;

varcitypop_90;

title'UnitedStates1990CensusofPopulationandHousing';

run;

例16.33配对数据的符号检验。

optionsnodatepageno=1linesize=80pagesize=60;

procunivariatedata=fdata.score1loccountmodes

alpha=.01cibasic（alpha=.05）cipctldf;

varscorechange;

labelscorechange='Change