Ich habe mich gefragt, ob jemand eine Funktion entwickelt hat, die Konfidenzintervalle erzeugen kann aus einem svyby-Objekt für Proportionen (in meinem Fall eine Kreuztabelle für einen binären Gegenstand im Umfragepaket). Ich vergleiche häufig Anteile über Gruppen, und es wäre sehr praktisch, eine Funktion zu haben, die Konfidenzintervalle extrahieren kann (mit der Umfragefunktion svyciprop anstatt confint). Das folgende Beispiel zeigt, was ich erreichen möchte.
Daten laden
%Vor%Erstellen Sie ein svyby-Objekt, das den Anteil der Variablen "both" über den Styp
hinweg vergleicht %Vor%Wäre es möglich, eine Funktion zu erstellen (zB durchCI (b, method="likelihood"), die dasselbe wie confint (b) erreicht, aber svyciprop verwendet? Es müsste im Grunde jede Ebene des svyby-Objekts durchlaufen und ein Vertrauensintervall erstellen Meine Versuche waren bisher erfolglos.
Es mag einen anderen Weg geben, aber ich benutze svyby (), weil es schnell und intuitiv ist.
Vielen Dank im Voraus.
svyby () hat ein Argument vartype =, um anzugeben, wie die Stichprobenunsicherheit angegeben werden soll. Verwenden Sie vartype="ci", um Konfidenzintervalle zu erhalten, zB
%Vor%Es ist leicht zu überprüfen, dass dies dasselbe ist, wie wenn man jedes Level von Hand macht, zB
%Vor%interessant .. diese zwei Befehle sollten nicht das gleiche Ergebnis geben .. der erste sollte wahrscheinlich einen Fehler oder eine Warnung werfen:
%Vor% Sie möchten Dr. Lumley auf dieses Problem aufmerksam machen -
Der Code in der Nähe von Zeile 80 von surveyby.R könnte wahrscheinlich leicht modifiziert werden, um svyciprop in svyby zu arbeiten.
aber ich kann etwas übersehen (und er mag es irgendwo in der Dokumentation notiert haben),
also lesen Sie alles sorgfältig durch, bevor Sie sich mit ihm in Verbindung setzen.
Wie auch immer, hier ist eine temporäre Lösung, die Ihr Problem lösen könnte
%Vor%