Una domanda da principianti: qualcuno sa come eseguire una regressione logistica con errori standard in cluster in R? In Stata è solo logit Y X1 X2 X3, vce(cluster Z)
, ma sfortunatamente non ho capito come fare la stessa analisi in R. Grazie in anticipo!Regressione logistica con robusti errori standard in cluster in R
risposta
Si potrebbe voler esaminare il pacchetto rms
(strategie di modellazione della regressione). Così, lrm
è modello di regressione logistica, e se fit
è il nome della vostra uscita, devi avere qualcosa di simile:
fit=lrm(disease ~ age + study + rcs(bmi,3), x=T, y=T, data=dataf)
fit
robcov(fit, cluster=dataf$id)
bootcov(fit,cluster=dataf$id)
È necessario specificare x=T
, y=T
nel modello di dichiarazione. rcs
indica spline cubiche limitate a 3 nodi.
Grazie mille! Ha funzionato a meraviglia! Leggerò più attentamente il manuale di rms e vedrò se esiste un modo per raggruppare i coefficienti per paese e anche per anno. Ancora una volta grazie! – danilofreire
Questa risposta è già molto buona ma potrebbe essere migliorata se fosse completamente replicabile. Non ho idea di dove vengano le variabili, quale sia l'output e perché 'rcs (bmi, 3)' sia necessario. – MERose
Ho battuto la testa contro questo problema negli ultimi due giorni; Ho trovato magicamente quello che sembra un nuovo pacchetto che sembra destinato a grandi cose - ad esempio, sto anche eseguendo nella mia analisi alcuni modelli Tobit robusti in cluster, e questo pacchetto ha anche quella funzionalità integrata. Per non parlare della sintassi è molto più pulito rispetto a tutte le altre soluzioni che ho visto (stiamo parlando di livelli di pulizia quasi-Stata).
Così, per il tuo esempio giocattolo, mi piacerebbe correre:
library(Zelig)
logit<-zelig(Y~X1+X2+X3,data=data,model="logit",robust=T,cluster="Z")
Et voilà!
Wow, sembra "funziona" in modi che il mio codice R non sembra mai. Questa nuova funzionalità? Se no, perché Zelig non è stato il modo canonico per risolvere questo in R? – Philip
Non so, ma spero che diventi così. [Il progetto] (http://zeligproject.org) sembra certamente ambizioso! Il [gruppo Google] (https://groups.google.com/forum/m/#!forum/zelig-statistical-software) non sembra però così attivo, quindi non sei sicuro di quanto siano rapidi i progressi. – MichaelChirico
Purtroppo, penso che il comando non funzioni nell'ultima versione di 'Zelig' (su CRAN). Ho appena eseguito alcuni modelli con e senza l'argomento 'cluster' e gli errori standard sono esattamente gli stessi. Credo che sia stato così dalla versione 4.0, l'ultima volta che ho usato il pacchetto. – danilofreire
- 1. Regressione dati pannello: errori standard robusti
- 2. R GBM regressione logistica
- 3. Errori standard robusti diversi della regressione della logica in Stata e R
- 4. Regressione logistica in Java
- 5. Regressione logistica con SciPy
- 6. Programmazione Regressione logistica con discesa del gradiente stocastico in R
- 7. R area di regressione logistica sotto curva
- 8. Regressione logistica ponderata in Python
- 9. Tabella di classificazione per la regressione logistica in R
- 10. errore standard robusto e in cluster in R per regressione probit e logit
- 11. regressione logistica binaria xgboost
- 12. Regressione logistica Wabbit Vowpal
- 13. R2WinBUGS - regressione logistica con dati simulati
- 14. regressione logistica sklearn con classi sbilanciate
- 15. sklearn Regressione logistica - caratteristiche importanti
- 16. Errori standard di Fama MacBeth in R
- 17. Come posso utilizzare la funzione di previsione in R in una regressione logistica installata anni fa?
- 18. tracciamento del limite decisionale della regressione logistica
- 19. Classificazione: PCA e regressione logistica usando sklearn
- 20. Tracciare i risultati di un modello di regressione logistica multivariata in R
- 21. Controllo della soglia nella regressione logistica in Scikit Learn
- 22. Statistiche di adattamento del modello per una regressione logistica
- 23. Regressione logistica quando la risposta è una proporzione (utilizzando JAGS)
- 24. Include robusti PHP per gestire gli errori?
- 25. Come calcolare gli errori standard di boot bootsrapped con R?
- 26. Tracciare la "linea di regressione" dalla regressione multipla in R
- 27. Python: Come utilizzare la regressione logistica multinomiale utilizzando SKlearn
- 28. come tracciare la regressione lineare in R?
- 29. Scikit Apprendimento: coefficienti del modello di regressione logistica: chiarimento
- 30. Errore standard nella regressione non lineare
La funzione 'vcovHC()' nel pacchetto 'sandwich' potrebbe anche essere utile (non è sicuro se si applica alle stime di regressione logistica) –
se si sta migrando da Stata si potrebbe trovare il pacchetto chiamato' 'plm'' utile. Inoltre, c'è il pacchetto chiamato '' pcse'' per implementare errori standard corretti del pannello manipolando la matrice di varianza covarianza dopo la stima –
Grazie mille per le risposte, Ben e Hubert. Verificherò anche i pacchetti che hai suggerito e vedrò se funzionano con stime logistiche. Grazie ancora! – danilofreire