���ǡ������ʤ��Ȥ��θ�������������Last modified: Aug 27, 2009

���㤨�С�2 ����ʿ���ͤκ��θ����Ԥ��Ȥ����ǡ�������ʿ���ͤ�ɸ���к������狼�äƤ��ʤ��Ȥ��ˡ�����ϤǤ��ʤ��������ʤ��뤤�ϡ��ʤ�Ȥ����ơ������Ԥ������ˤȹͤ���ͤ����롣

��2 ����ʿ���ͤκ��θ���ϡ�t.test �ؿ���Ȥ��йԤ��뤬�����Τ���ˤϳΤ��ˡ����ǡ�����ɬ�פǤ��롣

����������ʿ���ͤκ��θ���ʤ����� t ����ˤβ����򸫤�и���Τ����ɬ�פʾ���ϡ��Ʒ��Υǡ�������ʿ���͡�����ʬ���Ρ���� 6 �Ĥξ�������Ǥ褤���Ȥ��狼�롣�ɤ��ʬ�ۤ˽��äƤ��褦��̵�ط��ʤΤǤ��롣�ޤ�������ʬ����ɸ���к�����褹������Ǥ��뤫�顤��ɡ����󷲤ˤĤ��ơ��ǡ�������ʿ���ͤ�ɸ���к����狼�äƤ���С�ʿ���ͤκ��θ���ϹԤ���פȤ������Ȥˤʤ롣

�����Τ���ˡ����Ѥδؿ���񤯤��Ȥ�Ǥ��뤬����ä����Ѥ���ˡ�����롣����ϡ����뷲�Υǡ�����Ǥä��夲����ˡ�Ǥ��롣����ñ�ˤǤä��夲��ΤǤϤʤ����ǡ�������ʿ���͡�ɸ���к������ꤵ�줿��Τ�Ʊ���ˤʤ�ǡ�������ΤǤ��롣

���ǡ������ɤ��ʬ�ۤ˽����������ʤ��Ȥ������Ȥ��˽Ҥ٤������ޤ����������ʬ�ۤ�����ʬ�ۤ˽����ǡ�����������褦��R �ˤϡ������������������ؿ� rnorm �����롣ʿ���� mean��ɸ���к� sd ������ʬ�ۤ˽��� n �ĤΥǡ��� x ����������ˤ� x <- rnorm(n, mean, sd) �Ȥ�������Ǥ褤���������������ʳ��Ǥ� x ��ʿ���ͤ�ɸ���к��ϻ��ꤵ�줿��ΤȤϰۤʤ롣�ʤ�Ȥ��äƤ������פʤΤǤ��뤫�顣���Τ˻��ꤵ�줿ʿ���ͤ�ɸ���к�����������������ˤ�Ŭ������������������������ɸ�ಽ���롣ɸ�ಽ������������ϡ�ʿ���ͤ� 0��ɸ���к��� 1 �ˤʤ�Ȥ����Τ�ï�⤬�ΤäƤ��뤳�ȤǤ��롣��Ū�Ȥ�����������Ϥ���ɸ�ಽ���줿��������ˡ����ꤵ�줿ɸ���к���ݤ������ꤵ�줿ʿ���ͤ�ä���Ф褤�ΤǤ��롣��������С��Ǥ������������ʿ���ͤ�ɸ���к��ϡ����ꤵ�줿ʿ���ͤ�ɸ���к��ˤԤ���Ȱ��פ��롣�ʤԤ���Ȱ��פ���Ȥ����Ȥ��������դ�ɬ�ס��Ĥޤꡤ���ꤵ�줿ʿ���ͤ�ɸ���к��ϡ������줿�Ȥ���ο��ͤǤϤʤ����ݤ��줿���ͤǤ���Ȥ������ȡ������������Ǥϡ������ǡ��������Τ˸��Υǡ�����Ʊ���Ǥ���Ȥϸ����ʤ��Ȥ������ȤǤ����

�㡧
�ǡ������� 24 �ǡ�ʿ���ͤ� 123.45 ɸ���к��� 6.78 ������������롣
> options(digits=20)
> x <- rnorm(24)
> cat("mean =", mean(x), "   S.D.=", sd(x), "\n")
mean = -0.2130160848930967    S.D.= 0.950709790266453 # ����ʤΤ����ۤʤä��ͤˤʤ뤬�����ꤵ�줿�ͤǤʤ����Ȥ����餫
> x <- scale(x) # ɸ�ಽ��Ԥ��ؿ� scale ��Ŭ�Ѥ���
> cat("mean =", mean(x), "   S.D.=", sd(x), "\n")
mean = -2.313868136446147e-17    S.D.= 1 # ʿ���ͤ� -2.313868136446147��10��(-17)��Ȥ������Ȥǡ��ؤ� 0 �Ȥ�������
> x <- x*6.78+123.45 # ɸ�ಽ�ε�
> cat("mean =", mean(x), "   S.D.=", sd(x), "\n")
mean = 123.45    S.D.= 6.78 # ɽ������Ƥ���Ȥ����ʿ���ͤ�ɸ���к�����ĤȤ������ȤǤ���

�����ơ������ޤ����������Ϥ����������Ʒ����Ȥ˾�˽Ҥ٤���ˡ�ǥǡ�����Ǥä��夲�����ΤǤä��夲���ǡ����� t.test �ؿ���ȤäƸ��ꤹ��Ф褤����μ������Ȥ�����¾�Ǥ����ѤǤ������ʴؿ��Ǥ��뤫�顤�ʲ��Τ褦��������Ƥ����Ȥ褤��

gendat1 <- function(n, mean, sd) return(scale(rnorm(n))*sd+mean)
�����Ȥäơ��ʲ��Τ褦�˸����Ԥ����Ȥ��Ǥ��롣
�㡧
�������ǡ���������ʿ���͡���ɸ���к�
A ��    24        123.45       6.78
B ��    35        126.28       5.43

> (A <- gendat1(24, 123.45, 6.78)) # A ���Υǡ���
 [1] 126.4254447481493 119.2677005077077 109.9809264790594 125.2044488778086 118.5402943426712 126.7757599329591
 [7] 133.1987272298903 131.5923794221886 113.4644342955914 119.5887977878966 119.7263229865866 116.3732560196951
[13] 124.9521671141838 129.0740980418013 122.4442664845299 122.3567964723745 133.8439545000280 129.6773031494384
[19] 129.5823721125123 126.7703578980750 110.9939206252493 126.5980379505317 129.2924096537701 117.0758233673015

> (B <- gendat1(35, 126.28, 5.43)) # B ���Υǡ���
 [1] 139.9765197560334 132.1464449725871 136.6031115798909 125.9914781879804 123.5525306533870 130.6480113694447
 [7] 125.2960720595223 124.8400514520517 126.2640875155946 124.8183936737128 124.5685462179979 125.6605823013066
[13] 128.7106647045467 116.3830616911681 120.6424874454661 130.1845143553381 131.3087061533496 124.2051838899657
[19] 122.4152007594897 118.5603630994196 125.6499879613931 130.4590540768779 120.0483665135104 125.7390187089629
[25] 128.3522876051743 126.9468484651417 126.3746312355432 116.9733911064371 123.4602969373899 122.7968444006815
[31] 124.9700778037051 124.2663388696768 140.6244929523544 124.3357595465487 126.0265919783500
���������ʲ��η�������ͤ�Ĺ������ʿ���ͤ�ɸ���к��ϻ��ꤵ�줿��Τ�������
> cat("mean =", mean(A), "   S.D.=", sd(A), "\n")
mean = 123.45    S.D.= 6.78 

> cat("mean =", mean(B), "   S.D.=", sd(B), "\n")
mean = 126.28    S.D.= 5.43 
������ 2 �ĤΡ����ǡ����פ�Ȥäơ�t.test �ؿ��Ǹ����Ԥ����Ȥ��Ǥ��롣
> t.test(A, B)

	Welch Two Sample t-test

data:  gendat1(24, 123.45, 6.78) and gendat1(35, 126.28, 5.43) 
t = -1.7041, df = 42.164, p-value = 0.09571
alternative hypothesis: true difference in means is not equal to 0 
95 percent confidence interval:
 -6.180952255076989  0.520952255076993 
sample estimates:
mean of x mean of y 
   123.45    126.28 
�����Ѥδؿ���Ȥä���̤ϡ��ʲ��Τ褦�ˤʤ롣Ʊ����̤������뤳�Ȥ���ǧ�Ǥ��롣
> my.t.test(24, 123.45, 6.78^2, 35, 126.28, 5.43^2)

	�����������ˡ�ˤ�롤�󷲤�ʿ���ͤκ��θ���

data:  
n1 = 24, mean1 = 123.45, variance1 = 45.9684
n2 = 35, mean2 = 126.28, variance2 = 29.4849 
t = 1.7041, df = 42.164, p-value = 0.09571

��������ˡ�򸫤Ƥ���ȡ�����Ȥ����Τ����Υǡ����Τۤ�θ¤�줿���󤷤��ȤäƤ��ʤ����Ȥ˲���Ƶ��Ť���


�� ľ���Υڡ�������������� E-mail to Shigenobu AOKI

Made with Macintosh