similar to: RcppProgress: progress monitoring and interrupting c++ code, request for comments

Hello to the list I'm using a Rscript launched via the linux terminal as the title of this mail explains. The following is the code of the Rscript : myRscript.r -------------------------------- 1 #!/usr/bin/Rscript 2 3 #Read data 4 data=read.table("/home/morisseau/stage/recherche/gwas/data/CFHLC5/CFHLC_format_files_rm_07052010/output/windowstemp.txt") 5 6 data=t(data) 7

Hi all! I am using GenABEL on R for GWAS analysis. I am having a couple of issues: First, I am having a problem reading files (.map, & .ped, size 900Mb, using windows 32-bit) onto R in the "convert.snp.ped" statement. I am thinking this problem is likely due to the large size of the files & my version of R is not able to handle them, since I can read in smaller files.

Dear R gurus, I am interested in permutations-based cpu-intensive methods so I had to pay a little attention to Random Number Generators (RNG). For my needs, RNGs have to: 1) be fast. I profiled my algorithms, and for some the bottleneck was the RNG. 2) be scalable. Meaning that I want the RNG to remain fast as I add threads. 3) offer a long cycle length. Some basic generators have a

Hi R-help folks, I have been doing some single SNP association work using snpMatrix. This works well, but produces a lot of false positives, because of population structure in my data. I would like to correct the p-values (which snpMatrix gives me) for population structure, possibly using principle component analysis (PCA). My data is complicated, so here's a simple example of what

Hi Peter, Rui, Chrstophe and Gabriel, Thanks for your inputs -- the use of qnorm(., log=TRUE) is a good point in line with pnorm with which we devised log(p) as log(2) + pnorm(-abs(z), lower.tail = TRUE, log.p = TRUE) that could do really really well for large z compared to Rmpfr. Maybe I am asking too much since z <-20000 >

include/Rmath.h declares a set of 'logspace' functions for use at the C level. I don't think there are core R functions that call them. /* Compute the log of a sum or difference from logs of terms, i.e., * * log (exp (logx) + exp (logy)) * or log (exp (logx) - exp (logy)) * * without causing overflows or throwing away too much accuracy: */ double Rf_logspace_add(double

I am trying to compile R-2.4.1 in 64-bit on Solaris 10 running on AMD hardware. I am trying to do this with Sun Studio 11. My config.site looks like this: #! /bin/sh AR="/usr/ccs/bin/ar" TEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/tex" LATEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/latex" PDFTEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/pdftex"

I am trying to compile R-2.4.1 in 64-bit on Solaris 10 running on AMD hardware. I am trying to do this with Sun Studio 11. My config.site looks like this: #! /bin/sh AR="/usr/ccs/bin/ar" TEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/tex" LATEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/latex" PDFTEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/pdftex"

Hi All, Thanks for all your comments which allows me to appreciate more of these in Python and R. I just came across the matrixStats package, ## EXAMPLE #1 lx <- c(1000.01, 1000.02) y0 <- log(sum(exp(lx))) print(y0) ## Inf y1 <- logSumExp(lx) print(y1) ## 1000.708 and > ly <- lx*100000 > ly [1] 100001000 100002000 > y1 <- logSumExp(ly) > print(y1) [1] 100002000

I agree with many the sentiments about the wisdom of computing very small p-values (although the example below may win some kind of a prize: I've seen people talking about p-values of the order of 10^(-2000), but never 10^(-(10^8)) !). That said, there are a several tricks for getting more reasonable sums of very small probabilities. The first is to scale the p-values by dividing the

I am trying to compile R-2.4.1 in 64-bit on Solaris 10 running on AMD hardware. I am trying to do this with Sun Studio 11. My config.site looks like this: #! /bin/sh AR="/usr/ccs/bin/ar" TEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/tex" LATEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/latex" PDFTEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/pdftex"

I am trying to compile R-2.4.1 in 64-bit on Solaris 10 running on AMD hardware. I am trying to do this with Sun Studio 11. My config.site looks like this: #! /bin/sh AR="/usr/ccs/bin/ar" TEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/tex" LATEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/latex" PDFTEX="/usr/local/teTeX/bin/i386-pc-solaris2.10/pdftex"

>>>>> William Dunlap via R-devel >>>>> on Sun, 23 Jun 2019 10:34:47 -0700 writes: >>>>> William Dunlap via R-devel >>>>> on Sun, 23 Jun 2019 10:34:47 -0700 writes: > include/Rmath.h declares a set of 'logspace' functions for use at the C > level. I don't think there are core R functions that call

Dear R experts, recently I started developing a Rcpp package "OpenMPTest". Within that package I want to use OpenMP, as in the following code example: // header file #include <omp.h> using namespace Rcpp ; RcppExport SEXP testOpenMP( SEXP nThreads ) ; // cpp file SEXP testOpenMP( SEXP nThreads ) { BEGIN_RCPP NumericVector numberThreads = NumericVector( nThreads );

Hi, I was wondering if anybody knew how to trap SIGINTs (ie Ctrl-C) in backend C++ code for R extensions? I'm writing a package that uses the GPU for some hefty matrix operations in a tightly coupled parallel algorithm implemented in CUDA. The problem is that once running, the C++ module cannot apparently be interrupted by a SIGINT, leaving the user sat waiting even if they realise

You may want to look into using the log option to qnorm e.g., in round figures: > log(1e-300) [1] -690.7755 > qnorm(-691, log=TRUE) [1] -37.05315 > exp(37^2/2) [1] 1.881797e+297 > exp(-37^2/2) [1] 5.314068e-298 Notice that floating point representation cuts out at 1e+/-308 or so. If you want to go outside that range, you may need explicit manipulation of the log values. qnorm()

Hi! I am looking for suggestions how I can add rules (or change the configuration otherwise) without interrupting of services. Right now, when I need to do a ''shorewall restart'' all services are not available during this time. I tried playing with the ''routestopped'' parameter but without success. How do you get around this? Thanks in advance, Christian

Hi all, I am hoping someone can help me out on this. I want to be able to interrupt MOH every X seconds after the DIAL command is executed. The interrupt greeting is something like "please wait while we transfer your call". How can I do that? Within the DIAL options, I can't see any announce frequency or options that can help. Could anyone please tell me how that function can

Hello, if I interrupt Sweave while it's processing a file it seemingly leaves an open sink connection that hides printed output. Can this be changed to reset the sink on exit? I've been baffled by this for years. This is seen in Windows (R Under development (unstable) (2013-10-20 r64082)) and an older Linux (R version 3.0.0 (2013-04-03)). Run the code below in two parts with a manual

Dear fellow R users, is there a generic way to gracefully interrupt an R function without terminating the entire session ? I am mainly interested in this answer for Linux and MacOS. I found neither Esc nor Ctrl-C to work; it seems that R does not check for signals periodically? Also, an entirely unrelated question: I have been looking unsuccessfully for the R sources for the examples given in