Displaying 2 results from an estimated 2 matches for "mat50_s1".
Did you mean:
mat50_s2
2007 May 02
3
ED50 from logistic model with interactions
...it <- glm(Mature ~ Season * Size - 1, family = binomial, data=dat)
where Mature is a binary response, 0 for immature, 1 for mature. There
are 3 Seasons.
The Season * Size interaction is significant. I would like to compare the
size at 50% maturity between Seasons, which I have calculated as:
Mat50_S1 <- -fit$coef[1]/fit$coef[4]
Mat50_S2 <- -fit$coef[2]/(fit$coef[4] + fit$coef[5])
Mat50_S3 <- -fit$coef[3]/(fit$coef[4] + fit$coef[6])
But I am not sure how to calculate the standard error around each of
these estimates. The p.dose function from the MASS package does this
automatically, bu...
2007 May 03
0
unscrible pls
...- 1, family = binomial, data=dat)
>
> where Mature is a binary response, 0 for immature, 1 for mature. There
> are 3 Seasons.
>
> The Season * Size interaction is significant. I would like to compare the
> size at 50% maturity between Seasons, which I have calculated as:
>
> Mat50_S1 <- -fit$coef[1]/fit$coef[4]
> Mat50_S2 <- -fit$coef[2]/(fit$coef[4] + fit$coef[5])
> Mat50_S3 <- -fit$coef[3]/(fit$coef[4] + fit$coef[6])
>
> But I am not sure how to calculate the standard error around each of
> these estimates. The p.dose function from the MASS package doe...