2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Log rank test
Сообщение04.09.2013, 14:58 
Аватара пользователя
Даны две группы:
$Group 1: 143, 164, 188^*, 189, 190^*, 102, 206, 209^*, 213$
$Group 0: 142, 156^*, 163, 198^*, 205, 232, 233, 236^*, 241, 253^*, 261$
Нужно сделать лог-ранк тест.
Вот как я это делал:
$\frac{U(0)}{\sqrt{l(0)}}\sim N(0,1)$ где $U(0)=\frac{dl_p}{d\beta}|_{\beta=0}$, $l(0)=\frac{d^2l_p}{d^2\beta}|_{\beta=0}$
В моём случае $U(0)=\sum_{i=1}^n\delta_i [Z_i-\frac{N_1}{N_1+N_0}]=6-13\frac{1017}{1017+1477}=0.6988$
$N_i$ - сумма только тех, у кого не были "зацензурены" данные в группе $i$.
$Z_i=I(group=1)$
$\delta_i=I(observed)$

$l(0)=\sum_{i=1}^n\delta_i [\frac{N_1}{N_1+N_0}-{\frac{N_1}{N_1+N_0}}^2]=13[\frac{1017}{1017+1477}-{\frac{1017}{1017+1477}}^2]=3.139$

Соответственно: $\frac{U(0)}{\sqrt{l(0)}}=\frac{0,6988}{\sqrt{3,139}}=0,3944$
Проблема в том, что когда я использую Matlab для расчета log rank, у меня выходит совершенно другой результат.
Можете, пожалуйста, проверить и указать на ошибку?

 
 
 [ 1 сообщение ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group