# generate data
set.seed(5732)
test <- function(x)
        {.3*(1e6*(x^11*(1-x)^6)+1e4*(x^3*(1-x)^10))-2}
x <- (0:100)/100
p <- 1-1/(1+exp(test(x)))
y <- rbinom(x,3,p)

# fit the model and evaluate
logit.fit <- gssanova(cbind(y,3-y)~x,family="binomial")
logit.fit1 <- gssanova(y/3~x,"binomial",weights=rep(3,101))
est <- predict(logit.fit,data.frame(x=x),se=TRUE)

# plot the fit
plot(x,y/3,ylab="p(x)")
lines(x,p,col=6)
lines(x,1-1/(1+exp(est$fit)),col=4)
lines(x,1-1/(1+exp(est$fit+1.96*est$se)),col=5)
lines(x,1-1/(1+exp(est$fit-1.96*est$se)),col=5)