ld.ci: Confidence Intervals for the Relative Treatment Effects

Description Usage Arguments Details Author(s) References See Also Examples

Description

This function performs calculations of the two-sided confidence intervals for the relative treatment effects of the factors specified. The function performs calculations only if no observations are missing.

Usage

1
2
3
ld.ci(y, time, subject, group=NULL, alpha=0.05, time.name="Time", 
group.name="Group", description=TRUE, time.order=NULL, group.order=NULL,
rounds=4, plot.CI=TRUE, order.warning=TRUE)

Arguments

y

a vector of numeric variable of interest.

time

a vector of the sub-plot factor variable. See Details for more explanation.

subject

a vector of individual subjects.

group

a vector of the whole-plot factor variable; the default option is NULL. See Details for more explanation.

alpha

a numeric specifying the significance level of the confidence intervals; the default option is 0.05.

time.name

a vector of the sub-plot factor variable.

group.name

a vector of the whole-plot factor variable.

description

an indicator for whether a short description of the output should be shown; the default option is TRUE.

time.order

a character or numeric vector specifying the order of the time levels; the default option is NULL, in which case, the levels are in the order of the original data.

group.order

a character or numeric vector specifying the order of the group levels; the default option is NULL, in which case, the levels are in the order of the original data.

rounds

a numeric specifying the number of digits to be displayed.

plot.CI

an indicator for whether a plot of the confidence interval (CI) should be shown; the default option is TRUE.

order.warning

an indicator for whether a short description of the warning regarding the ordering of factors should be shown; the default option is TRUE.

Details

A whole-plot factor refers to a factor effective for each subject at all times. A sub-plot factor refers to a factor effective at a single time point for all time curves and all subjects. See Brunner et al. (2002) for more examples. Also, note that the interval for the relative treatment effects can only be interpreted as a confidence interval when the sample sample sizes are (approximately) the same (pp.60, Brunner et al., 2002).

The function returns a list with the following numeric components: the relative treatment effect (RTE), bias estimation (Bias), variance estimation (Variance), the lower and upper bound of the RTE (Lower bound, Upper bound, respectively), in the form of an n-by-9 data frame where n is the number of group factor levels times the number of time factor levels.

Author(s)

Kimihiro Noguchi, Karthinathan Thangavelu, Frank Konietschke, Yulia Gel, Edgar Brunner

References

Brunner, E., Domhof, S., and Langer, F. (2002). Nonparametric Analysis of Longitudinal Data in Factorial Experiments, Wiley, New York.

Brunner, E. and Langer, F. (1999). Nichtparametrische Analyse longitudinaler Daten, R. Oldenbourg Verlag, Munchen Wien.

Noguchi, K., Gel, Y.R., Brunner, E., and Konietschke, F. (2012). nparLD: An R Software Package for the Nonparametric Analysis of Longitudinal Data in Factorial Experiments. Journal of Statistical Software, 50(12), 1-23.

See Also

nparLD, ld.f1, ld.f2, f1.ld.f1, f1.ld.f2, f2.ld.f1, shoulder

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
## Example with the "Shoulder tip pain study" data ##
data(shoulder)
attach(shoulder)
group<-paste(group1, group2, sep="")
ex.ci<-ld.ci(y=resp, time=time, subject=subject, group=group, alpha=0.05, 
time.name="Time", group.name="Group", description=FALSE, time.order=c(1,2,3,4,5,6), 
group.order=c("YF","YM","NF","NM"))
# LD CI Calculations 
# ----------------------- 
# Order of the time and group levels.
# Time level:   1 2 3 4 5 6 
# Group level:   YF YM NF NM 
# The order may be specified in time.order or group.order.

## Summary of the output
ex.ci

