library(nlme)
library(ggplot2)
library(GGally)
library(splines)
library(nlme)
library(fields)
library(lattice)
require(ISLR)
library(grid)
library(dplyr)
library(MASS)
library(mgcv)
library(latticeExtra)
library(fields)

setwd("C:\Users\breni\Google Drive\Acadêmica\Mestrado\Splines")
dados <- read.table("leukemia.txt", header=TRUE)
dados$treatment  <- factor(dados$treatment)
dados$case  <- factor(dados$case)
head(dados)
attach(dados)


interaction <- dados %>%
  dplyr::select(height, age, case) %>%
  group_by(case, age) %>%
  summarise(Average = mean(height))

x11()
ggplot(interaction, aes(x=age, y=Average, colour=treatment, group=case)) +
  ggtitle("Evolução da altura em cm no tempo para cada individuo") +
  geom_line()