thr3ads.net - R help - [R] How can I extract part of the data in a panel dataset? [Apr 2013]

If this information is useful, please help other people find it:
Share via:

jpm miao

2013-Apr-08 14:29 UTC

[R] How can I extract part of the data in a panel dataset?

Taking the Grunfeld data, which is built-in in R, for example,

(1)How can I construct a dataset (or dataframe) that consists of the data
of all firms in 1951?

(2)How can I calculate the average capital in each form over the period
1951-1954?
What I can imagine is to categorize the data by firm, and then select the
data between 1951 and 1954 for each firm, but how can I do it?

Thanks,

Miao

Grunfeld data:
  invest value capital firm year  317.6 3078.5 2.8 General Motors 1935
391.8 4661.7 52.6 General Motors 1936  410.6 5387.1 156.9 General Motors
1937  257.7 2792.2 209.2 General Motors 1938  330.8 4313.2 203.4 General
Motors 1939  461.2 4643.9 207.2 General Motors 1940  512 4551.2 255.2 General
Motors 1941  448 3244.1 303.7 General Motors 1942  499.6 4053.7 264.1 General
Motors 1943  547.5 4379.3 201.6 General Motors 1944  561.2 4840.9 265 General
Motors 1945  688.1 4900.9 402.2 General Motors 1946  568.9 3526.5 761.5 General
Motors 1947  529.2 3254.7 922.4 General Motors 1948  555.1 3700.2
1020.1 General
Motors 1949  642.9 3755.6 1099 General Motors 1950  755.9 4833 1207.7 General
Motors 1951  891.2 4924.9 1430.5 General Motors 1952  1304.4 6241.7
1777.3 General
Motors 1953  1486.7 5593.6 2226.3 General Motors 1954  209.9 1362.4 53.8 US
Steel 1935  355.3 1807.1 50.5 US Steel 1936  469.9 2676.3 118.1 US Steel
1937  262.3 1801.9 260.2 US Steel 1938  230.4 1957.3 312.7 US Steel 1939
361.6 2202.9 254.2 US Steel 1940
...........

	[[alternative HTML version deleted]]

arun

2013-Apr-08 18:04 UTC

head link

[R] How can I extract part of the data in a panel dataset?

Hi,
1)
Grunfeld1951<-Grunfeld[Grunfeld$year==1951,]
?Grunfeld1951
#??? firm year??? inv?? value capital
#17???? 1 1951 755.90 4833.00 1207.70
#37???? 2 1951 588.20 2289.50? 342.10
#57???? 3 1951 135.20 1819.40? 671.30
#77???? 4 1951 160.62? 809.00? 203.50
#97???? 5 1951? 80.30? 327.30? 683.90
#117??? 6 1951? 95.30? 676.90? 177.20
#137??? 7 1951? 64.77? 179.00? 450.60
#157??? 8 1951? 54.38? 723.80? 129.70
#177??? 9 1951? 56.49? 339.90? 391.00
#197?? 10 1951?? 4.67?? 83.04??? 9.07

2)


?Grunfeld51_54<-Grunfeld[Grunfeld$year%in% 1951:1954,]
?head(Grunfeld51_54)
?#? firm year??? inv? value capital
#17??? 1 1951? 755.9 4833.0? 1207.7
#18??? 1 1952? 891.2 4924.9? 1430.5
#19??? 1 1953 1304.4 6241.7? 1777.3
#20??? 1 1954 1486.7 5593.6? 2226.3
#37??? 2 1951? 588.2 2289.5?? 342.1
#38??? 2 1952? 645.5 2159.4?? 444.2

?with(Grunfeld51_54,tapply(capital,list(firm),FUN=mean))
#??????? 1???????? 2???????? 3???????? 4???????? 5???????? 6???????? 7???????? 8
#1660.4500? 519.9000? 771.6500? 313.7750? 748.1000? 206.8500? 478.7500? 165.8750
?# ????? 9??????? 10 
?#429.0000?? 11.2525 


#or
?ddply(Grunfeld51_54,.(firm),summarise, Mean=mean(capital))
#?? firm????? Mean
#1???? 1 1660.4500
#2???? 2? 519.9000
#3???? 3? 771.6500
#4???? 4? 313.7750
#5???? 5? 748.1000
#6???? 6? 206.8500
#7???? 7? 478.7500
#8???? 8? 165.8750
#9???? 9? 429.0000
#10?? 10?? 11.2525


