similar to: horizontal grouped stacked plots and removing space between bars

Hi, I have a data frame like this: > dput(df) structure(list(ID = 1:8, Type = c("gmx mdrun -ntmpi 8 -ntomp 1 -s benchPEP.tpr -nsteps 10000 -resethway", "gmx mdrun -ntmpi 8 -ntomp 1 -s benchPEP.tpr -nsteps 10000 -resethway", "gmx mdrun -ntmpi 8 -s benchPEP.tpr -nsteps 4000 -resetstep 3000", "gmx mdrun -ntmpi 8 -s benchPEP.tpr -nsteps 4000 -resetstep

Hi; I have a R script that includes a call to genoud(); genoud process lasts about 4 seconds, what would be OK if I hadn't have to call it about 2000 times. This yields about 2 hours of processing. And I would like to use this script operationally; so that it should be run twice a day. It seems to me that the parallel processing option included in genoud() divides the task inside the function

Hi, I am trying to simulate animal movement in a gridded landscape made up of cells. At each time step (iteration), the animal moves from one cell to another in a random fashion. This is how I am simulating movement, where a and b are the x,y co-ordinates of the animal at the previous time step: for (i in 1:no.of.steps){ direction <- sample(1:8, 1) if(direction == 1){ a <- a b <- b -

Thanks! This is in regards to bug 3084, I am trying to find out why samba is crashing. I am using glibc-2.4 with nptl threads. It seems to be a memory or threading issue. Thanks for any help. Running smbd through gdb: Breakpoint 6, charset_name (ch=CH_DISPLAY) at lib/charcnv.c:57 57 const char *ret = NULL; (gdb) Continuing. Substituting charset 'ANSI_X3.4-1968' for

Dear wonderful R community, I've been creating color-coded concentric circles using the "points" function, but I just realized that what I would really like to do is draw color-coded concentric half-circles. (Because I want to communicate information about the diameters-- half-circles are sufficient to show diameters AND they leave ample room to one side of the figure to label

attempting to do multivariate modelling in R with known future conditions (in this case variable 'b') using MSBVAR and hc.forecast. The sample code (a paired down representation) does not give anywhere near the expected results - I am assuming that a forecast 8 steps out would approximate 'a' as the sequence 1.1,2.1,3.1,100.1 corresponding to the input set. I have varied the input

I saw the other topic posted 13 Sep 2010, about the fixme ... unimplemented error. I get that too. However, my system also crashes. I'm running Ubuntu 10.4 with Wine 1.3.1 on a AMD 1055t cpu (no overclocking) with a MSI 785G-e53 motherboard. The system is stable otherwise. This is executing the folding at home client, the Windows client beta 6.30. This is used by enough people and it

I am trying to use optim() to minimize a sum-of-squared deviations function based upon four parameters. The basic function is defined as ... SPsse <- function(par,B,CPE,SSE.only=TRUE) { n <- length(B) # get number of years of data B0 <- par["B0"] # isolate B0 parameter K <- par["K"]

GENERAL REFERENCE ON NONLINEAR OPTIMIZATION? What are your favorite references on nonlinear optimization? I like Bates and Watts (1988) Nonlinear Regression Analysis and Its Applications (Wiley), especially for its key insights regarding parameter effects vs. intrinsic curvature. Before I spent time and money on several of the refences cited on the help pages for "optim",

Hello list-mailers, I'm trying to use the 'Tps' function from the 'fields' package and I can getting this error message: > fishplot <- Tps(nodes.catch, res2) Error in if ((il > 1) & (il < nstep.cv)) { : missing value where TRUE/FALSE needed In addition: Warning message: no non-missing arguments to min; returning Inf nodes.catch is a 1600x2 matrix and

Saludos Cordiales Estoy tratando de obtener la descomposición de varianza para un modelo de vectores autorregresivos estructural pero el módulo vars no permite hacerlo de manera directa, sin embargo, encontré programaciones de funciones específicas para los vectores autorregresivos estructurales que se acoplan directamente al modelo VARs, alguien me puede decir en que parte debo incorporar la

Dear peRsons, I have a Sweave document which demonstrates the usage of step() function. With current R version 2.5.0 the step() function was changed so that the heading of trace=TRUE output for each model is printed using command message(): if (trace) message("\nStep: AIC=", format(round(bAIC, 2)), "\n",

Dear all, I have tried to modify the code of MASS:::stepAIC.default(), dropterm() and addterm() to use AICc instead of AIC for model selection. The code is appended below. Somehow the calculations are still not correct and I would be grateful if anyone could have a look at what might be wrong with this code... Here is a working example: ## require(nlme) model1=lme(distance ~ age + Sex, data =

Dear friends, are there enough R users in Cyprus to form a club? jason Dr. Iasonas Lamprianou Department of Education The University of Manchester Oxford Road, Manchester M13 9PL, UK Tel. 0044 161 275 3485 iasonas.lamprianou at manchester.ac.uk ----- Original Message ---- From: "r-help-request at r-project.org" <r-help-request at r-project.org> To: r-help at r-project.org