Graphs for Tezonteopan data

DATA <- read.csv("~/Documents/R/DATA.csv")
GandR.RUC.Age <- 0.06*DATA$RUC.P.Crown
GandR.RULI.Age <- (0.035*DATA$RULI.P.Crown)+1
GandR.RUCI.Age <- 0.045*DATA$RUCI.P.Crown
GandR.LUCI.Age <- 0.045*DATA$LUCI.P.Crown
GandR.LULI.Age <- (0.035*DATA$LUL.P.Crown)+1
GandR.LUC.Age <- 0.06*DATA$LUC.P.Crown
GandR.RLC.Age <- (0.06*DATA$RLC.P.Crown)+0.5
GandR.RLLI.Age <- 0.04*DATA$RLLI.P.Crown
GandR.RLCI.Age <- 0.04*DATA$RLCI.P.Crown
GandR.LLCI.Age <- 0.04*DATA$LLCI.P.Crown
GandR.LLLI.Age <- 0.04*DATA$LLLI.P.Crown
GandR.LLC.Age <- (0.06*DATA$LLC.P.Crown)+0.5
GoodmanRose <- cbind(GandR.RUC.Age, GandR.RULI.Age, GandR.RUCI.Age, GandR.LUCI.Age, GandR.LULI.Age, GandR.LUC.Age, GandR.RLC.Age, GandR.RLLI.Age, GandR.RLCI.Age, GandR.LLCI.Age, GandR.LLLI.Age, GandR.LLC.Age)

RDNE.RUC.Age <- 626.8*2.718^(0.011*DATA$RUC.P.Crown)/365
RDNE.RULI.Age <- 646.9*2.718^(0.010*DATA$RULI.P.Crown)/365
RDNE.RUCI.Age <- 418.4*2.718^(0.015*DATA$RUCI.P.Crown)/365
RDNE.LUCI.Age <- 418.4*2.718^(0.015*DATA$LUCI.P.Crown)/365
RDNE.LULI.Age <- 646.9*2.718^(0.010*DATA$LUL.P.Crown)/365
RDNE.LUC.Age <- 626.8*2.718^(0.011*DATA$LUC.P.Crown)/365
RDNE.RLC.Age <- 544.6*2.718^(0.014*DATA$RLC.P.Crown)/365
RDNE.RLLI.Age <- 358.3*2.718^(0.014*DATA$RLLI.P.Crown)/365
RDNE.RLCI.Age <- 353.3*2.718^(0.014*DATA$RLCI.P.Crown)/365
RDNE.LLCI.Age <- 353.3*2.718^(0.014*DATA$LLCI.P.Crown)/365
RDNE.LLLI.Age <- 358.3*2.718^(0.014*DATA$LLLI.P.Crown)/365
RDNE.LLC.Age <- 544.6*2.718^(0.014*DATA$LLC.P.Crown)/365
ReidDeanNE <- cbind(RDNE.RUC.Age, RDNE.RULI.Age, RDNE.RUCI.Age, RDNE.LUCI.Age, RDNE.LULI.Age, RDNE.LUC.Age, RDNE.RLC.Age, RDNE.RLLI.Age, RDNE.RLCI.Age, RDNE.LLCI.Age, RDNE.LLLI.Age, RDNE.LLC.Age)

GandRC.RUC.Age <- (0.0578*DATA$RUC.P.Crown)+0.22
GandRC.RULI.Age <- (0.035*DATA$RULI.P.Crown)+1
GandRC.RUCI.Age <- (0.0406*DATA$RUCI.P.Crown)+0.44
GandRC.LUCI.Age <- (0.0406*DATA$LUCI.P.Crown)+0.44
GandRC.LULI.Age <- (0.035*DATA$LUL.P.Crown)+1
GandRC.LUC.Age <- (0.0578*DATA$LUC.P.Crown)+0.22
GandRC.RLC.Age <- (0.0559*DATA$RLC.P.Crown)+0.91
GandRC.RLLI.Age <- (0.0382*DATA$RLLI.P.Crown)+0.18
GandRC.RLCI.Age <- (0.0354*DATA$RLCI.P.Crown)+0.46
GandRC.LLCI.Age <- (0.0354*DATA$LLCI.P.Crown)+0.46
GandRC.LLLI.Age <- (0.0382*DATA$LLLI.P.Crown)+0.18
GandRC.LLC.Age <- (0.0559*DATA$LLC.P.Crown)+0.91
GoodmanRoseC <- cbind(GandRC.RUC.Age, GandRC.RULI.Age, GandRC.RUCI.Age, GandRC.LUCI.Age, GandRC.LULI.Age, GandRC.LUC.Age, GandRC.RLC.Age, GandRC.RLLI.Age, GandRC.RLCI.Age, GandRC.LLCI.Age, GandRC.LLLI.Age, GandRC.LLC.Age)

