求助：累加如何进行编程

2014-3-27 22:08:04

复制代码

就行了，cumsum是前面所有元素累加，不是\[z_{i}=y_{i}+y_{i-1}\]

扫码加我拉你入群

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

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

2014-3-27 22:16:13

danica33 发表于 2014-3-27 22:08
就行了，cumsum是前面所有元素累加，不是

非常感谢详细的回答，这个得出来是一个数呢，还是一列数呢。比如在STATA中，我可以令变量b等于ai+bi-1
这样的出来的b，就是0，a2+a1, a3+a2+a1, a4+a3+a2+a1,等等

R中可以通过cumsum来得到这样的一列数据吗？

扫码加我拉你入群

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

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

2014-3-27 22:19:55

winniewang2222 发表于 2014-3-27 22:16
非常感谢详细的回答，这个得出来是一个数呢，还是一列数呢。比如在STATA中，我可以令变量b等于ai+bi-1
这 ...

得到的是一列数，cumsum就是向量a=(a1,...,an),,cumsum(a)得到a1,a1+a2,a1+a2+a3,....a1+...+an

扫码加我拉你入群

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

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

2014-3-27 22:26:33

danica33 发表于 2014-3-27 22:19
得到的是一列数，cumsum就是向量a=(a1,...,an),,cumsum(a)得到a1,a1+a2,a1+a2+a3,....a1+...+an

明白了~！非常感谢！！

再请教一个问题，我将dta的数据导入到了R, 用
>data<-read.dta("D:/test.dta")
然后输入data后，显示的是
a
1 0
2 3
3 4
4 6
...
这样的格式
我尝试对a进行加减乘除，比如a*d
我尝试
>c<-a*d
结果说找不到a
请问我是不是应该吧a转换一下格式还是怎样呢

扫码加我拉你入群

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

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

点击查看更多内容…

2014-3-27 22:33:02

winniewang2222 发表于 2014-3-27 22:26
明白了~！非常感谢！！

再请教一个问题，我将dta的数据导入到了R, 用

你attach(data)看行不行，还有d是什么

扫码加我拉你入群

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

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

2014-3-27 22:38:46

danica33 发表于 2014-3-27 22:33
你attach(data)看行不行，还有d是什么

谢谢建议，我试试
>d
[1] 1 4 5 7 8 ...
这个样子的，和a那样一列下来的格式不同

扫码加我拉你入群

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

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

2014-3-27 22:41:37

winniewang2222 发表于 2014-3-27 22:38
谢谢建议，我试试
>d
[1] 1 4 5 7 8 ...

现在可以加减吗

扫码加我拉你入群

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

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

2014-3-27 22:48:22

danica33 发表于 2014-3-27 22:41
现在可以加减吗

可以了！太好了！这似乎是把a给提出来了？

还有好几个问题，如果你要去休息，就再说啊，呵呵

之前有发帖问道排序的问题，回复建议order这个命令，因为我今天刚接触R，一下子很难很快弄明白怎么写code，不知道你熟不熟悉呢？

就是比如我将刚刚提取出来的a由大到小进行排序，然后产生一个新的变量x，x等于排序后a的序号。我尝试过这个命令
>cbind(year=as.numeric(rownames(data)),data)
可以得到year这个变量等于a的序号，但是问题在于，1，如何先对a排序，2,也是用attach(year)将year提取出来吗？

扫码加我拉你入群

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

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

2014-3-27 22:50:06

danica33 发表于 2014-3-27 22:41
现在可以加减吗

我试了试，不能将year提出来

> cbind(year=as.numeric(rownames(data)),data)
year  a
1    1  1
2    2  2
3    3  3
4    4  3
5    5  4
6    6  5
7    7  6
8    8 66
9    9  6
10 10  7
11 11 88
12 12 88
13 13 99
14 14 88
15 15  4
16 16  4
> attach(data)
下列对象被屏蔽了from data (position 3):

a
> a
[1]  1  2  3  3  4  5  6 66  6  7 88 88 99 88  4  4
> year
错误: 找不到对象'year'
>

扫码加我拉你入群

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

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

2014-3-27 22:55:48

winniewang2222 发表于 2014-3-27 22:48
可以了！太好了！这似乎是把a给提出来了？

还有好几个问题，如果你要去休息，就再说啊，呵呵

order(a)就可以得到a从小到大排列后对应序号，你rownames(data)得到什么东西？data是什么样子

扫码加我拉你入群

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

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

2014-3-27 22:56:09

cbind(year<-as.numeric(rownames(data)),data)

扫码加我拉你入群

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

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

2014-3-27 23:01:04

danica33 发表于 2014-3-27 22:55
order(a)就可以得到a从小到大排列后对应序号，你rownames(data)得到什么东西？data是什么样子

> order(a)
[1] 1 2 3 4 5 15 16 6 7 9 10 8 11 12 14 13
>
这是order后得出来的结果，并不是12345.。。16这样的序列

扫码加我拉你入群

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

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

2014-3-27 23:02:31

