Реализация статистического оператора Sweep на Python ⇐ Python

[Anonymous]

Я изучаю некоторые методы ведения статистики с отсутствующими данными из книги («Статистический анализ с отсутствующими данными» Литтла и Рубина). Одной из особенно полезных функций для работы с монотонными данными об отсутствии ответов является Оператор очистки (подробности на стр. 148–151). Я знаю, что модуль R gmm имеет функцию swp, которая делает это, но мне было интересно, реализовал ли кто-нибудь эту функцию в Python, в идеале для матриц Numpy для хранения входных данных. . Я искал StackOverflow, а также безуспешно провел несколько поисков в Интернете. Спасибо за любую помощь.
Вот определение.
Симметричная матрица G PxP называется просматриваемой по строке и столбцу k, если она заменена другой симметричной матрицей PxP H с элементами, определенными следующим образом:
h_kk = -1/g_kk
h_jk = h_kj = g_jk/g_kk for j != k
h_jl = g_jl - g_jk g_kl / g_kk j != k, l != k

G = [g11, g12, g13
g12, g22, g23
g13, g23, g33]
H = SWP(1,G) = [-1/g11, g12/g11, g13/g11
g12/g11, g22-g12^2/g11, g23-g13*g12/g11
g13/g11, g23-g13*g12/g11, g33-g13^2/g11]
kvec = [k1,k2,k3]
SWP[kvec,G] = SWP(k1,SWP(k2,SWP(k3,G)))

Inverse function
H = RSW(k,G)
h_kk = -1/g_kk
h_jk = h_kj = -g_jk/g_kk for j != k
h_jl = g_jk g_kl / g_kk j != k, l != k

G == SWP(k,RSW(k,G)) == RSW(k,SWP(k,G))


Подробнее здесь: https://stackoverflow.com/questions/157 ... p-operator