RUC <- list(RDNE.RUC.Age, GandR.RUC.Age, GandRC.RUC.Age)
names(RUC)[1] <- "RDNE"
names(RUC)[2] <- "GandR"
names(RUC)[3] <- "GandRC"

RULI <- list(RDNE.RULI.Age, GandR.RULI.Age, GandRC.RULI.Age)
names(RULI)[1] <- "RDNE"
names(RULI)[2] <- "GandR"
names(RULI)[3] <- "GandRC"

RUCI <- list(RDNE.RUCI.Age, GandR.RUCI.Age, GandRC.RUCI.Age)
names(RUCI)[1] <- "RDNE"
names(RUCI)[2] <- "GandR"
names(RUCI)[3] <- "GandRC"

LUCI <- list(RDNE.LUCI.Age, GandR.RUCI.Age, GandRC.RUCI.Age)
names(LUCI)[1] <- "RDNE"
names(LUCI)[2] <- "GandR"
names(LUCI)[3] <- "GandRC"

LULI <- list(RDNE.LULI.Age, GandR.LULI.Age, GandRC.LULI.Age)
names(LULI)[1] <- "RDNE"
names(LULI)[2] <- "GandR"
names(LULI)[3] <- "GandRC"

LUC <- list(RDNE.LUC.Age, GandR.LUC.Age, GandRC.LUC.Age)
names(LUC)[1] <- "RDNE"
names(LUC)[2] <- "GandR"
names(LUC)[3] <- "GandRC"

RLC <- list(RDNE.RLC.Age, GandR.RLC.Age, GandRC.RLC.Age)
names(RLC)[1] <- "RDNE"
names(RLC)[2] <- "GandR"
names(RLC)[3] <- "GandRC"

RLLI <- list(RDNE.RLLI.Age, GandR.RLLI.Age, GandRC.RLLI.Age)
names(RLLI)[1] <- "RDNE"
names(RLLI)[2] <- "GandR"
names(RLLI)[3] <- "GandRC"

RLCI <- list(RDNE.RLCI.Age, GandR.RLCI.Age, GandRC.RLCI.Age)
names(RLCI)[1] <- "RDNE"
names(RLCI)[2] <- "GandR"
names(RLCI)[3] <- "GandRC"

LLCI <- list(RDNE.LLCI.Age, GandR.LLCI.Age, GandRC.LLCI.Age)
names(LLCI)[1] <- "RDNE"
names(LLCI)[2] <- "GandR"
names(LLCI)[3] <- "GandRC"

LLLI <- list(RDNE.LLLI.Age, GandR.LLLI.Age, GandRC.LLLI.Age)
names(LLLI)[1] <- "RDNE"
names(LLLI)[2] <- "GandR"
names(LLLI)[3] <- "GandRC"

LLC <- list(RDNE.LLC.Age, GandR.LLC.Age, GandRC.LLC.Age)
names(LLC)[1] <- "RDNE"
names(LLC)[2] <- "GandR"
names(LLC)[3] <- "GandRC"

par(font.axis=2)
par(font.lab=2)

boxplot(RLC, xlim=c(0,8), ylim=c(0,6), main="Lower Canine", xlab="Right                               Left", ylab="Age at Formation", col=c("red","green","orange"))
boxplot(LLC, add=T, at=c(5,6,7), col=c("red","green","orange"))

plot of chunk unnamed-chunk-17

