similar to: Grofit

Hello, I'm trying to understand grofit's estimation of models, and fairly new to growth models generally. The data used by grofit consists of the vector of "experiments", that is the growth values for a vector of individuals measured at different times. Can I understand correctly that the program estimates parameters for the growth model based on a regression(linear or non

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

Does anyone know if it is possible to get the length of the log phase (the main growth phase) from a grofit logistic model? I can get the lag phase length, maximum growth and max rate of growth, but I also want to know the length of the log phase (ie, for how long do the bacteria grow exponentially?) Thank you for any help in advance [[alternative HTML version deleted]]

Dear all, I am analysing growth data - response variable - using GAM and GAMM models, and 4 covariates: mean size, mean capture year, growth interval, having tumors vs. not The models work fine, and fit the data well, however when I try to compare models using AIC I cannot get an AIC value. This is the code for the gam model:

Hello. I need to know how can R compute AIC when I study a regression model? For example, if I use these data: growth tannin 1 12 0 2 10 1 3 8 2 4 11 3 5 6 4 6 7 5 7 2 6 8 3 7 9 3 8 and I do model <- lm (growth ~ tannin) AIC(model) R responses: 38.75990 I know the following formula to compute AIC: AIC=

Hello, I would like to ask about the statistic used for initial values for models built with grofit. Is the mean (of the experiments or cases) at t1 used? Also at the other time points? regards, Russell [[alternative HTML version deleted]]

Dear all I have undertaken some phylogenetic and non-phylogenetic regressions with gls() in nlme with single preictor variables. A p value is associated with the intercept (upper p value) and another with the predictor variable (lower). Which p value is important? What does it mean if the intercept p value is insignificant but the predictor is still significant? Thanks a lot, and sorry for my

Hi there, I'm a PhD student investigating growth patterns in fish. I've been using the minpack.lm package to fit extended von Bertalanffy growth models that include explanatory covariates (temperature and density). I found the nls.lm comand a powerful tool to fit models with a lot of parameters. However, in order to select the best model over the possible candidates (without covariates,

To compare two different models, I've compared the result of using dropterm() on both. Single term deletions Model: growth ~ days + I(days^0.5) Df Sum of Sq RSS AIC <none> 2.8750 -0.2290 days 1 4.8594 7.7344 4.6984 I(days^0.5) 1 0.0234 2.8984 -2.1722 AND Single term deletions Model: growth ~ days + I(days^2) Df Sum

Hi, I'm a very new user of R and I hope not to be too "basic" (I tried to find the answer to my questions by other ways but I was not able to). I have 12 response variables (species growth rates) and two environmental factors that I want to test to find out a possible relation. The sample size is quite small: (7<n<12, depending on each species-case). I performed a

Hi I get the following error when I try and get the AICc for a gls regression using qpcR: > AICc(gls1) Loading required package: nlme Error in n/(n - p - 1) : 'n' is missing My gls is like this: > gls1 Generalized least squares fit by REML Model: thercarnmax ~ therherbmax Data: NULL Log-restricted-likelihood: 2.328125 Coefficients: (Intercept) therherbmax 1.6441405

Hello everybody. In my case, calculating multiple comparisons (Tukey) after lmer produced the following two errors: > sv.mc <- glht(model.sv,linfct=mcp(comp="Tukey")) Error in x$terms : $ operator not defined for this S4 class Error in factor_contrasts(model) : no 'model.matrix' method for 'model' found! What I have done before: > sv.growth <-

Hi,I have tumor growth curve data for a bunch of different mice in various groups. I want to compare the growth curves of the different groups to see if timing of drug delivery changed tumor growth.I am trying to run a mixed models with repeated measures over time with each mouse as a random effect with linear and quadratic terms for time.This took me a long time to figure out and I just wanted to

Hi I'm trying to replicate Smith et al.'s (http://www.sciencemag.org/content/330/6008/1216.abstract) findings by fitting their Gompertz and logistic models to their data (given in their supplement). I'm doing this as I want to then apply the equations to my own data. Try as a might, I can't quite replicate them. Any thoughts why are much appreciated. I've tried contacting the

Hi there, I am looking to compare nonlinear mixed effects models that have different nonlinear functions (different types of growth curve)embedded. Most of the literature I can find focuses on comparing nested models with likelihood ratios and AIC. Is there a way to compare model fits when models are not nested, i.e. when the nonlinear functions are not the same? Many thanks in advance! Lindsay

Dear all, I'm modeling growth curve of some ecosystems with respect to their rainfall-productivity relationship using a simple linear regression (ANPP(t)=a+b*Rain(t)) and a modified version of the Brody Model ANPP(t)=a*(1-exp(-b*rain(t))) I would like to know why the "best model" is function of the criteria that I use (maximizing the fit using R2 or testing the Null hypothesis with

Dear R users, I am having some problems with glm. The first is an error message "subscript out of bounds". The second is the fact that reasonable starting values are not accepted by the function. To be more specific, here is an example: > success <- c(13,12,11,14,14,11,13,11,12) > failure <- c(0,0,0,0,0,0,0,2,2) > predictor <- c(0,80*5^(0:7)) >

Dear R-list members, Some months ago I wrote a message on the usage of gnls (nlme library) and here I come again. Let me give an example: I have a 10 year length-at-age data set of 10 fishes (see growth.dat at the end of this message) and I want to fit a von Bertalanffy growth model, Li= Linf*(1-exp(-k*(ti-t0))) where Li = length at age i, Linf= asymptotic length, k= curvature parameter, ti=

Dear R users, Could somebody please help me to find a way of comparing nonlinear, non-nested models in R, where the number of parameters is not necessarily different? Here is a sample (growth rates, y, as a function of internal substrate concentration, x): x <- c(0.52, 1.21, 1.45, 1.64, 1.89, 2.14, 2.47, 3.20, 4.47, 5.31, 6.48) y <- c(0.00, 0.35, 0.41, 0.49, 0.58, 0.61, 0.71, 0.83, 0.98,

Dear help community, I'm a R-beginner and use it for my master thesis. I've got a mixed model and want to analyse it with lme. There are a lot Cofactors that coult be relevant. To extract the important ones I want to do the stepAIC, but always get an error warning. Structure of my data: data.frame': 72 obs. of 54 variables: $ Block : Factor w/ 3 levels