similar to: A naive lsoda question....

hi I'm puzzled as to the relation between the REML score computed by gam and the formula (4) on p.4 here: http://opus.bath.ac.uk/22707/1/Wood_JRSSB_2011_73_1_3.pdf I'm ok with this for poisson, or for quasipoisson when phi=1. However, when phi differs from 1, I'm stuck. #simulate some data library(mgcv) set.seed(1) x1<-runif(500) x2<-rnorm(500)

Dear all, I am trying to create a 6-plot layout - 3 rows and 2 columns - so that only the top two plots have variable widths, all else with their default setting. Using "layout(matrix(c(1,2,3,4,5,6),3,2,byrow=T),widths=c(5,2))" rescales column #2 of all three rows, whereas I would like to rescale that of only row #1. Is there a simple way of doing this? Thanks, Manojit

https://bugs.freedesktop.org/show_bug.cgi?id=78570 Priority: medium Bug ID: 78570 Assignee: nouveau at lists.freedesktop.org Summary: GT215: Freeze on X startup, 2014. QA Contact: xorg-team at lists.x.org Severity: blocker Classification: Unclassified OS: Linux (All) Reporter: phix.nay at gmail.com

All, If the post "Samba 3.6 -> 4.1 change that prevent sharing / - mount error(13)" makes it to the list, I apologize, but it has been over a day and I have no indication that the post ever made it to the list. The issue is that in versions prior to 4.x, I was able to create a root share in samba, but now, something prevents it from being mounted. The relevant smb.conf is:

I am running R 1.2.3 with ESS5.1.18 with Windows 98. I am trying to use lsoda in the odesolve apckage and am having problems. Question: The return value of the function of the system of ode's has to be a list that includes first, the ode's and second, "a vector (possibly with a `names' attribute) of global values that are required at each point in `times'." I

Hello there, Suppose you want to solve the following system of ODE's (a simple Lotka-Volterra predator prey model) dP/dt = beta*P*V - mu*P dV/dt = r*V - beta*P*V where P and V are the numbers of predators and prey. Now, this is easy to do, but suppose you have a system of equations like this, dP1/dt = beta1*P1*V1 - mu1*P1 dP2/dt = beta2*P2*V2 - mu2*P2 dV1/dt = r1*V1 - beta1*P1*V1

Does anyone know of example r-code/packages for carrying out analysis? Preferably this would have examples from real experimental data of two or more competing species... Thanks Chris Buddenhagen [[alternative HTML version deleted]]

Hi Greg, For quasi families I've used extended quasi-likelihood (see Mccullagh and Nelder, Generalized Linear Models 2nd ed, section 9.6) in place of the likelihood/quasi-likelihood in the expression for the (RE)ML score. I hadn't realised that this was possible before the paper was published. best, Simon ps. sorry for slow reply, the original message slipped through my filter for

Hello, I have a syntactic problem with ode. How do I specify vectors of equations in ordinary differential equation systems. (i.e. in my case I want to simulate an a priory undefined number of species that have different parameters but the same behaviour) I demonstrate this using the Lotka Volterra example. The code below does not work and I have not a good idea how to specify this right. ##

Dino, Alexey, there are a number of users suffering issues with your Powercom devices. Could you (Dino, and Keven if possible) please have a look at the below one, from Angela, and check for a fix? I've scheduled to release 2.6.1 next week, and having that fixed is part of the list. 2011/5/16 Angela Williams <angierfw at gmail.com> > Hi All > > On Friday 13 May 2011 at

Hello, I am working on a paper of ecological modelling in R. I have made a Lotka Volterra model for tree animal species. I was trying to get the model in R but there are some errors. Since i am not so experienced with R i can't find the bugs in my script. I hope if there is someone who can help me handle this problem. Greetz Roel This is the script i wrote: library(odesolve) LotVmod

dear all. Would someone be kind and willing to explain the code below for a person who has never used R? ( that is if one has enough time and inclination) It implements gillepsie's stochastic algorithm for Lotka Volterra model. What would help me tremendously is to see the breakdown of the line by line code into plain english. thanks for any insights or other comments. sean

I'm fairly new to R, and I'm trying to write out a population model that satisfies the following; the system consists of s species, i= 1, 2,...,s network of interactions between species is specified by a (s x s) real matrix, C[i,j] x[i] being the relative population of the "ith" species (0 =< x[i] =< 1, sum(x[i]=1) the evolution rule being considered is as follows;

Hi All! I only ever have complex setups. Customer site has a dedicated leased line from their ISP terminating on a Cisco router. Router is configuered with the first usable address on a /28 network - 196.x.y.73. The linux firewall is configured with the remaining 5 ip''s, 196.x.y.74 to 196.x.y.78 and 79 as the broadcast. Sounds normal but here is the twist. The primary or first ip

Hi R used with C-code experts, I had a look at the archives and did not find anything on this, so hopefully I am not doubling up. I have previously used the following approach where I needed some very small/large numbers (using Brobdingnag): surfacewithdiff <- function(t, y, p) { const=p["const"] kay =p["kay"] psii=p["psii"]

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 =

I have been trying to make use of the odesolve library on my university's Linux grid - currently R version 2.0.1 is installed and the system runs 64-bit Scientific Linux based on Redhat. I cannot seem to get lsoda working when I define the model as a shared C library. For example, the following snippet uses the mymod.c example bundled with the package: ### START rm(list=ls())

Hi All, I am running R 2.2.0 on Mac OS 10.4.2, dual G5 processors with 8 Gig RAM. I am running a simulation with lsoda that requires ~378 s to complete one set of time intervals. I need to optimize the parameters, and so need to considerably speed up the simulation. I have tried to figure out how to change the appropriate memory allocation and have search R help and Introductory

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

Dear Hank, Last question first: really, only you can say for sure if 4e-281 and 5e-11 are small enough; it depends on the units you measure your state variables in. However, this strategy cannot get the state variables to exactly 0. Obviously, you could get closer to 0.0 faster by setting the derivatives even larger in absolute value. You may run into problems with the solver when the