similar to: Fitting of Non-Linear Diff Equations and Parameter Estimation

Hello, I am developing an agent-based model to simulate the spread of infectious diseases in heterogeneous landscapes composed of habitat polygons (or clumps of connected cells). To simplify the model, I consider a habitat grid (or raster) containing the polygon ID of each cell. In addition, I have epidemiological parameters associated with each polygon ID. At each time step, the parameter values

Dear list - Does anyone have any ideas / comments about why I am receiving the following warning when I run lsoda: 1: lsoda-- at t (=r1), too much accuracy requested in: lsoda(start, times, model, parms) 2: for precision of machine.. see tolsf (=r2) in: lsoda(start, times, model, parms) I have tried changing both rtol and atol but without success. I saw the thread in the

It seems that the method `iteration` in the function `ode` can be useful in my case: > Method "iteration" is special in that here the function func should > return the new value of the state variables rather than the rate of > change. This can be used for individual based models, for difference > equations, or in those cases where the integration is performed within >

> On Jun 21, 2017, at 12:48 PM, Marine Regis <marine.regis at hotmail.fr> wrote: > > Hello, > > I am developing an agent-based model to simulate the spread of infectious diseases in heterogeneous landscapes composed of habitat polygons (or clumps of connected cells). To simplify the model, I consider a habitat grid (or raster) containing the polygon ID of each cell. In

I apologize that I've asked a similar question before, but being new to R I don't think I did a very good job of formating the question. I've included a text file since the date set is somewhat large. What I have is a huge string of numbers in a text file. The numbers are all separated by comma's and the groups are separated by a semicolon. What I would like to do is read each

