x <- 1:30
y <- 1/(1+exp(0.5*(15-x)))+rnorm(30)/20
TestRun <- gcBootSpline(x,y,"ID",grofit.control(nboot.gc=50))
plot(TestRun, pch=2, cex=3, colData="steelblue",
colSpline=c("orchid", "orchid1", "orchid2", "orchid3", "orchid4"))
Run the code above in your browser using DataLab