similar to: Alternative options: nonlinear model &autocorrelation?

Dear R community, I have a nonlinear model describing average daily soil temperature. What test should I use to compare differences in soil temperature of the two studied vegetation types depending upon month? Building linear contrasts for the developed nonlinear model does not help since this model does not include variable “Months” (only “Days”). 1) Just a Student’s test is not probably

Dear R community, I used a linear mixed model (named lm11) to model daily soil temperature depending upon vegetation cover and air temperature. I have almost 17,000 observations for six years. I can not account for autocorrelation in my model, since I receive the error message after applying the function: update(lm11, corr=corAR1()) Error: cannot allocate vector of size 220979 Kb Do

Dear R community, I used a linear mixed model (named lm11) to model daily soil temperature depending upon vegetation cover and air temperature. I have almost 17,000 observations for six years. I can not account for autocorrelation in my model, since I receive the error message after applying the function: update(lm11, corr=corAR1()) Error: cannot allocate vector of size 220979 Kb Do

Dear R-list members, Cynthia M. Jones wrote a paper (Fitting growth curves to retrospective size-at-age data, Fisheries Research 46(2000):123-129; abstract at http://www.elsevier.nl/gej-ng/10/19/44/70/24/37/abstract.html)where the SAS procedure MIXED, Macro NLINMIX (Littell et. al., 1996)was used to estimate the von Bertalanffy growth function parameters assuming that data from the same fish are

Hi, I have collected data on trees from 5 forest plots located within the same landscape. Data within the plots are spatially autocorrelated (calculated using Moran's I). I would like to do a ANCOVA type of analysis combining these five plots, but the assumption that there is no autocorrelation in the residuals is obviously violated. Does anyone have any ideas how to incorporate these spatial

Dear all, I would like to know whether any specification test for linear against nonlinear model hypothesis has been implemented in R using the quantreg package. I could read papers concerning this issue, but they haven't been implemented at R. As far as I know, we only have two specification tests in this line: anova.rq and Khmaladze.test. The first one test equality and significance of

Hi Version 0.7 of package tsDyn presented at useR! 2009 is now on CRAN, extended with several new features. The package tsDyn is aimed at estimating nonlinear time series models which exhibit regime specific properties. The regime switching dynamics can either be described by smooth transition (STAR and LSTAR) or threshold effects (SETAR). The package furthermore offers nonlinear models

Hi Version 0.7 of package tsDyn presented at useR! 2009 is now on CRAN, extended with several new features. The package tsDyn is aimed at estimating nonlinear time series models which exhibit regime specific properties. The regime switching dynamics can either be described by smooth transition (STAR and LSTAR) or threshold effects (SETAR). The package furthermore offers nonlinear models

Dear all, I am fitting a nonlinear mixed-effects model from a balanced panel of data using nlme. I would like to know whay would be the best options for formally testing for autocorrelation. Is it possible to carry out a Durbin-Watson test on a nlme object? As far as I've seen, I think the durbin.watson function from the car package just works on lm objects. Thank you very much, Antonio

Dear forum members, What is the formula to calculate denominator degrees of freedom (den df) for nonlinear mixed-effect models with covariates? My model is similar to a CO2 uptake example from Pinheiro and Bates (2000, page 376). In this CO2 dataset, there are two treatments and two types (84 observations in total), but den df for each parameter of the model is 64. Isn’t it too high? Your

Hi I implemented the age-structure model in Gove et al (2002) in R, which is a nonlinear statistical model. However running the model in R was very slow. So Dave Fournier suggested to use the AD Model Builder Software package and helped me implement the model there. ADMB was incredibly fast in running the model: While running the model in R took 5-10 minutes, depending on the

Dear all, I am trying to fit a mixed-effects non linear regression, but I have some trouble with it. My data are a balanced panel of 904 subjects with 8 observations (at regular periods) per subject. The functional form of my model is Y=Aexp(-BX1)X2 +e. I want to allow parameters A and B to vary among subjects and also include an autocorrelation term. I have already fitted a standard nonlinear

Hi, I have been brought back to the "R-Side" from MatLab. I have used R in graduate econometrics but only for statistics and regression (linear and nonlinear). But now I need to run general nonlinear optimization. I know about the add-in quadprog but my problem is not QP. My problem is a general nonlinear (obj funct) with linear constraints.I know about the "ms" and

I have the following set-up. 6 values of a continuous variable (let's say light intensity) are presented to a system. The input is presented as a random series of blocks lasting (say) 5 sec each. ---- ---- ---- etc ---- time -> The output is measured and sampled at say 10 samples/sec. Please ignore the fact that this is a time series and

Dear R community, Below you may find the details of my model (lm11). I receive the error message "Error: cannot allocate vector of size 220979 Kb" after applying the autocorrelation function update(lm11, corr=corAR1()). lm11<-lme(Soil.temp ~ Veg*M+Veg*year, data=a, random = list(Site=pdDiag(~Veg), Plot=pdDiag(~Veg))

I am using nonlinear regression to fit a couple of variables to a set of measurements. I would like to do some significance tests for the estimated parameters. I am able to check the confidence intervals using the Jacobian coming out of nonlinear regression. I do see in a paper which shows t-value (it says estimated by White method??), f-value, f-test, and j-test, are these available in matlab,

Dear R-Help, I am trying to use R to solve a nonlinear equation many times for different values. I am looking for a mathematical nonlinear equation solution which may not have a closed solution form. For example, I have equation: 2 = (t^2)/log(t) What is t? I am wondering how to solve it in R. Many thanks, Zhu Wang Statistical Science Department SMU.

Dear all, I'm interested in testing 2 nonlinear restrictions on coefficients of a nls object. Is there a package for doing this? Something in the lines of `test(nls object, res=c("res 1","res 2"),...)' I only found the function delta.method in the alr3 library that calculates the se of a singleton nonlinear restriction of a nls object using the delta method. Thanks in

Hi everyone, I've got a question concerning the function acf(.) in R for calculating the autocorrelation in my data. I have a table with daily returns of several stocks over time and I would like to calculate the autocorrelation for all the series (not only for one time series). How can I do this? After that I want to apply an autoregressive model based on the estimated lag in the

Hi All, Could any one tells me if R or S has the capacity to fit nonlinear regression with Huber's M estimation? Any suggestion is appreciated. I was aware of 'rlm' in MASS library for robust linear regression and 'nls' for nonlinear least squares regression, but did not seem to be able to find robust non-linear regression function. Thanks and regards, Ray Liu