Підвищення мінімального оцінювача


9

Припустимо, у мене є позитивних параметрів для оцінки та їх відповідних об'єктивних оцінок, отриманих оцінювачами , тобто , тощо.nμ1,μ2,...,μnnμ1^,μ2^,...,μn^E[μ1^]=μ1E[μ2^]=μ2

Я хотів би оцінити за допомогою оцінок, що знаходяться під рукою. Очевидно, що наївний оцінювач зміщений нижче, ніж min(μ1,μ2,...,μn)min(μ1^,μ2^,...,μn^)

E[min(μ1^,μ2^,...,μn^)]min(μ1,μ2,...,μn)

Припустимо, у мене також є матриця коваріації відповідних оцінювачів . Чи можливо отримати неупереджену (або менш упереджену) оцінку мінімуму, використовуючи дані оцінки та матрицю коваріації?Cov(μ1^,μ2^,...,μn^)=Σ


Чи готові ви використовувати Bayesian MCMC підхід чи вам потрібна формула закритої форми?
Мартін Модрак

Але простий підхід до вибірки в порядку? (також вам не потрібні пріори для аналізу Байєса, але це вже інша історія)
Martin Modrák

@ MartinModrák Я не маю досвіду вибірки підходів. Якщо я займаюся байезіаном, зазвичай роблю прості речі, що поєднуються. Але якщо ти думаєш, що це шлях, я продовжую вчитися.
Cagdas Ozgenc

Що ще вам відомо про ці оцінки? Чи знаєте ви вирази? Чи знаєте ви розподіл даних, які використовуються для оцінки цих параметрів?
wij

@wij Я можу спробувати оцінити деякі інші моменти оцінювачів, якщо це потрібно. У мене немає аналітичного вираження для розподілу оцінок. Рішення не повинно залежати від розповсюдження даних.
Cagdas Ozgenc

Відповіді:


4

Я не маю однозначної відповіді щодо існування об'єктивного оцінювача. Однак, з точки зору помилки оцінки, оцінка взагалі є непростою проблемою.min(μ1,,μn)

Наприклад, нехай , і . Нехай є цільовою кількістю, а - оцінка . Якщо ми використовуємо "наївний" оцінювач де , то помилка оцінки є верхньою межею до постійної. (Зверніть увагу, що похибка оцінки для кожного є ). Звичайно, якщоY1,,YNN(μ,σ2I)μ=(μ1,,μn)θ=miniμiθ^θθ^=mini(Y¯i)Yi¯=1Nj=1NYi,jL2

E[θ^θ]2σ2lognN
μiσ2Nμi'є далеко один від одного і дуже мала, похибка оцінки повинна бути зведена до . Однак, у гіршому випадку, немає оцінки працює краще, ніж наївний оцінювач. Ви можете точно показати, що де infimum переймає всі можливі оцінки на основі вибірки а приймає всю можливу конфігурацію .σσ2Nθ
infθ^supμ1,,μnE[θ^θ]2σ2lognN
θY1,,YNμi

Тому наївний оцінювач є мінімальним оптимальним до постійного, і кращої оцінки в цьому сенсі немає.θ


Надана додаткова інформація зовсім не допомагає? Яка додаткова статистика може бути корисною?
Cagdas Ozgenc

Вибачте за те, що ви заплуталися. Я не мав на увазі, що додаткова інформація (коваріація) не є корисною. Я просто хотів зазначити, що оцінити мінімум кількох засобів для населення є складним характером. Інформація про коваріацію має бути корисною. Наприклад, у звичайному випадку, якщо ми маємо ідеальні кореляції для всіх можливих пар, це означає, що випадкові спостереження походять від різного середнього значення + загального шумового терміна. У цьому випадку наївний оцінювач (мінімум зразкових засобів) є неупередженим.
JaeHyeok Shin

3

