Thanks everyone for the advice, you raise interesting points. Maybe the
best thing for me to do is do an ANOVA in R with binomial data (if
possible) and find the lowest dose that gives a significant difference from
the controls.
On Mon, Apr 2, 2012 at 2:45 PM, Danielle Duncan dldunc
I suppose I'll just report a LC10 using the dose.p function in the package
MASS using my glm fitted logistic regression on binomial data. Thanks
everyone for ideas input! The LOEC seems to be a flawed
calculation...I'll research it. Again, thanks all!
On Mon, Apr 2, 2012 at 2:45 PM, Danielle
be useful as
an example..
data(down)
with(down, plot(age, cases/births))
There is a paper of mine on R news 2008 discussing the package..
hope this helps you,
vito
On Mon, 2 Apr 2012 14:45:06 -0800, Danielle Duncan wrote
Hello, I used the glm function in R to fit a dose-response
:10, predict(glmfit,newdata=data.frame(X=1:10),type=response),
type=l,ylim=c(0,1),xlab=X,ylab=Y)
rug(jitter(X[Y==0]),side=1)
rug(jitter(X[Y==1]),side=3)
On Tue, Apr 3, 2012 at 3:19 PM, Danielle Duncan dldunc...@alaska.eduwrote:
Thanks, that is interesting, but what I'm really after is an easy
Hello, I used the glm function in R to fit a dose-response relationship and
then have been using dose.p to calculate the LC50, however I would like to
calculate the NOEL (no observed effect level), ie the lowest dose above
which responses start occurring. Does anyone know how to do this?
Greetings, I have a question that I'd like to get input on. I have a
classic toxicology study where I artificially fertilized and exposed
embryos to a chemical and counted defects. In addition, I kept track of
male-female pairs that I used to artificially fertilize and generate
embryos with. I
6 matches
Mail list logo