similar to: Package "demography" - calculating percentiles of survival probabilities distribution

Dear All, I wonderif there is a simple way to draw a regression line in the xyplot: more specifically, let: age <- c(20:30, 31:40 ) age.cut <- cut(age, breaks = 2 ) y<- rnorm(20) x <- rnorm(20,4,1) xyplot(y ~ x| age.cut, xlab="x", ylab="y") How to draw (in the plot given by xyplot) the two regression lines (y ~ x) corresponding to the two category

Good day all, I'm trying to plot a continuous line plot, which is followed by two forecast points eg. one forecast point is 12 months out, and another 24 months out from the last date of the line plot. In my attempts so far, the second plot (the forecast points) is scaled against a new axis scale, thus the two plots are not directly comparable (I need the forecast points to be scaled

1. This is an imaginary data on monthly outcomes of 2 years and I want to forecast the outcome for next 12 months of next year. data Data1; input Yr Jan Feb Mar Apr May June July Aug Sept Oct Nov Dec; datalines; 2008 12 13 12 14 13 12 11 15 10 12 12 12 2009 12 13 12 14 13 12 11 15 10 12 12 12 ; run; I converted the above data into the below format to use it in R as it was giving error: asking

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

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’,

Dear R users, I came up to a problem dealing with percentiles in R. >From my previous questions: I do have a big data.frame, with lots of columns and rows. The following command enables me to calculate means for all data frame. dat1$newID<-rep(1:(nrow(dat1)/12),each=12) #if nrow(dat1)/12 is integer dat2<-with(dat1,aggregate(cbind(dat1[,1:71]),by=list(newID),mean)) What I need is to

Dear List, I have some discrete data and want to calculate the percentiles and the percentile ranks for each of the unique scores. I can calculate the percentiles with quantile(). I know that "ecdf" can be used to calculate the empirical cumulative distribution. However, I don't know how to exact the cumulative probabilities for each unique element. The requirement is similar

Hi All: I am working with a dataset on Arsenic toxicity, and I am trying to calculate the 20th, 40th, 60th, 80th, and highest percentiles for a variable, dietary Moisture (variable name dMoist). The inbuilt function quantile(dMoist) would print 0, 25th, 50th, 75th, and 100th percentile. Does there exist a function that can calculate xth percentile (where x = 10th, 20th, ... etc) values?

Type e.g.: quantile(x,0.1) or Quantile(x,0.8) Which calculates the 10th and 80th quantile Matthias Templ -----Urspr??ngliche Nachricht----- Von: arinbasu at softhome.net [mailto:arinbasu at softhome.net] Gesendet: Montag, 26. April 2004 09:29 An: r-help at stat.math.ethz.ch Betreff: [R] Looking for help in calculating percentiles Hi All: I am working with a dataset on Arsenic toxicity,

Hello, I'm trying to get percentiles (PERCENTRANK for excel users) by factor in the following data.frame: myExample <- data.frame(Ret=seq(-2, 2.5, by=0.5),PE=seq(10,19),Sectors=rep(c("Financial","Industrial"),5)) myExample <- na.omit(myExample) Thanks to Patrick I I managed to put together the following lines which does it for the "Ret" column: myecdf

I was hoping to calculate the percentile for each date. So group all of one date together, calculating the value of the 2.5 and 97.5 percentile. -- View this message in context: http://r.789695.n4.nabble.com/Calculating-percentiles-for-multiple-dates-tp4638183p4638378.html Sent from the R help mailing list archive at Nabble.com.

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

Dear R-list, I try to calculate the 10- and 90-Percentile from grouped data. Using 'lappy' it works fine but I have no success with 'gapply': df <- data.frame(test = runif(100)) df <- df[order(df[,1]),] z <- rep(1:10, each = 10) lapply(split(df, z), function(x) quantile(x, probs=c(10,90)/100)) library(nlme) gapply(as.data.frame(df), which=1, function(x)

How do I calculate the 95th percentile when I know the 25th, the median and the 75th?? Thanks, Harmony Tenney [[alternative HTML version deleted]]

The standard R boxplot appears to use quartiles to determine the height of the rectangles and a range parameter - RNG - (default=1.5 I think) that determines the length of the whiskers as <= RNG x Interquartile Range. Is it possible to instead define the range as extending to the 95th percentile? If so, how would this be done? nb, I'm plotting multiple boxplots on a single chart so a

Hi there! Still struggling to translate Matlab code into R's tsDyn package. Here is my question: Is there in R an equivalent function to Matlab's prctile()? To the moment I thought it was quantile(), but I just realized I was wrong. The definition of the Matlab function: prctile Percentiles of a sample SyntaxY = prctile(X,p) Description Y = prctile(X,p) calculates a value that is

Hi All, I am looking at applications of percentiles to time sequenced data. I had just been using the quantile function to get percentiles over various periods, but am more interested in if there is an accepted (and/or R-implemented) method to apply weighting to the data so as to weigh recent data more heavily. I wrote the following function, but it seems quite inefficient, and not really very

Hi all, I'm just getting started with R and I would appreciate some help. I'm having trouble creating a boxplot with whiskers at the 95th and 5th percentiles instead of at 1.5 * IQR. I have read the relevant documentation, and checked existing mails on this topic. I found a small modification that should work : https://stat.ethz.ch/pipermail/r-help/2001-November/016817.html and tried to

I have a problem where I need to calculate the corresponding cohort percentile ranks for each of several variables. Essentially, what I need is a function that will calculate the distribution-free percentiles from each variable's data vector, returning a corresponding vector of percentiles: e.g.: percentile.my.data<-/function/(my.data) I tried to make ecdf() perform this task but