Bill tip change thru time

###all birds were juveniles

---

Stats

bill.tipX<-bill.tip
zz<-table(droplevels(bill.tip$id))
for(i in names(zz[zz!=3]))
bill.tipX<-bill.tipX[bill.tipX$id!=i,]

aov.out = aov(bill.tipX$bill.tip.length~bill.tipX$juvs + Error(bill.tipX$id/bill.tipX$juvs))

summary(aov.out)

## 
## Error: bill.tipX$id
##           Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 10   1.44   0.143               
## 
## Error: bill.tipX$id:bill.tipX$juvs
##                Df Sum Sq Mean Sq F value Pr(>F)  
## bill.tipX$juvs  1  0.201   0.201    5.43  0.042 *
## Residuals      10  0.370   0.037                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Error: Within
##           Df Sum Sq Mean Sq F value Pr(>F)
## Residuals 11  0.106  0.0096

---

The next one compares each year against the next one (0= first year)

with(bill.tipX, pairwise.t.test(bill.tip.length, juvs,
                   p.adjust.method="holm", paired=T,na.rm=T))

## 
##  Pairwise comparisons using paired t tests 
## 
## data:  bill.tip.length and juvs 
## 
##   0    1   
## 1 0.13 -   
## 2 0.13 0.18
## 
## P value adjustment method: holm