2016-10-29 6 views
0

Ich verwende das Umfragepaket in R für die Modellierung mit komplexen Vermessungsdaten. Ich versuche, das Grundniveau bestimmter Prädiktorvariablen, die in einer logistischen Regression verwendet werden, ohne Erfolg zurückzusetzen. Das Folgende ist ein reproduzierbares Beispiel:Zurücksetzen der Baseline-Ebene der Prädiktorvariablen in svyglm Funktion

library(RCurl) 
library(survey) 
data <- getURL("https://raw.githubusercontent.com/cbenjamin1821/careertech- ed/master/elsq1adj.csv") 
elsq1ch <- read.csv(text = data) 
#Specifying the svyrepdesign object which applies the BRR weights 
elsq1ch_brr<-svrepdesign(variables = elsq1ch[,1:16], repweights = elsq1ch[,18:217], weights = elsq1ch[,17], combined.weights = TRUE, type = "BRR") 
elsq1ch_brr 
#Log. Reg. model 
allCC <- svyglm(formula=F3ATTAINB~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH+F1RTRCC,family="binomial",design=elsq1ch_brr,subset=BYSCTRL==1&G10COHRT==1,na.action=na.omit) 
summary(allCC) 
##Attempting to reset baseline level for predictor variable to "PreAlg or Less"; the baseline level in the model is currently "Adv Alg" 
#Both attempts did not work 
elsq1ch$F1HIMATH <- C(elsq1ch$F1HIMATH,contr.treatment, base=1) 
elsq1ch$F1HIMATH <- relevel(elsq1ch$F1HIMATH,"PreAlg or Less") 
#Log. Reg. model with no changes in baseline levels for the predictor 
allCC <- svyglm(formula=F3ATTAINB~F1PARED+BYINCOME+F1RACE+F1SEX+F1RGPP2+F1HIMATH+F1RTRCC,family="binomial",design=elsq1ch_brr,subset=BYSCTRL==1&G10COHRT==1,na.action=na.omit) 
summary(allCC) 

Jede Führung wird sehr geschätzt.

Antwort

0

Anthony Damico versorgte mich mit diesem Ansatz, gearbeitet: elsq1ch_brr < - update (elsq1ch_brr, F1HIMATH = Nachholung (F1HIMATH "PreAlg or Less"))

Verwandte Themen