DATA Sample1 /**/ INPUT Price /* */ DATALINES

Size: px

Start display at page:

Download "DATA Sample1 /**/ INPUT Price /* */ DATALINES"

たけなりとみもと
5 years ago
Views:

1 3180, 3599, 3280, 2980, 3500, 3099, 3200, 2980, 3380, 3780, 3199, 2979, 3680, 2780, 2950, 3180, 3200, 3100, 3780, 3200

2 DATA Sample1 /**/ INPUT Price /* */ DATALINES /* TTEST*/ PROC TTEST DATA=Sample1 H0=3000 ALPHA=0.10 /* VAR*/ VAR Price TTEST Price N Price Price t t Pr > t Price /* SAS*/ DATA Stat_for_Sample1 INPUT _STAT_$ Price DATALINES N 20 MEAN STD /* TTEST*/ PROC TTEST DATA=Stat_for_Sample1 H0=3000 ALPHA=0.10 VAR Price

3 /* SAS*/ PROC UNIVARIATE DATA=Sample1 MU0=3000 VAR Price H0: Mu0= p Student t t Pr > t M 5 Pr >= M S 82 Pr >= S DATA Sample2 /* */ INPUT Before After DATALINES /* TTESTt*/ PROC TTEST DATA=Sample2 /* PAIRED2*/ PAIRED After*Before TTEST After - Before After - Before After - Before After - Before N t t Pr > t

4 DATA Sample2_Diff SET Sample2 Difference=Before-After /* */ /* UNIVARIATE*/ PROC UNIVARIATE DATA=Sample2_Diff VAR Difference DATA Sample3 /* SASSample3*/ INPUT Method$ DATALINES A 66 A 64 A 63 A 74 A 60 A 58 A 75 A 48 A 80 A 75 A 51 A 48 A 69 A 45 A 75 B 84 B 90 B 67 B 91 B 72 B 68 B 86 B 50 B 61 B 65 B 96 B 89 B 83 B 78 B 71 B 65 /* BOXPLOT*/ PROC BOXPLOT DATA=Sample3 PLOT *Method / BOXSTYLE=SCHEMATIC /* TTEST*/ PROC TTEST DATA=Sample3 /* CLASS*/ CLASS Method VAR /* VAR*/ Method N A B Diff (1-2) Method A B Diff(1-2) Method A B Diff(1-2)

t Pooled Equal 29 Satterthwaite Unequal 28.9 Folded F 15 14 t Pr > t -2.85 0.0079-2.87 0.0077 F Pr > F 1.29 0.

5 t Pooled Equal 29 Satterthwaite Unequal 28.9 Folded F t Pr > t F Pr > F DATA Stat_for_Sample3 INPUT Method$ _STAT_$ DATALINES A N 15 A MEAN 63.4 A STD B N 16 B MEAN 76 B STD /* TTEST*/ PROC TTEST DATA=Stat_for_Sample3 CLASS Method VAR

6 DATA Sample4 /* */ INPUT Event$ DATALINES Head Tail Head Head Tail Tail Tail Head Tail Tail Tail Tail Head Head Tail Head Tail Tail Head Tail PROC FREQ DATA=Sample4 /* FREQ*/ TABLES Event /BINOMIAL(LEVEL="Head" P=0.5) FREQ Event Head Tail Event = Head % 95% % % H0: = Z Pr < Z Pr > Z = 20 PROC FREQ DATA=Sample4 TABLES Event /BINOMIAL(LEVEL="Head" P=0.5) EXACT BINOMIAL /* EXACT*/

