我正在绘制一条河流的水平和流量之间的图表。然后,我想在散点上拟合一条直线。
我的数据框df1有4704行和2列 – 流量和水平。
fit <- lm(df1$Level ~ df1$Discharge)xx <- seq(0, max(df1$Discharge), length = 4704)plot(x = df1$Discharge, y = df1$Level, pch = 19)lines(xx, predict(fit, data.frame(x = xx)), col = "red")
我得到的不是一条直线,而是其他东西。
我想知道,有没有办法得到一条适合这些散点的直线?请指导我。
回答:
如果你能提供数据会更好,从我所看到的,没有提供包含正确变量名的数据框进行预测可能会导致问题。如果我用一些数据重新运行你的代码:
set.seed(100)df1 = data.frame(Level=runif(4704))df1$Discharge=20*df1$Level + runif(4704)fit <- lm(Level ~ Discharge, data = df1)xx <- seq(0, max(df1$Discharge), length = 4704)plot(x = df1$Discharge, y = df1$Level, pch = 19)
这一行会抛出错误:
lines(xx, predict(fit, data.frame(x = xx)), col = "red")Error in eval(predvars, data, env) : object 'Discharge' not found
如果你更改列名,它会工作:
lines(xx, predict(fit, data.frame(Discharge = xx)), col = "red")
请提供一个可复现的示例。