README.md

In baolinwu/IPAM: Inference problems in various ANOVA models

IPAM

• An R package for statistical tests in various ANOVA models
• References
• Ling,S., Johnson,C., Ellermann,J., Eberly,L. and Wu,B. (2019) A generalized Tukey method for weighted multiple pairwise comparisons in an ANOVA model. tech rep.
  • Bai,X., Bai,Y., and Wu,B. (2019) Generalized confidence interval calculation for the total variance in a single-factor random-effects ANOVA model with application to medical device comparison problems. tech rep.
• Sample R codes

 ## install the package
 devtools::install_github('baolinwu/IPAM')
 library(IPAM)
 ## Generalized Tukey test for weighted post hoc comparisons
 ## Compute post hoc comparison critical value given a FWER
 GTukey(c(15,10,5), alpha=0.05, W=3:1)
 ##  multiple testing adjusted p-value/FWER 
 PGTukey(rep(6,4), Qx=2.42, ik=1, W=c(5,rep(1,5)) )
 ## CI for total var
 A = rep(1:10, times=5:14)
 s2a = 0.05; s2 = 0.95
 Y = rnorm(10)[A]*sqrt(s2a) + rnorm(length(A))*sqrt(s2)
 BGPsv(Y,A)
 aa = GCIsv(Y,A); aa$GCI
 ## pulse-oximetry data
 data(poData)
 summary(poData$Y); table(poData$A)
 BGPsv(poData$Y,poData$A)
 aa = GCIsv(poData$Y,poData$A); aa$GCI
 hist(aa$Qt, prob=TRUE, breaks=20)
 d = mean(aa$Qt)^2/var(aa$Qt)*2
 a = mean(aa$Qt)/d
 d; a
 curve(dchisq(x/a,d)/a, add=TRUE, col=2)
 y = dchisq(aa$Qt/a,d)
 Q1 = aa$Qt[y>=quantile(y,0.1)]
 range(Q1)

baolinwu/IPAM documentation built on May 9, 2019, 2:21 a.m.