## [1] 164
## [1] 15633

Are samples equally distributed among sizes?

summary(pdb$Size.x) # 23% more samples from intacts than smalls
## Intact  Large  Small 
##     64     48     52

Community breakdown by groups

##### A sneak peak at comm. composition
#Mean, sd and total abundance per taxa 
total.abund.mean.table<- sapply(pdb[,c(10:15)], FUN = mean, na.rm = TRUE)
tamt<-sort(total.abund.mean.table,decreasing=T)
round(tamt,1)
## Flagellates       Algae     Amoebas    Ciliates    Rotifers     Diatoms 
##        63.0        11.3        10.3         6.8         2.2         1.6
#Quick barplot on mean abundances
barplot(tamt, cex.names=0.7,ylab="Mean abundance",col=colors()[550],yaxt="n", main="Protist abundance artificial treeholes") 
axis(2,cex.axis=0.7)

###Mean sneak peaks

##                                            
##               Leaves.x                     
##               Native         Standard      
##               mean     sd    mean     sd   
##  pdb$Tot.prot 73.79    114.9 119.4    243.3
##                                                   
##               Size.y                              
##               Intact       Large       Small      
##               mean   sd    mean  sd    mean  sd   
##  pdb$Tot.prot 44.86  53.66 109.8 238.5 144.1 225.3

Overview of sample Sorter effects

## Andrea  Angie     KC 
##     79     68     17

We keep Sorter ID as covariate in model/ Angie potentially undersampled flagellates. ###Factors

##                                                             
##                  Leaves.x                                   
##                  Native                Standard             
##                  mean     sd      sum  mean     sd      sum 
##  pdb$Rotifers     2.155     5.060  181  2.2564    4.235  176
##  pdb$Amoebas      9.750    10.317  819 10.0256   16.598  782
##  pdb$Diatoms      2.262     4.311  190  0.9359    2.116   73
##  pdb$Ciliates     6.369     9.146  535  7.3718   25.858  575
##  pdb$Flagellates 43.143   107.800 3624 85.9231  243.163 6702
##  pdb$Algae       10.107    15.061  849 12.8718   13.404 1004
##                                                              
##                  Size.x                                      
##                  Intact        Large          Small          
##                  mean   sd     mean   sd      mean    sd     
##  pdb$Rotifers     1.547  3.075  2.396   3.529   2.942   6.717
##  pdb$Amoebas      9.297 10.992 10.250  12.393  11.558  18.421
##  pdb$Diatoms      2.031  3.908  1.000   2.173   1.654   3.829
##  pdb$Ciliates     3.688  4.646  5.604  10.027  11.808  31.449
##  pdb$Flagellates 21.422 52.946 74.875 233.804 103.269 227.985
##  pdb$Algae        6.875  7.600 15.688  19.044  12.827  14.089

# STATS

##        chisq        ratio          rdf            p 
## 192.67494212   1.45965865 132.00000000   0.00045216
## Analysis of Deviance Table (Type II tests)
## 
## Response: Tot.prot
##                         Df   Chisq Pr(>Chisq)    
## Leaves.x                 1 12.3818  0.0004335 ***
## Size.x                   2 26.1606  2.086e-06 ***
## Sorter                   2  4.7928  0.0910463 .  
## Leaves.x:Size.x          2  7.1217  0.0284144 *  
## Leaves.x:Sorter          2  4.3267  0.1149387    
## Size.x:Sorter            4  8.3687  0.0789691 .  
## Leaves.x:Size.x:Sorter   4  5.3304  0.2550456    
## Residuals              131                       
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##   Sorter Size.x TH.length
## 1 Andrea Intact        31
## 2 Andrea  Large        24
## 3 Andrea  Small        24
## 4  Angie Intact        27
## 5  Angie  Large        17
## 6  Angie  Small        24
## 7     KC Intact         6
## 8     KC  Large         7
## 9     KC  Small         4

Do predators have an effect?

No apparent effect of predator p/a on protists.

Response to protist filter feeder culicids (all)

Response to Detritus type

Protist groups correlated with decomposition?

##             Rotifers Amoebas  Algae Diatoms Flagellates Ciliates Mass.loss
## Rotifers       1.000   0.351  0.076  -0.039      -0.025    0.051        NA
## Amoebas        0.351   1.000 -0.002  -0.025      -0.031   -0.014        NA
## Algae          0.076  -0.002  1.000   0.036       0.055   -0.031        NA
## Diatoms       -0.039  -0.025  0.036   1.000      -0.090   -0.011        NA
## Flagellates   -0.025  -0.031  0.055  -0.090       1.000    0.010        NA
## Ciliates       0.051  -0.014 -0.031  -0.011       0.010    1.000        NA
## Mass.loss         NA      NA     NA      NA          NA       NA         1

Among protists interactions