Lee Carter模型代码错误

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

data=read.csv(file="D://毕业论文相关//male.csv",header=TRUE)
newdata=read.csv(file="D://毕业论文相关//male-newyear.csv",header=TRUE)
break_year<-2018
T<-break_year-1998+1
n<-2035-break_year
D_all<-matrix(rep(0,19*38),19,38)
E_all<-matrix(rep(0,19*38),19,38)
m_all<-matrix(rep(0,19*38),19,38)
lnm_all<-matrix(rep(0,19*38),19,38)
for (t in 1:38) {
for (x in 1:19) {
   D_all[x,t]<-data[t+19*(x-1),3]
   E_all[x,t]<-data[t+19*(x-1),4]
   m_all[x,t]<-D_all[x,t]/E_all[x,t]
   lnm_all[x,t]<-log(m_all[x,t])  
}
}
D<-D_all[,1:T]
E<-E_all[,1:T]
m<-m_all[,1:T]
lnm_all<-lnm_all[,1:T]
m_original<-m_all[,(T+1):38]
age<-data[1:19,2]
year<-as.numeric(levels(factor(data[,1])))
oldyear<-year[1:T];newyear<-year[(T+1):(n+T)]
install.packages("demography")
install.packages("foreast")
library(forecast)
library(demography)
age_svd<-c(1:19)
std.m<-demogdata(m,E,ages=age_svd,years=oldyear,type=c("morality"),label="male",name="china",lambda=0)

std.m<-demogdata(m,E,ages=age_svd,years=oldyear,type=c("morality"),label="male",name="china",lambda=0)
Warning message:
In demogdata(m, E, ages = age_svd, years = oldyear, type = c("morality"),  :
  Unknown type


求问：这是怎么回事，怎么解决

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Carter CART CAR ART Lee