Re: [yourcast] Some questions related to YourCast - Yourcast

9 May 2017

Hi Antonio, I'm not sure of the issues with your analysis, but we're
working on a revision of yourcast that should make these easier.

Gary
--
*Gary King* - Albert J. Weatherhead III University Professor - Director,
IQSS <http://iq.harvard.edu/> - Harvard University
GaryKing.org - King(a)Harvard.edu - @KingGary <https://twitter.com/kinggary> -
617-500-7570 - Assistant &lt;king-assist(a)iq.harvard.edu&gt;du>: 617-495-9271

On Fri, Apr 28, 2017 at 9:41 PM, 林晓 &lt;linxiao6(a)mail2.sysu.edu.cn&gt; wrote:

...
  Dear Professor King,
 I was doing some literature reading and I came across your software. I try
 the package with my own data but I got some error.
 I modified my data and make it look like the example data in the file path
 of the package. yourprep runs smoothly with some data like below:

 dta$data$`204500`
            rspi2   popu2 time
 2004          NA 4858042 1994
 2005          NA 4785805 1995
 2006          NA 4704201 1996
 2007 0.006431800 4592693 1997
 2008 0.005454000 4522168 1998
 2009 0.005908178 4432729 1999
 2010 0.005238223 4317368 2000
 2011 0.005180368 4237330 2001
 2012 0.004544100 4176305 2002
 2013 0.004363600 4098927 2003
 2014 0.003821751 4056178 2004
 2015 0.003390900 3958900 2005
 2016          NA      NA 2006
 2017          NA      NA 2007
 2018          NA      NA 2008
 2019          NA      NA 2009
 2020          NA      NA 2010
 2021          NA      NA 2011
 2022          NA      NA 2012
 2023          NA      NA 2013
 2024          NA      NA 2014
 2025          NA      NA 2015

 But then I got error with yourcast() command.

 The scripts run below:
  dta <-
yourprep(dpath='./GARY_KING_CHINA_mortality',index.code="ggggaa",
 tag="csid",
 +                 year.var=TRUE, sample.frame=c(2004,2015,2016,2025),
 +                 G.names="cntry.codes.txt",
adjacency="adjacency.txt",
 +                 verbose=TRUE,summary=FALSE)
 Loading cross section files and checking for errors...
 csid204500.txt
 csid204501.txt
 csid204505.txt
 csid204510.txt
 csid204515.txt
 csid204520.txt
 csid204525.txt
 csid204530.txt
 csid204535.txt
 csid204540.txt
 csid204545.txt
 csid204550.txt
 csid204555.txt
 csid204560.txt
 csid204565.txt
 csid204570.txt
 csid204575.txt
 csid204580.txt
 csid204585.txt
 Loading auxiliary files...
 cntry.codes.txt
 adjacency.txt
 ...Finished
 Total number of cross sections: 19

 ylc <- yourcast(formula=log(rspi2/popu2) ~ time, dataobj=dta,
 model="LC",debug =TRUE)
 Error in seq.default(smpvec[ind], smpvec[ln], by = 1) :
   wrong sign in 'by' argument

 Thank you for your help!

 Antonio Lam