STATA: Loop문의 중요성..

* 지금 아래의 코드처럼, 간혹 이렇게 길고 지루한 코드를 짜기 위해서 머리가 아픈 적은 누구나 있겠죠...  그 밑에 축약된 코드문이 있는데,.쭉 한 번 내려가 보도록 하죠, 내려가면서 그 규칙성을 한 번 생각해보면 좋아요!! ...gogo

//28
replace income2031=0 if age_1999==28
//29
replace income2031=0 if age_1999==29
replace income2030=0 if age_1999==29
//30
replace income2031=0 if age_1999==30
replace income2030=0 if age_1999==30
replace income2029=0 if age_1999==30
//31
replace income2031=0 if age_1999==31
replace income2030=0 if age_1999==31
replace income2029=0 if age_1999==31
replace income2028=0 if age_1999==31
//32
replace income2031=0 if age_1999==32
replace income2030=0 if age_1999==32
replace income2029=0 if age_1999==32
replace income2028=0 if age_1999==32
replace income2027=0 if age_1999==32
//33
replace income2031=0 if age_1999==33
replace income2030=0 if age_1999==33
replace income2029=0 if age_1999==33
replace income2028=0 if age_1999==33
replace income2027=0 if age_1999==33
replace income2026=0 if age_1999==33
//34
replace income2031=0 if age_1999==34
replace income2030=0 if age_1999==34
replace income2029=0 if age_1999==34
replace income2028=0 if age_1999==34
replace income2027=0 if age_1999==34
replace income2026=0 if age_1999==34
replace income2025=0 if age_1999==34
//35
replace income2031=0 if age_1999==35
replace income2030=0 if age_1999==35
replace income2029=0 if age_1999==35
replace income2028=0 if age_1999==35
replace income2027=0 if age_1999==35
replace income2026=0 if age_1999==35
replace income2025=0 if age_1999==35
replace income2024=0 if age_1999==35
//36
replace income2031=0 if age_1999==36
replace income2030=0 if age_1999==36
replace income2029=0 if age_1999==36
replace income2028=0 if age_1999==36
replace income2027=0 if age_1999==36
replace income2026=0 if age_1999==36
replace income2025=0 if age_1999==36
replace income2024=0 if age_1999==36
replace income2023=0 if age_1999==36
//37
replace income2031=0 if age_1999==37
replace income2030=0 if age_1999==37
replace income2029=0 if age_1999==37
replace income2028=0 if age_1999==37
replace income2027=0 if age_1999==37
replace income2026=0 if age_1999==37
replace income2025=0 if age_1999==37
replace income2024=0 if age_1999==37
replace income2023=0 if age_1999==37
replace income2022=0 if age_1999==37
//38
replace income2031=0 if age_1999==38
replace income2030=0 if age_1999==38
replace income2029=0 if age_1999==38
replace income2028=0 if age_1999==38
replace income2027=0 if age_1999==38
replace income2026=0 if age_1999==38
replace income2025=0 if age_1999==38
replace income2024=0 if age_1999==38
replace income2023=0 if age_1999==38
replace income2022=0 if age_1999==38
replace income2021=0 if age_1999==38
//39
replace income2031=0 if age_1999==39
replace income2030=0 if age_1999==39
replace income2029=0 if age_1999==39
replace income2028=0 if age_1999==39
replace income2027=0 if age_1999==39
replace income2026=0 if age_1999==39
replace income2025=0 if age_1999==39
replace income2024=0 if age_1999==39
replace income2023=0 if age_1999==39
replace income2022=0 if age_1999==39
replace income2021=0 if age_1999==39
replace income2020=0 if age_1999==39
//40
replace income2031=0 if age_1999==40
replace income2030=0 if age_1999==40
replace income2029=0 if age_1999==40
replace income2028=0 if age_1999==40
replace income2027=0 if age_1999==40
replace income2026=0 if age_1999==40
replace income2025=0 if age_1999==40
replace income2024=0 if age_1999==40
replace income2023=0 if age_1999==40
replace income2022=0 if age_1999==40
replace income2021=0 if age_1999==40
replace income2020=0 if age_1999==40
replace income2019=0 if age_1999==40
//41
replace income2031=0 if age_1999==41
replace income2030=0 if age_1999==41
replace income2029=0 if age_1999==41
replace income2028=0 if age_1999==41
replace income2027=0 if age_1999==41
replace income2026=0 if age_1999==41
replace income2025=0 if age_1999==41
replace income2024=0 if age_1999==41
replace income2023=0 if age_1999==41
replace income2022=0 if age_1999==41
replace income2021=0 if age_1999==41
replace income2020=0 if age_1999==41
replace income2019=0 if age_1999==41
replace income2018=0 if age_1999==41
//42
replace income2031=0 if age_1999==42
replace income2030=0 if age_1999==42
replace income2029=0 if age_1999==42
replace income2028=0 if age_1999==42
replace income2027=0 if age_1999==42
replace income2026=0 if age_1999==42
replace income2025=0 if age_1999==42
replace income2024=0 if age_1999==42
replace income2023=0 if age_1999==42
replace income2022=0 if age_1999==42
replace income2021=0 if age_1999==42
replace income2020=0 if age_1999==42
replace income2019=0 if age_1999==42
replace income2018=0 if age_1999==42
replace income2017=0 if age_1999==42
//43
replace income2031=0 if age_1999==43
replace income2030=0 if age_1999==43
replace income2029=0 if age_1999==43
replace income2028=0 if age_1999==43
replace income2027=0 if age_1999==43
replace income2026=0 if age_1999==43
replace income2025=0 if age_1999==43
replace income2024=0 if age_1999==43
replace income2023=0 if age_1999==43
replace income2022=0 if age_1999==43
replace income2021=0 if age_1999==43
replace income2020=0 if age_1999==43
replace income2019=0 if age_1999==43
replace income2018=0 if age_1999==43
replace income2017=0 if age_1999==43
replace income2016=0 if age_1999==43
//44
replace income2031=0 if age_1999==44
replace income2030=0 if age_1999==44
replace income2029=0 if age_1999==44
replace income2028=0 if age_1999==44
replace income2027=0 if age_1999==44
replace income2026=0 if age_1999==44
replace income2025=0 if age_1999==44
replace income2024=0 if age_1999==44
replace income2023=0 if age_1999==44
replace income2022=0 if age_1999==44
replace income2021=0 if age_1999==44
replace income2020=0 if age_1999==44
replace income2019=0 if age_1999==44
replace income2018=0 if age_1999==44
replace income2017=0 if age_1999==44
replace income2016=0 if age_1999==44
replace income2015=0 if age_1999==44
//45
replace income2031=0 if age_1999==45
replace income2030=0 if age_1999==45
replace income2029=0 if age_1999==45
replace income2028=0 if age_1999==45
replace income2027=0 if age_1999==45
replace income2026=0 if age_1999==45
replace income2025=0 if age_1999==45
replace income2024=0 if age_1999==45
replace income2023=0 if age_1999==45
replace income2022=0 if age_1999==45
replace income2021=0 if age_1999==45
replace income2020=0 if age_1999==45
replace income2019=0 if age_1999==45
replace income2018=0 if age_1999==45
replace income2017=0 if age_1999==45
replace income2016=0 if age_1999==45
replace income2015=0 if age_1999==45
replace income2014=0 if age_1999==45
//46
replace income2031=0 if age_1999==46
replace income2030=0 if age_1999==46
replace income2029=0 if age_1999==46
replace income2028=0 if age_1999==46
replace income2027=0 if age_1999==46
replace income2026=0 if age_1999==46
replace income2025=0 if age_1999==46
replace income2024=0 if age_1999==46
replace income2023=0 if age_1999==46
replace income2022=0 if age_1999==46
replace income2021=0 if age_1999==46
replace income2020=0 if age_1999==46
replace income2019=0 if age_1999==46
replace income2018=0 if age_1999==46
replace income2017=0 if age_1999==46
replace income2016=0 if age_1999==46
replace income2015=0 if age_1999==46
replace income2014=0 if age_1999==46
replace income2013=0 if age_1999==46
//47
replace income2031=0 if age_1999==47
replace income2030=0 if age_1999==47
replace income2029=0 if age_1999==47
replace income2028=0 if age_1999==47
replace income2027=0 if age_1999==47
replace income2026=0 if age_1999==47
replace income2025=0 if age_1999==47
replace income2024=0 if age_1999==47
replace income2023=0 if age_1999==47
replace income2022=0 if age_1999==47
replace income2021=0 if age_1999==47
replace income2020=0 if age_1999==47
replace income2019=0 if age_1999==47
replace income2018=0 if age_1999==47
replace income2017=0 if age_1999==47
replace income2016=0 if age_1999==47
replace income2015=0 if age_1999==47
replace income2014=0 if age_1999==47
replace income2013=0 if age_1999==47
replace income2012=0 if age_1999==47
//48
replace income2031=0 if age_1999==48
replace income2030=0 if age_1999==48
replace income2029=0 if age_1999==48
replace income2028=0 if age_1999==48
replace income2027=0 if age_1999==48
replace income2026=0 if age_1999==48
replace income2025=0 if age_1999==48
replace income2024=0 if age_1999==48
replace income2023=0 if age_1999==48
replace income2022=0 if age_1999==48
replace income2021=0 if age_1999==48
replace income2020=0 if age_1999==48
replace income2019=0 if age_1999==48
replace income2018=0 if age_1999==48
replace income2017=0 if age_1999==48
replace income2016=0 if age_1999==48
replace income2015=0 if age_1999==48
replace income2014=0 if age_1999==48
replace income2013=0 if age_1999==48
replace income2012=0 if age_1999==48
replace income2011=0 if age_1999==48
//49
replace income2031=0 if age_1999==49
replace income2030=0 if age_1999==49
replace income2029=0 if age_1999==49
replace income2028=0 if age_1999==49
replace income2027=0 if age_1999==49
replace income2026=0 if age_1999==49
replace income2025=0 if age_1999==49
replace income2024=0 if age_1999==49
replace income2023=0 if age_1999==49
replace income2022=0 if age_1999==49
replace income2021=0 if age_1999==49
replace income2020=0 if age_1999==49
replace income2019=0 if age_1999==49
replace income2018=0 if age_1999==49
replace income2017=0 if age_1999==49
replace income2016=0 if age_1999==49
replace income2015=0 if age_1999==49
replace income2014=0 if age_1999==49
replace income2013=0 if age_1999==49
replace income2012=0 if age_1999==49
replace income2011=0 if age_1999==49
replace income2010=0 if age_1999==49


