search for: demographi

Hi, I am using the package "demography" from Rob Hyndman for the Lee-Carter-Model. It is an amazing powerful tool but I am struggling with one issue: *I want to compute different quintiles for the cumulative survival probability derived from the Lee-Carter-Forecast (e.g. the 50%-quintile, 75%-quintile and 99%-quintile) for the next 10 years. * I am sure the package possess this

Hi, I am using the package "demography" from Rob Hyndman for the Lee-Carter-Model. It is an amazing powerful tool but I am struggling with one issue: I want to compute different percentiles of the survival probability distribution derived from the Lee-Carter-Forecast (e.g. the 50%tile, 60%tile, 75%tile and 99%tile) for each of the next 10 years. Is there any possibility to retrieve

Dear list: I have tried demography package in my Ubuntu(10.04). What I interested in was lifetable analysis. Chiang(1984)'s reference have been copied to steeper my learning curve. step by step, the example presented by the textbook also have been examined. In the course I find that the paramter "fraction of last age intervl of life" could not be take as a option in the module

Buenos días, Necesito calcular la vida media de cada uno de nuestros modelos de producto pero tengo un histórico incompleto (para los productos más nuevos tengo la vida completa de las cohortes, pero para los más antiguos tengo sólo la parte final de la vida de las cohortes). He decidido calcular la esperanza de vida dinámica con un modelo actuarial, concretamente la librería 'demography'

Dear lists: I could not calculate what is wrong with my installing. Could you help me to find the solution? Thanks. Following are the installing screen massages: install.packages("demography") 将程序包安装入‘/home/zzk/R/i686-pc-linux-gnu-library/2.12’ (因为‘lib’没有被指定) --- 在此連線階段时请选用CRAN的鏡子 --- 载入Tcl/Tk接口... 完成 also installing the dependencies ‘akima’, ‘rgl’, ‘misc3d’, ‘quadprog’, ‘zoo’, ‘locfit’,

I've attached a perl script that I wrote to build smbpasswd from /etc/passwd. Please note that this writes a samba-2.0 format smbpasswd file. This script will: # - create accounts that exist in /etc/passwd and not in smbpasswd # - remove accounts that exist in smbpasswd and not in /etc/passwd # - preserve existing password, flags and LCT fields in smbpasswd # - sets new account passwords to

The generic stats::median method is defined as median <- function (x, na.rm = FALSE) {UseMethod("median")} I suggest that this should become median <- function (x, na.rm = FALSE, ...) {UseMethod("median")} This would allow additional S3 methods to be developed with additional arguments. Currently I have to over-ride this generic definition in the demography

I have some problem in running R-cran's Demography package. The hmd.mx function need Rcurl. I tried to install RCurl, but meet the following error: ********************************************************************* ... * installing *source* package ?RCurl? ... checking for curl-config... no Cannot find curl-config ERROR: configuration failed for package ?RCurl? * removing

I am trying to read a SAS dataset into R that is in a .SD2 format, which is the Version 6 standard format from SAS. I see the routines that read the SAS XPORT format (foreign, Hmisc), but is there any way to read this one? Any help would be greatly appreciated. Corey Sparks Corey Sparks Assistant Professor Department of Demography and Organization Studies University of Texas-San Antonio One UTSA

Greetings- We use smbclient -M to send a Message of the Day to clients every once in a while. The clients are NT4 (SP3), and the server is Samba 1.9.19-prealpha. However, the same problem appears with version 2.0.0. The issue is: once the file to be sent exceeds approximately 675 bytes (smbclient reports sending 975 bytes, presumably including some overhead), the message turns up blank on the

Hi, I have a data.frame with 100 variables and I have assigned a "label", "units" and "category" attribute to each variable. I would like to reorder the variables in the data.frame by the "category" attributes but can't find a way. For example, the first variable is: > attributes(hh$aez) $levels [1] "coastal" "forest"

I have trouble using function hmd.mx I have all packages required. I call function: slovenia <- hmd.mx("SVN", username, password, "Slovenia") : and I get this error: NAs introduced by coercion And if I then call slovenia, I get: > slovenia Mortality data for Slovenia Series: female male total Years: 1983 - 2009 Ages: 0 - 110 But I do not get any value

Dear list, I am trying to use gregexpr to see if entries in a dataframe have either of two possible values for a string. here's an example text<-c("fat", "rat", "cat", "dog", "log", "fish") If I just wanted to find if any one of the elements in text match the pattern "at" I would do gregexpr("\\at", text)

Dear All, I am fairly new to R. I work mainly in SAS. Now, I know that SAS is approved by the FDA for submissions. My question is, does the FDA approve {R} for clinical trial submissions. Also has anyone ever tried to produce TFL's using R. I would like to know how difficult it to produce the TFL's in R as compared to SAS. I know that in SAS it is not difficult once you know

> Dear list, > > I have some data on a geneaology, here is a subset: > warmerge[1:11,c(1,6,25)] > Warrior SibID birth.year > 1100 3793 2013 1926 > 4 2013 2024 1934 > 1094 3769 2024 1918 > 632 2747 2037 1928 > 176 2083 2039 1944 > 187 2085 2039 1949 > 192 2086 2039 NA > 495

Update (for me): * On HPUX (cc/f77/make), make check passes the Examples but fails in arith-true.R because is.infinite(.Machine$double.base ^ .Machine$double.max.exp) gives FALSE. Interestingly > .Machine$double.base [1] 2 > .Machine$double.max.exp [1] 1024 > 2 ^ 1024 [1] 1.797693e+308 > 2 ^ 1025 [1] 1.797693e+308 > 2 ^ 2024 [1] 1.797693e+308 so I guess we need help here. *

Dear R Users, I'm working on a problem where I have a multivariate response vector of counts and a continuous predictor. I've thought about doing this the same way you would do a Multvariate regression model with normally distributed data, but since these data are counts, they are probably better modeled with a Poisson distribution. For example y1<-rpois(100,3.5) y2<-rpois(100,1.5)

...rgl.viewpoint(coo[i]) #filename <- paste("pic", formatC(i,digits=1,flag="00"), ".png", sep="") #rgl.snapshot(filename) } ----------------------------------------------------- Dr. Carlo Giovanni Camarda Research Scientist Max Planck Institute for Demographic Research Laboratory of Statistical Demography Konrad-Zuse-Straße 1 18057 Rostock - Germany Phone: +49 (0)381 2081 172 Fax: +49 (0)381 2081 472 camarda@demogr.mpg.de http://www.demogr.mpg.de/en/staff/camarda/default.htm ----------------------------------------------------- ---------- This mail ha...

...lots persp3d(x,y,z1, color="#CCCCFF", front="lines") persp3d(x,y,z2, col="red", add=TRUE) play3d(spin3d(axis=c(0,0,1), rpm=10), duration=5) ----------------------------------------------------- Dr. Carlo Giovanni Camarda Research Scientist Max Planck Institute for Demographic Research Laboratory of Statistical Demography Konrad-Zuse-Straße 1 18057 Rostock - Germany Phone: +49 (0)381 2081 172 Fax: +49 (0)381 2081 472 camarda@demogr.mpg.de http://www.demogr.mpg.de/en/staff/camarda/default.htm ----------------------------------------------------- ---------- This mail...

Dear R users, I am interested in taking the columns from multiple dataframes, the problem is that the different dataframes have different combinations of the same variable names, here's a simple example: a<-rep(1:10) b<-rep(1:10) c<-rep(21:30) d<-rep(31:40) dat.a<-data.frame(a,b,c,d) names(dat.a)<-c("a", "b", "c", "d")