R语言分组计算 - 经管之家

12235

收藏 2015-12-17

请问用R语言如何给某一变量从小到大等分成十组，然后再对每一组分别进行其他计算？
对R语言某一变量从小到大等分成十组后，可否分别保存这十组数据？
例如以下数据，可否按V1从小到大分成十组，然后在每一组中，对V2和V3进行计算？还有对V1分成十组后，能不能保存下来每一组数据。
谢谢，拜托有好心人可以帮忙。

V1	V2	V3
16.50694	-0.13895	14710415
18.6736	-0.11111	13075925
18.4233	0.075893	14068294
15.50773	0.128631	15877909
17.2779	0.170343	18582602
17.4706	0.324607	24614650
9.254718	0.143874	28156047
14.8777	0.322046	37223578
14.7202	-0.00418	37067912
7.830189	-0.00892	36737122
11.6707	0.37447	50494085
12.2747	0.105588	55825637
5.46873	0.055417	57427590
9.3111	0.316497	75603255
9.6813	0.048855	79296817
4.081657	0.052105	83428599
7.1357	0.201851	1E+08
7.5524	-0.24912	75290241
3.982347	0.069845	88525516
5.688	-0.13731	76370458
6.0896	-0.00451	76026447
6.331345	-0.14932	64674081
5.5228	0.050355	67930720
5.9079	-0.14821	57862662
9.568371	-0.23385	46175411
6.8326	0.076048	49686940
6.9579	-0.02885	48253663

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

jiangbeilu

2015-12-17 16:46:47

1.新建一个v4变量

2.利用order排序v1,然后分成10组，用cut函数就可以，赋值给v4

3.根据v4,用tapply，分组计算v2,v3

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-17 21:50:49

jiangbeilu 发表于 2015-12-17 16:46
1.新建一个v4变量

2.利用order排序v1,然后分成10组，用cut函数就可以，赋值给v4

请问cut的代码该怎么打啊？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-17 22:05:36

jiangbeilu 发表于 2015-12-17 16:46
1.新建一个v4变量

2.利用order排序v1,然后分成10组，用cut函数就可以，赋值给v4

请问cut和tapply具体的代码能写一下吗？万分感谢。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jiangbeilu

2015-12-17 22:59:32

复制代码

这个就是把1:10这10个数分成了两组。
?tapply
可以查看具体用法的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-17 23:14:10

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

我隐身啦

2015-12-18 00:30:42

suimong 发表于 2015-12-17 23:14

谢谢你的回复。我试了你的代码，确实分成了10组，但是每一组的数量并不一样啊。这是按照什么分类的？我已经排好序了，我想等分成10组，每组数量一样，请问怎么分？就是像如果用stata的话类似于group（10）这样的。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-18 00:41:12

我隐身啦发表于 2015-12-18 00:30
谢谢你的回复。我试了你的代码，确实分成了10组，但是每一组的数量并不一样啊。这是按照什么分类的？我已 ...

你是要按数量等分啊，我以为你是要按区间大小等分。那你就连cut函数都不需要了，直接生成一个和你的data frame行数相等的vector，用等量的1:10填充就好了。
比如你的data frame有200行，那么你可以写V4=sort(rep(1:10,20))

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-18 02:08:51

suimong 发表于 2015-12-18 00:41
你是要按数量等分啊，我以为你是要按区间大小等分。那你就连cut函数都不需要了，直接生成一个和你的data ...

我试了确实很好用。可是如果行数不能整除10怎么办？比如有199行或者201行就会报错误
Error: wrong result size (200), expected 199 or 1，
怎么办？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-18 08:43:49

那要看你准备把多出来的一行或几行放到第几组去了？这个是你的分析逻辑决定的

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-18 11:18:07

suimong 发表于 2015-12-18 08:43
那要看你准备把多出来的一行或几行放到第几组去了？这个是你的分析逻辑决定的

可是要怎么决定啊？好像没有选择的地方啊。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-18 12:30:22

我隐身啦发表于 2015-12-18 11:18
可是要怎么决定啊？好像没有选择的地方啊。

这个是看你做数据分析的目的是什么，和R没有关系了。只有你知道为什么要对数据进行分组计算，所以你也自然要考虑不能等分的情况下怎么处理最合适。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-18 12:36:21

suimong 发表于 2015-12-18 12:30
这个是看你做数据分析的目的是什么，和R没有关系了。只有你知道为什么要对数据进行分组计算，所以你也自然 ...

可是用stata中的group（10），就可以分10组，如果多一行好像就随机有一组多一行。但是用R这样做就报错进行不下去了啊。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-18 12:48:04

随机分的话对于你数据分析来说不见得是好事情吧？你如果多了9行呢？这9行是自成一组呢，还是全扔进最后一组呢，还是平均分给之前的组呢？平均分是怎么分？是开头的九组还是最后九组每组多一行呢？不同的分析目的会导致不同的分法，最后产生的结果也是不尽相同的。R和Stata并不能帮你决定怎么分析数据。
你可以写清楚你的数据有几行，多出来的行怎么分，然后我可以把代码怎么写告诉你

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-18 22:35:53

