On Jul 9, 2012, at 20:23 , Lawrence, Adaku wrote: > Hello, > > I am trying to determine LD50 and LD95 using dose.p in MASS however some of > the Residual variance is larger than the degrees of freedom. Please can > anyone help with any advice as to how i can correct for this?
Er, in what sense is that a problem? Your code is not reproducible, at least some output to look at might help. -pd > > Here is the model as inputted into R > > > > y<-cbind(dead,n-dead) > > model<-glm(y~log(conc),binomial) > summary(model) > > xv<-seq(min(log(conc)-1),max(log(conc)+1),0.01) > lines(xv,predict(model,list(conc=exp(xv)),type="response")) > > dose.p(model,p=c(0.10,0.25,0.5,0.75,0.90,0.95)) > > > > Thanks > > Adaku > > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd....@cbs.dk Priv: pda...@gmail.com ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
