Single Adult 1

Seriamente Yourdreamslut Ro Gb Cs Your Dream Slut µÚ4Õ SASµÄ»ù±¾Í³¼Æ·ÖÎö¹¦ÄÜ

Seriamente Yourdreamslut Ro Gb Cs Your Dream Slut

Seriamente Seriamente à Yourdreamslut P Yourdreamslut R Seriamente W Seriamente Y Seriamente ø Seriamente Ð Yourdreamslut ÇÎ0ý Yourdreamslut wsearch%searchEsearchD1%C% Seriamente 6searchDsearch%3 Seriamente C0%search8C Seriamente % Seriamente E Seriamente D Seriamente %2searchEsearch% Seriamente 41A Yourdreamslut %BRAZZERS%B3%F6%C6%B74 Seriamente B Seriamente %search7 Seriamente F8aSe1I Seriamente S Yourdreamslut x Yourdreamslut 2 Yourdreamslut B Seriamente a Yourdreamslut dosearch2 Seriamente 8%C3%AB%D4%F3%B6%AB%B5%C4%C7%E9%C8%CB%D7%D4%B0%D7%C2%BC% Seriamente E Seriamente o Seriamente ¥searchò Seriamente Ø Seriamente ½ Yourdreamslut î Seriamente Ö Yourdreamslut ö

µ±·½²î·ÖÎöµÄÕý̬·Ö²¼¼Ù¶¨»ò·½²îÏàµÈ¼Ù¶¨²»ÄÜÂú×ãʱ£¬¶Ôµ¥ÒòËØÎÊÌ⣬¿ÉÒÔʹ Ó÷DzÎÊý·½²î·ÖÎöµÄKruskal-Wallis¼ìÑé·½·¨¡£ÕâÖÖ¼ìÑé²»ÒªÇó¹Û²âÀ´×ÔÕý̬·Ö²¼×ÜÌ壬²» ÒªÇó¸÷×éµÄ·½²îÏàµÈ£¬ÉõÖÁÖ¸±ê¿ÉÒÔÊÇÓÐÐò±äÁ¿£¨±äÁ¿È¡ÖµÖ»ÓдóС֮·Ö¶øûÓвî¾àµÄ¸ÅÄî £¬±ÈÈçÄ¥ËðÁ¿¿ÉÒÔ·ÖΪ´ó¡¢ÖС¢Ð¡Èýµµ£¬µÃ²¡µÄ³Ì¶È¿ÉÒÔ·ÖΪÖØ¡¢Çá¡¢ÎÞ£¬µÈµÈ£©¡£

NPAR1WAY¹ý³ÌµÄµ÷ÓÃÓëANOVA¹ý³Ì²»Í¬£¬ÒòΪËüÊǵ¥ÒòËØ·½²î·ÖÎö¹ý³Ì£¬ËùÒÔÖ»ÒªÓÃCLASS Óï¾ä¸ø³ö·ÖÀà±äÁ¿£¨ÒòËØ£©£¬ÓÃVARÓï¾ä¸ø³öÖ¸±ê¾Í¿ÉÒÔÁË£¬Ò»°ã¸ñʽΪ£º

			PROC NPAR1WAY DATA=Êý¾Ý¼¯ WILCOXON;
				CLASS ÒòËØ;
				VAR Ö¸±ê;
			RUN;

×¢ÒâÕâÑùµÄÓï¾ä¸ñʽÓë1.4.2ÖÐÁ½¶ÀÁ¢Ñù±¾±È½ÏµÄ×ö·¨ÍêÈ«Ïàͬ¡£µ±¡°ÒòËØ¡±ÓÐÁ½¸öˮƽ ʱ£¬Ö´ÐÐWilcoxonÖȺͼìÑ飬¶à¸öˮƽʱִÐÐKruskal-Wallis¼ìÑé¡£

