# vout1 = VPA�̌v�Z����
# yprs.year.shortr = RPS�̎Q�ƔN

# Frec�����ōĕ]������ꍇ
# 
#  ������C�ʏ��ABC�v�Z�iABC�Z��N��2018�N�C�ڕW�N��2028�N�j  
# ���̏ꍇ�C2017�N��F��current F, 2018�N�ȍ~��F��Frec�ƂȂ� 

fout1<-future.vpa(vout1,multi=0.5,waa.year=2016,Frec=list(stochastic=FALSE,future.year=2028,Blimit=64000),ABC.year=2018,nyear=20,rec.arg=list(rps.year=rps.year.shortr))

# �����čĕ]���iABC�Z��N��2018�N�C�ڕW�N��2028�N�j  
# �ĕ]���ł́C2017�N�ȍ~��F��Frec�Ƃ���
# RVPA�ł� ABC.year�ȍ~�N��F��Frec  �ɒu�������̂�ABC.year��1�N���点�΂悢

fout1r<-future.vpa(vout1,multi=0.5,waa.year=2016,Frec=list(stochastic=FALSE,future.year=2027,Blimit=64000),ABC.year=2017,nyear=20,rec.arg=list(rps.year=rps.year.shortr))

# �����čāX�]���́H=> 2016�N��F��Frec�ɒu�������D
# �����C2016�N��F�͂��ł�VPA�̌��ʂƂ��ďo�͂���Ă���̂ŁC���ʂɂ���2016�N��F��VPA���ʂ��g���Ă��܂�
# ���̂��߁CVPA���ʂ���ŏI�N��F�̒l���O�ɒu��������K�v������

vout1r <- vout1 # VPA���ʂ��āX�]���p��vpa1r�ɃR�s�[
vout1r$faa["2016"] <- 0 # 2016�N��F��0�ɒu��������D0�����Ă����ƁC�����\����F���g����d�l�ɂȂ��Ă���
fout1rr <- future.vpa(vout1r,multi=0.5,waa.year=2016,Frec=list(stochastic=FALSE,future.year=2026,Blimit=64000),ABC.year=2016,nyear=20,rec.arg=list(rps.year=rps.year.shortr),start.year=2016) # �����\���̊J�n�N��2016�N�Ƃ���i���������2015�N��numbers at age��F at age���g���đO�i�v�Z�����C2016�N��F��Frec���g�����ƂɂȂ�j�D�܂��CABC.year=2016, future.year=2026�ɕύX�D