2017-09-15 7 views
0

プールされた回帰における特定の変数の個々の係数をすべて抽出したいという問題が発生しました。プールされたOLS推定から個々のスロープ係数を抽出するR

私のデータは次のようになり、私はplmパッケージと要約を使用してX.

Observation name date Y X 
1   A 1 Y1 X1 
2   A 2 Y2 X2 
3   B 1 Y3 X3 
4   B 2 Y4 X4 

にYを退行、Rだけしかし私にはXの1つの係数を与え、私はでX変数の係数を持つようにしたいですそれぞれの回帰。誰もこれで私を助けることができますか?

明確にするために、私が望むのは、X_nに関連するすべてのベータ版、下の図の1です。混乱させて申し訳ありません。 Pictures for clarification

答えて

0

、あなたは、コード例を参照して、個体あたりOLSモデルに合わせてパッケージplmから機能pvcm()を使用したり、手でそれを行うことができます以下。ヘルプページ(?pvcm)から例の際に

ビル:

library(plm) 
data("Produc", package = "plm") 

form <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp 

pvcm(form, data = Produc, model = "within") # individual coefficients 
## Coefficients: 
##    (Intercept) log(pcap) log(pc) log(emp)  unemp 
## ALABAMA   8.49604 -1.4426440 0.2795010 1.835250 0.00735450 
## ARIZONA   4.66528 -0.1627084 -0.0052207 1.075828 -0.00365798 
## ARKANSAS   3.24565 -0.5056503 0.3212473 1.234017 0.00149221 
## CALIFORNIA   0.27935 0.2639377 0.2484033 0.699135 -0.01074510 
## [...] 

###### same using OLS on splitted data for first entity (ALABAMA): 
l <- split(Produc, Produc$state) 
plm(form, data = l[[1]], model = "pooling") 
## Coefficients: 
## (Intercept) log(pcap)  log(pc) log(emp)  unemp 
## 8.4960384 -1.4426440 0.2795010 1.8352498 0.0073545 
+0

あなたの優雅なソリューションをありがとう! –

0

あなたはplmから結果を取得し、データフレームにそれらを置くためにほうきパッケージからtidyaugmentを使用することができます。あなたは(あなたの例ではnameによって分割)エンティティの異なる係数を使用する場合は

library(plm) 
library(broom) 

data("Grunfeld", package = "plm") 

p <- plm(inv ~ value + capital, data = Grunfeld, model = "pooling") 

df <- tidy(p) 

#   term estimate std.error statistic  p.value 
# 1 (Intercept) -42.7143694 9.51167603 -4.490730 1.207357e-05 
# 2  value 0.1155622 0.00583571 19.802589 9.542703e-49 
# 3  capital 0.2306785 0.02547580 9.054808 1.347370e-16 

df$estimate 

# [1] -42.7143694 0.1155622 0.2306785 

augment(p) 

#   inv value capital  .fitted  .resid 
# 1 317.60 3078.50 2.80 313.689629 3.910371 
# 2 391.80 4661.70 52.60 508.135423 -116.335423 
# 3 410.60 5387.10 156.90 616.023978 -205.423978 
# 4 257.70 2792.20 209.20 328.216223 -70.516223 
# 5 330.80 4313.20 203.40 502.648328 -171.848328 
# 6 461.20 4643.90 207.20 541.741311 -80.541311 
# 7 512.00 4551.20 255.20 542.101267 -30.101267 
# 8 448.00 3244.10 303.70 402.237879 45.762121 
# 9 499.60 4053.70 264.10 486.662133 12.937867 
# 10 547.50 4379.30 201.60 509.871765 37.628235 
+0

こんにちは。まず、あなたの返事をありがとうが、これは私が欲しいものではありません。私はこの結果を得ることができます。私が望むのは、変数 "値"の推定係数です。たとえば、最初に推定された資本係数、資本の第2の推定係数など1つ –

+0

私はあなたの助けをまず感謝します。しかし、これは推定ベータではなく私のデータセットの初期値を与えるだけでしょうか? –

+0

あなたの質問は不明です。おそらくあなたは[再現可能な例](https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example)を提供して、さらに助けてください。 –

関連する問題