Halabi, Anan
2010-Sep-21 05:12 UTC
[R] Estimating Weibull parameters by maximum liklihood - with censored and non censored data?
Below my code: sampleSize <- 20 shape.true <- 1.82 scale.true <- 987 sampWB <- rweibull(sampleSize, shape=shape.true, scale=scale.true) print(sampWB) censidx <- sample(1:length(sampWB), length(sampWB)*0.3) Censored.data <- sampWB[censidx] print(Censored.data) NonCensored.data <- sampWB[-censidx] print(NonCensored.data) how can I estimte weibull parameters by maximum liklihood with censored and non censored?
Bill.Venables at csiro.au
2010-Sep-21 05:30 UTC
[R] Estimating Weibull parameters by maximum liklihood - with censored and non censored data?
library(help = survival)
help("survreg", package = "survival")
-----Original Message-----
From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
On Behalf Of Halabi, Anan
Sent: Tuesday, 21 September 2010 3:12 PM
To: r-help at r-project.org
Subject: [R] Estimating Weibull parameters by maximum liklihood - with censored
and non censored data?
Below my code:
sampleSize <- 20
shape.true <- 1.82
scale.true <- 987
sampWB <- rweibull(sampleSize, shape=shape.true, scale=scale.true)
print(sampWB)
censidx <- sample(1:length(sampWB), length(sampWB)*0.3)
Censored.data <- sampWB[censidx]
print(Censored.data)
NonCensored.data <- sampWB[-censidx]
print(NonCensored.data)
how can I estimte weibull parameters by maximum liklihood with censored and non
censored?
______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.