作者chien533 (谣指部总司令)
看板Statistics
标题[问题] R的lag function
时间Sun Jun 15 12:07:55 2008
这几天用R的lag function来做某个变数的lag effect
原本以为很简单的一个函式居然产生无解的问题
就是我照着reference manual上的范例去做
但是另存出来却丝毫没有任何效果
举简单的例子来看
假设随便用个 a<- c(1,2,3,4,5)
而在使用 b<-lag(a,k=1) 之後
另存新资料出来的结果应该是 {.,1,2,3,4} 才对
可是居然还是 {1,2,3,4,5}
最後我直接使用手册上的范例:
lag(ldeaths, 12)
如果是列印在萤幕上:
[原始资料]
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1973 3035 2552 2704 2554 2014 1655 1721 1524 1596 2074 2199 2512
1974 2933 2889 2938 2497 1870 1726 1607 1545 1396 1787 2076 2837
1975 2787 3891 3179 2011 1636 1580 1489 1300 1356 1653 2013 2823
1976 3102 2294 2385 2444 1748 1554 1498 1361 1346 1564 1640 2293
1977 2815 3137 2679 1969 1870 1633 1529 1366 1357 1570 1535 2491
1978 3084 2605 2573 2143 1693 1504 1461 1354 1333 1492 1781 1915
[lag後]
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
1972 3035 2552 2704 2554 2014 1655 1721 1524 1596 2074 2199 2512
1973 2933 2889 2938 2497 1870 1726 1607 1545 1396 1787 2076 2837
1974 2787 3891 3179 2011 1636 1580 1489 1300 1356 1653 2013 2823
1975 3102 2294 2385 2444 1748 1554 1498 1361 1346 1564 1640 2293
1976 2815 3137 2679 1969 1870 1633 1529 1366 1357 1570 1535 2491
1977 3084 2605 2573 2143 1693 1504 1461 1354 1333 1492 1781 1915
这是没问题的。
但随便用个 View() 函示来看,使用前和使用後的资料是完全一样的。
请问该怎样让使用 lag function 过後所另存出来的资料真的有 lag 後的效果?
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 152.2.149.157
1F:推 lin15:好像物件要是time series 物件才可以使用lag 06/15 12:26
2F:→ bcs:手动…> n=length(xvar), xvar[-n] or x=ts(xvar),lag(x,1) 06/15 12:41
3F:推 clickhere:a<-ts(1:5);a;b<-lag(a,k=1);b;is.ts(a);is.ts(b) 06/15 12:46
4F:→ clickhere:str(a);str(b); 06/15 12:50
5F:→ clickhere:Use object, save life. View() is only worked on MS. 06/15 12:54
6F:→ chien533:搞定了,我用的是bcs的方法,谢谢! 06/15 13:37