similar to: FW: Random seed problem in MCMC coupling of chains

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,

Folks, I am finding problems with using "dse": > library(dse1) Loading required package: tframe Error: c("package '%s' required by '%s' could not be found", "setRNG", "dse1") > library(dse2) Loading required package: setRNG Error: package 'setRNG' could not be loaded In addition: Warning message: there is no package called

I would like to request help with the installation of dse in raqua in mac os x 10.3. I get the following error message after the messages indicating that parts were successfully installed. I would be most grateful for a solution. ----------------------------------------- * Installing *source* package 'setRNG' ... ** R ** inst ** help >>> Building/Updating help pages for

https://bugzilla.netfilter.org/show_bug.cgi?id=855 Summary: Can create user chains with names unparsable by iptables-restore Product: iptables Version: 1.4.x Platform: All OS/Version: Debian GNU/Linux Status: NEW Severity: minor Priority: P5 Component: iptables AssignedTo:

In Solaris with R-1.7.0beta from 2003-04-01 I still have a problem with "R CMD check" after I build a package: * creating setRNG-manual.tex ... OK * checking setRNG-manual.tex ... ERROR Could not create DVI version. This typically indicates Rd problems. I can manually gererate the dvi from setRNG-manual.tex, so I don't think my latex is defective as I speculated last week. The file

Hello! I am crossposting this to R-help and BioC, since it is relevant to both groups. I wrote a wrapper for Entrez search utility (link for this is provided bellow), which can add some new search functionality to existing code in Bioconductor's package 'annotate'*. http://eutils.ncbi.nlm.nih.gov/entrez/query/static/esearch_help.html Entrez search utuility returns a XML document

I unintentionally missed to cc to r-devel. -------- Original Message -------- Subject: Re: [Rd] levels for list and data.frame Date: Tue, 21 Mar 2006 20:50:21 +0100 From: Gregor Gorjanc <gregor.gorjanc at gmail.com> Reply-To: gregor.gorjanc at gmail.com To: Martin Maechler <maechler at stat.math.ethz.ch> References: <441F2C49.8080703 at gmail.com> <17439.50351.172900.987055

Hello to Rusers! I am puzzled with R and I really do not know where to look in for my problem. I am moving from SAS and I have difficulties in translating SAS to R world. I hope I will get some hints or pointers so I can study from there on. I would like to do something like this. In SAS I can write a macro as example bellow, which is afcourse a silly one but shows what I don't know how

James, thanks for the response. I understand now my puzzle. tmp <- data.frame(y1=1:4, f1=factor(c("A", "B", "C", "D"))) tmp$y2 <- NA tmp[1:2, "y2"] <- 2 Does the job. I see that I should add a full column. In my case adding full column of NAs and then adding values, solves my problems. Thanks to all. james.holtman at convergys.com

Duncan, you are a king! Thanks a lot for this cookie. It really helped me. Thanks for the code as well as detailed explanation at the end. >Hi Gregor. > >Here is a function that will collect all of the nodes in the >XML document whose names are in the vector elementNames > >getElements = >function(elementNames) >{ > els = list() > > startElement = function(node,

I am having problems with the package dse. I just installed R 2.2.1 and reinstalled all packages. I am running Windows XP Pro with all updates. Below there are two examples of error messages generated when trying to execute some simple programs. The code was taken directly from the package documentation. Any help on this will be greatly appreciated. Merry Christmas Fernando

Dear R developers, I needed to align the title of legend for some of my plots. If there is interest to include this into base R, bellow are the examples and the relevent diffs. Regards, Gregor x <- 1:10 y1 <- 1.5 * x y2 <- 2.0 * x y <- cbind(y1, y2) source(file=url("http://gregor.gorjanc.googlepages.com/legend.R")) matplot(x, y, type="l", lty=c(1, 2)) ##

-----Original Message----- From: Duncan Murdoch [mailto:murdoch at stats.uwo.ca] Sent: ned 2004-10-17 15:34 To: Gorjanc Gregor Cc: r-help at stat.math.ethz.ch Subject: Re: [R] Plotcorr: colour the ellipses to emphasize the differences On Sun, 17 Oct 2004 02:51:58 +0200, "Gorjanc Gregor" <Gregor.Gorjanc at bfro.uni-lj.si> wrote: [removed old stuff] >>I would also like to

Hi All, I have k identical parallel pieces of code running, each using n.rand random numbers.? I would like to use the same RNG (for now), and set the seeds so that I can guarantee that there are no overlaps in the random numbers sampled by the k pieces of code.? Another side goal is to have reproducibility of my results.? In?the past I have used C with SPRNG for this task, but I'm hoping

Jun hi again! Sorry for bothering. Although I had no problem with building Windows package under Linux, I noticed that man pages are not processed as they should or maybe I have a problem with my setup. My package is called 'GeneticsPed'. Instead of HTML or text files I just got a GeneticsPed.Rd.gz in man directory of a package (this is within GeneticsPed_0.1.zip file!). The file

G?ran, thanks. Happy easter! -----Original Message----- From: G?ran Brostr?m [mailto:gb at stat.umu.se] Sent: ned 2005-03-27 07:38 To: Gorjanc Gregor Cc: r-help at stat.math.ethz.ch Subject: Re: [Rd] F90 On Sun, Mar 27, 2005 at 03:26:39AM +0200, Gorjanc Gregor wrote: > Hello! > > Has anyone successfully compiled F90 sources in R-package? I found the > same question on r-devel list

On 12 October 2006 at 13:02, Gregor Gorjanc wrote: | Hello! | | I encounter the following problem when I try to load package foreign | | library(foreign) | Error in dyn.load(x, as.logical(local), as.logical(now)) : | unable to load shared library '/usr/lib/R/library/foreign/libs/foreign.so': | /usr/lib/R/library/foreign/libs/foreign.so: undefined symbol: | Rf_allocString | Error:

Hi! If I use format() on numeric vector, na.encode argument does not have any e= ffect. This was reported before: - https://stat.ethz.ch/pipermail/r-help/2007-October/143881.html - http://tolstoy.newcastle.edu.au/R/e2/devel/06/09/0360.html It works for other (say character) classes! > format(c("a", NA), na.encode=3DTRUE) [1] "a " "NA" >

Hello! I am producing a set of images and I would like them to be sorted by names I give. I was able to produce my names and add integer to them. That is easy. But my problem lies in sort of file from this process: figure_10.png figure_11.png figure_12.png ... figure_1.png figure_20.png ... So I would like to convert integers to something like 01 if upper limit for this conert is 10 or 001

Hello! Is it possible to set na.rm=TRUE in a global way? I'am constantly forgeting on this when performing analyses. I agree that one should be carefull with this when developing some code, but not necesarilly so in data analysis. Lep pozdrav / With regards, Gregor Gorjanc ---------------------------------------------------------------------- University of Ljubljana PhD student