similar to: interpretation of flexsurvreg output from flexsurv package

I'm working on a survival analysis project and, based on previous work of my colleagues on similar projects, am trying to estimate an AFT model with a generalized gamma distribution. When I execute the following code: >flexsurvreg(timesurv ~ cycles + rcycle + cyctime, dist="gengamma") I get the following error message: Error in optim(optpars, minusloglik.flexsurv, t = Y[,

Dear R I am trying to understand and use the flexible parametric survival model suggested by Royston and Parmar. However I am stuck trying to plot the adjusted survival curves for different covariates in the following code: library(flexsurv) library(graphics) spl <- flexsurvspline(Surv(futime, fustat) ~ rx+ecog.ps+resid.ds+age, data = ovarian, k=2, scale="odds") spl the code

Hi all, Sorry if this has been answered already, but I couldn't find it in the archives or general internet. Is it possible to implement the gompertz distribution as survreg.distribution to use with survreg of the survival library? I haven't found anything and recent attempts from my side weren't succefull so far. I know that other packages like 'eha' and

Hi, After loading R 4.0.0 on Ubuntu 20.04. R version 4.0.0 (2020-04-24) -- "Arbor Day" Copyright (C) 2020 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()'

List gurus, I'm trying to code a Gompertz growth curve function as part of a larger project and have run across a problem due to my ignorance of environments. Some sample data and the function are as follows: growth <- data.frame(age = c(1.92, 3, 5.83, 3.17, 15.5, 1.17, 5.58, 13.33, 14.29, 5.83, 13.79, 6.33, 13.75, 16.83, 13, 11.67, 0.25, 1.73, 9.46, 5.67), length = c(157, 165, 179,

Hello, all. I have co-opted a number of functions that can be used to plot the hazard/survival functions and associated density distribution for a Gompertz mortality model, given known parameters. The Gompertz hazard model has been shown to fit relatively well to the human adult lifespan. For example, if I wanted to plot the hazard (i.e., mortality) functions: pop1 <- function (t) {

Full_Name: Alberto Gaidys Version: 2.9.0 OS: Mac OS X 10.5.7 Submission from: (NULL) (201.81.185.155) When calling LM or GLM from within a function, R gives a message error that it can't find the specified weights object "Erro em eval(expr, envir, enclos) : objeto 'W' n?o encontrado" (Error in eval(expr, envir, enclos) : object 'W' not found). The call from within

Dear list: I am attempting to use what I thought would be a pretty straightforward practical application of Cox regression. I figure users of the survival package must have come across this problem before, so I would like to ask you how you dealt with it. I have set up an illustrative example and included it at the end of this post. I took a sample of 100 data points from each of two populations

Dear Sir/ Madam, I'm having trouble de-bugging the following - which works perfectly well with optim or optimx - but not with mle2. I'd be really grateful if someone could show me what is wrong. Many thanks in advance. JSC: gompertz<- function (x,t=data) { a3<-x[1] b3<-x[2] shift<-data[1] h.t<-a3*exp(b3*(t-shift))

Hi, r-help members. I have a question about summing two density distributions. I have two samples from which I've estimated hazard parameters for a Gompertz mortality model. With those parameters, I can calculate the PDF (survival function times hazard function) of ages-at-death in a birth cohort subject to the hazard function at each age. I'd like to combine these two density

Hi, all. I'm working with published paleodemographic data (counts of skeletons that have been assigned into an age-range category, based upon observed morphological characteristics). For example, the following is the age distribution from a prehistoric cemetery in Egypt: naga <-

Dear R-users, I would like to show you a simple example that gives an overview of one of my current issue. Although my working setting implies a different parametric model (which cannot be framed in the glm), I guess that what I'll get from the following example it would help for the next steps. Anyway here it is. Firstly I simulated from a series of exposures, a series of deaths (given a

> From: Thomas Lumley > > On Sun, 30 Oct 2005, Jonathan Rougier wrote: > > > I'm not sure about this. Perhaps I am a dinosaur, but my feeling is > > that if people are writing functions in R that might be subject to > > simple operations like outer products, then they ought to be writing > > vectorised functions! > > I would agree. How about an

Dear All, I want to use grofit package for biological growth curves. My dataset only includes "age" variable and "size" variable. I want to use logistic, gompertz and richards growth curves to predict age from size. How can I implement this data set to the function in grofit package? Best regards, Deniz [[alternative HTML version deleted]]

Hello, I want to fit a Gompertz model for tree diameter growth that depends on a 4 levels edaphic factor (?Drain?) and I don?t manage to introduce the factor variable in the formula. Dinc is the annual diameter increment and D is the Diameter. >treestab > Dinc D Drain [1,] 0.03 26.10 2 [2,] 0.04 13.05 1 [3,] 0.00 24.83 1 [4,] 0.00 15.92 4

Dear all i want to fit the self start model in nls. i have two question. i have a function, (asfr ~ I(((a*b)/c))+ ((c/age)^3/2)+ exp((-b^2)*(c/age)+(age/c)-2) i am wondering how to build the selfstart model. there is lost of example, (i.e. SSgompertz, SSmicman, SSweibull, etc). my question is, how to derive the function of parameters. and also which model to use for get the initials values. In the

Hi, I have a problem with fitting a nonlinear mixed-effects model to my data. I was able to fit quite simple exponential model but now I?m trying to fit the following Gompertz model: > gomp<-nlme(paino~b0+b1*exp(-exp(-b2*(daydeg.scale-b3))), data=group1, fixed=b0+b1+b2+b3~1, start=c(b0=150, b1=3000, b2=0.6, b3=2500)) I?m getting the following error: Error in chol((value + t(value))/2) :

Hi All I'm analysing growth rates using a gompertz (logistic) curve and am attempting to fit parameters for all of my study birds using nlsList. I've been looking for an option in nlsList to set min and max limits for parameter values during estimation, although I have failed to find any under 'control'. Other users seem to have reported that using control in the port algorithm

I haven't used r in quite a while but would like to get back into it. I have a problem that I would like to solve with r. I have some multinomial data that looks to follow an asymmetric sigmoidal growth pattern. Solving a multinomial gompertz partial differential equation in r is what I’m after. Would anyone be able to provide me the code and packages to do something like this? Regards,

Hi, all. I've written a function that returns the survival function for a Gompertz mortality model. I've specified the two model parameters. Using a simple integration, I can calculate the life expectancy at any age. Is there a way I can use a loop with the integration that will quickly return life expectancy over a range of ages, say 0 to 80, so that I don't have to manually type