我有一个数据框，其中有几列和多行。例如，下面的数据框显示第1列到第4列是我的土壤湿度数据，5到8是蒸发蒸腾（ET），第9列到第12列是降水数据（p）。我有156 Colum of soil moisture，ET，P。现在我想用土壤水分列进行多元线性回归到相应的降水和蒸发蒸发柱，其中我的土壤水分是一个因变量，ET和P是自变量。我想在运行多元回归时使用以下等式作为我的自变量，并计算a1，a2和b1，b2等系数。
ΔΦ(j.t+1) = a1P(j.t) + a2P(J.t-1) + b1ET(j.t) + b2ET(J.t-1)+εt+1
其中Δθ(j,t+1)是t和t+1之间的土壤湿度差异，a和b是要确定降水的系数（{{1} }）和蒸发蒸腾（P）和ET是模型的误差项。
我的数据框架如下：
1 2 3 4 5 6 7 8 9 10 11 12
NA NA NA NA楠楠楠楠楠楠楠楠NA NA NA NA 645 732 1010 1050 169.7152863 264.919281 174.5597839 117.1220169 1.788969846 1.769300308 1.881795708 1.798497618楠楠楠楠183.5778961 272.1306458 164.4124298 94.27450562 3.41084 3.41084 NA NA 576 657 840 817 166.3976135 263.9378357 161.8202972 97.60928345 1.995551665 2.037976549 2.144578129 2.065179551 533 597 747 744 151.7550201 265.6147766 155.722702 91.79171753
到目前为止，我使用以下代码以R编程语言创建数据框和矩阵，并希望有更多帮助来完成任务。
εt+1答案 0 :(得分：0)
线性回归函数lm()采用以下参数：
lm(dependentVariable ~ independentVariable1 + independentVariable2, data = df)
其他人可能会指出解决方法，但我认为这意味着所有变量（列）必须位于一个数据框中。所以，你的第一行很好：
data <- read.table("MasterPiece.csv", header = TRUE, sep = ",")
（注意：建议在分配时使用<-代替=）。您不需要将变量放在单独的数据框中。
例如，如果您的变量名称为：smdi1 ...，modet1 ...，prcp1 ...您可以使用以下命令运行模型：
m1 <- lm(smdi1 ~ modet1 + prcp1, data = data)
summary(m1)