1
Bagaimana menghitung matriks topi untuk regresi logistik di R?
Saya ingin menghitung matriks topi langsung di R untuk model logit. Menurut Long (1997) matriks topi untuk model logit didefinisikan sebagai: H=VX(X′VX)−1X′VH=VX(X′VX)−1X′VH = VX(X'VX)^{-1} X'V X adalah vektor variabel independen, dan V adalah matriks diagonal dengan pada diagonal.π(1−π)−−−−−−−√π(1−π)\sqrt{\pi(1-\pi)} Saya menggunakan optimfungsi ini untuk memaksimalkan kemungkinan dan menurunkan hessian. Jadi saya …