similar to: Filling points in a trellis object

Hi everyone, I am using eigen() to extract the 2 major eigenpairs from a large real square symmetric matrix. The procedure is already rather efficient, but becomes somehow slow for real time needs with moderately large matrices (few thousand lines). The R implementation statically extracts all eigenvalues (and optionally associated eigenvectors). I heard about optimizations of the eigen

I have a stupid problem that is currently driving me crazy... Let us suppose that I want to draw a big red square in the middle of my current device (say X11) I tried the following code : pushViewport(viewport(xscale=c(0,1), yscale=c(0,1)), just=c("center", "center")) vp1 <- viewport(x=unit(0.5, "native"), y=unit(0.5, "native"), width=unit(0.4,

Hi everyone, >From the currently available version of the package VBmix, I would like to retrieve the intermediate .tex file that generates the VBmix-manual.pdf file. Formerly, running R CMD check with --no-clean allowed to get this tex source from a hidden directory : this feature was removed, but can apparently still be accessed through R CMD Rd2pdf --no-clean. Surprisingly, while the

(sorry for repetition: the previous mail resulted from a weird manipulation in the forum) With the following code : dat1 <- matrix(nrow=4, ncol=2) dat1[1,] <- c(-2, 1) dat1[2,] <- c(-1.7, 0.9) dat1[3,] <- c(0.1, 0.6) dat1[4,] <- c(0.5, 0.5) theplot <- xyplot(V2 ~ V1, as.data.frame(dat1), pch=c(4,1,5,4)) plot(theplot, prefix="theplot") # for a predictable name

Hello everyone, Considering the following code sample : ---- indexes <- function(vec) { vec <- which(vec==TRUE) return(vec) } mat <- matrix(FALSE, nrow=10, ncol=10) mat[1,3] <- mat[3,1] <- TRUE ---- Issuing apply(mat, 1, indexes) returns a 10-cell list, as expected. Now if I do: ---- mat[1,3] <- mat[3,1] <- FALSE apply(mat, 1, indexes) ---- I would expect a

Dear R contributors, I'm running into trouble when trying to allocate some large (but in theory viable) vector in the context of C code bound to R through .Call(). Here is some sample code summarizing the problem: SEXP test() { int size = 10000000; double largevec[size]; memset(largevec, 0, size*sizeof(double)); return(R_NilValue); } If size if small enough (up to 10^6), everything is

Dear All,   any thoughts on how I can do the following:   let us say we have:   a <-c(2,4,16,28,48) b <-c(10,4,2,0.4,0.03) d <-cbind(a,b)   what I would like to do is to extract values of column b in the matrix based on the corresponding values of column a. For example: I would like to extract all b values that have a corresponding a value that is less than 24 into a numeric vector, so

Dear R contributors, Considering the following sample C code, that illustrates two possible uses of a Cholesky decomp for inverting a matrix, equally valid at least in theory: SEXP test() { int d = 2; int info = 0; double mat[4] = {2.5, 0.4, 0.4, 1.7}; double id[4] = {1.0, 0.0, 0.0, 1.0}; double lmat[3]; F77_CALL(dpotrf)("L", &d, mat, &d, &info); lmat[0] = mat[0]; lmat[1]

Good day. I create a data frame like this: > data <- data.frame(a=1:10,b=11:20,c=21:30) I can subset this data.frame by saying: > data[data$a>7,] and I get this result a b c 8 8 18 28 9 9 19 29 10 10 20 30 I understand I should get the same result by saying

Hi all. I would like to highlight overlapping regions of two densities and I could not find a way to do it. Here is the sample code: myd <- c(2,4,5, 4,3,2,2,3,3,3,2,3,3,4,2,4,3,3,3,2,2.5, 2, 3,3, 2.3, 3, 3, 2, 3) myd1 <- myd-2 plot(range(density(myd)$x, density(myd1)$x), range(density(myd)$y, density(myd1)$y), type = "n") lines(density(myd), col=1, lwd=4)

Dear R contributors, Say I want to debug some C code invoked through .Call() - say "varbayes" in the VBmix package. following the instructions in "Writing R Extensions", I perform the following actions : R -d gdb run library(VBmix) CTRL+C break varbayes signal 0 mod <- varbayes(as.matrix(iris)[,1:4], 2) The breakpoint is indeed activated, seemingly at the correct position

Hi Pierrick, You're storing largevec on the stack, which is probably causing a stack overflow. Allocate largvec on the heap with malloc or one of the R memory allocation routines instead and it should work fine. Karl On Thu, Dec 18, 2014 at 12:00 AM, Pierrick Bruneau <pbruneau at gmail.com> wrote: > > Dear R contributors, > > I'm running into trouble when trying to

This isn't the help list for LAPACK, but as far as I can tell, dppsv expects a symmetric matrix input compacted as triangular, not a Choleski decomposed one. So try assigning lmat before the call to dpotrf. -pd > On 20 Dec 2014, at 22:06 , Pierrick Bruneau <pbruneau at gmail.com> wrote: > > Dear R contributors, > > Considering the following sample C code, that

Problem: I am trying to get the histogram and box plot x axis to match. I?ve tried using the expand_limits function to make the axis match but that didn?t make the axis match. The histogram?s axis are still consistently larger than the ones for the box plot (though the function did help). Does anyone have a suggestion as to what I should do instead? Background: I am building a Shiny app that

Hi, I feel that i'm having a small problem with my upsmon configuration (maybe even simply a permissions problem) since I upgrade today to 2.0.2 I have apcsmart set up and talk to my ups no problem (upsc get answer). I have upsd running no problem! But when I try to start upsmon, it's unable to talk to my upsd. Having try a few option I found that saying RUN_AS_USER=root makes my upsmon

Hello, I'm trying to tranform a numeric vector into a character vector. > x=c(2.00,1.20,5.00,6.56) > y= as.character(x) > y [1] "2" "1.2" "5" "6.56" What I want is : [1] "2.00" "1.20" "5.00" "6.56" Does someone know how to do this please ? Benoit Bruneau

On Thu, 2014-12-11 at 14:00 +0100, Pierrick Bruneau wrote: > Dear R contributors, > > Say I want to debug some C code invoked through .Call() - say > "varbayes" in the VBmix package. following the instructions in > "Writing R Extensions", I perform the following actions : > > R -d gdb > run > library(VBmix) > CTRL+C > break varbayes >

Hi! My boss asked me to be able to share some directories on a FreeBSD Samba server with users already created (and used) on a NT 4.0 PDC server. So I followed http://www.sugoi.org/bits/index.php?bit_id=10 I just replaced the Windows 2k part by adding the NetBIOS name of the machine to the NT PDC. When I did the smbpasswd trick it worked like charm: testsmb# smbpasswd -j CH-DOMAIN -r PDC

I've been reading this list for a few weeks now and I've given advice on questions that look challenging but I've deleted MANY questions like these: "How do I (easy question found in the documents)?" Though I don't count myself an expert, I've known enough experts to see that they _HATE_ it when you don't invest some time before asking a question. I too have

hello, I'm looking for a R function proceeding MRPP (Multi-Response Permutation Procedures). Is it available? Mielke,P.W., Jr. 1984. Meteorological applications of permutation techniques based on distance functions. Pages 813-830. In P.R. Krishnaiah and P.K, Sen, eds, Handbook of Statistics, Vol.4. Elsevier Science Publisher. I know PC-Ord can do it, but I'd rather a R function