* 규칙성이 대충 어떤지 감은 오시죠... 그 감을....

forvalues ii = 0(1)21 {
    loc yy = 2031 - `ii'
  forvalues jj = 2031(-1)`yy' {
    replace income`jj' = 0 if age_1999==28 + `ii'
  }
}

* 맞는지 확인하고 싶다면....ㅋ
forvalues ii = 0(1)21 {
    loc yy = 2031 - `ii'
  forvalues jj = 2031(-1)`yy' {
    di " 나이는 `=28+`ii''  "
    di " 연도는 `jj' "
  }
}




댓글

댓글 쓰기

이 블로그의 인기 게시물

pdf 문자 인식(OCR)

acrobat을 이용하면 대부분 pdf문서의 문자들을 인식하여(OCR기능), 한글, 워드 및 엑셀에서 유용하게 사용할 수 있다. 그러나, 종종 렌더링(rendering)문제가 생겨 문자를 인식하지 못하는 경우가 있다. 이때, 해결방법은 두 가지가 같다. 첫째, 1) pdf문서 인쇄 시 프린터를  microsoft XPS document writer로 인쇄해서 xps문서 만들기 2) xps문서를 pdf로 전환 3) 전환된 pdf문서에서 OCR 기능 실행하면  ~~ 끝!! 둘째, 1) pdf문서를 다른 이름으로 저장할 때, 그 형식으로 tif 혹은 tiff로 저장 2) 페이지별로 생성된 그림파일들을 pdf문서로 결합 3)결합된 pdf문서에서 OCR 기능 실행하면 ~~끝!! (다른 그림 파일들은 안됨!! 사실, 해보지는 않았음 ~ㅋ)
자세한 내용 보기

