Nothing
het.perm <- function(x1, x2, B = 999) {
n1 <- dim(x1)[1] ; n2 <- dim(x2)[1]
x <- rbind(x1, x2)
ina <- c( rep(1, n1), rep(2, n2) )
ni <- c(n1, n2)
p <- dim(x)[2] ## dimensionality of the data
n <- n1 + n2 ## sample size of the data
S <- rowsum(x, ina)
kapa <- numeric(2)
mi <- rowsum(x, ina) / ni
kapa[1] <- Directional::vmf.mle( x1, fast = TRUE )$kappa
kapa[2] <- Directional::vmf.mle( x2, fast = TRUE )$kappa
tw <- Rfast::colsums(kapa * ni * mi)
Tt <- sum( kapa * ni * sqrt( Rfast::rowsums(mi^2) ) ) - sqrt( sum(tw^2) )
ptt <- numeric(B)
for (i in 1:B) {
ind <- Rfast2::Sample(ina, n)
S <- rowsum(x, ind)
kapa <- numeric(2)
mi <- rowsum(x, ind) / ni
kapa[1] <- Directional::vmf.mle( x[ind == 1, ], fast = TRUE )$kappa
kapa[2] <- Directional::vmf.mle( x[ind == 2, ], fast = TRUE )$kappa
tw <- Rfast::colsums(kapa * ni * mi)
ptt[i] <- sum( kapa * ni * sqrt( Rfast::rowsums(mi^2) ) ) - sqrt( sum(tw^2) )
}
p.value <- ( sum(ptt > Tt) + 1 ) / (B + 1)
statistic <- 2 * Tt ; names(statistic) <- "Permutation het test statistic"
parameter <- "NA" ; names(parameter) <- "df"
alternative <- "The 2 directional mean vectors differ"
method <- "Permutation ANOVA for 2 directional mean vectors using the heterogeneous approach"
data.name <- c("data ", " groups")
result <- list( statistic = statistic, parameter = parameter, p.value = p.value,
alternative = alternative, method = method, data.name = data.name )
class(result) <- "htest"
return(result)
}
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.