similar to: generating loglogistic distribution in R

I am trying to use R to do loglogistic hazard estimation. My plan is to generate a loglogistic hazard sample data and then use survreg to estimate it. If everything is correct, survreg should return the parameters I have used to generate the sample data. I have written the following code to do a time invariant hazard estimation. The output of summary(modloglog) shows the factor loading of

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 every one, I am trying to generate a normally distributed random variable with the following descriptive statistics, min=1, max=99, variance=125, mean=38.32, 1st quartile=38, median=40, 3rd quartile=40, skewness=-0.274. I know the "rnorm" will allow me to simulate random numbers with mean 38.32 and Sd=11.18(sqrt(125)). But I need to have the above mentioned descriptive

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

Dear All, Can someone please tell me if there is a provision in R to fit a random coefficient multinomial logistic regression. Thanks in advance for your help and suggestions. Regards Murthy.N.M

Hi, Can Someone inform me how to generate RV's using the below CDF, by inverse technique. Thanks for your help and time. My CDF is as follows \[ F(x)=0 \ \text{if} \ x < 0\]\[ F(x)=\{\frac{x-x_i}{x_{i+1}-x_{i}}*(p_{i+1}-p_{i})\}+p_{i}\ \forall \ x_{i}\leq x < x_{i+1} \] \[ F(x)=1 \ \text{if} \ x > x_{i+1} \] Regards Murthy

Dear List: We are generating data such that students are clustered in schools for some item response data for a simulation study. One component of our simulation is to generate measurement error from a logistic distribution with a mean of 0 and standard deviation of 1.7 to match the logistic curve of the Rasch model. We are generating an error term for each of the 40 hypothetical test items a

We need help to get Micro Planner X-Pert 3.5.1 running on Wine on Ubuntu and are willing to pay for this help. If you can help there is a free trial download of Micro Planner X-Pert available for download from www.microplanning.com We are running Ubuntu 11.04 if we can get this sorted we will be able to kiss windows good bye Hope someone can help. Jim Radford

Hi! I was wondering if PERT or CPM was implemented in R. I looked in the search engines but didn't find anything. Since there are so many packages, I thought I'd double check via the discussion. Thanks, Laura. [[alternative HTML version deleted]]

Hi any, Can some please detail me the createX command in bayesm package? To make things easy for you to help me, let me put forward my problem Suppose I have 3 covariates (say X matrix) and my Y has 3 categories say (1,2,3). Now from the CreateX I understand that the data matrix say 'Xa' must be of dimension n* (naxp), where 'na' is the number of variables and 'p' is

Hello, is it possible to create two uncorrelated random vectors for a given distribution. In fact, I would like to have something like the function "rnorm" or "rlogis" with the extra property that they are uncorrelated. Thanks for your help, Luba [[alternative HTML version deleted]]

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.