[yourcast] About preparing data for YourCast

Dear Professors, I am using YourCast in forecasting disease-specific mortalities in China. The data I input is in following format in .csv format. As showed, there are missings in between observed years. YourCast works quite well for LC and OLS model. However, when I try model "map" to smooth over age and cities, I met two problem. year sex2 city2 age r2 time 1987 0 0 10 0.05 0 1988 0 0 10 0.05 1 1989 0 0 10 NA 2 1990 0 0 10 NA 3 1991 0 0 10 0 4 1992 0 0 10 0 5 1993 0 0 10 NA 6 1994 0 0 10 NA 7 1995 0 0 10 NA 8 1996 0 0 10 NA 9 1997 0 0 10 NA 10 1998 0 0 10 NA 11 1999 0 0 10 NA 12 2000 0 0 10 NA 13 2001 0 0 10 NA 14 2002 0 0 10 0 15 2003 0 0 10 0.12 16 2004 0 0 10 0.14 17 2005 0 0 10 0.27 18 2006 0 0 10 0.08 19 2007 0 0 10 0 20 2008 0 0 10 0 21 2009 0 0 10 NA 22 2010 0 0 10 NA 23 2011 0 0 10 0.07 24 2012 0 0 10 NA 25 2013 0 0 10 NA 26 2014 0 0 10 NA 27 2015 0 0 10 NA 28 2016 0 0 10 NA 29 2017 0 0 10 NA 30 2018 0 0 10 NA 31 2019 0 0 10 NA 32 2020 0 0 10 NA 33 1. length of age; the file name is in csidggggaa.csv format and I have age group 0, 1, 5, 10, 15,20,...,85, and file names are db100100.csv db100101.csv db100105.csv db100110.csv...etc. (1001 represents cities) YourCast returns " 'names' attribute [19] must be the same length as the vector [18] " and I guess it maybe related to the name of age group? + G.names="G.name.txt",sample.frame=c(1987,2011,2012,2030)) Loading cross section files and checking for errors... ...Finished Total number of cross sections: 76 model="map", sample.frame=c + (1987,2011,2012,2030)) ... Running MAP model Preparing for smoothing over age groups Preparing for smoothing of time trend over age groups Preparing for smoothing over time Error in names(whoyg) <- age.char : 'names' attribute [19] must be the same length as the vector [18] Data is inconsistent. Error in m$coeff : $ operator is invalid for atomic vectors 2. A second issue happend after I use only age group 10 to 85 to bypass the 1st problem. Here is what I get + G.names="G.name.txt",sample.frame=c(1987,2011,2012,2030)) Loading cross section files and checking for errors... ...Finished Total number of cross sections: 64 model="map", sample.frame=c + (1987,2011,2012,2030)) ... Running MAP model Preparing for smoothing over age groups Preparing for smoothing of time trend over age groups Preparing for smoothing over time Choleski for symmetric, positive definite matrix sample.improper.normal, warning:ranks (chol and code) do not agree [1] 48 48 32 Error in backsolve(CH, diag(dim)) : singular matrix in 'backsolve'. First zero in diagonal [33] In addition: Warning message: In chol.default(wmat, pivot = T) : the matrix is either rank-deficient or indefinite Data is inconsistent. Error in m$coeff : $ operator is invalid for atomic vectors Thank you so much for your comments! Weiming Zhu MD. PhD. China Center for Health Development Studies Peking University