Data.frame和不同規格的線性迴歸模型

這是我的data.frame：Data.frame和不同規格的線性迴歸模型

data <- matrix(rnorm(50*5),nrow=50) 
m <- data.frame(data) 
m 

      X1   X2   X3   X4 
1 -0.47903358 1.92799699 -0.584364168 -1.475276350 
2 0.05464517 1.90064618 1.721449935 1.215188405 
3 2.23782855 0.16839367 -0.074212064 -0.719745849 
4 0.44476924 0.37247676 0.300423948 1.722222410 
5 -0.66889333 0.26094142 -0.001254037 0.617630295 
6 -1.17663634 -1.93020031 -2.813433555 0.350943310 
7 2.16609882 -0.06594556 0.133278642 1.693078802 
8 -0.16090564 -1.14454452 -1.340407149 0.574943347 
9 0.46691283 0.98394614 -0.277273498 0.030499576 
10 1.19049154 1.03160103 0.626364442 0.429348535 
11 0.72593476 0.89130574 -0.599052058 2.
12 -0.04748923 -1.38044706 0.058364936 0.825752213 
13 -0.90407306 0.14790680 0.002608154 -0.991979288 
14 -1.50113350 0.78807179 -2.371160466 0.559174890 
15 0.04075635 -1.14486991 1.585769771 -0.495937216 
16 1.29724449 0.47729721 -0.871033202 1.524746619 
17 1.12384656 -0.07589516 -0.463837527 0.206108262 
18 0.65092671 0.29660588 -1.096981115 -0.326524160 
19 -0.54953490 -0.06792527 -0.543922865 0.211205138 
20 0.19187346 -0.65190412 -1.674653668 0.712704645 
21 0.97108788 0.91087005 0.182424518 -0.296413534 
22 1.64975949 0.22140448 -0.075355571 -0.472299544 
23 1.45613296 -1.29101443 1.231839302 0.774475499 
24 -0.15324579 -1.31596417 -0.343467527 0.087122108 
25 0.20305822 -0.96867715 -0.311745324 -1.483964874 
26 -2.22916775 1.46982603 -0.254705650 -0.003573105 
27 0.26135095 -0.44053524 0.208284330 0.459095264 
28 0.55020951 -0.37742141 0.346186077 0.825811292 
29 -0.78695097 -0.21014858 0.003055816 2.142121347 
30 -0.33506969 0.31572240 -0.139970552 1.344697176 
31 -0.94164226 -0.88366990 -0.906886069 0.774590884 
32 3.23764133 -0.47556802 0.166917302 0.959307700 
33 -0.64114926 0.16089400 1.942563410 0.142073261 
34 0.03848502 -0.83066179 0.934403254 1.730448144 
35 -0.34621698 0.65682342 1.601365990 0.160012282 
36 0.62848314 0.33126945 0.188823549 2.380319403 
37 1.02807329 1.07299920 0.398083643 -0.208241920 
38 -0.50616984 1.68998408 0.427342040 -0.345264684 
39 1.31721397 -1.17111247 1.754036030 0.755414080 
40 1.23714888 -0.12253370 0.887065669 0.001770354 
41 0.59541135 0.66438057 0.880643498 1.441513269 
42 0.41100680 0.52757460 1.219644246 1.450262235 
43 0.08984815 1.44830794 0.113193523 -0.620546777 
44 -1.51411951 -0.17346499 -0.698924584 -0.767921061 
45 -1.03652530 1.27296507 1.740954529 -0.467568299 
46 1.14533381 -0.67304957 -0.384.637582075 
47 0.37996214 0.53566774 -1.413930797 1.222728459 
48 -0.53705734 0.77591264 0.375236588 -1.430924067 
49 0.64051195 -0.22562951 -1.667050192 0.345009059 
50 -1.19521122 1.18436746 -0.905004199 -0.005221906

好了，我想現在估計3個不同的模型。

X1 = X2 + X3 + X4 
X1 = X2^2 + X3(Lag1) + X4^2 
ln(X1)=X2^2 + X3(Lag1) + X4^2

作爲即時通訊的一個先行者，我不能實現這個功能。我正在閱讀lm函數，但仍無法實現這些模型。我能夠實現的方式非常基本，即每個模型都有一個data.frame。

謝謝

來源

2015-10-16 Laorie

要使用'lm（）'，你需要一個公式來回歸。您的結果數據必須是數據框中的一列。公式中的「t」值應該是數據框中要使用的常量。一旦你有了這個，你可以創建一個字符串你的公式和退步。所以，Y需要在您的data.frame中，並且需要將您作爲data.frame以及您的x值。一個例子公式可以是字符串：'Y〜X1 + X2 + X3'。 '〜'圖標可以等效爲等號。 – giraffehere

感謝您解答長頸鹿。我編輯了這個問題。 – Laorie

您還需要列爲X3的滯後變量，但我會寫出如何執行此操作的答案。請注意，您需要爲迴歸中的每個變量添加一列。 – giraffehere

要使用lm（），您至少需要兩件事。你需要一個公式，你需要數據來應用該公式。在你的情況下，你將如何編寫代碼來對你的數據運行一個非常簡單的，基本的線性迴歸。再次，您需要在數據中使用「X3_Lag1」變量（實際上是一個稱爲滯後應用的列），以對滯後數據進行迴歸（創建一個X3向下按一行的新列），然後剪切數據以找到該缺失值在第一排。

# Getting the first model. 
formula1 <- "X1 ~ X2 + X3 + X4" 
lin1 <- lm(formula1, data) 
summary(lin1) # This will give you the coefficients for each value. 

# Getting the second model. 
formula2 <- "X1 ~ X2^2 + X3_Lag1 + X4^2" 
lin2 <- lm(formula2, data) 
summary(lin2) 

# Getting the third model. 
formula3 <- "log(X1) ~ X2^2 + X3_Lag1 + X4^2" 
lin3 <- lm(formula3, data) 
summary(lin3)

警告：我相信你會得到第三個公式的錯誤，因爲你正在取一個負數的自然對數，這是不可能的。請注意這種情況。

來源

2015-10-16 21:22:07 giraffehere

我在想如何在沒有包含列的情況下做到這一點。我做了這個'lm（X1〜X2 + I（X3^2），m）'。我不需要包含x3^2列。和lm（log（X1）〜X2 + I（X3^2），m）'。滯後的變量是問題。 – Laorie

您不需要該運算符的「X3^2」列，該運算符將起作用，從而爲您提供X3列中每個數字的平方。然而，對於滯後的變量，就我所知，您將需要一個「X3_Lag1」列，您需要在數據中創建這一列。據我所知，在公式本身中沒有創建滯後的方法，您必須事先在公式之外和數據中進行。我對log（X1）的關注是你的X1列中有負值，這在對數中是不允許的。 – giraffehere

Data.frame和不同規格的線性迴歸模型

回答

相關問題