簡體   English   中英

#obs每月的中位數圖。 在r-studio中每月在第二軸上

[英]Median plot per month with #obs. on second axis per month in r-studio

我有一個由2個變量組成的data.frame,每個250萬個obs。

str(values)
data.frame':    2529905 obs. of  2 variables:
 $ Date : Factor w/ 498 levels "1977-11","1978-06",..: 108 60 12 108 58 108 132 188 51 60     ...
$ Value: num  223000 171528 110269 426000 172436 ...
> head(values)
 Date    Value
1 2003-01 223000.0
2 1999-01 171528.0
3 1992-01 110268.6
4 2003-01 426000.0
5 1998-11 172436.5
6 2003-01 334000.0

我想用每個日期的中位數制作一個data.frame:

library(plyr)
medianperdate = ddply(values, .(Date), summarize, median_value = median(Value))

> str(medianperdate)
'data.frame':   498 obs. of  2 variables:
 $ Date        : Factor w/ 498 levels "1977-11","1978-06",..: 1 2 3 4 5 6 7 8 9 10 ...
 $ median_value: num  106638 84948 85084 75725 88487 ...
> head(medianperdate)
     Date median_value
1 1977-11    106638.35
2 1978-06     84947.65
3 1985-07     85083.79
4 1986-05     75724.58
5 1986-11     88487.14
6 1986-12     98697.20

但是我想要的是一個額外的列,它每月統計觀察值(例如2003-01,使用的數據將是對象“值”

在另一列中,我定義了它的類別是:

a : < 200 000 
b : < 300 000 & > 200 000
c : < 300 000 & > 2000000

我將繼續嘗試此操作,但由於我已經停滯了幾個小時,因此非常感謝您的幫助!!

如果不清楚,我能理解什么。 以下testdataframe展示了我希望數據框如何顯示

> testdf
Year MedianValue HouseClass #Observations
1 1999-1      200000          B           501
2 1999-2      150000          A           664
3 1999-3      250000          C           555

就像我對您先前的問題的回答0

library(data.table)
dt <- data.table(df)


dt2 <- dt[,list(
   medianvalue = median(value),
   obs = .N
   ),
   by = "Date"
]

dt2[,HouseClass := "c"]
dt2[obs < 300000,HouseClass := "b"]
dt2[obs < 200000,HouseClass := "a"]

您可以在apply中編寫函數,並應用類似函數(包括plyr函數)。 它看起來像這樣:

ddply(values, .(Date), .fun = function(x) {
  median <- median(x)
  value <- ifelse(median < 200000, 'A', ifelse(median < 300000, 'B', 'C'))
  n <- length(x)
  return(c(median, value, n))
})

暫無
暫無

聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM