V <- 0.02 # in litre
m <- 2 # in gram
MW.Ad <- 75
M.Ad <- 55.5
Non_T_trt <- c(0,0,0,0,0,0,1,1,1,1,1,1)
T_trt <- c(1,1,1,2,2,2,1,1,1,2,2,2)
Rep <- c(1,2,3,1,2,3,1,2,3,1,2,3)
IGC <- c(2,4)
Ce2 <- c(0.030, 0.031, 0.032, 0.033, 0.034, 0.035, 0.030, 0.031, 0.032, 0.033, 0.034, 0.035)
Ce4 <- c(0.030, 0.031, 0.032, 0.033, 0.034, 0.035, 0.030, 0.031, 0.032, 0.033, 0.034, 0.035)
Ce <- data.frame(Ce2, Ce4)
my.LnKd<- LnKd(V, m, MW.Ad, M.Ad, Non_T_trt, T_trt, Rep, IGC, Ce)
Run the code above in your browser using DataLab