# ESCI 197C  Tree movement data 10/23/2026
# This section contains R commands only; see below for R output

twig.move <- c(3.5, 4, 2, 7.5, 4.5, 1.5, 3, 9.5, 5, 2, 8.5, 7.5, 20, 
  4.6, 5.5, 10, 6.5, 6.5, 10, 6, 3, 6.5, 12, 3, 2.5)  # Ren's data
twig.move
length(twig.move)  # what is sample size?
mean(twig.move)
st.dev(twig.move)
sd(twig.move)
min.year <- 30*24*365

# 25 twigs/branch
# 75 branches/tree
25*75  # twigs/tree

mean(twig.move) * min.year * 25 * 75  # ave tree movement/year (cm)
# should sum all twig movements recorded within 2-minute interval to estimate tree movement
sum(twig.move) * min.year * 25 * 75  # tree movement/year (cm)
sum(twig.move) * min.year * 25 * 75 / (100*1000) # tree movement/year (km)
# Compare with Nadkarni's estimate of 300,143 km/yr
sd(twig.move) * min.year * 25 * 75 / (100*1000) # st.dev tree movement/year (km)

#################################################################################

> # ESCI 197C  Tree movement data 10/23/2026
> # This section contains R commands and R output

> twig.move <- c(3.5, 4, 2, 7.5, 4.5, 1.5, 3, 9.5, 5, 2, 8.5, 7.5, 20, 
+   4.6, 5.5, 10, 6.5, 6.5, 10, 6, 3, 6.5, 12, 3, 2.5)  # Ren's data
> twig.move
 [1]  3.5  4.0  2.0  7.5  4.5  1.5  3.0  9.5  5.0  2.0  8.5  7.5 20.0  4.6  5.5
[16] 10.0  6.5  6.5 10.0  6.0  3.0  6.5 12.0  3.0  2.5
> length(twig.move)  # what is sample size?
[1] 25
> mean(twig.move)  # what is average twig movement distance?
[1] 6.184
> sd(twig.move)  # standard deviation in twig measurements
[1] 4.052333
> min.year <- 30*24*365  # Number of 2-minute intervals in one year
> 
> # 25 twigs/branch
> # 75 branches/tree
> 25*75  # twigs/tree
[1] 1875
> 
> mean(twig.move) * min.year * 25 * 75  # ave tree movement/year (cm)
[1] 3047166000
> # should sum all twig movements recorded within 2-minute interval to estimate tree movement
> sum(twig.move) * min.year * 25 * 75  # tree movement/year (cm)
[1] 76179150000
> sum(twig.move) * min.year * 25 * 75 / (100*1000) # ave tree movement/year (km)
[1] 761791.5
> # Compare with Nadkarni's estimate of 300,143 km/yr
> sd(twig.move) * min.year * 25 * 75 / (100*1000) # st.dev tree movement/year (km)
[1] 19967.87