我尝试使用lm对这里提到的对数-对数图数据使用mammal dataset来拟合幂函数。我根据描述使用了以下代码,但未能重现结果。
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")]我不明白为什么我的预测值与实测值有很大不同。我做错了什么?
发布于 2021-06-24 16:26:54
博客中有一个错误,正确的答案是:
代谢率= 10^0.5758 *质量^ 0.7063 mammals$MetRate_predict <- 10^0.57584 * (mammals$AdultBodyMass_g ^ 0.70630)
https://stackoverflow.com/questions/68036242
复制相似问题