Displaying 20 results from an estimated 1000 matches similar to: "Computing Mean Lifetime from Hazard"
2008 Jul 31
2
Help with hazard plots
Hello. ?I am hoping someone will be willing to help me understand something about hazard plots created with muhaz(...). ?I have some background in statistics (minor in grad school), but I haven't been able to figure one thing about hazard plots. ?I am using hazard plots to track customer cancellations. ?I figure I can treat a cancellation as a "death", and if someone is still a
2008 Sep 05
1
Confidence Intervals on Hazard Plots
Hello, Is it possible to create confidence intervals for hazard rates? ?I'm creating two muhaz objects:?
haz1 <- muhaz(NumDaysCustomer[cRV=="true"],status[cRV=="true"])
haz2 <- muhaz(NumDaysCustomer[cRV=="false"],status[cRV=="false"])
and plotting them. ?
There are many, many more observations in the cohort cRV=="false" than
2007 Jan 23
1
Estimate and plot hazard function using "muhaz" package
Dear R users,
I am trying to use "muhaz" and "plot.muhaz" functions in "muhaz" package to
estimate and plot hazard funciton. However function "muhaz" always gives
error message "Error in Surv(times, delta) : object "times" not found". I
could not even run their sample codes in the user's manual as follows:
data(ovarian)
2006 Jun 12
0
non parametric estimates of the hazard with right censored data
Hi,
I want to plot non parametric estimates of the empirical hazard function for
right censored data. I've tried many functions from different packages
(muhaz, Design, survival, eha, event), but none of them gave me what I
wanted. Am I missing something?
Here's what I want. The data below is the same used by Kiefer (J. Economic
Literature, 1988), which in turns use a subset of the data
2003 Nov 05
1
Estimate hazard function from right-censored data only
Dear All,
I would like to ask if it is possible to estimate a hazard function
using the muhaz command when all the data is right-censored. My data
has information of the number of weeks people has been unemployed but
all of them are unemployed at the date of the survey, that is, I cannot
observed when the individuals leave the unemployment state.
I appreciate your help,
Best,
Monica L.
2009 Jun 20
1
Plotting Cumulative Hazard Functions with Strata
Hello:
So i've fit a hazard function to a set of data using
kmfit<-survfit(Surv(int, event)~factor(cohort))
this factor variable, "cohort" has four levels so naturally the strata
variable has 4 values.
I can use this data to estimate the hazard rate
haz<-n.event/n.risk
and calculate the cumulative hazard function by
H<--log(haz)
Now, I would like to plot this
2018 Apr 04
1
parfm unable to fit models when hazard rate is small
Hello, I would like to use the parfm package: https://cran.r-project.org/web/packages/parfm/parfm.pdfhttps://cran.r-project.org/web/packages/parfm/parfm.pdf in my work. This package fits parametric frailty models to survival data. To ensure I was using it properly, I started by running some small simulations to generate some survival data (without any random effects), and analyse the data using
2010 May 25
2
Relative Risk/Hazard Ratio plots for continuous variables
Dear all,
I am using Windows and R 2.9.2 for my analyses. I have a large dataset and
I am particularly interested in looking at time to an event for a continuous
variable. I would like to produce a plot of log(relative risk) or relative
risk (also known as hazard ratio) against the continuous variable.
I have spent a long time looking for advice on how to do this but my search
has proved
Looking for new maintainer of orphans R2HTML SemiPar cghseg hexbin lgtdl monreg muhaz operators pamr
2014 Aug 08
2
Looking for new maintainer of orphans R2HTML SemiPar cghseg hexbin lgtdl monreg muhaz operators pamr
Dear maintainers and R-devel,
Several orphaned CRAN packages are about to be archived due to
outstanding QC problems, but have CRAN and BioC packages depending on
them which would be broken by the archival (and hence need archiving
alongside).
Therefore we are looking for new maintainers taking over maintainership
for one or more of the following packages:
R2HTML SemiPar cghseg hexbin lgtdl
2011 Jan 26
0
baseline hazard function
Dear colleagues, I have the following dataset. It is modelled on the data included in Box-Seteffenheiser and Jones "Event History Modelling"
Using the following code, I try to find the baseline hazard function
haz_1<-muhaz(bpa$time, bpa$censored, subset=(bpa$year=="2010" | bpa$ban=="1"), min.time=1, max.time=3)
I think I'm doing everything right, but what I
2010 Mar 16
1
Problems compiling a simple source package on MacOS X
Hello. I wrote a package (that contains C source) that I've been
compiling and running on both Linux and Windows for about a year.
However, that same source package fails to compile on MacOS (10.4.11,
PowerPC G4, Xcode installed, gcc version 4.0.1, make version 3.80, ld
version cctools-590.23.2.obj~17).
There is nothing platform-specific in the code-- just numerical
functions that I
2009 Jun 03
2
Create a time interval from a single time variable
I am trying to set up a data set for a survival analysis with time-varying covariates. The data is already in a long format, but does not have a variable to signify the stopping point for the interval. The variable DaysEnrolled is the variable I would like to use to form this interval. This is what I have now:
ID Age DaysEnrolled HAZ WAZ WHZ Food onARV
2004 May 10
1
Explaining Survival difference between Stata and R
Dear Everybody:
I'm doing my usual "how does that work in R" thing with some Stata
projects. I find a gross gap between the Stata and R in Cox PH models,
and I hope you can give me some pointers about what goes wrong. I'm
getting signals from R/Survival that the model just can't be estimated,
but Stata spits out numbers just fine.
I wonder if I should specify initial
2001 Oct 11
2
Where's MVA?
Hi All:
Package TSERIES is stated to depend on MVA. However, there is no MVA package to be found under the list of package sources.
Best wishes,
ANDREW
tseries: Package for time series analysis
Package for time series analysis with emphasis on non-linear and non-stationary modelling Version: 0.7-6
Depends: ts, mva, quadprog
Date: 2001-08-27
Author: Compiled by Adrian
2014 Jan 30
2
objecto ggplot
Daniel,
Disculpa si no entiendo, pero, entonces un objecto ggplot en en realidad
una función
y no un gráfico
El 30 de enero de 2014, 11:09, Marta Garcia <marta000garcia@gmail.com>escribió:
> Sí, devuelve Error: No layers in plot, no le he agregado capas
> por eso el error, la funcion me dice que es un gráfico
>
>
> El 30 de enero de 2014, 10:31, José María Mateos
2012 Jan 05
1
Fwd: WHO Anthro growth curve macros and R&In-Reply-To=<CAAOCNNZawGtKkWpgFMYADSyxWGTeWEDxqVVHv7=Azo=1G+H9gg@mail.gmail.com>
Dear Gustaf,
I wish you a happy new year!
I am a PhD Staff at LSHTM and I have been trying to use WHO anthro macros in R with the file you posted on internet but when I double check the WHZ, HAZ, WAZ output with output from ENA or STATA, it differs slightly. Do you have any idea why that might be the case?
Many thanks
Severine
2007 Oct 23
2
2-D numerical integration over odd region
Hello all,
I'm hoping to find a way to evaluate the following sort of integral in R.
\int_a^b \int_{g(y)}^Inf f(x,y) dx dy.
The integral has no closed form and so must be evaluated numerically. The "adapt" package provides
for multidimensional integration but does not appear to allow the limits of integration to be a
function. I need to evaluate a number of integrals of this
2007 Jan 19
1
Error in basehaz function ?
Hello R-users.
I believe that the way basehaz (in the survival package) compute the
baseline hazard function is false.
I come to question this function when it gives me hazard probabilities
greater than 1.
Looking at the code I think I've localised the error :
hazard probability is computed as :
H <- -log(surv)
but it seems to me that hazard probabilities is rather an instantaneous
2011 Nov 30
1
Nomogram with stratified cph in rms package, how to get failure probability
Hello,
I am using Dr. Harrell's rms package to make a nomogram. I was able to make
a beautiful one. However, I want to change 5-year survival probability to
5-year failure probability.
I couldn?t get hazard rate from Hazard(f1) because I used cph for the model.
Here is my code:
library(rms)
f1 <- cph(Surv(retime,dfs) ~
age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11,
2011 May 18
1
Integral Symbol
Dear All,
I am documenting a R package. That means writing the *.Rd files inside the \man folder of the package structure
I was wondering how to write the symbol for an integral function in a formula.
Similar to this one in LaTeX:
\int_{0}^{10} \Omega(t)dt
I already tried
\deqn{\int_{0}^{10} \Omega(t)dt}
but it does not work. Any idea? Which math symbols does R-help recognise?
Regards,