how to convert survival rate

martinig / Fitness-and-dispersal-MA

https://martinig.github.io/Fitness-and-dispersal-MA/

0 stars 0 forks source link

how to convert survival rate #2

Closed martinig closed 9 months ago

martinig commented 1 year ago

[x] code to convert survival rate to effect size @itchyshin

itchyshin commented 1 year ago

@martinig - please give me a couple of examples here so I can work it out

martinig commented 9 months ago

Treat it like a proportion with error.

h <- n/n1 + n/n2
    p <- n1/n # prop for n1
    q <- n2/n # prop for n2
    s1 <- 1/sqrt(sd1^2/(4*m1*(1-m1)))
    s2 <- 1/sqrt(sd2^2/(4*m2*(1-m2)))
    m1 <- asin(sqrt(m1/100))
    m2 <- asin(sqrt(m2/100))
    s_pool <- sqrt( ((n1-1)*s1^2 + (n2-1)*s2^2) / (n - 2) )
    j <- 1 - (3 / (4*n - 9))
    d <- ((m2 - m1) / s_pool) * j
    r_pb <-  d / sqrt(d^2 + h)
    r_b <- r_pb*(sqrt(p*q)/dnorm(qnorm(p)))
    r <- r_b #r_b = r