corMat1 <- cor.mat(K=2, M = 2, n=107, n0=198, n0t = 43)$cormat
fwer_critical(ntrt=4, fwer=0.025, corMat=corMat1)
#
#$ntrt
#[1] 4
#
#$fwer
#[1] 0.025
#
#$corMat
# [,1] [,2] [,3] [,4]
#[1,] 1.0000000 0.3508197 0.2746316 0.2746316
#[2,] 0.3508197 1.0000000 0.2746316 0.2746316
#[3,] 0.2746316 0.2746316 1.0000000 0.3508197
#[4,] 0.2746316 0.2746316 0.3508197 1.0000000
#
#$pairwise_alpha
#[1] 0.006657461
#
#$critical_val
#[1] 2.475233
Run the code above in your browser using DataLab