Hello R experts, what follows is my reproducible example: mydata<-structure(list(ped.avg = c(335.9, 110.8, 645.7, 638.9, 1468.1, 126.4, 4811.1, 88.5, 868.5, 656.6, 723.6, 654, 2.8, 15, 14.2, 17.5, 15.4, 112.1, 424.7, 18.3, 19.9, 28.6, 25.6, 23.5, 15.4, 27, 62.1, 15.6, 74.6), ped.erst = c(96, 53.2, 615.2, 616.5, 512.9, 56.2, 1851.8, 57.1, 579.5, 613.2, 601.1, 613.6, 1.3, 6.3, 6.5, 6.1,

Hi I have a doubt regarding Mini-OS and Xen. If it is true that Xen reserves in the virtual address space the top 64MB on 32bit systems and 168MB in PAE systems, doesn''t it exist with a Mini-OS kernel which has only 32MB in size i.e. going by its default config file? Obviously I''m missing something. Regards, Bhaskar. _______________________________________________ Xen-devel

Hi, Am I hitting some limit in match? Consider the following example: > tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1) > sum(tst) [1] 76 > seq(100,125,by=.2) [1] 100.0 100.2 100.4 100.6 100.8 101.0 101.2 101.4 101.6 101.8 102.0 102.2 [13] 102.4 102.6 102.8 103.0 103.2 103.4 103.6 103.8 104.0 104.2 104.4 104.6 [25] 104.8 105.0 105.2 105.4 105.6 105.8 106.0 106.2 106.4 106.6 106.8

Hi, Am I hitting some limit in match? Consider the following example: > tst<-seq(100,125,by=.2)%in%seq(0,800,by=.1) > sum(tst) [1] 76 > seq(100,125,by=.2) [1] 100.0 100.2 100.4 100.6 100.8 101.0 101.2 101.4 101.6 101.8 102.0 102.2 [13] 102.4 102.6 102.8 103.0 103.2 103.4 103.6 103.8 104.0 104.2 104.4 104.6 [25] 104.8 105.0 105.2 105.4 105.6 105.8 106.0 106.2 106.4 106.6 106.8

On Xen-3.3.0, domain0 Megaraid SAS (SAS 1068 controller) driver is not loading correctly if vtd support in Xen is enabled. It fails at the point of initializing firmware. I wasn''t seeing this error with Xen-3.2.2-rc3 (Unstable version), though with vtd disabled in Xen-3.3.0, it is working. Looks like a degrade problem. Any clues? Thx, Venkat

Hi thanks for your posts so far. Since what I had assumed was a 32 bit kernel which I compiled didn''t work, I religiously set out on compiling a 64 bit version as my system is currently booted in 64 bit mode. I installed SUNStudio12 on the system and set the SPRO_VROOT to point towards it and when I do a make all in usr/src directory, I get the following error when make cwds and runs in

Hello friends!!! I have a list of values called y. The list is y=c(221.0, 212.0, 206.0, 202.7, 198.4, 195.1, 192.2, 189.7, 187.6, 185.8); y is f(x) and x=c(10,20,30,40,50,60,70,80,90,100). I only have x and y. I don?t know f(x). I would like interpolate f(x) to obtain other values as f(15), f(25), f(35) and if it was possible obtain f(110), f(120). is there any function that allow me obtain

Dear list members, how can multiple imputation realized for anova tables in R? Concretely, how to combine F-values and R^2, R^2_adjusted from multiple imputations in R? Of course, the point estimates can be averaged, but how to get standarderrors for F-values/R^2 etc. in R? For linear models, lm.mids() works well, but according to Rubins rules, standard errors have to be used together with

Hi R-users, I'm not sure what's wrong and how do I check the error? I have 3 sets of data > aug_dt1  [1]  22.8  49.6  77.9  67.6  91.2  48.9 124.6  54.6  85.8  62.8  63.3 174.7  32.3  22.7  99.9 123.2   4.8 149.0  66.0 [20]  73.9  37.4  22.4  69.5  87.3  66.9  87.2  81.5 109.1  47.4  22.8  92.4  60.9  77.5  91.3  91.9  78.2  39.4 107.5 [39]  57.9  65.5  76.0  47.2 145.5  95.9 106.0

I'm sure this is probably simple, but I can't find an answer... I have a data frame (Galton's data on heights of parents and children), in grouped form, parent child frequency 73.5 72.2 1 73.5 73.2 3 72.5 68.2 1 72.5 69.2 2 72.5 70.2 1 72.5 71.2 2 72.5 72.2 7 72.5 73.2 2 72.5 74.2 4 71.5 65.2 1 71.5 66.2 3 71.5 67.2 4 71.5 68.2 3 71.5 69.2 5 71.5 70.2 10 ... and need the ungrouped

Hello all, I have not been able to find an answer to this problem. I feel like it might be so simple though that it might not get a response. Suppose I have a dataframe like the one I have copied below (minus the 'calib' column). I wish to create a column like calib where I am subtracting the 'Count' when 'stain' is 'none' from all other 'Count' data for

hello, > subset(swiss, Agriculture > 60 & Examination != c(14,16), select = c(Agriculture,Examination,Catholic)) Agriculture Examination Catholic Broye 70.2 16 3.30 Glane 67.8 14 4.20 Aigle 62.0 21 5.16 Avenches 60.7 19 5.23 Cossonay 69.3 22 5.62

hi Michael, am sorry for the incomplete reply . csv file data havinmg like this , >getSymbols("IBM") >weekly_data = to.weekly(week_name) >dataframe=data.frame(weekly_data) >outputfile_name="F:\\R-programs\\Outputfile.csv" >write.table(dataframe, file =outputfile_name,sep = ",",col.names =TRUE,row.names = T) > datafrom_table <-

Hello, I'm new to R, and I'm trying to set the number of tickmarks for a plot using xaxp, but R seems to ignore this setting completely. Maybe you've an idea what I'm doing wrong. I'm using R 1.4.1, and here's hat I'm doing: > d[,'seq'] [1] 28913 16323 13922 6237 4257 3881 4100 3781 2694 2064 1769 1550 [13] 1539 480 >

Hi, I am trying to use the Tps procedure to visualize 3 variables in XYZ space, but keep getting the error messages below: Warning messages: 1: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV, Krig.fgcv, : GCV search gives a minimum at the endpoints of the grid search 2: In Krig.find.gcvmin(info, lambda.grid, gcv.grid$GCV.one, Krig.fgcv.one, : GCV search gives a minimum at the endpoints