作者developme223 (allen)
看板R_Language
標題[問題]時間序列資料 求中位數與標準差
時間Tue Nov 17 15:36:18 2015
文章分類提示:
- 問題: 當你想要問問題時,請使用這個類別
[問題類型]:
程式諮詢(我想用R 做某件事情,但是我不知道要怎麼用R 寫出來)
[軟體熟悉度]:
請把以下不需要的部份刪除
入門(寫過其他程式,只是對語法不熟悉)
[問題敘述]:
請簡略描述你所要做的事情,或是這個程式的目的
Code date close_price
3052 2012/5/21 50
3052 2012/5/22 43
3052 2012/5/23 44
3052 2012/5/24 51
3041 2012/8/21 17
3041 2012/8/22 16
3041 2012/8/23 21
.
.
.
5521 2014/5/22 22
5521 2014/5/23 21
5521 2014/5/24 20
8011 2015/5/01 102
8011 2015/5/02 103
我的資料大概長成上面的型態,Code為公司代碼(每間公司資料長短皆不同)
想請問各位如果要算每個Code的Close_price之中位數及標準差,要怎麼做呢?
[程式範例]:
[環境敘述]:
版本:R i386 3.2.2
[關鍵字]:
median、tapply
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 120.126.58.61
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1447745780.A.67E.html
1F:→ celestialgod: 關鍵字tapply, median 11/17 15:46
2F:推 Edster: aggregate(as.matrix(X[,3]),by=list(X[,1]), quantile, 11/17 15:46
3F:→ Edster: probs = 0.5) 11/17 15:47
4F:→ Edster: 標準差的結構同上, quantile 換成 sd) 11/17 15:48
※ 編輯: developme223 (120.126.58.61), 11/17/2015 15:52:59
5F:→ celestialgod: dat %>% group_by(me=median(close_price), 11/17 15:53
6F:→ celestialgod: sd = sd(close_price)) 11/17 15:53
7F:→ celestialgod: 先library(dplyr),以及dat是你的資要 11/17 15:53
8F:→ celestialgod: 我推關鍵字是指你可以透過那兩個函數處理XD 11/17 15:54
9F:→ developme223: 謝謝兩位高手解答。 11/17 16:44