22、温度数据与生物数据的分析与验证

温度数据与生物数据的分析与验证

1. 室内外温度分析

室内地下室温度通常在14°C至23°C之间波动。不过，部分数据会超出这个范围，比如夏季高温表明房屋有供暖系统但无空调，冬季的低温峰值可能是开窗导致的。室外最低温度可达 -20°C，而室内最低温度一般高于10°C。2003年9月和10月左右似乎出现了异常情况，可能是供暖系统在此期间发生了变化。

1.1 室外温度的季节性趋势处理

为了识别室外温度中意外的冷暖测量值，我们可以利用年度温度的周期性。由于年度温度大致呈周期性变化，我们可以用二次多项式来拟合数据，以去除季节性趋势。

# Model the data as a second order polynomial 
year.model <- lm(outside ~ poly(timestamp, 2), data = glarp) 
# Display the regression and the data 
ggplot(glarp, aes(x=timestamp)) + clean_theme + 
  geom_line(aes(y = outside), color = "gray") + 
  geom_line(aes(y = predict(year.model)),  
            color = "darkred", size = 2) + 
  ggtitle("Outside temperature versus polynomial fit") 

通过上述代码，我们可以得到室外温度与多项式拟合的图形。从图中可以看出，年度