Nothing
Sloc <-
function (reference = "F2", i = NULL, genZ = NULL)
{
ans <- NULL
if (reference == "F2") {
ans <- matrix(c(1, -1, -0.5, 1, 0, 0.5, 1, 1, -0.5),
byrow = TRUE, ncol = 3)
}
else if (reference == "P1") {
ans <- matrix(c(1, 0, 0, 1, 1, 1, 1, 2, 0), byrow = TRUE,
ncol = 3)
}
else if (reference == "P2") {
ans <- matrix(c(1, -2, 0, 1, -1, 1, 1, 0, 0), byrow = TRUE,
ncol = 3)
}
else if (reference == "F1") {
ans <- matrix(c(1, -1, -1, 1, 0, 0, 1, 1, -1), byrow = TRUE,
ncol = 3)
}
else if (reference == "Finf") {
ans <- matrix(c(1, -1, 0, 1, 0, 1, 1, 1, 0), byrow = TRUE,
ncol = 3)
}
else if (reference == "UWR") {
ans <- matrix(c(1, -1, -1/3, 1, 0, 2/3, 1, 1, -1/3),
byrow = TRUE, ncol = 3)
}
else if (reference == "noia") {
if (is.null(genZ) || is.null(i)) {
stop("zmat is necessary for statistical model computation")
}
f <- Z2freq(genZ[, (3 * i - 2):(3 * i)])
ans <- matrix(c(1, -f[2] - 2 * f[3], (-2 * f[2] * f[3])/(f[1] +
f[3] - (f[1] - f[3])^2), 1, 1 - f[2] - 2 * f[3],
(4 * f[1] * f[3])/(f[1] + f[3] - (f[1] - f[3])^2),
1, 2 - f[2] - 2 * f[3], (-2 * f[1] * f[2])/(f[1] +
f[3] - (f[1] - f[3])^2)), byrow = TRUE, ncol = 3)
}
else if (reference == "G2A") {
if (is.null(genZ) || is.null(i)) {
stop("zmat is necessary for statistical model computation")
}
f <- Z2freq(genZ[, (3 * i - 2):(3 * i)])
p <- f[1] + f[2]/2
ans <- matrix(c(1, -2 * (1 - p), -2 * (1 - p) * (1 -
p), 1, -1 + 2 * p, 2 * p * (1 - p), 1, 2 * p, -2 *
p * p), byrow = TRUE, ncol = 3)
}
else {
stop("unknown reference ", reference, ".")
}
return(ans)
}
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.