用于cuzick测试的软件包安装

Question

我应该为cuzick测试安装哪个软件包？我已经尝试过PMCMRplus和install_github（＆＃34; raredd / rawr＆＃34;）但是我仍然无法找到任何cuzick测试包

Answer 1

cuzick = function(x,z,test.type=c("two.sided", "upper", "lower")) {

    N = length(z)
    n = unique(z)

    ranks=rank(x)

    T = sum(ranks*z)

    p = (table(z)/N)
    E_Z = sum(unique(z)*p)
    E_T = 0.5*N*(N+1)*E_Z

    Var_Z = sum(unique(z)^2*p) - E_Z^2
    Var_T = N^2*(N+1)/12*Var_Z

    Zscore = (T-E_T)/sqrt(Var_T)

    if(test.type == "two.sided") {
        pval = 2*pnorm(-abs(Zscore))
    } else if(test.type == "upper") {
        pval = pnorm(Zscore,lower.tail=F)
    } else  pval = pnorm(Zscore,lower.tail=T)

    out = data.frame(cbind(Zscore,pval,test.type))
    colnames(out) = c("Z","p","testType")
    return(out)
}