boxplot(RLLI, xlim=c(0,8), ylim=c(0,6), main="Lower Lateral Incisor", xlab="Right                               Left", ylab="Age at Formation", col=c("red","green","orange"))
boxplot(LLLI, add=T, at=c(5,6,7), col=c("red","green","orange"))

plot of chunk unnamed-chunk-18

boxplot(RLCI, xlim=c(0,8), ylim=c(0,6), main="Lower Central Incisor", xlab="Right                               Left", ylab="Age at Formation", col=c("red","green","orange"))
boxplot(LLCI, add=T, at=c(5,6,7), col=c("red","green","orange"))

plot of chunk unnamed-chunk-19

boxplot(RUC, xlim=c(0,8), ylim=c(0,6), main="Upper Canine", xlab="Right                               Left", ylab="Age at Formation", col=c("red","green","orange"))
boxplot(LUC, add=T, at=c(5,6,7), col=c("red","green","orange"))

plot of chunk unnamed-chunk-20

boxplot(RULI, xlim=c(0,8), ylim=c(0,6), main="Upper Lateral Incisor", xlab="Right                               Left", ylab="Age at Formation", col=c("red","green","orange"))
boxplot(LULI, add=T, at=c(5,6,7), col=c("red","green","orange"))

plot of chunk unnamed-chunk-21

boxplot(RUCI, xlim=c(0,8), ylim=c(0,6), main="Upper Central Incisor", xlab="Right                               Left", ylab="Age at Formation", col=c("red","green","orange"))
boxplot(LUCI, add=T, at=c(5,6,7), col=c("red","green","orange"))

plot of chunk unnamed-chunk-22

RDNE.estimates <- read.csv("~/Documents/R/RDNE estimates.csv")

