Nothing
`tableHWE` <-
function (x, strata, ...)
{
if (!inherits(x, "setupSNP")) {
stop("x must be an object of class 'setupSNP'")
}
colSNPs <- attr(x, "colSNPs")
# VM seleccion incorrecta de datos!! data.SNPs <- x[colSNPs,,drop=FALSE]
data.SNPs <- x[ , colSNPs, drop=FALSE]
tt <- mclapply( data.SNPs, table, ...)
ans <- cbind( "HWE (p value)" = unlist(mclapply(tt, SNPHWE, ...)) )
if (!missing(strata)) {
# VM buscar en x si existe
strata.name <- deparse(substitute(strata))
if(!exists(strata.name) & strata.name %in% names(x)) {
strata<-x[,strata.name]
}
if (length(table(strata)) > 5) {
stop("strata looks numeric")
}
strates <- names(table(strata) > 0)
n.strata <- length(strates)
i <- 1
while (i <= n.strata) {
data.SNPs.temp <- subset(data.SNPs, strata == strates[i])
# tt <- apply(data.SNPs.temp, 2, table)
## VM fallaba si todos los SNPs devuelven 3 genotipos porque el resultado es una matriz, no lista
tt <- mclapply(data.SNPs.temp, table, ...)
temp <- cbind("HWE (p value)" = unlist(mclapply(tt, SNPHWE, ...)))
ans <- cbind(ans, temp)
i <- i + 1
}
dimnames(ans)[[2]] <- c("all groups", strates)
}
class(ans) <- c("tableHWE", "matrix")
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.