[논문 리뷰] Exploring recursion for optimal estimators under cascade rotation

Jan Kowalski, Jacek Wesołowski|arXiv (Cornell University)|2014. 12. 02.

Statistical and numerical algorithms참고 문헌 14인용 수 3

한 줄 요약

이 논문은 반복 설문 조사에서 샘플링 단위가 일시적으로 결석한 후 복귀하는 '공백'이 포함된 로테이션 패턴을 가진 인구 평균의 최적 선형 추정기의 분석적 재귀 공식을 제시한다. 이는 설문 조사 표본 추출 분야에서 오랫동안 남아있던 과제를 해결한다. 대수적 연산자 접근법을 사용하여 문제를 체비셰프 다항식과 관련된 다항식의 근을 찾고, 행렬의 질량을 분석하는 것으로 환원함으로써, 재귀의 차수는 로테이션 패턴에서 가장 큰 공백의 크기 + 1과 일치함을 보여준다.

ABSTRACT

We are concerned with optimal linear estimation of means on subsequent occasions under sample rotation where evolution of samples in time is designed through a cascade pattern. It has been known since the seminal paper of Patterson (1950) that when the units are not allowed to return to the sample after leaving it for certain period (there are no gaps in the rotation pattern), one step recursion for optimal estimator holds. However, in some important real surveys, e.g. Current Population Survey in the US or Labour Force Survey in many countries in Europe, units return to the sample after being absent in the sample for several occasions (there are gaps in rotation patterns). In such situations difficulty of the question of the form of the recurrence for optimal estimator increases drastically. This issue has not been resolved yet. Instead alternative sub-optimal approaches were developed, as K-composite estimation (see e.g. Hansen et al. (1955)), AK-composite estimation (see e.g. Gurney and Daly (1965) or time series approach (see e.g. Binder and Hidiroglou (1988)). In the present paper we overcome this long-standing difficulty, that is, we present analytical recursion formulas for the optimal linear estimator of the mean for schemes with gaps in rotation patterns. It is achieved under some technical conditions: ASSUMPTION I and ASSUMPTION II (numerical experiments suggest that these assumptions might be universally satisfied). To attain the goal we develop an algebraic operator approach which allows to reduce the problem of recursion for the optimal linear estimator to two issues: (1) localization of roots (possibly complex) of a polynomial Q_p defined in terms of the rotation pattern (Q_p happens to be conveniently expressed through Chebyshev polynomials of the first kind), (2) rank of a matrix S defined in terms of the rotation pattern and the roots of the polynomial Q_p.

연구 동기 및 목표

• 반복 설문 조사에서 단위가 일시적으로 결석한 후 복귀하는 로테이션 패턴을 가진 경우, 최적의 선형 무편향 추정기(BLUE)에 대한 재귀 공식을 유도하는 데 오랫동안 열려 있던 문제를 해결하는 것.
• 원래는 복귀하지 않는 단위에만 유효했던 패트릭슨(1950)의 한 단계 재귀를, 결석 후 복귀하는 경우로 확장하는 것.
• 재귀 공식의 문제를 다항식 근의 국소화와 행렬 질량 분석으로 환원하는 일반적인 대수적 프레임워크를 개발하는 것.
• 재귀의 차수가 로테이션 패턴에서 가장 큰 공백의 크기에 의해 결정됨을 규명하는 것.
• 기술적 가정 하에 재귀 계수의 정확한 공식을 제공하며, 수치적 증거는 이러한 가정이 실제 설문 조사 설계에서 널리 만족됨을 시사한다.

제안 방법

• 최적의 선형 추정기를 과거 관측치와 이전 추정치의 재귀적 선형 조합으로 공식화하여, 패트릭슨의 한 단계 재귀를 일반화한다.
• 관측치와 가중치의 시간적 변화를 모델링하기 위해 이동 및 반사 연산자로 구성된 대수적 연산자 프레임워크를 도입한다.
• 첫 번째 종류의 체비셰프 다항식을 통해 표현할 수 있는 다항식 $ Q_p(x) $ 를 정의하며, 이 다항식의 근이 재귀의 구조를 결정한다.
• 재귀 계수를 찾는 문제를 두 가지 핵심 과제로 환원한다: (1) $ Q_p $ 의 근을 국소화하고, (2) 로테이션 패턴과 이 근들로부터 유도된 행렬 $ S $ 의 질량을 계산한다.
• 대각 행렬과 상관 행렬에 대한 행렬 역행 및 추적 항등식을 적용하여 추정기의 분산을 표현한다.
• 연산자 방정식에서 유도된 제약 조건을 사용하여, 도출된 추정기가 편향 없고 최적임을 검증한다.

실험 결과

연구 질문

• RQ1 공백이 있는 로테이션 체계에서, 단위가 결석 후 복귀하는 경우 최적의 선형 무편향 추정기(BLUE)에 대한 재귀 공식을 도출할 수 있는가?
• RQ2 이러한 체계에서 재귀의 차수는 무엇에 의해 결정되는가—특히 로테이션 패턴의 가장 큰 공백의 크기와 관련이 있는가?
• RQ3 재귀 계수는 로테이션 패턴과 상관 구조에 대해 어떻게 분석적으로 표현할 수 있는가?
• RQ4 대수적 연산자 접근법이 유효하고 최적의 재귀를 도출하는 데에는 어떤 조건이 필요한가?
• RQ5 로테이션 패턴에서 유도된 다항식 $ Q_p $ 는 체비셰프 다항식과 체계적으로 연결될 수 있는가—이를 통해 근의 국소화를 단순화할 수 있는가?

주요 결과

• 최적 선형 추정기의 재귀 차수는 정확히 로테이션 패턴에서 가장 큰 공백의 크기 + 1과 일치한다.
• 재귀 계수는 첫 번째 종류의 체비셰프 다항식으로 표현 가능한 다항식 $ Q_p $ 의 근을 통해 정확히 계산할 수 있다.
• 추정기의 최적성은 두 가지 기술적 가정(ASSUMPTION I 및 II) 하에서 보장되며, 수치 실험 결과 이 가정들은 실제 설문 조사 설계에서 널리 만족됨을 시사한다.
• 최적 추정기의 분산은 $ Q_p $ 의 근과 행렬 $ S $ 를 포함하는 연립방정식의 해로부터 유도된 $ heta_{0,0} $ 와 동일함을 보여주며, 이는 추정기의 최소 분산을 확인한다.
• 이 방법은 패트릭슨(1950)의 한 단계 재귀를 공백이 있는 경우로 일반화하여, 이전에는 존재하지 않았던 완전한 분석적 해를 제공한다.
• 실제 설문 조사 체계의 예시를 통해 이론적 프레임워크가 커런트 인구 조사 및 고용력 조사와 같은 설문 조사에 적용 가능함을 확인한다.