±ÈÈ磬ΪÁË·ÖÎöÉÏÃæµÄ½ººÏ°åÀý×ÓÖи÷ÅÆ×ÓµÄÄÍÄ¥ÐÔÓÐÎÞÏÔÖø²îÒ죬ȡ¶¨0.10µÄ¼ìÑé ˮƽ£¬¿ÉÒÔÓÃÈçϵÄNPAR1WAY¹ý³Ì£º

proc npar1way data=sasuser.veneer wilcoxon;
 class brand;
 var wear;
run;

  N P A R 1 W A Y P R O C E D U R E
 Wilcoxon Scores (Rank Sums) for Variable WEAR
  Classified by Variable BRAND

 Sum of  Expected Std Dev Mean
 BRAND  N Scores  Under H0  Under H0  Score

 ACME 4 40.0 42.0 10.4830691 10.0000000
 CHAMP  4 44.0 42.0 10.4830691 11.0000000
 AJAX 4 12.0 42.0 10.4830691  3.0000000
 TUFFY  4 69.0 42.0 10.4830691 17.2500000
 XTRA 4 45.0 42.0 10.4830691 11.2500000

 Average Scores Were Used for Ties
¡¡
 Kruskal-Wallis Test (Chi-Square Approximation)
 CHISQ = 11.982  DF = 4  Prob > CHISQ = 0.0175

½á¹û·ÖΪÁ½¸ö²¿·Ö£¬µÚÒ»²¿·ÖÊǸ÷×éµÄÖȺ͵ÄÇé¿ö£¬°üÀ¨¹Û²â¸öÊý£¨N£©¡¢ÖȺͣ¨Sum of Scores£©¡¢ÔÚ¸÷×éÎÞÏÔÖø²îÒìµÄÁã¼ÙÉèϵÄÆÚÍûÖȺͣ¨Expected Under H0£©¡¢ÔÚÁã¼ÙÉèÏ µÄ±ê×¼²î£¨Std Dev Under H0£©¡¢Æ½¾ùÖȺͣ¨Mean Score£¬ÎªÖȺͳýÒÔ×éÄÚ¹Û²âÊý£©¡£ÏÂÃæ µÄ¡°Average Scores Were Used for Ties¡±ÊÇ˵µ±Ãû´ÎÏàͬʱ£¨ÈçÁ½¸öµÚ2£©ÓÃÃû´ÎµÄƽ¾ù Öµ£¨(2+3)/2=2.5£©¡£µÚ¶þ²¿·ÖΪKruskal-Wallis¼ìÑéµÄ½á¹û£¬°üÀ¨½üËÆµÄ Í³¼ÆÁ¿£¬×ÔÓɶȣ¬¼ìÑéµÄpÖµ£¨Prob > CHISQ £©¡£ÏÖÔÚpÖµ0.0175СÓÚÔ¤¶¨µÄˮƽ0.10ËùÒÔ½áÂÛÊǸ÷ÖÖÅÆ×ӵĽººÏ°åµÄÄÍÄ¥ÐÔÄÜÓÐÏÔÖø²îÒì ¡£×¢Ò⣬Kruskal-Wallis¼ìÑéÊǷDzÎÊý¼ìÑ飬ÔÚͬµÈÌõ¼þÏ·DzÎÊý¼ìÑéÒ»°ã±È²ÎÊý¼ìÑéµÄ¹¦ ЧµÍ£¬ËùÒÔÕâÀïµÄpÖµ0.0175±ÈÓÃANOVA¹ý³ÌµÃµ½µÄpÖµ0.0017Òª´ó¡£

¶àÖرȽÏ

