x1<-rnorm(120,20,2)
x2<-rnorm(120,100,8)
x3<-rpois(120,10)
x4<-rnbinom(120,mu=10, size=10)
A<-rep(c("a1","a2","a3"), c(40,40,40))
B<-rep(rep(c("b1","b2","b3","b4"), c(10,10,10,10)), times=3)
dat<-data.frame(x1=x1, x2=x2, x3=x3, x4=x4, A=A, B=B)
test<-pairwiseMEP(x=dat,
ep=c("x1","x2","x3", "x4"),
f="A", by="B",
method=c("Param.ratio","Param.ratio","Negbin.ratio","Negbin.ratio"))
test
as.data.frame(test)
as.data.frame(test, whichep=c(1,2))
as.data.frame(test, whichep=c(3,4))
as.data.frame(test, whichep=c("x1","x2"))Run the code above in your browser using DataLab