�I�[�v���\�[�X�̓��v��������E���́uR�v���g���Ď��H�I�ȓ��v��͂�g�ɕt���悤�B�A��2��ڂ̍���́A���v�I����Ƃ͉����Ƃ����������肵�����������ɁA��͂�R�Ŏ��ۂɌv�Z�����Ă݂܂��B
�@��������v��������́uR�v���g���āA��������̓I�ȃf�[�^����v�Z�����Ă݂����Ǝv���܂��B�����������v�I����Ƃ͉����A�ǂ��������Ƃ�����Ă���̂��Ƃ����T���ƁA���v�I����ɑ���ᔻ�ɂ��Ă����Љ�����Ǝv���܂��B
�@�܂����O���̂����炢�ł��B
�@�O��̏I���̂ق����A�������i���ނ����j�ł���uJ���[�K�[�̌��ʏo�������z�͓��{�l�̌��ʏo�������z�Ɠ����ł���v���ǂ����v�w�I�Ɍ������邽�߂ɁA�J�C��挟������܂����B
�@�܂��͑O��̌��ʂ��Č����邱�Ƃ���n�߂܂��傤�B�K�v�ȃf�[�^��ǂݍ���ŁA�ϐ���ݒ肵�܂��B���������łɈȉ��̃R�[�h�𑖂点�Ă���ꍇ�́A�������Ƃ�2�x����K�v�͂���܂���B
> conn <- url("http://euler.bakfoo.com/public/jleagers.rda") > load(conn) > conn <- url("http://euler.bakfoo.com/public/japanpop.rda") > load(conn) > slice <- 7:10 > japrecent <- apply(data.japan[slice,], 2, sum)[-c(1,2)]/length(slice) > ls() [1] "conn" "data.japan" "data.jleagers" "japrecent" [5] "slice"
�@data.jleagers�ɂ�J���[�K�[�̌��ʏo�������z�������Ă��āAdata.japan�ɂ͓��{�l�̌��ʏo�������z�������Ă��܂��Bjaprecent�ɂ�1975�N����1990�N�܂ł̓��{�l�̌��ʂ̏o�����������Ă��܂��B�����āAR�ɂ����āA�J�C��挟������s����菇�͎��̂Ƃ���ł����B
> prob <- japrecent/sum(japrecent) > chisq.test(data.jleagers,p=prob) Chi-squared test for given probabilities data: data.jleagers X-squared = 89.5635, df = 11, p-value = 2.03e-14
�@���̌��ʁAp-value���\���ɏ������Ƃ������ƂŁA�A�������ł���uJ���[�K�[�̌��ʏo�������z�͓��{�l�̌��ʏo�������z�Ɠ����ł���v�Ƃ��������͍̗p�ł��Ȃ��i���m�ɂ͍̗p���Â炢�j�A���p�ł���Ƃ������Ƃł����B���̈�A�̌���́A�J�C���K���x����ƌĂ�܂��B�W�{�f�[�^�̕��z�i���̏ꍇJ���[�K�[�̌��ʏo�������z�j���A���ꂪ��������W�c�̕��z�i���̏ꍇ�A���{�l�̌��ʏo�������z�j�Ɠ����ł��邩�ǂ���������������@�ł��B
�@�����œ��R���̂悤�ȋ^�₪�����Ă���Ǝv���܂��B�Ȃ��A�������ׂ�̂��H �J�C���Ƃ����͉̂����H p-value�͂ǂ������v�Z�ɂ�蓱�������ꂽ�̂��H�����āA���̃J�C��挟����܂��v�I����͊��S�����̃c�[���Ȃ̂��H
�@���̈�A�̌v�Z�ɂ��Ă̌����������Ő������܂��B
�@���̐��w�藝�����藧���Ă��邱�Ƃ��A�J�C���K���x���肪�������鍪���ɂȂ��Ă��܂��B
�藝
�@�J�e�S������k�A����J�e�S��i�ɊY������ϑ��x����ni���̃J�e�S���̊��ғx����ei�Ƃ���Ƃ��ɁA�ȉ��̊m���ϐ�x2
�̕��z���A�T���v�������\���ɑ傫�������Ƃ��Ɏ��R�xk-1�̃J�C��敪�z�ɋ߂Â��B
�@�����Ȃ萔�����o���Ă��܂��܂������S�z���Ȃ��ł��������B���̌シ���ɈӖ���������܂��B
�@�����Łu�J�e�S���v�Ƃ́A�T���v�����Ƃ�Ƃ��̋敪���̂��Ƃł��B�Ⴆ���܂܂Řb��ɂ��Ă���J���[�K�[�̌��ʏo�����̏ꍇ�Adata.jleagers�ɂ���Ajan,�Afeb�Amar�c�c�A�Ƃ����a�����ɋ敪�������ďo�����𐔂��܂��̂ŁA�a�������J�e�S���ɂȂ�܂��B
> data.jleagers jan feb mar apr may jun jul aug sep oct nov dec [1,] 58 52 60 125 137 110 101 98 108 88 67 69
�@���̏����ł́A��L�̒藝�̏ؖ��ɂ܂ŗ�������܂���B����������ǎ҂́A�荠�Ȑ������v�w�̋��ȏ����������Ă��������B�M�҂��悭���p����̂́A�w�����ォ�炨���b�ɂȂ��Ă���ȉ��̖{�ł��B
�@�ؖ��͂��Ȃ�����ɁA�藝�̈Ӗ����Ȍ��ɐ����������Ǝv���܂��B
�@�܂��A�藝�ɂ���x2�Ƃ́A�J�e�S��i���Ƃ̊ϑ��x���Ɗ��ғx���̍��̑��Βl�A
���A���ׂẴJ�e�S���ɂ��đ������킹�����̂ł��B�܂�A����x2�l�͊ϑ��������z���A���҂��Ă��镪�z�ɔ�ׂĂǂꂾ���ގ����Ă��邩��\���l�ł��B
�@���̂��߂ɁAx2���A���R�xk-1�̃J�C��敪�z�����肵���Ƃ��ɁA�J�C��敪�z�̊m���Ƃ��ď\���ɏ����Ȓl���Ƃ邩�A�܂�ϑ��������z�Ɗ��҂��Ă��镪�z���ǂ̒��x�̊m���œ����ł��邩�ǂ����A�Ƃ������Ƃ������̊�ɂȂ�܂��B
�@�����ŃJ�C��敪�z�Ƃ́A����m�����x�������A���m�����z�̂��Ƃł��B����܂Ŋm���ϐ��̐��������K���z�̐������m�����x���̐��������Ă��Ȃ��̂ŁA�����Ȃ艽�̂�������Ƃ����J���W�ɂȂ�̂ł͂Ȃ����Ǝv���܂��B�����ŊȒP�ɐ������܂��ƁA�m���ⓝ�v�̐��E�ł́A������ϐ����m��I�Ȓl�����킯�łȂ��A����m���́u���v�����̂ł��B������m���ϐ��Ƃ����܂��B
�@�Ⴆ�������H��100km/h�ʼn^�]����Ƃ��܂��B���̂Ƃ��Ƀh���C�o�[�͂ł��邾��100km/h��ۂƂ��Ƃ��Ă��A���ۂɂ͘H�ʏ⑼�̎����ԂƂ̊W�A�����ĕ���g���l���Ȃǂ̊O���v���ɂ���āA100km/h�𒆐S�ɁA���镝���������l�ɂȂ�͂��ł��B
�@���̍������H�𑖂�Ԃ̑��x�̂悤�ɁA�u���v���������ϐ����m���ϐ��ŁA���̊m���ϐ��������̌`���m�����x���Ƃ����܂��B�m�����x���́u���x�v�Ƃ������t�̂Ƃ���A�ϕ�����ƑS�m���A�܂�u�m��=1�v�ɂȂ�悤�Ȋ��ł��B
�@�藝�ɂ���x2�Ƃ����m���ϐ������A�J�C��敪�z�̊m�����x���Ƃ́A���R�x���p�����[�^�Ƃ��āA�ȉ��̂悤�Ȍ`���������̂ł��B
�@���̐}�������ɂ́AR�ɂ���g���݂̃J�C���m�����x��dchisq���������`���curve�Ƒg�݂킹�āA���ɂ̂悤�ɂ��܂��B
> curve(dchisq(x,1), from=0, to=10)
�@���Ȃ݂ɁA�J�C���m�����x���̎��R�x��1����12�܂œ������Ǝ��̂悤�ɂȂ�܂��B
Copyright © ITmedia, Inc. All Rights Reserved.