#or
aggregate(capital~firm,data=Grunfeld51_54,mean)
#?? firm?? capital
#1???? 1 1660.4500
#2???? 2? 519.9000
#3???? 3? 771.6500
#4???? 4? 313.7750
#5???? 5? 748.1000
#6???? 6? 206.8500
#7???? 7? 478.7500
#8???? 8? 165.8750
#9???? 9? 429.0000
#10?? 10?? 11.2525


A.K.




----- Original Message -----
From: jpm miao <miaojpm at gmail.com>
To: r-help <r-help at r-project.org>
Cc: 
Sent: Monday, April 8, 2013 10:29 AM
Subject: [R] How can I extract part of the data in a panel dataset?

Taking the Grunfeld data, which is built-in in R, for example,

(1)How can I construct a dataset (or dataframe) that consists of the data
of all firms in 1951?

(2)How can I calculate the average capital in each form over the period
1951-1954?
What I can imagine is to categorize the data by firm, and then select the
data between 1951 and 1954 for each firm, but how can I do it?

Thanks,

Miao

Grunfeld data:
? invest value capital firm year? 317.6 3078.5 2.8 General Motors 1935
391.8 4661.7 52.6 General Motors 1936? 410.6 5387.1 156.9 General Motors
1937? 257.7 2792.2 209.2 General Motors 1938? 330.8 4313.2 203.4 General
Motors 1939? 461.2 4643.9 207.2 General Motors 1940? 512 4551.2 255.2 General
Motors 1941? 448 3244.1 303.7 General Motors 1942? 499.6 4053.7 264.1 General
Motors 1943? 547.5 4379.3 201.6 General Motors 1944? 561.2 4840.9 265 General
Motors 1945? 688.1 4900.9 402.2 General Motors 1946? 568.9 3526.5 761.5 General
Motors 1947? 529.2 3254.7 922.4 General Motors 1948? 555.1 3700.2
1020.1 General
Motors 1949? 642.9 3755.6 1099 General Motors 1950? 755.9 4833 1207.7 General
Motors 1951? 891.2 4924.9 1430.5 General Motors 1952? 1304.4 6241.7
1777.3 General
Motors 1953? 1486.7 5593.6 2226.3 General Motors 1954? 209.9 1362.4 53.8 US
Steel 1935? 355.3 1807.1 50.5 US Steel 1936? 469.9 2676.3 118.1 US Steel
1937? 262.3 1801.9 260.2 US Steel 1938? 230.4 1957.3 312.7 US Steel 1939
361.6 2202.9 254.2 US Steel 1940
...........

??? [[alternative HTML version deleted]]

______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

arun

2013-Apr-08 18:16 UTC

head link

[R] How can I extract part of the data in a panel dataset?

Hi,
Just to add:
load:

library(plyr) #forgot

if you are using ddply()


#you can directly use subset() in aggregate()
?aggregate(capital~firm,data=subset(Grunfeld,year%in% 1951:1954),mean)
#?? firm?? capital
#1???? 1 1660.4500
#2???? 2? 519.9000
#3???? 3? 771.6500
#4???? 4? 313.7750
#5???? 5? 748.1000
#6???? 6? 206.8500
#7???? 7? 478.7500
#8???? 8? 165.8750
#9???? 9? 429.0000
#10?? 10?? 11.2525
A.K.



----- Original Message -----
From: arun <smartpink111 at yahoo.com>
To: jpm miao <miaojpm at gmail.com>
Cc: R help <r-help at r-project.org>
Sent: Monday, April 8, 2013 2:04 PM
Subject: Re: [R] How can I extract part of the data in a panel dataset?

Hi,
1)
Grunfeld1951<-Grunfeld[Grunfeld$year==1951,]
?Grunfeld1951
#??? firm year??? inv?? value capital
#17???? 1 1951 755.90 4833.00 1207.70
#37???? 2 1951 588.20 2289.50? 342.10
#57???? 3 1951 135.20 1819.40? 671.30
#77???? 4 1951 160.62? 809.00? 203.50
#97???? 5 1951? 80.30? 327.30? 683.90
#117??? 6 1951? 95.30? 676.90? 177.20
#137??? 7 1951? 64.77? 179.00? 450.60
#157??? 8 1951? 54.38? 723.80? 129.70
#177??? 9 1951? 56.49? 339.90? 391.00
#197?? 10 1951?? 4.67?? 83.04??? 9.07

