Displaying 20 results from an estimated 300 matches similar to: "DeSolve Package and Moving Average"
2007 Feb 20
2
scaler plugin fixes
Hi
I think this is mostly a question of preference for my part, but I
think it would give an overall better quality feel. What I'm talking
about is mostly the movement it uses. As I've come to understand, it
uses velocity and direction, and a target point. This means that it
updates it velocity and direction every X seconds (timestep), and
gives it the possibility to be
2012 Feb 03
3
[fields] image.plot abends with NAs in image.plot.info
summary: image.plot-ing two sets of netCDF data, with the second
derived from the first. First plots to PDF as expected (title, data,
legend). Second plots the data and title, but abends before drawing
the legend, with
> Error in if (del == 0 && to == 0) return(to) :
> missing value where TRUE/FALSE needed
> Calls: plot.layers.for.timestep -> image.plot -> seq ->
2009 Jun 12
2
External signal in ODE written in C (using deSolve and approx1?)
Dear list
The deSolve package allows you to specify the model code in C or Fortran.
Thanks to the excellent vignette this works fine. However I have not yet
managed to use forcing functions in C code.
In pure R code this works very well with approxfun() specified outside the
model:
###############################################
#Model
lvml <- function(t, x, parms) {
2000 Aug 29
1
Newbie question: Linear regression with error bars.
Hello guys,
I am a total newbie on R, having downloaded it, read the documentation and
started playing with it right now.
My general question is what 'lr' model can be used for doing a linear
regression on points that have a variance associated with them (ie. Monte
Carlo simulation results).
Actually my Data sets look like:
Timestep Energy Variance_of_the_Energy
0.0005 -14.876840
2007 Jun 24
1
movement backend
Hi
I'm just wondering what the status is on the new movement
architecture that I think David mentioned a few months ago. The point
was to get more predictable and accurate movement in the plugins.
Since the current method recomputes speed and direction every
"timestep". So the movemnt is actually computed on the fly, instead
of having a predetermined path which is
2010 Mar 12
1
Exporting animation of a series of plots
Hi All,
Thank you in advance for your help.
I made an agent-based model in R and I plot it using image() each
timestep creating a nice animation. I can't seem to find in any of the
R codes a means to export the series of images so that they can be
played as a movie or animation.
Thank you again.
Flana
2008 May 28
1
Writing intermediate analysis to disk
Is there a way to write and analysis to disk and then reconstruct the whole
thing back into an object.
wavCWT() #wmtsa package
I am running out of memory on my computer and I was wondering if there was a
way to iterate through this process (as it is an iterative process anyway-
it just stores the whole thing to memory). Or is there a way to set the
scale that I want to look at so that wavCWT
2011 Oct 17
1
compressing/reducing data for plot
Hello,
I have simulation results in the form of
Time V I
0.000000000000e+000 7.218354344368e-001 5.224478627497e-006
1.000000000000e-009 7.218354344368e-001 5.224477718002e-006
2.000000000000e-009 7.218354344368e-001 5.224477718002e-006
4.000108361244e-009 7.218354344368e-001 5.224478627497e-006
8.000325083733e-009
2010 Jun 07
1
Desolve package: How to pass thousand of parameters to C compiled code?
Hi,
I have used DeSolve package for my ODE problem regarding infectious disease transmission and currently am trying to pass lots (roughly a thousand) of model parameters to the C compiled model (I have to use C compiled code instead of R code purely because of the speed).
I can't go define it one by one as it gonna take ages to finish and also quite difficult to revise. I have read the
2011 Aug 12
1
deSolve output
Hi,
I've solved a simple differential equation describing the degradation of
amino acid carbon (THAA-C) using deSolve.
Code is a follows:
# Input of model parameters, a and b describes form of curve, i is apparent
initial age of Org. C.
parameters <- c(a = a, b = b, i=i)
# Initial value of the model, G
state = c(G = G)
#specifies the function degradation as a function of
2011 Jan 25
1
deSolve: Problem solving ODE including modulo-operator
I have a problem integrating the 'standard map' (
http://en.wikipedia.org/wiki/Standard_map
http://en.wikipedia.org/wiki/Standard_map ) with deSolve:
By using the modulo-operator '%%' with 2*pi in the ODEs (standardmap1), the
resulting values of P and Theta, should not be greater than 2pi. Because
this was not the case, i was thinking that the function 'ode' has some
2012 Feb 03
1
Logistic population growth and deSolve
Hello, I am new to R and I am having problems trying to model logistic
population growth with the deSolve package. I would like to run the model
for four populations with the same initial population and carrying capacity
but with different growth rates and put the results into a data frame. When
I run the following lines of code I get unexpected results from "output" but
the format is
2013 Feb 21
1
using and event in deSolve
Hi All
Having been pointed the use of events and roots in deSolve, I was able to
implement the Izchikevich model of spiking neurons. However, I'm not too
sure of defining the event. The deSolve documentation says:
An event is triggered when the ball hits the ground (height = 0) Then
velocity (y2) is reversed
and reduced by 10 percent. The root function, y[1] = 0, triggers the event:
>
2011 Oct 03
0
deSolve - Function daspk on DAE system - Error (Vince)
Vince,
When that happens, one possible reason is that your DAE is of index > 1,
which cannot be solved by daspk.
The solver radau, also from deSolve can handle DAEs up to index 3, but
you need to rewrite the problem in the form M*y' = f(x,y), where M is a
mass matrix.
If you do that for your problem, and solve it with radau, then radau
complains that the "matrix is repeatedly
2009 Jun 29
0
new version of package deSolve on CRAN
Dear R users,
an improved version of package deSolve (version 1.3) is now available on
CRAN. deSolve, the successor of R package odesolve, is a package to
solve initial value problems (IVP) of:
- ordinary differential equations (ODE),
- differential algebraic equations (DAE) and
- partial differential equations (PDE).
The implementation includes stiff integration routines based on the
ODEPACK
2009 Jun 29
0
new version of package deSolve on CRAN
Dear R users,
an improved version of package deSolve (version 1.3) is now available on
CRAN. deSolve, the successor of R package odesolve, is a package to
solve initial value problems (IVP) of:
- ordinary differential equations (ODE),
- differential algebraic equations (DAE) and
- partial differential equations (PDE).
The implementation includes stiff integration routines based on the
ODEPACK
2018 May 31
0
Help in dynamic simulation using deSolve
Hi R help,
I wanted to simulate two pool model (A&B) using deSolve package for time 0 to 12 by 1.? Initial values of the state variables are A=5, B=3. The fluxes are as follows1) Flux into A= 5 units per unit time?2) Flux from A to B= 0.33) Flux out of A=0.1?4) Flux from B to A=0.35) Flux out of B=0.3
Here is the R code I compiled to estimate the size of A and B and graph the output
2011 Oct 02
0
deSolve - Function daspk on DAE system - Error
I'm getting this error on the attached code and breaking my head but can't
figure it out. Any help is much appreciated. Thanks, Vince
CODE:
library(deSolve)
Res_DAE=function(t, y, dy, pars) {
with(as.list(c(y, dy, pars)), {
res1 = -dS -dES-k2*ES
res2 = -dP + k2*ES
eq1 = Eo-E -ES
eq2 = So-S -ES -P
return(list(c(res1, res2, eq1, eq2)))
})
}
pars <- c(Eo=0.02,
2011 Aug 16
1
deSolve, extracting variable values from inside ode function
I'm just getting to grips with using ode function and have used the
examples and vignettes to produce a small model of a one-pool,
michaelis-menten, enzyme kinetic reaction. The rate of flux of substrate
into pool A is constant (fluxoa) however the rate of flux out of pool A
is controlled by the HMM equation (v = Vmax/ ( 1 + (Km / Concentration A
)) ). This function works fine and
2013 Apr 21
1
lsoda question from deSolve package
Dear List,
Wonder if you have some thoughts on the following question using lsoda in desolve:
I have the following data and function:
require(deSolve)
times <- c(0:24)
tin <- 0.5
D <- 400
V <- 26.3
k <-0.056
k12 <- 0.197118
k21 <- 0.022665
yini <- c(dy1 = 0,dy2 = 0)
events <- data.frame(var = "dy1",time = c(10,15),value = c(200,100),method =