0

我尝试使用 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")]

我不明白为什么我的预测值与测量值有很大不同。我究竟做错了什么?

4

1 回答 1

0

博客中有错误正确答案是:

代谢率 = 10^0.5758 * 质量 ^ 0.7063 哺乳动物$MetRate_predict <- 10^0.57584 * (哺乳动物$AdultBodyMass_g ^ 0.70630)

于 2021-06-24T08:26:54.177 回答