#エアコンのデータ 1961年以降の普及率
#出所)内閣府経済社会総合研究所「消費動向調査」
#http://www.esri.cao.go.jp/jp/stat/menu.html#shohi-z
N<-c(0.4,0.7,1.3,1.7,2.0,2.0,2.8,3.9,4.7,5.9,7.7,9.3,12.9,12.4,17.2,19.5,25.7,29.9,35.5,39.2,41.2,42.2,49.6,49.3,52.3,54.6,57.0,59.3,63.3,63.7,68.1,69.8,72.3,74.2,77.2,77.2,79.3,81.9,84.4,86.2,86.2,87.2,88.8,87.1)
#便宜上、1961年を導入1年目とする。
t<- 1:length(N)
plot(t,N)
#初期値を与えて推定
result<-nls(N~m*(1-exp(-(p+q)*t))/(1+q/p*exp(-(p+q)*t)), start = list(p=0.01,q=0.5,m=80))
summary(result)
#推定値をもちいて 普及率を計算し実測値と比較
Nhat<-predict(result)
lines(t,Nhat)