Hodges-Lehmann estimation for a variable of a 2x2 study

It performs Hodges-Lehmann estimation for 2x2 study. This is usually for Tmax variable.

hodges(bedata, Var)

Arguments

bedata

Data table name. This should have at least the following columns and a variable column to be tested.

  
 GRP : Group or Sequence, 'RT' or 'TR'
 PRD : Period, 1 or 2
 SUBJ : Subject ID
 TRT : Treatment or Drug, 'R' or 'T'

Var

Variable to be estimated. This should be one of the column names in bedata table. Usually 'Tmax'

Details

It nonparametrically tests Var variable equivalency from a 2x2 study. This is done for a variable which we cannot assume log-normal distribution.

Value

Wilcoxon Signed-Rank Test

A kind of nonparametric test

Hodges-Lehmann Estimate

90% confidence interval in the original scale and the percent scale

Author

Kyun-Seop Bae <k@acr.kr>

Examples

hodges(NCAResult4BE, "Tmax") 
#> $`Wilcoxon Signed-Rank Test`
#>   p-value 
#> 0.2326894 
#> 
#> $`Hodges-Lehmann Estimate`
#>                            Lower Limit Point Estimate Upper Limit
#> 90% Confidence Interval       -0.33000       -0.03500      0.1050
#> 90% Confidence Interval(%)    74.37661       97.28237    108.1529
#>