suimong 发表于 2015-12-18 12:48
随机分的话对于你数据分析来说不见得是好事情吧？你如果多了9行呢？这9行是自成一组呢，还是全扔进最后一组 ...

非常感谢你。我之前以为随机分的，以为这不会有太大影响，没想到每个分法还有不同的影响，还有不同的代码。那我再看看，如果有什么问题希望能再向你请教，非常感谢你。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-21 20:55:07

suimong 发表于 2015-12-18 00:41
你是要按数量等分啊，我以为你是要按区间大小等分。那你就连cut函数都不需要了，直接生成一个和你的data ...

你好，最近经常遇到分组的问题，比如我的data frame有21649行，我想分成3组，多出来的一行放到最后一组，代码该怎么写？还是用sort吗？
如果我还想把这21649行分成5组，多出来的四行平均分到后面四组，代码又该怎么写？
希望能得到你的解答，谢谢。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-22 01:15:48

我隐身啦发表于 2015-12-21 20:55
你好，最近经常遇到分组的问题，比如我的data frame有21649行，我想分成3组，多出来的一行放到最后一组， ...

复制代码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-22 14:12:11

suimong 发表于 2015-12-22 01:15

太谢谢了

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-23 11:27:31

suimong 发表于 2015-12-22 01:15

你好，不好意思打扰了，还有个问题向你请教。
我要求实际资本结构和目标资本结构的偏离程度。现在我有了实际的资本结构，我想知道目标的资本结构，影响目标资本结构的变量我都有了。要求出目标资本结构的拟合值，我才能知道实际资本结构与目标资本结构偏离的程度。怎么求目标资本结构的拟合值？
我看有的用广义线性模型( GLM) 的ogit 最大似然估计( logit QMLE) 方法进行参数估计，有的用先求初始值再用非线性模型进行估计，有的用广义矩估计（GMM）。
这些是都可以吗？该用哪种？R语言的代码方便告知一下吗？非常感谢。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-23 11:27:38

suimong 发表于 2015-12-22 01:15

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-23 15:21:16

我隐身啦发表于 2015-12-23 11:27
你好，不好意思打扰了，还有个问题向你请教。
我要求实际资本结构和目标资本结构的偏离程度。现在我有了 ...

你如果有具体的代码问题我可以帮你，研究问题帮不了你了。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-23 22:21:28

suimong 发表于 2015-12-23 15:21
你如果有具体的代码问题我可以帮你，研究问题帮不了你了。

请问GMM估计的方程有什么要求吗？
gmm(g,x,t0=NULL,gradv=NULL, type=c("twoStep","cue","iterative"),
wmatrix = c("optimal","ident"), vcov=c("HAC","iid","TrueFixed"),
kernel=c("Quadratic Spectral","Truncated", "Bartlett", "Parzen", "Tukey-Hanning"),
crit=10e-7,bw = bwAndrews, prewhite = FALSE, ar.method = "ols", approx="AR(1)",
tol = 1e-7, itermax=100,optfct=c("optim","optimize","nlminb", "constrOptim"),
model=TRUE, X=FALSE, Y=FALSE, TypeGmm = "baseGmm", centeredVcov = TRUE,
weightsMatrix = NULL, traceIter = FALSE, data, eqConst = NULL,
eqConstFullVcov = FALSE, ...)
help里面有这么多代码？该怎么用啊？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

suimong

2015-12-23 23:22:40

我隐身啦发表于 2015-12-23 22:21
请问GMM估计的方程有什么要求吗？
gmm(g,x,t0=NULL,gradv=NULL, type=c("twoStep","cue","iterative"), ...

这么多都只是参数，而且绝大部分都是带默认值的，不用怕。
每个函数的帮助文档都分为几个section，例如Usage, Arguments, Details, Value, Example这几个基本上所有包里的主要函数都会有。想了解原理的话从Details开始看；如果觉得看着头大的话那么就从Examples开始看，自己一行一行对着Examples敲一遍，边敲边运行，敲完一遍基本就能明白个大概了。一般包的作者会把函数的基本用法放在Examples的最前面，所以如果你秉承实用主义原则的话甚至后面一些例子不看都行。当你初步明白函数的用法之后再回头去看Arguments里对每个参数的解释和Details部分就会更清楚一些。
gmm包的作者把帮助文档写的非常规范，详尽和清晰，从Examples部分开始看，不会看不懂的。
另外，gmm的帮助文档里，尤其是Details部分有不少的数学公式，因此建议看pdf版。百度："gmm cran"，在r-project.org的页面上找reference manual

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

我隐身啦

2015-12-24 00:04:41

suimong 发表于 2015-12-23 23:22
这么多都只是参数，而且绝大部分都是带默认值的，不用怕。
每个函数的帮助文档都分为几个section，例如U ...

恩恩，好的。我仔细看看。谢谢了啊。太感谢了。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群