2)


?Grunfeld51_54<-Grunfeld[Grunfeld$year%in% 1951:1954,]
?head(Grunfeld51_54)
?#? firm year??? inv? value capital
#17??? 1 1951? 755.9 4833.0? 1207.7
#18??? 1 1952? 891.2 4924.9? 1430.5
#19??? 1 1953 1304.4 6241.7? 1777.3
#20??? 1 1954 1486.7 5593.6? 2226.3
#37??? 2 1951? 588.2 2289.5?? 342.1
#38??? 2 1952? 645.5 2159.4?? 444.2

?with(Grunfeld51_54,tapply(capital,list(firm),FUN=mean))
#??????? 1???????? 2???????? 3???????? 4???????? 5???????? 6???????? 7???????? 8
#1660.4500? 519.9000? 771.6500? 313.7750? 748.1000? 206.8500? 478.7500? 165.8750
?# ????? 9??????? 10 
?#429.0000?? 11.2525 


#or
?ddply(Grunfeld51_54,.(firm),summarise, Mean=mean(capital))
#?? firm????? Mean
#1???? 1 1660.4500
#2???? 2? 519.9000
#3???? 3? 771.6500
#4???? 4? 313.7750
#5???? 5? 748.1000
#6???? 6? 206.8500
#7???? 7? 478.7500
#8???? 8? 165.8750
#9???? 9? 429.0000
#10?? 10?? 11.2525


#or
aggregate(capital~firm,data=Grunfeld51_54,mean)
#?? firm?? capital
#1???? 1 1660.4500
#2???? 2? 519.9000
#3???? 3? 771.6500
#4???? 4? 313.7750
#5???? 5? 748.1000
#6???? 6? 206.8500
#7???? 7? 478.7500
#8???? 8? 165.8750
#9???? 9? 429.0000
#10?? 10?? 11.2525


A.K.




----- Original Message -----
From: jpm miao <miaojpm at gmail.com>
To: r-help <r-help at r-project.org>
Cc: 
Sent: Monday, April 8, 2013 10:29 AM
Subject: [R] How can I extract part of the data in a panel dataset?

Taking the Grunfeld data, which is built-in in R, for example,

(1)How can I construct a dataset (or dataframe) that consists of the data
of all firms in 1951?

(2)How can I calculate the average capital in each form over the period
1951-1954?
What I can imagine is to categorize the data by firm, and then select the
data between 1951 and 1954 for each firm, but how can I do it?

Thanks,

Miao

Grunfeld data:
? invest value capital firm year? 317.6 3078.5 2.8 General Motors 1935
391.8 4661.7 52.6 General Motors 1936? 410.6 5387.1 156.9 General Motors
1937? 257.7 2792.2 209.2 General Motors 1938? 330.8 4313.2 203.4 General
Motors 1939? 461.2 4643.9 207.2 General Motors 1940? 512 4551.2 255.2 General
Motors 1941? 448 3244.1 303.7 General Motors 1942? 499.6 4053.7 264.1 General
Motors 1943? 547.5 4379.3 201.6 General Motors 1944? 561.2 4840.9 265 General
Motors 1945? 688.1 4900.9 402.2 General Motors 1946? 568.9 3526.5 761.5 General
Motors 1947? 529.2 3254.7 922.4 General Motors 1948? 555.1 3700.2
1020.1 General
Motors 1949? 642.9 3755.6 1099 General Motors 1950? 755.9 4833 1207.7 General
Motors 1951? 891.2 4924.9 1430.5 General Motors 1952? 1304.4 6241.7
1777.3 General
Motors 1953? 1486.7 5593.6 2226.3 General Motors 1954? 209.9 1362.4 53.8 US
Steel 1935? 355.3 1807.1 50.5 US Steel 1936? 469.9 2676.3 118.1 US Steel
1937? 262.3 1801.9 260.2 US Steel 1938? 230.4 1957.3 312.7 US Steel 1939
361.6 2202.9 254.2 US Steel 1940
...........

??? [[alternative HTML version deleted]]

______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reasonably Related Threads

Search for more reasonably related threads

R help - Apr 2013 - How can I extract part of the data in a panel dataset?

[R] How can I extract part of the data in a panel dataset?

[R] How can I extract part of the data in a panel dataset?

[R] How can I extract part of the data in a panel dataset?

Reasonably Related Threads