房地产抑制政策何时会出？

发表于 | 分类于

房价快速上涨的时候，上车的欢喜，没上车的愁，ZF为了照顾劳苦大众的情绪，总会象征性的出点抑制政策。以下是2007年以来的主要的抑制政策：

  • 2007年9月，二套房首付提升至50%
  • 2009年12月，个人住房转让营业税免征时限由2年恢复到5年
  • 2010年4月，“新国十条”，差别化信贷政策；“京十二条”，北京家庭只能新购一套商品房
  • 2013年3月，“京十九条”，二手房20%个税
  • 2013年4月，二套首付7成
  • 2016年9月，提高首付比例，重启“70/90”

下面用不同的方法来拟合 是否出抑制政策，1代表是，0代表否

logistic回归

logistic回归学习，请参见R语言与二分类logistic回归

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
library(readxl)
library(lubridate)
library(ggplot2)
setwd("/Users/lixinyao/Desktop")
mydata = read_excel("月度全市场量价07年到16年.xlsx",sheet = 1,
                    col_types = rep("text",3))
head(mydata)
mydata$月 = ymd(20061201) + months(1:length(mydata$月))
mydata$网签量 = as.numeric(mydata$网签量)
mydata$成交均价 = round(as.numeric(mydata$成交均价),0)
mydata$月 = str_trim(mydata$月)
# 是否出政策
mydata$政策 = NULL
month = c("2007-09-01","2009-12-01",
          "2010-04-01","2013-03-01",
          "2013-04-01","2016-09-01")
for(i in 1:length(mydata$月)){
  if(mydata$月[i] %in% month) mydata$政策[i] = 1
  else mydata$政策[i] = 0
}
# logistic
logistic.fit = glm(政策~网签量+成交均价,
              binomial(link = "logit"),
              data = mydata)
mydata$logitpre = predict(logistic.fit)
mydata$logitpre = ifelse(mydata$logitpre>-1,1,0)

查看拟合结果，分类准确率为 $\frac{109+4}{109+2+2+4} = 96.58\%$ ，实际无政策，预测有政策的时点为 2009年11月、2016年3月；实际有政策，预测无政策的时点为 2007年9月、2013年4月

1
2
3
4
> table(mydata$政策,mydata$logitpre)
      0   1
  0 109   2
  1   2   4

svm

1
2
3
4
5
library(e1071)
mydata = mydata %>% select(-logitpre)
svm.fit = svm(政策~.,data = mydata)
mydata$svmpre = predict(svm.fit)
mydata$svmpre = ifelse(mydata$svmpre>0.025,1,0)

分类准确率为100%

1
2
3
4
> table(mydata$政策,mydata$svmpre)
      0   1
  0 111   0
  1   0   6

当svmpre值大于0.025时，政府就会出抑制政策

喂他一颗糖