그래프 그리기 05:= EU-15 국가들의 GDP 대비 사회지출 비중과 사회지출 대비 사회보장기여금의 비중:= 4사분면 만들기

이미지
// 이번에 간단하게 그래프에 4사분면을 구획하는 방법입니다. // 이를 위해 사용한 자료는 2010년 기준  "EU-15 국가들의 GDP 대비 사회지출 비중과 사회지출 대비 사회보장기여금의 비중"입니다. // 이것은 국가별로 사회지출에 필요한 재원조달에 있어서 사회보장기여금(=사회보험료)이 어느 정도 역할하는가에 대한 국가간 비교를 통해  복지모델 유형화를 하거나 혹은 유형을 확인하는 것입니다. // 1999년 Bonoli가 처음 시도했다고 해서, 일명 "Bonoli Matrix"라고도 합니다. // 자, 보면, 중간에 선이 들어가 있는데, 이 선은 15개 국가들만의 평균값입니다.. 즉, 각 평균값을 기준으로 선을 그으면 4사분면이 생기죠....^^ // 변수: ser  (사회지출비중),  sc_se (사회지출 대비 사회보장기여금의 비중) #delimit ;  twoway (scatter ser sc_se , mlabel(wid) mlabposition(3)) if yr == 2010   , xlabel(3 61 , grid )   ylabel(23 33  , angle(horizontal))  xline(41.6 ,  lwidth(thick) lpattern(solid) lcolor(gs13) )  yline(26.9 ,  lwidth(thick) lpattern(solid) lcolor(gs13) )    ; // 보시면, 이해하실거란 생각해서 별 다른 설명은 하지 않겠습니다. // 다만, 2010년이 아니라 여러 연도에 대해서 동일한 그래프를 반복해서 그리고 싶죠?!^^ // 그러러면, 횡축과 종축의 최대값과 최소값, 그리고 각각의 평균값들을  로컬변수로 받아내서 반복문을 만들면 되겠죠!!
자세한 내용 보기

