APC模型本质上是一个广义线性模型,只不过外加了年龄时期队列的影响,把响应变量发病率和年龄时期等相关变量合成列表结构,指定响应变量分布类型和连接函数,调用apc包里的apc.fit.mode函数就可以了。例如:
>data.list<- data.Italian.bladder.cancer()
>fit.apc<- apc.fit.model(data.list,"poisson.dose.response","APC")
>fit.ac<- apc.fit.model(data.list,"poisson.dose.response","AC")
另外Epi包是一个传染病模型的统计分析包,主要还是APC模型,当然外加了许多数据集和许多对象操作函数,比如画图等,而apc包则是专门做APC模型分析的包,也就是说不仅可以应用与传染病,也可以用于别的领域。可以输入以下命令查看相应的介绍。
> help(package="Epi")
> help(package="apc")