X <- read.csv("mice.csv", header = TRUE)     # replace by your own path here
boxplot(X[,1] ~ X[,2])

lm1 <- lm(X[,1] ~ X[,2])     

 summary(lm1)              # "milking" information from lm1
 anova(lm1)

  plot(lm1$resid)

  plot(lm1$fit, lm1$resid)   # residuals vs means

   qqnorm(lm1$res)

   # experiments with different shapes


    #1 
    n <- 1000
    X <- runif(n, 2, 12)
    hist(X)
    
    qqnorm(X)

    #2 
    X <- rt(n, df = 3)
    hist(X)

    qqnorm(X)