similar to: Survreg function for loglogistic hazard estimation

Hi Dear, Can someone please inform me on genreating random variables of Loglogistic, and PERT beta distributions? Thanks for your time and help, in advance. Regards, Murthy.

Hi, I'm estimating a loglogistic aft (accelerated failure time) model, just a simple plain vanilla one (without time dependent covariates), I'm comparing the results that I obtain between aftreg (eha package) and survreg(surv package). If I don't use any covariate the results are identical , if I add covariates all the coefficients are the same until a precision of 10^4 or 10^-5 except

Hi, is there any similar function in R to S function qq.loglogistic, which produces a Q-Q plot? Thanks a lot Pete

Dear all, What is the cumulative distribution (with parameterization) used within survreg with respect to the log-logistic distribution? That is, how are the parameters linked to the survivor function? Best regards, Mario [[alternative HTML version deleted]]

Hola a tod@s, Hace poco en la lista se discutieron algunas aproximaciones para determinar la distribución de probabilidad que potencialmente podrían haber generado los datos (ver [1]) y una de ellas fué el AIC. Haciendo uso del programa enviado por Pablo Verde (ver [1]), estos son los resultados para mis datos: weibull , AIC = 69839.44 exponential , AIC = 79488.77 gaussian , AIC = 69413.03

I a trying to use survreg to fit a Weibull distribution. From the last few messages I understand how to interpret the parameters. However, how do I get the covariance of lambda and alpha? Is there a predict command for that? Thanks -- View this message in context: http://n4.nabble.com/Survreg-function-for-loglogistic-hazard-estimation-tp893503p1604918.html Sent from the R help mailing list

Hi list! I'm getting an error message when trying to fit an accelerated failure time parametric model using the aftreg() function from package eha: > Error in optim(beta, Fmin, method = "BFGS", control = list(trace = > as.integer(printlevel)), : > non-finite finite-difference value [2] This only happens when adding four specific covariates at the same time in the

Dear all, I am trying to find a suitable R-function for parametric proportional hazard regressions. The package survival contains the coxph() function which performs a Cox regression which leaves the base hazard unspecified, i.e. it is a semi-parametric method. The package Design contains the function pphsm() which is good for parametric proportional hazard regressions when the underlying base

Hi I'm trying to predict the values for a survreg object called loglogistic_na. Here is the definition of loglogistic_na and following that the syntax used for the predict function. But upon execution I don't get any output. Not sure what I'm doing wrong: loglogistic_na <- survreg(Surv(time_na,event_na) ~ t_na, dist="loglogistic") summary(loglogistic_na)

Dear all, I am looking for routines which allow multi-variate sampling from non-normal distributions (loglogistic) given correlations among the variables. Unfortunately, I could not find a suitable package for R. Does anybody know one? Many thanks and best regards, Stefan Albrecht [[alternative HTML version deleted]]

I am used to using the replace command to quickly replace strings in file, but it's not available on some of my fresh CentOS 5.4 servers. "yum info replace", "yum whatprovides replace", and "yum search replace" doesn't show me which package(s) to install to get it. So, does anyone know which package to install to get the "replace" command? Google

Greetings, Version 0.1.4-0 of the ChainLadder package is available via CRAN and http://code.google.com/p/chainladder/ The ChainLadder package, which is focused on claims reserving methods typically carried out by property/casualty insurance actuaries, has recently been enhanced to implement the methods in David Clark's 2003 CAS (Casualty Actuarial Society) *Forum* paper "LDF

On Wed, Aug 1, 2012 at 3:08 PM, <fra.meucci@hotmail.it> wrote: > Dear Göran Broström, > I am trying to use AFTREG function for R to estimate a loglogistic > survival function, including time dependent covariates. > Actually, my Subset includes some partial events; the idea is to model > this kind of events using something similar to “weights” in the SURVREG > function.

This kind of error seems to surprise R users. It surprises me that it doesn't happen much more frequently. The "BFGS" method of optim() from the 1990 Pascal version of my book was called the Variable Metric method as per Fletcher's 1970 paper it was drawn from. It really works much better with analytic gradients, and the Rvmmin package which is an all-R version that adds bounds

Hello. I am not certain even how to search the archives for this particular question, so if there is an obvious answer, please smack me with a large halibut and send me to the URLs. I have been experimenting with fitting curves by using both maximum likelihood and maximum spacing estimation techniques. Originally, I have been writing distribution-specific functions in 'R' which work

Hi, I'm doing a stats project using R to work out the size of a t-test and wilcoxon test depending on the distribution and sample size. I just can't get it to work - I want to put my results from the function size() into an array.At the moment I keep getting the error message:Error in res[distribution, test, samplesize] <- results : subscript out of boundsCan anyone tell me where

Hi, I' am doing a stats project using R to work out the size of a t-test and wilcoxon test depending on the distribution and sample size. I just can't get it to work - I want to put my results from the function size() into an array.At the moment I keep getting the error message:Error in res[distribution, test, samplesize] <- results : subscript out of boundsCan anyone tell me where

Hi, I'm trying to get better at things like lapply but it still stumps me. I have a function I've written, tested and debugged using individual calls to the function, ala: ResultList5 = DoAvgCalcs(IndexData, Lookback=5, SampleSize=TestSamples , Iterations=TestIterations ) ResultList8 = DoAvgCalcs(IndexData, Lookback=8, SampleSize=TestSamples , Iterations=TestIterations ) ResultList13

Here is the function that makes the data.frames in the list: funweek <- function(df) if (length(df$elapsed_time) > 5) { res = fitdist(df$elapsed_time,"exp") year = df$sale_year[1] sample = df$sale_week[1] mid = df$m_id[1] estimate = res$estimate sd = res$sd samplesize = res$n loglik = res$loglik aic = res$aic bic = res$bic chisq =

Hi Hugo, I've been able to replicate your bug, including for other distributions (runif, rexp, rgamma, etc) which shouldn't be surprising since they're probably all drawing from the same pseudo-random number generator. ?Interestingly, it does not seem to depend on the choice of seed, I am not sure why that is the case. I'll point out first of all that the R-devel mailing list is