jmpamao 发表于 2014-3-27 22:56
cbind(year

谢谢！~ 我之前尝试的是这个，但是无法将year分离出来为一个独立的变量，与a进行运算

扫码加我拉你入群

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

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

2014-3-27 23:03:46

winniewang2222 发表于 2014-3-27 23:01
> order(a)
[1] 1 2 3 4 5 15 16 6 7 9 10 8 11 12 14 13
>

对啊，a从小到大排列后对应原来的序号，你的nrownames(data)得到的是行号

扫码加我拉你入群

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

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

2014-3-27 23:07:08

看好了是 year <- 不是=
<-是创立一个对象year，如果是=号，只是cbind 函数中数据的列名参数为year

扫码加我拉你入群

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

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

2014-3-27 23:07:50

danica33 发表于 2014-3-27 23:03
对啊，a从小到大排列后对应原来的序号，你的nrownames(data)得到的是行号

这是原来的序号？原来比如最后两个值是4，那序号也不应该是14,13吧...

> a
[1] 1 2 3 3 4 5 6 66 6 7 88 88 99 88 4 4
> order(a)
[1] 1 2 3 4 5 15 16 6 7 9 10 8 11 12 14 13
>

扫码加我拉你入群

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

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

2014-3-27 23:11:44

jmpamao 发表于 2014-3-27 22:56
cbind(year

太棒了！原来换成<-就可以了~~！
再用attach(data）是分离出了year和a

这之前有个问题就是，我怎样可以先将a从小到大进行排序呢，然后再产生排序后的序号呢

扫码加我拉你入群

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

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

2014-3-27 23:14:50

winniewang2222 发表于 2014-3-27 23:07
这是原来的序号？原来比如最后两个值是4，那序号也不应该是14,13吧...

> a

order(a)是a从小到大排列后第i大的元素在原来的向量a中的位置，比如a排列后
sort(a)##a从小到大排列
[1] 1 2 3 3 4 4 4 5 6 6 7 66 88 88 88 99
>
第5大元素为4，在原来a中位置5

扫码加我拉你入群

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

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

2014-3-27 23:18:55

danica33 发表于 2014-3-27 23:14
order(a)是a从小到大排列后第i大的元素在原来的向量a中的位置，比如a排列后
sort(a)##a从小到大排列
[ ...

找到方法了，用rank
> a
[1]  1  2  3  3  4  5  6 66  6  7 88 88 99 88  4  4
> rank(a)
[1]  1.0  2.0  3.5  3.5  6.0  8.0  9.5 12.0  9.5 11.0 14.0 14.0 16.0 14.0  6.0  6.0
> d<-rank(a)
> d
[1]  1.0  2.0  3.5  3.5  6.0  8.0  9.5 12.0  9.5 11.0 14.0 14.0 16.0 14.0  6.0  6.0
>

不知道可否对3.5 3.5, 9.5 9.5这样的数字取整呢，因为序号代表一个人的话，不能说第三个半人吧

扫码加我拉你入群

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

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

2014-3-27 23:20:47

winniewang2222 发表于 2014-3-27 23:07
这是原来的序号？原来比如最后两个值是4，那序号也不应该是14,13吧...

> a

最后面的两个4， 4 在a中的序号为 15,16 最后两个，对吧

order 后，看看 15,16 在什么位置，就是把4,4,按照大小放前面了。
运行下面看看
a
order(a)
a[order(a)]

扫码加我拉你入群

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

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

2014-3-27 23:25:56

winniewang2222 发表于 2014-3-27 23:18
找到方法了，用rank
> a
[1] 1 2 3 3 4 5 6 66 6 7 88 88 99 88 4 4

你要的year是秩啊，看你要怎么处理相等的数的秩，可以看一下rank的帮助，?rank，比如
> rank(a,ties.method="fir")
[1] 1 2 3 4 5 8 9 12 10 11 13 14 16 15 6 7

扫码加我拉你入群

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

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

2014-3-27 23:28:33

jmpamao 发表于 2014-3-27 23:20
最后面的两个4， 4 在a中的序号为 15,16 最后两个，对吧

order 后，看看 15,16 在什么位置，就 ...

这个比rank更加好，我尝试了order(a[order(a)]),总算得到想要的了！谢谢大家帮忙！

> a
[1]  1  2  3  3  4  5  6 66  6  7 88 88 99 88  4  4
> order(a)
[1]  1  2  3  4  5 15 16  6  7  9 10  8 11 12 14 13
> a[order(a)]
[1]  1  2  3  3  4  4  4  5  6  6  7 66 88 88 88 99
> order(a[order(a)])
[1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16
>

还有个问题想请教就是，不知道各位是否熟悉核密度估计呢，我查到命令是density(), 我尝试density(a),但是也有问题，我想得到的是新生成一个变量，比如y, y是等于x对应的密度函数值

扫码加我拉你入群

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

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

2014-3-27 23:30:59

danica33 发表于 2014-3-27 23:25
你要的year是秩啊，看你要怎么处理相等的数的秩，可以看一下rank的帮助，?rank，比如
> rank(a,ties.met ...

原来这个方法也可以，谢谢~！记下了^^

还有个问题是关于核密度估计，之前也回复jmpa了，不知道你熟悉吗？

扫码加我拉你入群

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

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

2014-3-27 23:39:33

density 会重新模拟数据的。如果你的数据长度小于 512，模拟到512个。kernel = c("gaussian",
"epanechnikov", "rectangular", "triangular", "biweight",
"cosine", "optcosine")

你要查看的话

复制代码

扫码加我拉你入群

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

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