search for: mortality

i am trying to replicate the following graph using xyplot : attach(x) plot ( jitter(type), mortality, pch=16, xlim = c(0.25, 3.75)) lines ( c(1-0.375,1.375) , c ( median(mortality[type==1]), median(mortality[type==1])), lwd=5,col=2) lines ( c(2-0.375,2.375) , c ( median(mortality[type==2]), median(mortality[type==2])), lwd=5,col=2) lines ( c(3-0.375,3.375) , c ( median(mortality[type==3]), median...

...ened for the sake of example). Gender Age M 70 F 65 M 70 Each row represents a person with an age/gender combination. We could put this data into a data frame. Now, I would like to do some actuarial analysis on this data set. To do so, I need to create and store a mortality curve for each person in the table (a mortality curve is a matrix with 2 columns: date and survival probability). I can write a function that returns a mortality curve given gender and age. The question is the following: In what data format should I store all these mortality curve objects? Should...

...interesante, sin embargo cuando mandamos a leer desde la tabla csv original de los autores no hace el análisis porque algo falta y no nos percatamos de que es. Adjunto tabla original Saludos cordiales #ESTO ES LO QUE VIENE EN LA AYUDA Y TRABAJA BIEN rm(list = ls()) library(FactoMineR) data(mortality) library(epicalc) des(mortality) #Agregue esto para comparar las dos datas res<-MFA(mortality,group=c(9,9),type=c("f","f"), # REALIZA EL ANÁLISIS BIEN name.group=c("1979","2006")) plot(res,choix="freq",invisible="ind",hab...

Hola. Como te dijo Carlos, el problema está en los nombres de las columnas y en los nombres de las filas. Cuando hice la importación (con dd<-read.csv('mortality.csv'), tuve problemas con las filas de nombre: - Malignant tumour of the larynx trachea bronchus and lungs - Malignant tumour of the lip pharynx and mouth - Other endocrinological metabolic and nutritional conditions que se me corrían a la derecha creándome una columna más. Corregido...

Dear all, I'm trying reproduce an example of bioassays Yielding Concentration-Mortality Data particularly control - adjustment model from book Bioassay of Entomopathogenic Microbes and Nematodes chapter 7 with R. I used glm with family=binomial and link=probit, but I do not know how to implement parameter gamma (control mortality - mortality of the untreated control insect in this ex...

Dear R Users: I have the individual mortality rate and 95% CI of 100 hospitals, how to do the plot with the individual hospital in the Yaxis, and the mortality rate and 95% CI in the Xais and a overall mean as a reference line? Thanks and regards, Xin [[alternative HTML version deleted]]

...re..06." > head(mort$X) [1] Accidental poisoning [2] Addiction to prescription medication [3] Alcohol abuse and alcohol-related psychosis [4] Asthma [5] Blood and hematopoietic disorders [6] Cerebrovascular disease 65 Levels: 1109 1869 36 ... Viral hepatitis En cambio mira el dataset "mortality": > names(mortality) [1] "15-24 (79)" "25-34 (79)" "35-44 (79)" [4] "45-54 (79)" "55-64 (79)" "65-74 (79)" [7] "75-84 (79)" "85-94 (79)" "95 and more (79)" [10] &...

Hi all, I am not sure how to handle interactions with categorical predictors in the GAM models. For example what is the different between these bellow two models. Tests are indicating that they are different but their predictions are essentially the same. Thanks a bunch, > gam.1 <- gam(mortality.under.2~ maternal_age_c+ I(maternal_age_c^2)+ + s(birth_year,by=wealth) + + + wealth + sex + + residence+ maternal_educ + birth_order, + ,data=rwanda2,family="binomial") > > gam.2 <- gam(mortality.under.2~ maternal_age_c+...

Hello, all. I have co-opted a number of functions that can be used to plot the hazard/survival functions and associated density distribution for a Gompertz mortality model, given known parameters. The Gompertz hazard model has been shown to fit relatively well to the human adult lifespan. For example, if I wanted to plot the hazard (i.e., mortality) functions: pop1 <- function (t) { x=c(0.03286343, 0.04271132) a3<-x[1] b3<-x[2] shift<-15...

Hi i am doing an environmental research The equation is as follow: gam(y1 ~ x1 + s(x2) + s(x3) + s(x4), family = gaussian, fit = true) I would like to obtain the beta coefficient and 95CI of x4 (or s(x4)), what should I do? Thanks, Lung -- View this message in context: http://r.789695.n4.nabble.com/mgcv-beta-coefficient-and-95-CI-tp3320491p3320491.html Sent from the R help mailing list

...;red"), lty = c(1,2), lwd=c(1,2), xlab="Time", ylab="Survival probability") legend(10, 0.2, c("males","females"), col = c("black","red"), lty = c(1,2), lwd=c(1,2),bty='n') I would like to make a similar plot of the cumulative mortality hazard over time for both males and females with associated confidence intervals. I would like to make sure I obtained proper estimates of the cumulative hazard for both males and females with associated confidence intervals but I have doubts regarding my calculations # I used basehaz ?H0? to get...

Hi All, I am relatively new to R (took a class 2 yrs ago ...), and am hoping someone can point me in the right direction for a problem I'd like to solve in R : I would like to calculate the relative mortality of a particular impairment, relative to the standard population. I.e. I'm trying to find S_relative(t) in the eqn below : S_impar(t) = S_standard(t) * S_relative(t) where S_standard(t) is the mortality of the general population, and S_impar is the mortality of the the impaired group under...

...Xkids = 0.232, Xteens = 0.375, Xadults = 0.235, Xaged = 0.52, cases = sample(x=100:500,size=10), conc = sample(x=20:125,size=10) ) ## Two of the published dose-response relationships adult.CP.mortality = list(end.point = "Cardiopulmanory mortality in adults over 30", pollutant = "PM10", relationship = "log-linear", beta = c(0.0562,0.1551,0.2541), Xpop = c("Xba...

Hi I am trying to carry out some maximum likelihood estimation and I'm not making much headway, and I'm hoping that someone will be able to point me in the right direction. I am modelling mortality statistics. One way to do this is to model the mortality rate (or, more accurately, log of the mortality rate, log_m) as (say) a constant plus a proportion of age, plus time, so: r_1 <- lm(formula=log_m ~ age + time) summary(r_1) However, an alternative approach is to use try and estimate the...

Hi, In my analysis of impacts of insecticide-treated bednets on malaria, I look at the relationship between malaria incidence and mosquito behaviors. The condensed data set is copied here. Ordinary regression (lm) shows that Incidence was negatively related to Mortality. This makes sense because the latter reflected the strength of killing mosquitoes by insecticide-treated nets. Since the original data set has a complex structure with more parameters and scenarios. I guess a tree model would help explore the structure of the data. However, regression tree (rpart(...

Hello, I have placed a package on CRAN about two surveys of mortality in Iraq that were published in the Lancet. http://cran.at.r-project.org/src/contrib/Descriptions/lancet.iraqmortality.html > install.packages("lancet.iraqmortality") ... > library(lancet.iraqmortality) Loading required package: foreign > ?lancet.iraqmortality > vignette(&q...

Hello, I have placed a package on CRAN about two surveys of mortality in Iraq that were published in the Lancet. http://cran.at.r-project.org/src/contrib/Descriptions/lancet.iraqmortality.html > install.packages("lancet.iraqmortality") ... > library(lancet.iraqmortality) Loading required package: foreign > ?lancet.iraqmortality > vignette(&q...

Hi everyone, We are looking for some data sets working with relative risk mortality. so, someone know where can I find the data.mgus dataset and the data.mgus? Using 1384 records from Minnesota. This data set are used in the : Robert A. Kyle, Terry M. Therneau, S. Vincent Rajkumar, Janice R. Offord, Dirk R. Larson, Matthew F. Plevak, and L. Joseph Melton III. A long-term study...

I have a mortality data over many years and I wish to plot the data and also add some smoother to clearly highlight the trend. How could I do that in R with base graphics or ggplot? I have the following sample data: require(lubridate) mdate<-seq(ymd('2000-01-01'),ymd('2010-12-31'), by = '1 day...

V&R describes binomial GLMs with mortality out of 20 budworms. Is it appropriate to use the same approach with mortality out of numbers as low as 3? I feel reticent to do so with data that is not very continuous. There are one continuous and one categorical independent variables. Would it be more appropriate to treat the response as an...