·½²î·ÖÎöÖ»¼ìÑé¸÷×éÊÇ·ñûÓÐÈκÎÁ½Á½Ö®¼äµÄ²îÒ죬µ«²»¼ìÑéµ½µ×ÊÇÄÄÁ½×éÖ®¼äÓÐÏÔ Öø²îÒì¡£ÔÚÈý¸ö»ò¶à¸ö×éÖ®¼ä½øÐÐÁ½¸ö»ò¶à¸ö±È½ÏµÄ¼ìÑé½Ð×ö¶àÖرȽϡ£¶àÖرȽÏÔÚͳ¼Æѧ ÖÐûÓÐÒ»¸ö¹«ÈϵĽâ¾ö·½·¨£¬¶øÊÇÌṩÁËÈô¸ÉÖÖ¼ìÑé·½·¨¡£ÒòΪ¶àÖرȽÏÒª½øÐв»Ö»Ò»´ÎµÄ ±È½Ï£¬ËùÒÔÔÚ¶àÖرȽϵļìÑéˮƽÓÐÁ½ÖÖ£º×Ü´íÎóÂÊ£¨experimentwise error rate£©ºÍµ¥´Î ±È½Ï´íÎóÂÊ¡£×Ü´íÎóÂÊÊÇÖ¸ËùÓбȽϣ¨±ÈÈ磬Îå¸ö×éÁ½Á½Ö®¼ä±È½ÏÓÐ10´Î£©µÄ×ܵÚÒ»Àà´íÎó ¸ÅÂÊ£¬µ¥´Î±È½Ï´íÎóÂÊÊÇָÿһ´Î±È½ÏµÄµÚÒ»Àà´íÎó¸ÅÂÊ¡£ÏÔÈ»£¬×Ü´íÎóÂÊÒª±Èµ¥´Î±È½Ï´í ÎóÂʸߡ£

ÔÚANOVA¹ý³ÌÖÐʹÓÃMEANSÓï¾ä¿ÉÒÔ½øÐжàÖرȽϡ£¸ñʽÈçÏ£º

		MEANS ÒòËØ / Ñ¡Ïî;

Èç¹û²»Ê¹ÓÃÑ¡ÏÔòÖ»¶ÔÒòËصĸ÷ˮƽ¼ÆËãÖ¸±êµÄƽ¾ùÖµºÍ±ê×¼²î£¬±ÈÈ磺

proc anova data=sasuser.veneer;
 class brand;
 model wear=brand;
 means brand;
run;

 Level of  -------------WEAR------------  
 BRAND  N  Mean  SD  

 ACME  4 2.32500000  0.17078251  
 AJAX  4 2.05000000  0.12909944  
 CHAMP  4 2.37500000  0.17078251  
 TUFFY  4 2.60000000  0.14142136  
 XTRA  4 2.37500000  0.09574271  

Òª½øÐÐÁ½Á½±È½Ï£¬ÓжàÖÖ·½·¨£¬¿ÉÒÔÔÚMEANSÓï¾äµÄÑ¡ÏîÖÐÖ¸¶¨¼ìÑé·½·¨¡£

Öظ´t¼ìÑéµÄÏë·¨ºÜ¼òµ¥£ºÔÚÊʵ±µÄ¼ìÑéˮƽ϶ÔÁ½×éÖ®¼ä½øÐÐÁ½Ñù±¾t¼ìÑé²¢¶Ô ËùÓÐ×éÁ½Á½Ö®¼ä¼ìÑé¡£¿ØÖƵÄÊÇÿ´Î±È½ÏµÄµÚÒ»Àà´íÎó¸ÅÂÊ¡£È±Ê¡Ê¹ÓÃ0.05ˮƽ¡£×¢ÒâÕâÑù ¼ìÑéµÄ×Ü´íÎóÂʽ«´ó´ó¸ßÓÚÿ´Î±È½ÏµÄ´íÎóÂÊ¡£±ÈÈ磬ÔÚÉÏÃæ³ÌÐòºó¼ÓÈ루ANOVAÊǽ»»¥Ê½¹ý ³Ì£©

 means brand / t;
run;

 T tests (LSD) for variable: WEAR 

  NOTE: This test controls the type I comparisonwise error rate not the 
 experimentwise error rate. 

  Alpha= 0.05 df= 15 MSE= 0.020833  
 Critical Value of T= 2.13  
 Least Significant Difference= 0.2175 

  Means with the same letter are not significantly different. 

  T Grouping  Mean  N BRAND  

 A 2.6000  4 TUFFY  

 B 2.3750  4 XTRA 
 B 
 B 2.3750  4 CHAMP  
 B 
 B 2.3250  4 ACME 

 C 2.0500  4 AJAX 

