使用逻辑增长曲线拟合模型

Question

我正在尝试在https://www.icao.int/MID/Documents/2014/Aviation%20Data%20Analyses%20Seminar/8991_Forecasting_en.pdf中复制模型（第40页）使用此语法

air<-data.frame (Año=seq (1960,1996,1),
     Y=C (2.535,2.629,2.778,2.988,3.276,3.593,3.941,4.535,4.874,5.247,5.52,5.639,6.126,6.451,6.691,6.953,7.293,7.627,9.434,9.295,9.290,9.331,9.538,9.769,10.35,10.936,11.426,12.305,12.907,13.054,13.571,13.273,14.415,14.683,15.648,16.248,16.941),
     Tiempo=seq (1,37,1)

对于第一个值，我使用logit

a<-coef(lm(logit(Y)~Año,data=air));a

并且对于模型，使用

fit<-nls(Y~phi1/(1+phi2*exp(-phi3*Año)),
            start=list(phi1=25,phi2=a[1],phi3=a[2]),data=air)

我正在尝试到达HTML中的同一张表，但我有

“ nlsModel（formula，mf，start，wts）中的错误：初始参数估计时的奇异梯度矩阵“

Answer 1

似乎您的问题是如何在R中实现逻辑回归。您可能要考虑使用另一个函数。通常，逻辑回归是通过glm函数实现的。如果将“家庭”参数设置为“二项式”，它将实现逻辑回归。GLM documentation in R