# Load 'dataset' (n = 2000, 700 = support, 700 = oppose, 600 = don't know)
load("dataset.rda")

#Subset respondents
datasetSupOpp <- subset(dataset, !is.na(dataset$Q1SupOppBinary))
datasetDontKnow <- subset(dataset, is.na(dataset$Q1SupOppBinary))

#Create weighted version of SupOpp subset
datasetSupOpp.W <- svydesign(ids = ~1, data = datasetSupOpp, weights =
datasetSupOpp$Weight)

#Setup prediction model
#SupOppBinary describes support or opposition (1 = support, 0 = opposition)
predictionModel <- svyglm(SupOppBinary~Age + Social_Grade + Education +
Gender, design=datasetSupOpp.W)

#Attempt predict
predictions <- predict.glm(predictionModel, newdata = datasetDontKnow, family = binomial)