我尝试使用 lm 在此处提到的对数日志图数据上拟合幂函数https://statisticsbyjim.com/regression/log-log-plots/使用哺乳动物数据集。我根据描述使用了以下代码,但无法重现结果。
mammals <- read.csv("mammals.csv", header=TRUE, stringsAsFactors = FALSE)
plot(log10(mammals$AdultBodyMass_g), log10(mammals$BasalMetRate_mLO2hr)) # log-log plot
lmMammals <- lm(log10(BasalMetRate_mLO2hr) ~ log10(AdultBodyMass_g), data=mammals)
summary(lmMammals)
# Metabolic Rate = 0.5758 Mass ^ 0.7063
mammals$MetRate_predict <- 0.57584 * (mammals$AdultBodyMass_g ^ 0.70630)
mammals[1:5, c("BasalMetRate_mLO2hr", "MetRate_predict")]
我不明白为什么我的预测值与测量值有很大不同。我究竟做错了什么?