그래프 그리기 06:= 등탄력적 효용함수 그래프 그리기

이미지
// 스타타를 이용하여 일반 함수, 즉 y=x ,  y=x^2 + x + 3 와 같은 함수를 그리고 싶을때 // 오늘은 경제학에서 자주 사용되는 등탄력적 효용함수를 그려보도록 하겠습니다. 이때, 탄력도의 차이에 따른 함수 형태의 차이를 보고자 합니다. // 그리고 탄력도의 차이에 따른 함수에 따라 소득수준에 따른 한계효용의 변화 역시 비교합니다. // 등탄력적 효용함수의 일반적 형태...         // 이러한 효용함수의 특징은 소득수준에 따라 탄력도가 항상  ε 로 같다는 것입니다. // 즉,  ε 값에 따라 그 함수를 그리면 됩니다. // 함수를 그리는 기본적인 코드는 다음과 같습니다. twoway (function y = (x^(1-0) - 1)/(1-0),     range(0 5)     // ε=0일때, x값 범위를 0과 5로하여 그래프를 그리시오....    // 즉,  ε=0일때, 위 함수는 y=x-1이 됩니다...    // 이건 스타타 명령어 줄에서 아무런 자료가 없어도 그려집니다... 위의 function이라는 부분땜시...ㅋ // 여하튼,  ε 가 0 , 0.5, 1, 2, 5일때 각 함수를 하나의 그래프로 그리려면.... #delimit ; twoway (function y = (x^(1-0) - 1)/(1-0),     ra(0 5) lcolor(black) lwidth(medium) lpattern(solid))        (function y = (x^(1-0.5) - 1)/(1-0.5), ra(0 5) lcolor(blue) lwidth(medium) lpattern(solid))        (function y = (x^(1-1...
자세한 내용 보기