```r
oneway.test(d$Lshell~d$Site) 
```

```
## 
## 	One-way analysis of means (not assuming equal variances)
## 
## data:  d$Lshell and d$Site
## F = 9.6559, num df = 5.000, denom df = 31.194, p-value = 1.207e-05
```