公积金房贷优惠比较（R）

假设首套房总价464975，贷款期限20年月付，首付和贷款比例2:8，且贷款金额向下取整万数，可选优惠985折。根据央行2022年10月1日发布的五年期以上首套个人住房公积金贷款利率3.1%计算是否选择优惠如下：

i=0.031/12
v=1/(1+i)
a=(1-v**(20*12))/i
b=464975
x=b*0.985-floor(b*0.985*0.8/10000)*10000
y=b-floor(b*0.8/10000)*10000
z=x-y
k1=(floor(b*0.985*0.8/10000)*10000)/a
k2=(floor(b*0.8/10000)*10000)/a
m=data.frame(b,x,y,z,k1,k2)
names(m)=c("总价","优惠首付","一般首付","优惠后-前","优惠月供","一般月供")

month=1:(12*20)
m1=data.frame(month)
m1["优惠月供"]=k1
w=floor(b*0.985*0.8/10000)*10000
m1["优惠月供本金"]=0
m1["优惠月供利息"]=0
m1["优惠剩余本金"]=0
m1[1,4]=w*i
m1[1,3]=k1-m1[1,4]
m1[1,5]=w-m1[1,3]
for (k in 2:(12*20)){m1[k,4]=m1[k-1,5]*im1[k,3]=k1-m1[k,4]m1[k,5]=m1[k-1,5]-m1[k,3]
}
m1["月供"]=k2
m1["月供本金"]=0
m1["月供利息"]=0
m1["剩余本金"]=0
w1=floor(b*0.8/10000)*10000
m1[1,8]=w1*i
m1[1,7]=k2-m1[1,8]
m1[1,9]=w1-m1[1,7]
for (k in 2:(12*20)){m1[k,8]=m1[k-1,9]*im1[k,7]=k2-m1[k,8]m1[k,9]=m1[k-1,9]-m1[k,7]
}