fit = lm(time~cyan)
histogram(~residuals(fit)|cyan,type = "density") #lattice library

qqnorm(residuals(fit))
qqline(residuals(fit))

t.test(time~cyan,mu=0,alt = "two.sided",conf = 0.95,var.eq = T,paired = F)