q<-read.table("shibouritu07M.txt",header=T)
q[,2]
px<-1-q[,2]
lx<-dx<-Lx<-ex<-rep(0,108)
lx[1]<-100000
for (i in 1:107){
lx[i+1]<-px[i]*lx[i]
dx[i]<-lx[i]-lx[i+1]
}
for (i in 1:108){
Lx[108-i]<-lx[109-i]+Lx[109-i]
ex[108-i]<-1/2+Lx[109-i]/lx[109-i]
}
data.frame(x=q[,1],S=q[,2],=px,=round(lx,digits=0),S=round(dx,digits=0),ϗ]=round(ex,digits=2))
Dx<-Cx<-Nx<-Mx<-rep(0,108)
v<-1/1.015
Dx<-v^(0:107)*lx #ׂ(0-107)̃xNg쐬
Cx<-v^(0:107+0.5)*dx
Nx<-rev(cumsum(rev(Dx)))@#ŏINݐς
Mx<-rev(cumsum(rev(Cx)))
KISU<-data.frame(DX=round(Dx,digits=0),CX=round(Cx,digits=0),
NX=round(Nx,digits=0),MX=round(Mx,digits=0))
KISU

TEIKI<-function(x,n){
(Mx[x+1]-Mx[x+n+1])/(Nx[x+1]-Nx[x+n+1])
}
SEIZON<-function(x,n){
Dx[x+n+1]/(Nx[x+1]-Nx[x+n+1])
}

TEIKI(30,30)
SEIZON(30,30)
TEIKI(30,30)+SEIZON(30,30)

V<-function(x,n,t){
(Mx[x+t+1]-Mx[x+n+1]+Dx[x+n+1])/Dx[x+t+1]-(Mx[x+1]-Mx[x+n+1]+Dx[x+n+1])/(Nx[x+1]-Nx[x+n+1])*(Nx[x+t+1]-Nx[x+n+1])/Dx[x+t+1]
}
V(30,10,0:10)
V(40,10,0:10)

cipois<-function(x){
ppois(mu,x)-alpha  # ́AKz2W΍ɑΉ2.27̗p.
}
alpha<-0.0227
d<-c(57,50,56,83,139)@#10-14ΎS
n<-c(443870,522463,597643,670188,730207)@#10-14Όoߌ_
for(i in 1:5){
mu<-d[i]
result<-uniroot(cipois,c(0,200))@#
print(result$root,digits=1)@#S̏l
print((1-exp(-result$root/n[i]))*1000,digits=2)@#S̎S
}