EDIT: Наступні відповіді на інше питання, ніж те, що було задано - воно оформлено так, ніби вважається випадковим, але не працює, коли вважається фіксованим, що, мабуть, було на увазі ОП. Якщо виправлено, у мене немає кращої відповіді, ніжμμμmin(μ^1,...,μ^n)


Якщо ми розглядаємо лише оцінки середнього значення та коваріації, ми можемо розглядати як єдиний зразок із багатоваріантного нормального розподілу. Простий спосіб отримати мінімальну оцінку - це потім набрати велику кількість зразків з , обчислити мінімум кожного зразка, а потім взяти середнє значення цих мінімумів.(μ1,...,μn)MVN(μ^,Σ)

Описану вище процедуру та її обмеження можна зрозуміти в байесівських термінах - взяття позначень із Вікіпедії на MVN , якщо є відомою коваріацією оцінювачів і у нас є одне спостереження, спільний задній розподіл є де і виникають з попереднього місця, де перед спостереженням будь-яких даних ми беремо попередній ). Оскільки ви, ймовірно, не бажаєте ставити пріори на , ми можемо прийняти обмеження як , в результаті чого передній стає рівним, а задній стаєΣμMVN(μ^+mλ01+m,1n+mΣ)λ0mμMVN(λ0,m1Σμm0μMVN(μ^,Σ). Однак, враховуючи плоскість до цього, ми неявно робимо припущення, що елементи сильно різняться (якщо всі дійсні числа однаково ймовірні, отримання подібних значень малоймовірне).μ

А швидке моделювання показує , що оцінка цієї процедури дещо завищеною , коли елементи відрізняються багато і недооцінює , коли елементи подібні. Можна стверджувати, що без будь-якого попереднього знання це правильна поведінка. Якщо ви бажаєте вказати хоча б деяку попередню інформацію (наприклад, ), результати можуть трохи поводитися для вашого випадку використання.min(μ)μmin(μ)m=0.1

Якщо ви готові взяти на себе більшу структуру, можливо, ви зможете вибрати кращий розподіл, ніж звичайний багатоваріантний. Також може мати сенс використовувати Стен або інший пробовідбірник MCMC, щоб відповідати оцінкам в першу чергу. Ви отримаєте набір зразків які відображають невизначеність у самих оцінниках, включаючи їх структуру коваріації (можливо, багатшу, ніж може надати MVN). Ще раз ви можете обчислити мінімум для кожного зразка, щоб отримати задній розподіл за мінімумами, і взяти середнє значення цього розподілу, якщо вам потрібна бальна оцінка.μ(μ1,...,μn)


Зауважте, що я не намагаюся оцінити мінімум N випадкових величин. Я намагаюся оцінити мінімум N параметрів. Здається, ваша пропозиція є оцінкою для тоді як мені потрібна оцінкаE[min(μ1^,μ2^,...,μn^)]min(μ1,μ2,...,μn)
Cagdas Ozgenc

Я спробував відредагувати відповідь, щоб пояснити обґрунтування, сподіваюся, що це допомагає.
Мартін Модрак

Таким чином, цей метод вибірки дає кращі результати порівняно з простим оцінкою , який також добре працює, коли далеко один від одного і недооцінюється, коли вони поруч. Щоб це було корисно, воно повинно працювати, коли вони поруч. min(μ1^,μ2^,...,μn^)μi
Cagdas Ozgenc

Також зауважте, що всі це додатні числа, тому вам не потрібна від'ємна частина реальної лінії. мкi
Cagdas Ozgenc

1
Ви праві, що я ігнорую знаки і не бачу простого способу їх розміщення. Також оцінювач, який я запропонував, працює краще, коли вважається випадковим, але він гірший, ніж для фіксованого . Я не думаю, що я можу це врятувати, і я не впевнений, який найкращий шлях вперед - я схильний спробувати видалити відповідь, оскільки це насправді не відповідає на питання, але (сподіваюся) відповідь також містить деякі ідеї, які може комусь бути корисним. мкмiн(мк^)мк
Мартін Модрак
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.