similar to: how to start R script editor by default

Hello! I am performing coupling of chains in MCMC and I need the same value of seed for two chains. I will show demo of what I want: R code, which might show my example is: niter <- 3 nchain <- 2 tmpSeed <- 123 for (i in 1:niter) { # iterations for (j in 1:nchain) { # chains set.seed(tmpSeed) a <- runif(1) cat("iter:", i, "chain:", j,

Devang Patel wrote: > On Jul 25, 2008, at 12:58 PM, Marc de Kruijf wrote: > > I'd like to write a pass that does both: implements a code transformation and saves information that can be accessed by subsequent passes. > > Ideally, we want to use two separate pass. However, it is quiet possible that your requirement is unique. Would it be possible to provide more info. on what

Thanks for the fast response John. On Thu, Oct 1, 2015, at 04:51 PM, John Criswell wrote: > Dear Ed, > > First, someone has updated the DSA code in the poolalloc project to LLVM > 3.7, and a Master's student worked for me over the summer to update a > large chunk of SAFECode to LLVM 3.7. However, the update to LLVM 3.7 > isn't finished (we need to finish integrating

Dear R users, I'm pretty new on using lmer package. My response is binary and I have fixed treatment effect (2 treatments) and random center effect (7 centers). I want to test the effect of treatment by fitting 2 models: Model 1: center effect (random) only Model 2: trt (fixed) + center (random) + trt*center interaction. Then, I want to compare these 2 models with Likelihood Ratio Test.

I don't know the details of pls (in the pls.pcr package, I assume), but if you use validation="CV", that says you want to use CV to select the best number of components. Then why would you specify ncomp as well? Andy > From: ryszard.czerminski at pharma.novartis.com > > When I try to use ncomp parameter in pls procedure I get > following error: > > >

Hi all, I am currently doing a consistency analysis using an lmer model and trying to use parametric bootstrapping for the confidence intervals. My model is like this: model<-lmer(y~A+B+(1|C/D)+(1|E),binomial) where E is the individual level for consistency analysis, A-D are other fixed and random effects that I have to control for. Following Nakagawa and Scheilzeth I can work out the

Hello, I'm working on a problem using nested for-loops and I don't know if it's a problem with the order of the loops or something within the loop so any help with the problem would be appreciated. To briefly set up the problem. I have 259 trees (from 11 different species, of unequal count for each species) of which I am trying to predict biomass. For each tree species I have 10000

Hello, I would appreciate *any* ideas on this problem. I'm the maintainer of a package ("subselect"), which on CRAN's Daily Package Checks is OK on all flavours of R, except r-devel Linux x86_64, where there is a "memory not mapped" segfault with the very first example that is tried out (output below). Additionally, a user with an AMD64 machine has just reported a

Hi, I'm trying to get the pool allocator and safe code building against llvm trunk. I've run into a build error, and I see that in the past another user was told just not to build the pool allocator for use with safecode [1]. However, I really want the pool allocator transforms, so I just wanted to check why the suggestion was not to use it. Has it been superseded in some way by something

The relaying-configuration seems to be straight-forward, unfortunately it does not work. That's my configuration: <relays-on-demand>1</relays-on-demand> <relay> <server>http://mp3stream7.apasf.apa.at</server> <port>8000</port> <mount>/</mount> <local-mount>/oe3.mp3</local-mount> <on-demand>1</on-demand>

Dear R-help list readers, I am fitting mixed models with the lme function of the nlme package. If I get convergence depends on how the method (ML/REM) and which (and how much) parameters will depend randomly on the cluster-variable. How get the bist fit without convergence? I set the parameters msVerbose and returnObject to TRUE: lmeControl(maxIter=50000, msMaxIter=200, tolerance=1e-4,

Hello, My problem is, that I have to build hundreds of GARCH models to obtain volatility forecasts. I would like to run a loop, that would build those forecasts for me. There is no problem, with writing only the results of the forecasts, but I'd like to have stored results of the models in some file, that I could check later, what are the models like, to be able to compare if I should use

I would like to build a model in R to simulate the seed dispersal by one plant. The plant produced 5 seeds and the probability of falling inside the eight closest space was 0.8 and in the next space 0.2 and in the rest space 0: 0 0 0 0 0 0 0.2 0.2 0.2 0.2 0.2 0 0.2 0.8 0.8 0.8 0.2 0 0.2 0.8 1 0.8 0.2 0 0.2 0.8 0.8 0.8 0.2 0 0.2 0.2 0.2 0.2 0.2 0 0 0 0 0 0 0

## I would like help in using variable values in plotmath expressions ## in lattice x <- 1:10 y <- 1:10 pval <- .95 plot(y ~ x, ## works as I want in base graphics main=substitute(list(alpha * " = " * group("",list(p),"")), list(p=pval))) plot(y ~ x, ## doesn't work as intended: "pval" is displayed main=substitute(list(alpha *

On Jul 25, 2008, at 12:58 PM, Marc de Kruijf wrote: > I'd like to write a pass that does both: implements a code > transformation and saves information that can be accessed by > subsequent passes. Ideally, we want to use two separate pass. However, it is quiet possible that your requirement is unique. Would it be possible to provide more info. on what your pass does ? >

Hi , This Is 'Suresh M' from AMD compiler team. 1) I obtained the source code of llvm and pool-alloc (relase_32 versions ) using the following commands, svn co http://llvm.org/svn/llvm-project/llvm/branches/release_32 llvm cd llvm/projects svn co http://llvm.org/svn/llvm-project/poolalloc/branches/release_32 poolalloc 2) Was able to build the binaries and libraries

moments::anscombe.test(x) does give errors when x has too few values or if all the values in x are the same > moments::anscombe.test(c(1,2,6)) Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed > moments::anscombe.test(c(2,2,2,2,2,2,2,2)) Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed You can use tryCatch() to

Hi Bill, Thanks for prompt reply and letting me know a way around it. I have more than 1200 observations and not all the values are the same. However, my data points are quite similar, for example, 0.079275, 0.078867, 0.070716 in millions and etc. I have run the data without converting it to millions and I still get the same error message. As I have kurtosis value, it should be fine for the

Hi all, While running the anscombe.test in R, I'm getting an error of *Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed* for a few time series columns whereas for most of the series the function is working fine. I have checked for those specific columns for missing values. However, there is no NA/NAN value in the dataset. I have also run kurtosis for

Hi everybody, I would like to apply cor.test to a matrix with m rows and n columns and get the results in a list of matrices , one matrix for p.val, one for the statistic, one for the correlation and 2 for upper and lower confidence intervals, something analog with cor() applied to a matrix. I have done my own function to get a matrix of p.values and i suppose i can build similar functions for