library(survey)

data(api)

# example data set
head( apiclus2 )

# instead of var1 - var4, use these four variables:
ftable( apiclus2[ , c( 'sch.wide' , 'comp.imp' , 'both' , 'awards' ) ] )

# move it over to x for faster typing
x <- apiclus2


# also give x a column of all ones
x$one <- 1

# run the glm() function specified.
logit.1 <-
    glm(
        comp.imp ~ target + cnum + growth ,
        data = x ,
        family = binomial( link = 'logit' )
    )

summary( logit.1 )

# now create the survey object you've described
dclus <-
    svydesign(
        id = ~dnum + snum , # cluster variable(s)
        strata = ~stype ,   # stratum variable
        weights = ~pw ,     # weight variable
        data = x ,
        nest = TRUE
    )

# weighted counts
svyby(
    ~one ,
    ~ sch.wide + comp.imp + both + awards ,
    dclus ,
    svytotal
)


# weighted counts formatted differently
ftable(
    svyby(
        ~one ,
        ~ sch.wide + comp.imp + both + awards ,
        dclus ,
        svytotal ,
        keep.var = FALSE
    )
)


# run the svyglm() function specified.
logit.2 <-
    svyglm(
        comp.imp ~ target + cnum + growth ,
        design = dclus ,
        family = binomial( link = 'logit' )
    )

summary( logit.2 )