similar to: range of group variables

Dear list: I am attempting to use what I thought would be a pretty straightforward practical application of Cox regression. I figure users of the survival package must have come across this problem before, so I would like to ask you how you dealt with it. I have set up an illustrative example and included it at the end of this post. I took a sample of 100 data points from each of two populations

I am a relative newbie to survival analysis and R in general, but would like to use the coxme package to analyse some data I currently have. The data is relative to survival times of drosophila melanogaster populations to infection with pathogens, and has the variables: Time, Status, Treatment (4 treatments + 2 controls) Population Replicate ?and I'm currently using the following call

Dear List: I have a data frame prepared in the couting process style for including a binary time-dependent covariate. The first few rows look like this. PtNo Start End Status Imp 1 1 0 608.0 0 0 2 2 0 513.0 0 0 3 2 513 887.0 0 1 4 3 0 57.0 0 0 5 3 57 604.0 0 1 6 4 0 150.0 1 0 The outcome

Dear All, We need CentOS 5 for a software, and installed CentOS 5.4. But when we try to install Developments Tools like gcc it gives error that repositories are not found. Can you please help us in install gcc on CentOS 5.4? -bash-3.2$ cat /etc/redhat-release CentOS release 5.4 (Final) -bash-3.2$ uname -a Linux i80pc57 2.6.18-164.el5PAE #1 SMP Thu Sep 3 04:10:44 EDT 2009 i686 i686 i386 GNU/Linux

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]),

I am wondering if there is an R function that could estimate a generalized nonlinear mixed model. >From my reading it seems that nlme from the nlme package can fit nonlinear mixed models, while lmer from the lme4 package can fit generalized linear mixed models. One alternative I?ve found is gnlmix from the repeated package, although this only allows for a single random effect. Is there

Hello R Community, I have the following design question. I have a data set that looks like this (shortened 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

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

Estimados Queremos con el paquete FactoMineR hacer este tipo de tabla de mortalidad que lea los datos desde de una tabla csv Realizamos lo que viene en la ayuda y es muy 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

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) + +

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) {

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 exaple)

Hello all, I am hoping to use survival analysis to examine whether parasite attack increases nest death in a species of social wasp. I therefore have data for 1. Whether the nest "died" in the 6 week census period ("Status", where 1=died, 0=survived) 2. The day number of death/last recorded day it was observed alive. 3. Whether the nest was attacked by the parasite (0/1 as

Hallo Everybody How do you specify arguments for a function used within another function? Here is my problem: I am reconstructing a calculator for the burden of disease due to air pollution from publications and tools published by the WHO. The calculations make use of published dose-response relationships for particular health end-points. This is then applied to populations with known or

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]]

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

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

Hi guys, How could I specify minimum and maximum x-axis values in logi.hist.plot? My code is something like this: plot(mydata$Temperature,mydata$Mortality, ,xlab="Temperature",ylab="Probability of mortality") curve(predict(temp.glm,data.frame(Temperature=x),type="resp"),add=TRUE, col="red") points(mydata$Temperature,fitted(temp.glm),pch=20)

Hi... I'm trying to understand the following syntax: cor.test(~mortality + hardness,data=water,method="pearson") which is the same as: cor.test(water$mortality,water$hardness,data=water,method="pearson") Can anyone point me to the correct doc or explain to me how to interpret "~mortality + hardness"? Thanks, Joe

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 ordered factor with four levels? If so, what family