保護中: Principal component analysis (PCA)

Data is here. password: ecoepi_2008 Now we are going to divide households into 4 categories using Asset factor Score defined by PCA. 1. Advanced (simple) program. *--- Start use house_pcf, clear pca car-toiletcat3, comp(1) ** Asset factor score calculation using comp1 predict pc1, score ** quintile category xtile pc1_qtl = pc1 , nquantiles (4) 2. Basic program: calculate one by one *--- Start use house_pcf, clear ** "#delimit ; " means that stata program reads through until #delimit cr regardless line change. #delimit ; pca car fridge tv stereo ironbox cell motorcycle fan radio bike sofa xtable electricity bed livestock poultry sewing gas wallcat floorcat roofcat lightingcat1 lightingcat2 lightingcat3 lightingcat4 cookingcat1 cookingcat2 cookingcat3 landcat1 landcat2 landcat3 phasecat1 phasecat2 phasecat3 phasecat4 toiletcat1 toiletcat2 toiletcat3 ; #delimit cr ** example: but we use only comp1 /* ---------------------------------------------------------------------------------- Variable | Comp1 Comp2 Comp3 Comp4 Comp5 Comp6 Comp7 -------------+---------------------------------------------------------------------- car | 0.0484 0.1166 0.1899 -0.0162 -0.1068 -0.0005 -0.0118 fridge | 0.0630 0.1170 0.2875 0.0072 -0.1014 -0.0140 -0.0538 tv | 0.1114 0.2800 0.1743 -0.0609 -0.0646 0.0230 0.0542 stereo | 0.0378 0.1512 0.1058 -0.0119 -0.1139 0.0122 0.0374 ironbox | 0.0443 0.2753 0.0302 -0.0946 0.1463 0.0456 0.2213 cell | 0.1504 0.2398 -0.0780 -0.1382 0.2594 0.0826 0.0319 motorcycle | 0.0752 0.1501 0.0955 0.0236 -0.0052 0.0136 0.0363 fan | 0.0278 0.0882 0.1593 0.0588 -0.0546 -0.0129 -0.0955 radio | 0.0727 0.2379 -0.1276 -0.0005 0.2104 0.0342 0.1148 bike | -0.0227 0.2422 -0.0949 -0.1491 0.0726 0.0584 0.3133 sofa | -0.0243 -0.2807 0.1077 0.0116 -0.1581 0.0114 0.0554 xtable | -0.0374 0.2102 -0.0843 0.1268 0.1532 -0.0676 -0.0574 electricity | 0.1245 0.1532 0.4449 0.0123 -0.1094 -0.0389 -0.0901 bed | -0.0367 0.1733 -0.0606 0.1104 0.1477 -0.0528 -0.0666 livestock | -0.1738 0.2678 -0.0419 -0.0202 0.1412 0.0390 0.1746 poultry | -0.2278 0.2375 -0.0266 -0.0280 0.0724 0.0042 0.0438 sewing | 0.0157 0.0844 0.0067 -0.0664 -0.0084 0.0273 0.1451 gas | 0.1967 0.0632 0.0852 0.0041 0.1287 -0.0039 0.0878 wallcat | -0.3002 0.0480 0.0537 -0.0203 -0.0097 -0.0141 0.0192 floorcat | -0.3202 -0.0476 0.0832 -0.0209 -0.0423 -0.0220 -0.0318 roofcat | -0.1288 -0.0770 0.1457 -0.0249 0.0550 -0.0223 0.0518 lightingcat1 | 0.1333 0.2375 -0.3911 -0.0137 -0.3880 -0.0530 -0.2810 lightingcat2 | 0.1435 0.1542 0.4534 0.0115 -0.0984 -0.0403 -0.0995 lightingcat3 | -0.1688 -0.2873 0.2728 0.0418 0.4163 0.0488 0.2695 lightingcat4 | 0.0063 0.0550 -0.0142 -0.1126 -0.0389 0.0480 0.1124 cookingcat1 | 0.3583 -0.0767 -0.0426 -0.0480 0.0930 0.0162 0.0649 cookingcat2 | -0.3636 0.0682 0.0014 0.0614 -0.0806 -0.0182 -0.0564 cookingcat3 | 0.0488 0.0294 0.1680 -0.0590 -0.0446 0.0094 -0.0304 landcat1 | -0.3460 0.1957 0.0407 0.0126 -0.0798 -0.0122 -0.0758 landcat2 | 0.3343 -0.1475 -0.0478 0.0218 0.1540 0.0002 0.0770 landcat3 | 0.0683 -0.1206 0.0085 -0.0693 -0.1356 0.0250 0.0075 phasecat1 | 0.0177 -0.0308 -0.0123 -0.5195 0.0827 -0.5759 -0.0527 phasecat2 | 0.0007 -0.0243 -0.0025 -0.2098 -0.1159 0.7744 -0.0302 phasecat3 | -0.0065 0.0078 -0.0263 0.4027 -0.3536 -0.1536 0.5573 phasecat4 | -0.0143 0.0514 0.0427 0.3967 0.3754 0.0344 -0.4673 toiletcat1 | 0.0122 0.0314 0.0504 0.0065 -0.0486 -0.0065 0.0302 toiletcat2 | 0.1522 0.0735 -0.0838 0.4908 0.0246 -0.0765 0.0873 toiletcat3 | 0.0639 0.0626 0.1799 0.0425 -0.0335 -0.0166 -0.0375 ------------------------------------------------------------------------------------ */ ***** Lots of components!! but we use only first one, so we add the option " comp(1)" #delimit ; pca car fridge tv stereo ironbox cell motorcycle fan radio bike sofa xtable electricity bed livestock poultry sewing gas wallcat floorcat roofcat lightingcat1 lightingcat2 lightingcat3 lightingcat4 cookingcat1 cookingcat2 cookingcat3 landcat1 landcat2 landcat3 phasecat1 phasecat2 phasecat3 phasecat4 toiletcat1 toiletcat2 toiletcat3, comp(1) ; #delimit cr * Only comp1 shown. ***後ほど使う平均と分散のために。以下・・・・。 #delimit ; sum car fridge tv stereo ironbox cell motorcycle fan radio bike sofa xtable electricity bed livestock poultry sewing gas wallcat floorcat roofcat lightingcat1 lightingcat2 lightingcat3 lightingcat4 cookingcat1 cookingcat2 cookingcat3 landcat1 landcat2 landcat3 phasecat1 phasecat2 phasecat3 phasecat4 toiletcat1 toiletcat2 toiletcat3 ; #delimit cr ** あと平均分散も必要か。。。面倒だけどそれをもとめる。comp1をもちいる。 **再度 use house_pcf, clear #delimit; order car fridge tv stereo ironbox cell motorcycle fan radio bike sofa xtable electricity bed livestock poultry sewing gas wallcat floorcat roofcat lightingcat1 lightingcat2 lightingcat3 lightingcat4 cookingcat1 cookingcat2 cookingcat3 landcat1 landcat2 landcat3 phasecat1 phasecat2 phasecat3 phasecat4 toiletcat1 toiletcat2 toiletcat3 ; #delimit cr *面倒だが代入。Comp1の値を。 gen car_afs= 0.0484 gen fridge_afs= 0.063 gen tv_afs= 0.1114 gen stereo_afs= 0.0378 gen ironbox_afs= 0.0443 gen cell_afs= 0.1504 gen motorcycle_afs= 0.0752 gen fan_afs= 0.0278 gen radio_afs= 0.0727 gen bike_afs= -0.0227 gen sofa_afs= 0.0243 gen xtable_afs= -0.0374 gen electricity_afs= 0.1245 gen bed_afs= -0.0367 gen livestock_afs= -0.1738 gen poultry_afs= -0.2278 gen sewing_afs= 0.0157 gen gas_afs= 0.1967 gen wallcat_afs= -0.3002 gen floorcat_afs= -0.3202 gen roofcat_afs= -0.1288 gen lightingcat1_afs= 0.1333 gen lightingcat2_afs= 0.1435 gen lightingcat3_afs= -0.1688 gen lightingcat4_afs= 0.0063 gen cookingcat1_afs= 0.3583 gen cookingcat2_afs= -0.3636 gen cookingcat3_afs= 0.0488 gen landcat1_afs= -0.346 gen landcat2_afs= 0.3343 gen landcat3_afs= 0.0683 gen phasecat1_afs= 0.0177 gen phasecat2_afs= 0.0007 gen phasecat3_afs= -0.0065 gen phasecat4_afs= -0.0143 gen toiletcat1_afs= 0.0122 gen toiletcat2_afs= 0.1522 gen toiletcat3_afs= 0.0639 gen car_mean= 0.0031589 gen fridge_mean= 0.0021607 gen tv_mean= 0.0729575 gen stereo_mean= 0.0166662 gen ironbox_mean= 0.2204266 gen cell_mean= 0.3996993 gen motorcycle_mean= 0.0199768 gen fan_mean= 0.0069874 gen radio_mean= 0.7710571 gen bike_mean= 0.2815446 gen sofa_mean= 0.6721842 gen xtable_mean= 0.9374795 gen electricity_mean= 0.0145561 gen bed_mean= 0.9452124 gen livestock_mean= 0.3873165 gen poultry_mean= 0.6353896 gen sewing_mean= 0.0152763 gen gas_mean= 0.0467261 gen wallcat_mean= 0.6090255 gen floorcat_mean= 0.6730036 gen roofcat_mean= 0.1183271 gen lightingcat1_mean= 0.34535 gen lightingcat2_mean= 0.0151502 gen lightingcat3_mean= 0.6039341 gen lightingcat4_mean= 0.0291082 gen cookingcat1_mean= 0.1901226 gen cookingcat2_mean= 0.783377 gen cookingcat3_mean= 0.020043 gen landcat1_mean= 0.7180325 gen landcat2_mean= 0.2187217 gen landcat3_mean= 0.0567884 gen phasecat1_mean= 0.331156 gen phasecat2_mean= 0.2246824 gen phasecat3_mean= 0.2217765 gen phasecat4_mean= 0.222385 gen toiletcat1_mean= 0.0002111 gen toiletcat2_mean= 0.1383635 gen toiletcat3_mean= 0.0040607 gen car_sd= 0.0561154 gen fridge_sd= 0.0464331 gen tv_sd= 0.2600683 gen stereo_sd= 0.1280183 gen ironbox_sd= 0.4145369 gen cell_sd= 0.4898396 gen motorcycle_sd= 0.1399211 gen fan_sd= 0.0832989 gen radio_sd= 0.4201551 gen bike_sd= 0.4497552 gen sofa_sd= 0.4694203 gen xtable_sd= 0.2421001 gen electricity_sd= 0.1197682 gen bed_sd= 0.2275666 gen livestock_sd= 0.48714 gen poultry_sd= 0.4813238 gen sewing_sd= 0.1226506 gen gas_sd= 0.211053 gen wallcat_sd= 0.4879717 gen floorcat_sd= 0.4691189 gen roofcat_sd= 0.3229971 gen lightingcat1_sd= 0.4754852 gen lightingcat2_sd= 0.122151 gen lightingcat3_sd= 0.4890815 gen lightingcat4_sd= 0.168111 gen cookingcat1_sd= 0.3924002 gen cookingcat2_sd= 0.4119461 gen cookingcat3_sd= 0.1401481 gen landcat1_sd= 0.4499604 gen landcat2_sd= 0.4133819 gen landcat3_sd= 0.2314393 gen phasecat1_sd= 0.470632 gen phasecat2_sd= 0.4173756 gen phasecat3_sd= 0.4154442 gen phasecat4_sd= 0.415851 gen toiletcat1_sd= 0.0145282 gen toiletcat2_sd= 0.3452833 gen toiletcat3_sd= 0.0635949 save house_pcf2, replace use house_pcf2, clear ***ここでAsset factor scoreを計算 foreach var of varlist car-toiletcat3 { replace `var'=`var'_mean if `var'==. gen `var'_score= (`var'-`var'_mean)*`var'_afs/`var'_sd } egen household_score=rsum(car_score-toiletcat3_score) sum household_score, detail ** quartileにわける gen hscorecat=1 if household_score &lt;=r(p25) replace hscorecat=2 if household_score &gt;r(p25) &amp; household_score&lt;=r(p50) replace hscorecat=3 if household_score &gt;r(p50) &amp; household_score&lt;=r(p75) replace hscorecat=4 if household_score &gt;r(p75) keep houseid phase household_score hscorecat grid subgrid sort houseid phase save house_score, replace

2024年4月
日	月	火	水	木	金	土
	1	2	3	4	5	6
7	8	9	10	11	12	13
14	15	16	17	18	19	20
21	22	23	24	25	26	27
28	29	30