##############################################
# Analysis of Longitudinal Data (SoSe 2017)  #
#                                            #
# Solution for tutorial 2                    #  
# Date: 18.05.2017                          #
##############################################

###############
# Exercise 1 #
###############

library(nlme)
load("rats.long.RData")
str(rats.long)

# Transformation of TIME:
rats.long$logT     <- log(1+(rats.long$TIME-45)/10) 

# take out NAs
rats.long          <- rats.long[complete.cases(rats.long),]