#     Group Time Nobs RankMeans    RTE    Bias Variance  Lower  Upper
#1  GroupYF    1   14  123.9643 0.5019  0.0008   0.1680 0.3792 0.6243
#2  GroupYF    2   14  100.2857 0.4056  0.0003   0.1366 0.3000 0.5222
#3  GroupYF    3   14   89.2500 0.3608 -0.0005   0.0979 0.2722 0.4615
#4  GroupYF    4   14  101.4643 0.4104  0.0000   0.1136 0.3131 0.5164
#5  GroupYF    5   14   72.1071 0.2911 -0.0004   0.0389 0.2351 0.3554
#6  GroupYF    6   14   84.3214 0.3407 -0.0002   0.0808 0.2605 0.4329
#7  GroupYM    1    8  107.4375 0.4347 -0.0002   0.3093 0.2786 0.6060
#8  GroupYM    2    8  113.4375 0.4591  0.0002   0.2607 0.3120 0.6142
#9  GroupYM    3    8   87.3750 0.3532  0.0001   0.1464 0.2471 0.4774
#10 GroupYM    4    8   76.6875 0.3097  0.0002   0.0873 0.2277 0.4070
#11 GroupYM    5    8   92.0625 0.3722 -0.0001   0.2011 0.2487 0.5167
#12 GroupYM    6    8   92.0625 0.3722 -0.0001   0.2011 0.2487 0.5167
#13 GroupNF    1   11  154.3636 0.6255  0.0004   0.2954 0.4500 0.7708
#14 GroupNF    2   11  174.3636 0.7068  0.0007   0.1890 0.5573 0.8188
#15 GroupNF    3   11  162.4545 0.6584 -0.0001   0.1941 0.5130 0.7768
#16 GroupNF    4   11  182.2273 0.7387  0.0007   0.1704 0.5926 0.8425
#17 GroupNF    5   11  146.8182 0.5948  0.0001   0.2383 0.4408 0.7307
#18 GroupNF    6   11  133.5000 0.5407 -0.0018   0.1531 0.4208 0.6555
#19 GroupNM    1    8  126.7500 0.5132 -0.0005   0.3339 0.3422 0.6810
#20 GroupNM    2    8  168.6250 0.6834  0.0005   0.2855 0.5021 0.8196
#21 GroupNM    3    8  176.5000 0.7154  0.0016   0.3400 0.5090 0.8555
#22 GroupNM    4    8  172.5625 0.6994  0.0013   0.3102 0.5066 0.8376
#23 GroupNM    5    8  150.0000 0.6077 -0.0003   0.3585 0.4175 0.7683
#24 GroupNM    6    8  122.8125 0.4972 -0.0025   0.2469 0.3502 0.6447

Example output

Loading required package: MASS
 LD CI Calculations 
 ----------------------- 
 Order of the time and group levels.
 Time level:   1 2 3 4 5 6 
 Group level:   YF YM NF NM 
 The order may be specified in time.order or group.order (does not affect the calculation).

     Group Time Nobs RankMeans    RTE    Bias Variance  Lower  Upper
1  GroupYF    1   14  123.9643 0.5019  0.0008   0.1680 0.3792 0.6243
2  GroupYF    2   14  100.2857 0.4056  0.0003   0.1366 0.3000 0.5222
3  GroupYF    3   14   89.2500 0.3608 -0.0005   0.0979 0.2722 0.4615
4  GroupYF    4   14  101.4643 0.4104  0.0000   0.1136 0.3131 0.5164
5  GroupYF    5   14   72.1071 0.2911 -0.0004   0.0389 0.2351 0.3554
6  GroupYF    6   14   84.3214 0.3407 -0.0002   0.0808 0.2605 0.4329
7  GroupYM    1    8  107.4375 0.4347 -0.0002   0.3093 0.2786 0.6060
8  GroupYM    2    8  113.4375 0.4591  0.0002   0.2607 0.3120 0.6142
9  GroupYM    3    8   87.3750 0.3532  0.0001   0.1464 0.2471 0.4774
10 GroupYM    4    8   76.6875 0.3097  0.0002   0.0873 0.2277 0.4070
11 GroupYM    5    8   92.0625 0.3722 -0.0001   0.2011 0.2487 0.5167
12 GroupYM    6    8   92.0625 0.3722 -0.0001   0.2011 0.2487 0.5167
13 GroupNF    1   11  154.3636 0.6255  0.0004   0.2954 0.4500 0.7708
14 GroupNF    2   11  174.3636 0.7068  0.0007   0.1890 0.5573 0.8188
15 GroupNF    3   11  162.4545 0.6584 -0.0001   0.1941 0.5130 0.7768
16 GroupNF    4   11  182.2273 0.7387  0.0007   0.1704 0.5926 0.8425
17 GroupNF    5   11  146.8182 0.5948  0.0001   0.2383 0.4408 0.7307
18 GroupNF    6   11  133.5000 0.5407 -0.0018   0.1531 0.4208 0.6555
19 GroupNM    1    8  126.7500 0.5132 -0.0005   0.3339 0.3422 0.6810
20 GroupNM    2    8  168.6250 0.6834  0.0005   0.2855 0.5021 0.8196
21 GroupNM    3    8  176.5000 0.7154  0.0016   0.3400 0.5090 0.8555
22 GroupNM    4    8  172.5625 0.6994  0.0013   0.3102 0.5066 0.8376
23 GroupNM    5    8  150.0000 0.6077 -0.0003   0.3585 0.4175 0.7683
24 GroupNM    6    8  122.8125 0.4972 -0.0025   0.2469 0.3502 0.6447

nparLD documentation built on May 2, 2019, 1:09 p.m.

Related to ld.ci in nparLD...