7 H0: = Z Pr < Z Pr > Z Pr <= P = 2 * = 20 DATA Sample5 /* */ INPUT Event$ Frequency DATALINES Head 8 Tail 12 PROC FREQ DATA=Sample5 TABLES Event /BINOMIAL(LEVEL="Head" P=0.5) EXACT BINOMIAL WEIGHT Frequency /* WEIGHT*/ PROC FORMAT /* FORMAT*/ VALUE Approve 1="" 2="" DATA Sample6 INPUT Last_Year This_Year Freq FORMAT Last_Year This_Year Approve. DATALINES PROC FREQ DATA=Sample6 /* AGREE*/ TABLES Last_Year*This_Year / AGREE EXACT MCNEM WEIGHT Freq Last_Year * This_Year McNemar (S) Pr > S Pr >= S

Sex * Answer ----------------------------------------------- 2 2 2 1 1 1 1.0714 1.0747 0.

8 DATA Sample7 /* */ INPUT Sex$ Answer$ Cnt DATALINES Female Yes 18 Female _No 12 Male Yes 14 Male _No 16 PROC FREQ DATA=Sample7 TABLES Sex*Answer / CHISQ /* CHISQ*/ WEIGHT Cnt Sex * Answer Mantel-Haenszel Cramer V

/* MAUTOLOCDISPLAY */ OPTIONS MAUTOLOCDISPLAY %LET x=%trim(abc) 1 OPTIONS MAUTOLOCDISPLAY 2 %LET x=%trim(abc) MAUTOLOCDISPLAY(TRIM): C: Program Files SAS SAS 9.1 core sasmacro trim.

10 /* MAUTOLOCDISPLAY */ OPTIONS MAUTOLOCDISPLAY %LET x=%trim(abc) 1 OPTIONS MAUTOLOCDISPLAY 2 %LET x=%trim(abc) MAUTOLOCDISPLAY(TRIM): C: Program Files SAS SAS 9.1 core sasmacro trim.sas =0 >0 <0 FILENAME myfile "C: myfile.txt" %MACRO chk(filenm) %IF %SYSFUNC(FILEREF(&filenm))=0 %THEN %PUT "FILENAME &filenm " %ELSE %IF %sysfunc(fileref(&filenm))>0 %THEN %PUT "FILENAME &filenm " %ELSE %IF %SYSFUNC(FILEREF(&filenm))<0 %THEN %PUT "FILENAME &filenm " %MEND %chk(myfile) LEFTMARGIN : RIGHTMARGIN : TOPMARGIN : BOTTOMMARGIN : PROC TEMPLATE DEFINE STYLE styles.test PARENT=styles.rtf STYLE BODY FROM DOCUMENT/ LEFTMARGIN=3in/**/ RIGHTMARGIN=2.5in /**/ TOPMARGIN=2.5in/**/ BOTTOMMARGIN=2.5in/**/ END /*TEMPLATE*/ ODS RTF FILE="c: temp test.rtf" STYLE=styles.test

11 OPTIONS NONOTES ODS LISTING CLOSE /* */ < > ODS LISTING /* */ ODS NORESULTS /* GDEVICE0 */ LIBNAME GDEVICE0 'D: sasdev gdev' /* */ PROC CATALOG C=SASHELP.DGDEVICE /* CGMMW6C */ COPY OUT=GDEVICE0.DEVICES SELECT CGMMW6C / ET=DEV QUIT

12 PROC NLP TECH=NONE PHESSIAN PARMS x= , y= f=2*sin(x)*(cos(y)+1) MAX f PROC IML START func(x) f=2*sin(x[1])*(cos(x[2])+1) RETURN(f) FINISH x0={ } CALL NLPFDD(f,g,h,"func",x0)/* f*/ /* gh 2*/ PRINT g h QUIT

H22 BioS (i) I treat1 II treat2 data d1; input group patno treat1 treat2; cards; ; run; I

H22 BioS (i) I treat1 II treat2 data d1; input group patno treat1 treat2; cards; ; run; I H BioS (i) I treat II treat data d; input group patno treat treat; cards; 8 7 4 8 8 5 5 6 ; run; I II sum data d; set d; sum treat + treat; run; sum proc gplot data d; plot sum * group ; symbol c black