normal_data <- matrix(rnorm(200*6), 200, 6)
mydesign <- c(0,0,0,1,1,1)
norm_result <- RBM_T(normal_data,mydesign,50,0.05)
unif_data <- matrix(runif(200*7, 0.10, 0.95), 200, 7)
mydesign2 <- c(0,0,0, 1,1,1,1)
unif_result <- RBM_T(unif_data,mydesign2,100,0.05)
normdata_F <- matrix(rnorm(200*9, 0, 2), 200, 9)
mydesign_F <- c(0, 0, 0, 1, 1, 1, 2, 2, 2)
aContrast <- c("X1-X0", "X2-X1", "X2-X0")
normresult_F <- RBM_F(normdata_F, mydesign_F, aContrast, 100, 0.05)
unifdata_F <- matrix(runif(200*18, 0.15, 0.98), 200, 18)
mydesign2_F <- c(rep(0, 6), rep(1, 6), rep(2, 6))
aContrast <- c("X1-X0", "X2-X1", "X2-X0")
unifresult_F <- RBM_F(unifdata_F, mydesign2_F, aContrast, 100, 0.05)
Run the code above in your browser using DataLab