Re: [R] GSL in C code for R

2010-04-30 Thread alexander russell
Hello, I'm following this, since it seems I'm having to try to get C code going in order to estimate dynamic models from biological data.Without saying that gsl and related packages are proving to be obdurate, and knowing that R is proving to be an excellent solution to the widest range of

[R] grofit

2010-04-15 Thread alexander russell
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]] __

Re: [R] estimation of parameters with grofit

2010-03-23 Thread alexander russell
to the approximately 12 that we read will give us a basic regression? regards, shfets On Mon, Mar 22, 2010 at 2:43 PM, alexander russell ssv...@gmail.com wrote: Hello, I'm trying to understand grofit's estimation of models, and fairly new to growth models generally. The data used by grofit consists

[R] estimation of parameters with grofit

2010-03-22 Thread alexander russell
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

[R] standard errors in bbmle

2009-10-23 Thread alexander russell
Hello, Mle2 is a little unforthcoming in the matter of standard errors? Is there a way to ask the program to supply standard errors along with estimates in cases when it doesn't print them 'voluntarily'? regards, s [[alternative HTML version deleted]]

[R] Fwd: When modeling with negbin from the aod package...

2009-10-21 Thread alexander russell
-- Forwarded message -- From: alexander russell ssv...@gmail.com Date: Tue, Oct 20, 2009 at 4:34 PM Subject: Re: [R] When modeling with negbin from the aod package... To: Matthieu Lesnoff matthieu.lesn...@gmail.com Hello again, It seems that, though we have a simple estimate

[R] When modeling with negbin from the aod package...

2009-10-14 Thread alexander russell
Hi, When modeling with negbin from the aod package, parameters for a given count y | lambda~Poisson(lambda) with lambda following a Gamma distribution Gamma(r, theta) are estimated. The intercept is called phi. Some other parameters may be also be estimated from factors in the data: the

[R] interactions and stall or memory shortage

2009-08-31 Thread alexander russell
Hello, After putting together interaction code that worked for a single pair of interactions, when I try to evaluate two pairs of interactions( flowers*gopher, flowers*rockiness) my computer runs out of memory, and the larger desktop I use just doesn't go anywhere after about 20 minutes. Is it

[R] coeftab

2009-08-23 Thread alexander russell
Hello, I'd like to ask: going on from Richard McElreath's coursework in which he says(my paraphrase) basically that the constancy with which a parameter appears satisfactory with respect to its standard error over a series of different models is a good measure of its usefulness, where is the

[R] parameter asterisks

2009-08-04 Thread alexander russell
Hello, Is there a clearcut answer as to why R prints 'NA' sometimes instead of standard errors? mle2(minuslogl = nlikfun4, start = list(a = 1, c = 1, d = 0.2, b = 0.1, b1 = 0.1), method = Nelder-Mead) Coefficients: Estimate Std. Error z value Pr(z) a 3.83845751 0.47320236 8.1117 4.993e-16

[R] is AIC always 100% in evaluating a model?

2009-07-03 Thread alexander russell
Hello, I'd like to say that it's clear when an independent variable can be ruled out generally speaking; on the other hand in R's AIC with bbmle, if one finds a better AIC value for a model without the given independent variable, versus the same model with, can we say that the independent variable

[R] mle2 for poisson data

2009-05-25 Thread alexander russell
Hello, I'd like to obtain parameters for my poisson data. Having made the function poisregfun = function(a,b) {Y.pred = exp(a + b*x) - sum(dpois(y, lambda = Y.pred, log = TRUE))} where Y is a vector of data to be regressed on another vector x. I'm asking for mle2( poisregfun, start = list(a =