similar to: estimation of lambda and gamma with std errors for a weibull model

Dear R users, Please can you help me with a relatively straightforward problem that I am struggling with? I am simply trying to plot a baseline survivor and hazard function for a simple data set of lung cancer survival where `futime' is follow up time in months and status is 1=dead and 0=alive. Using the survival package: lung.wbs <- survreg( Surv(futime, status)~ 1, data=lung,

I was wondering if someone familiar with survival analysis can help me with the following. I would like to fit a Weibull curve, that may be dependent on a covariate, my dataframe "labdata" that has the fields "cov", "time", and "censor". Do I do the following? wieb<-survreg(Surv(labdata$time, labadata$censor)~labdata$cov,

Dear All, I would like to fit some parametric survival models using left truncated, right censored data in R. However I am having problems finding a function to fit parametric survival models which can handle left truncated data. I have tested both the survreg function in package survival: fit1 <- survreg(Surv(start, stop, status) ~ X + Y + Z, data=data1) and the psm function in package

Hello I would like to compare the results obtained with a classical non parametric proportionnal hazard model with a parametric proportionnal hazard model using a Weibull. How can we obtain the equivalence of the parameters using coxph(non parametric model) and survreg(parametric model) ? Thanks Virginie

Dear R-users, I need to use the aftreg function in package 'eha' to estimate failure times for left truncated survival data. Apparently, survreg still cannot fit such models. Both functions should be fitting the accelerated failure time (Weibull) model. However, as G?ran Brostr?m points out in the help file for aftreg, the parameterisation is different giving rise to different

Dear R-users I would like to fit weibull parameters using "Method of moments" in order to provide the inital values of the parameter to de function 'fitdistr' . I don`t have much experience with maths and I don't know how to do it. Can anyone please put me in the rigth direction? Borja [[alternative HTML version deleted]]

Hi All, I have conducted the following survival analysis which appears to be OK (thanks BRipley for solving my earlier problem). > surv.mod1 <- survreg( Surv(timep1, relall6)~randgrpc, data=Dataset, dist="weibull", scale = 1) > summary(surv.mod1) Call: survreg(formula = Surv(timep1, relall6) ~ randgrpc, data = Dataset, dist = "weibull", scale = 1)

Hi, I was wondering if someone can help me interpret the results of running weibreg. I run the following and get the following R output. > weibreg(Surv(time, censor)~covar) fit$fail = 0 Call: weibreg(formula = Surv(time, censor)~covar) Covariate Mean Coef Rel.Risk L-R p Wald p covar 319.880 -0.002 0.998 0.000 log(scale) 0.000 8.239

Hello, I'm quite new to R and want to make a Weibull-regression with the survival package. I know how to build my "Surv"-object and how to make a standard-weibull regression with "survreg". However, I want to fit a translated or 3-parametric weibull dist to account for a failure-free time. I think I would need a new object in survreg.distributions, but I don't know how

Hi I have been asked to provide Weibull parameters from a paper using Kaplan Meir survival analysis. This is something I am not familiar with. The survival analysis in R works nicely and is the same as commercial software (only the graphs are superior in R). The Weibull does not and produces an error (see below). Any ideas why this error should occur? My approach may be spurious.

Dear All, I have two questions regarding distribution fitting. I have several datasets, all left-truncated at x=1, that I am attempting to fit distributions to (lognormal, weibull and exponential). I had been using fitdistr in the MASS package as follows: fitdistr<-(x,"weibull") However, this does not take into consideration the truncation at x=1. I read another posting in this

Hi I am using R to fit a survival function to my data (with a weibull distribution). Data: Survival of individuals in relation to 4 treatments ('a','b','c','g') syntax: ---- > survreg(Surv(date2)~males2, dist='weibull') But I have some problems interpreting the outcome and getting the parameters for each curve. --------- Value Std.

Dear R-users I have some datasets, all left-censoring, and I would like to fit distributions to (weibull,exponential, etc..). I read one solution using the function survreg in the survival package. i.e survreg(Surv(...)~1, dist="weibull") but it returns only the scale parameter. Does anyone know how to successfully fit the exponential, weibull etc... distributions to left-censoring

I am trying to generate predictions from a weibull survival curve but it seems that the predictions assume that the shape(scale for survfit) parameter is one(Exponential but with a strange rate estimate?). Here is an examle of the problem, the smaller the shape is the worse the discrepancy. ### Set Parameters scale<-10 shape<-.85 ### Find Mean scale*gamma(1 + 1/shape) ### Simulate Data

Hi all I try to make a weibull survival analysis on R. I know make this on GLIM, and now I try to make the GLIM exercice GLEX8 on R to learning and compare the test. The variables are: time censor group bodymass In GLIM I make: $calc %s=1 $ to fit weibull rather than exponential $input %pcl weibull $ $macro model group*bodymass $endmac$ $use weibull t w %s $ Then, GLIM estimate an alpha for the

Dear r-help subscribers, A couple of weeks ago I sent the following message to the r-help mail list. It hasn't generated any response, and I could really use some help on this. Anyone able to help? Thanks again, Roger Dungan >> I am working on some survival analysis of some interval censored failure time data in R. I have done similar analysis before using PROC LIFEREG in SAS. In

Hello! I have a matrix with data and a column indicating whether it is censored or not. Is there a way to apply weibull and exponential maximum likelihood estimation directly on the censored data, like in the paper: Backtesting Value-at-Risk: A Duration-Based Approach, P Chrisoffersen and D Pelletier (October 2003) page 8? The problem is that if I type out the code as below the likelihood

Hi all, I'm trying to generate a Weibull distribution including four covariates in the model. Here is the code I used: T = rweibull(200, shape=1.3, scale=0.004*exp(-(-2.5*b1+2.5*b2+0.9*x1-1.3*x2)/1.3)) C = rweibull(n, shape=1.5, scale=0.008) #censoring time time = pmin(T,C) #observed time is min of censored and true event = time==T # set to 1 if event is observed

Dear R-users! I would like to fit exponential, Weibull and log-logistic via glm() like functions. Does anyone know a way to do this? Bellow is a bit longer description of my problem. Hm, could family() be adjusted/improved/added to allow for these distributions? SAS procedure GENMOD alows to specify deviance and variance functions to help in such cases. I have not tried that option and I do not

Hi, I am trying to fit a non-linear model for a parasite dataset. Initially, I tried log-transforming the data and conducting a 2-way ANCOVA, and found that the equal variance of populations and normality assumptions were violated. Gaba et al. (2005) suggests that the Weibull Distribution is best for highly aggregated parasite distributions, and performs better (lower type 1 and 2 error rates)