boxplot(LLC.age.est ~ DATA.N.S, data=RDNE.estimates, main="Left Lower Canine", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-24

boxplot(LLLI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Left Lower Lateral Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-25

boxplot(LLCI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Left Lower Central Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-26

boxplot(RLCI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Right Lower Central Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-27

boxplot(RLLI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Right Lower Lateral Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-28

boxplot(RLC.age.est ~ DATA.N.S, data=RDNE.estimates, main="Right Lower Canine", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-29

boxplot(RUC.age.est ~ DATA.N.S, data=RDNE.estimates, main="Right Upper Canine", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-30

boxplot(RULI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Right Upper Lateral Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-31

boxplot(RUCI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Right Upper Central Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-32

boxplot(LUCI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Left Upper Central Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-33

boxplot(LULI.age.est ~ DATA.N.S, data=RDNE.estimates, main="Left Upper Lateral Incisor", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-34

boxplot(LUC.age.est ~ DATA.N.S, data=RDNE.estimates, main="Left Upper Canine", xlab="Supplementation Status", ylab="Age at Formation", ylim=c(0,5), col=c("springgreen1", "tan1"))

plot of chunk unnamed-chunk-35

boxplot(LLC.age.est ~ sex, data=RDNE.estimates, main="Left Lower Canine", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-36

boxplot(LLLI.age.est ~ sex, data=RDNE.estimates, main="Left Lower Lateral Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-37

boxplot(LLCI.age.est ~ sex, data=RDNE.estimates, main="Left Lower Cental Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-38

boxplot(RLCI.age.est ~ sex, data=RDNE.estimates, main="Right Lower Central Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-39

boxplot(RLLI.age.est ~ sex, data=RDNE.estimates, main="Right Lower Lateral Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-40

boxplot(RLC.age.est ~ sex, data=RDNE.estimates, main="Right Lower Canine", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-41

boxplot(RUC.age.est ~ sex, data=RDNE.estimates, main="Right Upper Canine", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-42

boxplot(RULI.age.est ~ sex, data=RDNE.estimates, main="Right Upper Lateral Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-43

boxplot(RUCI.age.est ~ sex, data=RDNE.estimates, main="Right Upper Central Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-44

boxplot(LUCI.age.est ~ sex, data=RDNE.estimates, main="Left Upper Central Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-45

boxplot(LULI.age.est ~ sex, data=RDNE.estimates, main="Left Upper Lateral Incisor", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-46

boxplot(LUC.age.est ~ sex, data=RDNE.estimates, main="Left Upper Canine", xlab="Sex", ylab="Age at Formation", ylim=c(0,5.5), col=c("darkorchid3", "forestgreen"))

plot of chunk unnamed-chunk-47

quantile(RDNE.estimates$Height.Percentile, na.rm=TRUE)

##   0%  25%  50%  75% 100% 
##    1    2    3   15   94

height.quantile <- subset(RDNE.estimates, Height.Percentile <= 2 | Height.Percentile >=15)
height.quantile$Height.Category[height.quantile$Height.Percentile <=2] <- "lower"
height.quantile$Height.Category[height.quantile$Height.Percentile >=15] <- "upper"

boxplot(LLC.age.est ~ Height.Category, data=height.quantile, main="Left Lower Canine", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-49

boxplot(LLLI.age.est ~ Height.Category, data=height.quantile, main="Left Lower Lateral Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-50

boxplot(LLCI.age.est ~ Height.Category, data=height.quantile, main="Left Lower Central Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-51

boxplot(RLCI.age.est ~ Height.Category, data=height.quantile, main="Right Lower Central Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-52

boxplot(RLLI.age.est ~ Height.Category, data=height.quantile, main="Right Lower Lateral Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-53

boxplot(RLC.age.est ~ Height.Category, data=height.quantile, main="Right Lower Canine", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-54

boxplot(RUC.age.est ~ Height.Category, data=height.quantile, main="Right Upper Canine", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-55

boxplot(RULI.age.est ~ Height.Category, data=height.quantile, main="Right Upper Lateral Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-56

boxplot(RUCI.age.est ~ Height.Category, data=height.quantile, main="Right Upper Central Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-57

boxplot(LUCI.age.est ~ Height.Category, data=height.quantile, main="Left Upper Central Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-58

boxplot(LULI.age.est ~ Height.Category, data=height.quantile, main="Left Upper Lateral Incisor", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-59

boxplot(LUC.age.est ~ Height.Category, data=height.quantile, main="Left Upper Canine", sub="Height", xlab="Quartile", ylab="Age at Formation", col=c("dodgerblue", "firebrick"))

plot of chunk unnamed-chunk-60

quantile(RDNE.estimates$Weight.Percentile, na.rm=TRUE)

##   0%  25%  50%  75% 100% 
##  1.0 12.0 29.5 52.0 90.0

weight.quantile <- subset(RDNE.estimates, Weight.Percentile <=12 | Weight.Percentile >=52)
weight.quantile$Weight.Category[weight.quantile$Weight.Percentile <=12] <- "lower"
weight.quantile$Weight.Category[weight.quantile$Weight.Percentile >=52] <- "upper"

boxplot(LLC.age.est ~ Weight.Category, data=weight.quantile, main="Left Lower Canine", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-62

boxplot(LLLI.age.est ~ Weight.Category, data=weight.quantile, main="Left Lower Lateral Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-63

boxplot(LLCI.age.est ~ Weight.Category, data=weight.quantile, main="Left Lower Central Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-64

boxplot(RLCI.age.est ~ Weight.Category, data=weight.quantile, main="Right Lower Central Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-65

boxplot(RLLI.age.est ~ Weight.Category, data=weight.quantile, main="Right Lower Lateral Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-66

boxplot(RLC.age.est ~ Weight.Category, data=weight.quantile, main="Right Lower Canine", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-67

boxplot(RUC.age.est ~ Weight.Category, data=weight.quantile, main="Right Upper Canine", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-68

boxplot(RULI.age.est ~ Weight.Category, data=weight.quantile, main="Right Upper Lateral Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-69

boxplot(RUCI.age.est ~ Weight.Category, data=weight.quantile, main="Right Upper Central Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-70

boxplot(LUCI.age.est ~ Weight.Category, data=weight.quantile, main="Left Upper Central Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-71

boxplot(LULI.age.est ~ Weight.Category, data=weight.quantile, main="Left Upper Lateral Incisor", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-72

boxplot(LUC.age.est ~ Weight.Category, data=weight.quantile, main="Left Upper Canine", sub="Weight", xlab="Quartile", ylab="Age at Formation", ylim=c(0,5), col=c("forestgreen", "gold"))

plot of chunk unnamed-chunk-73