½á¹ûÏÈ˵Ã÷Á˼ìÑéµÄÖ¸±êÊDZäÁ¿WEAR£¬È»ºó˵Ã÷ÁËÕâÖÖ¼ìÑé¿ØÖƵ¥´Î±È½ÏµÄµÚÒ»Àà´íÎó¸Å Âʶø²»ÊÇ×ܵĵÚÒ»Àà´íÎó¸ÅÂÊ¡£ÏÂÃæ¸ø³öÁ˼ìÑéµÄһЩָ±ê£¬±ÈÈçˮƽ£¨Alpha£©Îª0.05£¨¿Ø ÖƵ¥´Î±È½ÏµÄµÚÒ»Àà´íÎó¸ÅÂÊ£©£¬×ÔÓɶȣ¨df£©Îª15£¬Îó²îµÄ¾ù·½£¨MSE£¬ÊÇ·½²î·ÖÎö±íÖÐÎó ²îµÄ¾ù·½£©Îª0.020833£¬Á½Ñù±¾t¼ìÑéµÄtͳ¼ÆÁ¿µÄÁÙ½çÖµ£¨Critical Value of T£©Îª2.13£¬ Èç¹ûÁ½Ñù±¾t¼ìÑéµÄtͳ¼ÆÁ¿Öµ¾ø¶ÔÖµ³¬¹ýÁÙ½çÖµÔòÈÏΪÁ½×éÓÐÏÔÖø²îÒ죬»òÕߵȼ۵أ¬Èç¹û Á½×éµÄ¾ùÖµÖ®²î¾ø¶ÔÖµ´óÓÚ×îСÏÔÖø²î±ð£¨Least Significant Difference£©0.2175Ò²ÊÇÓÐ ÏÔÖø²îÒì¡£ÏÂÃæÁгöÁ˼ìÑéµÄ½á¹û£¬°ÑÒòËظ÷ˮƽµÄÖ¸±êƽ¾ùÖµÓÉ´óµ½Ð¡ÅÅÁУ¬È»ºó°ÑÁ½Á½ ±È½ÏµÄ½á¹ûÓõÚÒ»ÁеÄ×ÖĸÀ´±íʾ£¬×ÖĸÏàͬµÄˮƽûÓÐÏÔÖø²îÒ죬×Öĸ²»Í¬µÄˮƽÓÐÏÔÖø ²îÒì¡£ËùÒÔÎÒÃÇ¿´µ½£¬Öظ´t¼ìÑéµÄ½á¹û°ÑÎåÖÖÅÆ×Ó·Ö³ÉÁËA¡¢B¡¢CÈý¸ö×飬TUFFYµ¥¶ÀÊÇÒ»×é £¬ËüµÄÄ¥ËðÁ¿×î´ó£»XTRA¡¢CHAMP¡¢ACMEÊÇÒ»×飬ÕâÈýÖÖÁ½Á½Ö®¼äûÓÐÏÔÖø²îÒ죻AJAXµ¥¶ÀÊÇ Ò»×飬ÆäÄ¥ËðÁ¿×îС¡£

Bonferroni t¼ìÑéͨ¹ý°Ñÿ´Î±È½ÏµÄ´íÎóÂÊÈ¡µÃºÜСÀ´¿ØÖÆ×ÜÎó²îÂÊ¡£±ÈÈ磬¹² ÓÐ10´Î±È½Ïʱ£¬°Ñÿ´Î±È½ÏµÄ´íÎóÂÊ¿ØÖÆÔÚ0.005¾Í¿ÉÒÔ±£Ö¤×Ü´íÎóÂʲ»³¬¹ý0.05£¬µ«ÊÇ£¬Õâ ÑùµÃµ½µÄʵ¼Ê×ܵÚÒ»Àà´íÎóÂÊ¿ÉÄÜÒª±ÈÔ¤¶¨µÄˮƽСµÃ¶à¡£ÔÚMEANSÓï¾äÖÐʹÓÃBONÓï¾ä¿ÉÒÔ Ö´ÐÐBonferroni t¼ìÑ飬ȱʡ×Ü´íÎóÂÊ¿ØÖÆˮƽΪ0.05¡£¶ÔÉÏÃæµÄ½ººÏ°åÊý¾ÝÔö¼ÓÈçÏÂÓï¾ä £º

 means brand / bon;
