'ds.glm' code amended as 'c()' cannot be used for aggregate function like 'glm.ds' which was passed on the starting values as a vector.