[英]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.