run;

 Bonferroni (Dunn) T tests for variable: WEAR  

  NOTE: This test controls the type I experimentwise error rate, but
 generally has a higher type II error rate than REGWQ. 

  Alpha= 0.05 df= 15 MSE= 0.020833  
 Critical Value of T= 3.29  
 Minimum Significant Difference= 0.3354 

  Means with the same letter are not significantly different. 

 Bon Grouping  Mean  N BRAND  

 A 2.6000  4 TUFFY  
 A  
 B  A 2.3750  4 XTRA  
 B  A  
 B  A 2.3750  4 CHAMP  
 B  A  
 B  A 2.3250  4 ACME  
 B  
 B 2.0500  4 AJAX  

½á¹ûÏÈ˵Ã÷Á˼ìÑéÀàÐͺÍÖ¸±ê£¨±äÁ¿WEAR£©£¬È»ºó˵Ã÷Á˼ìÑé¿ØÖÆ×ܵÚÒ»Àà´íÎóÂÊ£¬µ«Ò» °ã±ÈREGWQ·½·¨µÄµÚ¶þÀà´íÎó¸ÅÂʸߣ¨¼ìÑ鹦Ч½ÏµÍ£©¡£ÏÂÃæ¸ø³öÁ˼¸¸ö¼ìÑéÓõÄÖµ¡£×îºó¸ø ³öÁËBonferroni t¼ìÑéµÄ½á¹û£¬ÓÐÏàͬ·Ö×é×ÖĸµÄÒòËØˮƽ¼äÎÞÏÔÖø²îÒ죬·ñÔòÓÐÏÔÖø²îÒì ¡£ÎÒÃÇ¿´µ½£¬TUFFYÓëXTRA¡¢CHAMP¡¢ACMEûÓÐÏÔÖø²îÒ죬ÓëAJAXÓÐÏÔÖø²îÒ죻XTRA¡¢CHAMP¡¢ACME Á½Á½Ö®¼äûÓÐÏÔÖø²îÒ죬¶øÇÒÓëÆäËüÁ½¸öÒ²¶¼Ã»ÓÐÏÔÖø²îÒ죻AJAXÓëTUFFYÓÐÏÔÖø²îÒ죬ÓëÆä ËüÈý¸öûÓÐÏÔÖø²îÒì¡£Æä·Ö×éÊÇÓн»²æµÄ¡£

ÓÃBonferroni t¼ìÑé¿ØÖÆ×Ü´íÎóÂʹýÓÚ±£ÊØ£¬¹¦Ð§½ÏµÍ£¬²»Ò×·¢ÏÖʵ¼Ê´æÔÚµÄÏÔ Öø²îÒì¡£REGWQ·½·¨¿ÉÒÔ¿ØÖÆ×Ü´íÎóÂʲ¢ÇÒÒ»°ã±ÈBonferroni t¼ìÑéÒªºÃ¡£ÕâÖÖ·½·¨Ö´Ðжà½× ¶ÎµÄ¼ìÑ飬Ëü¶ÔÒòËØˮƽµÄ¸÷ÖÖ×Ó¼¯½øÐмìÑé¡£ÔÚMEANSÓï¾äÖÐÓÃREGWQÑ¡Ïî¿ÉÒÔ½øÐÐREGWQ¼ì Ñé¡£ÀýÈ磬ÔÚÇ°ÃæµÄÀý×ÓºóÔÙÔËÐÐ

 means brand/ regwq;
run;

  Ryan-Einot-Gabriel-Welsch Multiple Range Test for variable: WEAR  
  
  NOTE: This test controls the type I experimentwise error rate. 
  
  Alpha= 0.05 df= 15 MSE= 0.020833  
  
 Number of Means 2 3 4 5  
 Critical Range  0.264828 0.2917322 0.2941581  0.31516  
  
  Means with the same letter are not significantly different. 
  
 REGWQ Grouping  Mean  N BRAND 
  
 A 2.6000  4 TUFFY 
 A  
 A 2.3750  4 XTRA  
 A  
 A 2.3750  4 CHAMP 
 A  
 A 2.3250  4 ACME  
  
 B 2.0500  4 AJAX  

¿É¼ûËü±ÈBonferroni·½·¨·¢ÏÖÁ˽϶àµÄÏÔÖø²îÒ죬³ýÁËTUFFYºÍAJAXÈÔÓÐÏÔÖø²îÒìÒÔÍ⣬ »¹·¢ÏÖXTRA¡¢CHAMP¡¢ACMEÒ²¶¼ÓëAJAXÓÐÏÔÖø²îÒì¡£

MEANSÓï¾äµÄÑ¡Ïî¿ÉÒÔͬʱʹÓá£ÔÚMEANSÓï¾äÖпÉÒÔÓÃALPHA=ˮƽֵÀ´Ö¸¶¨¼ìÑéµÄË® ƽ¡£ANOVA¹ý³ÌÖл¹ÌṩÁËÆäËüµÄ¶àÖرȽϷ½·¨£¬Çë×Ô¼º²Î¿¼ÓйØ×ÊÁÏ¡£

¶àÒòËØ·½²î·ÖÎö

SASÌṩÁËÈô¸É¸ö·½²î·ÖÎö¹ý³Ì£¬¿ÉÒÔ¿¼ÂǶà¸öÒòËØ¡¢Óн»»¥×÷Óá¢ÓÐǶÌ×µÈÇé¿öµÄ ·½²î·ÖÎö¡£ÓÃGLM¹ý³Ì»¹¿ÉÒÔÓÃÒ»°ãÏßÐÔÄ£ÐÍÀ´´¦Àí·½²î·ÖÎöÎÊÌâ¡£ÔÚÕâÀïÎÒÃÇÖ»½éÉÜÈçºÎÓÃANOVA ¹ý³Ì½øÐоùºâÉè¼ÆµÄ¶àÒòËØ·½²î·ÖÎö¡£

ÀýÈ磬ΪÁËÌá¸ßÒ»ÖÖÏ𽺵Ķ¨Ç¿£¬¿¼ÂÇÈýÖÖ²»Í¬µÄ´Ù½ø¼Á£¨ÒòËØA£©¡¢ËÄÖÖ²»Í¬·ÖÁ¿µÄÑõ »¯Ð¿£¨ÒòËØB£©¶Ô¶¨Ç¿µÄÓ°Ï죬¶ÔÅä·½µÄÿÖÖ×éºÏÖظ´ÊÔÑéÁ½´Î£¬×ܹ²ÊÔÑéÁË24´Î£¬µÃµ½ÈçÏ µÄ½á¹û£º

±í¸ñ 1 Ïð½ºÅä·½ÊÔÑéÊý¾Ý

ÎÒÃÇÊ×ÏÈ°ÑÊý¾ÝÊäÈëΪSASÊý¾Ý¼¯¡£ÊäÈëµÄ°ì·¨¿ÉÒÔÊÇÖ±½ÓÊäÈë¸÷¸ö¹Û²â£¬ÀýÈ磺

data rubber;
 input A B STREN;
 cards;
1 1 31
1 1 33
1 2 34
1 2 36
1 3 35
1 3 36
1 4 39
1 4 38
2 1 33
¡­¡­¡­¡­¡­
;
run;