ВсОА 2023, Сбор 1 · анализ на данни — Задача 1. D11 Одиннадцатый круг Оушена

Q23S1D

Условие

В файле delSct.dat представлены данные о переменных типа δ\delta Sct в затменных двойных системах:

  • TIC — номер объекта в каталоге,
  • RA и DEC — экваториальные координаты звезды в градусах,
  • TeffTESS — эффективная температура по данным TESS,
  • TeffGAIA — эффективная температура по данным GAIA,
  • LGAIA — светимость по данным GAIA,
  • Etype — тип затменной переменной.

a. Определите коэффициент линейной корреляции и его значимость для температур TeffTESS и TeffGAIA для всей выборки объектов, в диапазоне температур до 6 000 K, в диапазоне температур от 6 000 K до 9 000 K.

Коэффициент линейной корреляции для набора величин (xi,yi)(x_i, y_i), i=1,…,Ni = 1, \ldots, N вычисляется как

r=∑i=1N(xi−xˉ)(yi−yˉ)∑i=1N(xi−xˉ)2∑i=1N(yi−yˉ)2.r = \frac{\sum_{i=1}^{N}(x_i-\bar{x})(y_i-\bar{y})}{\sqrt{\sum_{i=1}^{N}(x_i-\bar{x})^2}\sqrt{\sum_{i=1}^{N}(y_i-\bar{y})^2}}.

Здесь xˉ\bar{x} и yˉ\bar{y} — выборочные средние значения.

Значимость корреляции можно определить следующим образом. Вероятность получить ∣r∣|r| больше измеренного при нулевой корреляции случайных величин xx и yy при небольшом объеме выборки в рамках нуль-гипотезы об отсутствии корреляции можно рассчитать по формуле (Press et al., 1997):

P=1−SN−2(t),t=rN−21−r2,P = 1 - S_{N-2}(t), \qquad t = r\sqrt{\frac{N-2}{1-r^2}},

где SN−2(t)S_{N-2}(t) — распределение Стьюдента с N−2N-2 степенями свободы:

Sκ(t)=Γ(κ/2)π(κ−1)⋅Γ((κ−1)/2)∫−tt(1+u2κ−1)−κ/2du.S_{\kappa}(t) = \frac{\Gamma(\kappa/2)}{\sqrt{\pi(\kappa-1)}\cdot\Gamma((\kappa-1)/2)}\int_{-t}^{t}\left(1+\frac{u^2}{\kappa-1}\right)^{-\kappa/2}du.

В таблице приведены значения SκS_{\kappa} (в процентах в верхней строке) для некоторых величин κ\kappa (левый столбец) при разных значениях tt (клетки самой таблицы):

50%60%70%80%90%95%97%99%99.5%99.9%
310.680.851.061.311.702.042.272.753.023.63
500.680.851.051.301.682.012.232.682.943.49
850.680.851.041.291.661.992.212.642.883.42
1270.680.841.041.291.661.982.192.612.853.37
1600.680.841.041.291.651.982.192.612.853.35

b. На диаграмме рассеяния (TeffTESS, TeffGAIA) в диапазоне температур от 6 000 K до 9 000 K проведите аппроксимационную прямую вида TeffGAIA = k⋅k\cdot TeffTESS, определите её угловой коэффициент.

Решение

Внимание: Непълно решение

В архива няма официално решение на тази задача.


Оригинал в Архива: Q23S1D.pdf