樹木年輪分析

# 樹木年齡與樹圍資料集 # Orange # 指定精準度 options(digits=4) # 清除殘餘變數 rm(list=ls()) # 前 10 筆資料 head(Orange, n=10) # 改變欄位名稱 colnames(Orange) <- c("Tree", "Age", "Size") # 計算 Orange$Ratio.LINEAR <- Orange$Size / Orange$Age Orange$Ratio.LOG <- log(Orange$Size) / log(Orange$Age) # 分析 Orange.Fit = lm(Size ~ Age, Orange) anova(Orange.Fit) # 繪圖 qplot(x=Age, y=Size, data=Orange, main="樹圍與樹齡關係", xlab="Age", ylab="Size", colour=Size, size=Age)