Viewed the code, but it did not help me much.
Maybe depending on the function argument for Wright that you specify.
vrtest functions below:
> Wright
function (y, kvec)
{
y <- as.matrix(y)
n <- nrow(y)
W_mat <- matrix(NA, nrow = length(kvec), ncol = 3)
for (i in 1:length(kvec)) {
k <- kvec[i]
W <- Wright_stat(y, k)
W_mat[i, ] <- cbind(W$WR1, W$WR2, W$WS1)
}
VR <- W_mat
rownames(VR) <- paste("k=", kvec, sep = "")
colnames(VR) <- c("R1", "R2", "R3")
return(list(Stats = VR))
}
<bytecode: 0x0000024df212f038>
<environment: namespace:vrtest>
> vrtest:::Wright_stat
function (y, k)
{
y <- as.matrix(y)
n <- nrow(y)
ranking <- as.matrix(rank(y))
r1 <- (ranking - 0.5 * (n + 1))/sqrt((n - 1) * (n + 1)/12)
r2 <- qnorm(ranking/(n + 1))
s <- sign(y)
s[s == 0] <- -1
R1 <- stat(r1, k)
R2 <- stat(r2, k)
S1 <- stat(s, k)
return(list(WR1 = R1, WR2 = R2, WS1 = S1))
}
<bytecode: 0x0000024df2127498>
<environment: namespace:vrtest>
> vrtest:::stat
function (x, k)
{
x <- as.matrix(x)
n <- nrow(x)
index <- 1:k
summ <- 0
for (i in k:n) {
summ <- summ + sum(x[index])^2
index <- index + 1
}
vr1 <- sum(x^2)/n
vr2 <- summ/(n * k)
vr <- vr2/vr1
tem1 <- 2 * (2 * k - 1) * (k - 1)
tem2 <- 3 * k * n
vrstat <- (vr - 1)/sqrt(tem1/tem2)
}
<bytecode: 0x0000024df211b770>
<environment: namespace:vrtest>
>