dat <- read.table(text ="con.1 con.2 treat.1 treat.2
                  8100        9599        8144        7000",
                  header = T)

disp <- 0.01776080