Nothing
"nullinfo" <-
function(sel.frac)
{
z <- -qnorm(sel.frac/2)
2*( z*dnorm(z) + sel.frac/2 )
}
"info" <- function(sel.frac,theta=0,cross)
{
if(cross=="bc")
info.bc(sel.frac,theta)
else if(cross=="f2")
info.f2(sel.frac,theta)
else
stop("Unknown cross ", cross, ".")
}
"info.bc" <-
function(sel.frac,theta=0)
{
nullinfo(sel.frac)*deflate.bc(theta)
}
"info.f2" <-
function(sel.frac,theta=0)
{
defl <- deflate.f2(theta)
list( add=nullinfo(sel.frac)*defl$add, dom=nullinfo(sel.frac)*defl$dom )
}
"deflate" <- function(theta,cross)
{
if(cross=="bc")
deflate.bc(theta)
else if(cross=="f2")
deflate.f2(theta)
else
stop("Unknown cross ", cross, ".")
}
"deflate.bc" <-
function(theta)
{
theta1 <- recomb(0.5*genetic.dist(theta))
q <- ((1-theta1)^2)/(theta1^2 + (1-theta1)^2)
A <- (1-4*q*(1-q))*(1-theta)
A
}
"deflate.f2" <-
function(theta)
{
t <- recomb(0.5*genetic.dist(theta))
num <- 6*t^4 - 12*t^3 +10*t^2 - 4*t + 1
den <- 8*t^4 - 16*t^3 +12*t^2 - 4*t + 1
add.factor <- deflate.bc(theta)
dom.factor <- (add.factor^2)*num/den
list(add=add.factor,dom=dom.factor)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.