similar to: counting row repetitions without loop

Dear list, I would like to generate random numbers from a Voigt distribution, hopefully in a way as simple as getting random numbers from a normal distribution with 'rnorm'. Is there any package to do this? Speed is an issue in this application. Or, as the Voigt distribution is a convolution of a Gaussian and a Lorentzian, can I simply combine random numbers from rnorm and rcauchy in some

Dear R community, I would be greatful if somebody could shed light on the following. I have created a set of 6 points to check how centroid agglomeration works in cluster analysis: > Y <- data.frame(x=c(-1,1,1,-1,10,12),y=c(1,1,-1,-1,0,0)) It is quite intuitive to understand that the last clusters to be joined will be {1,2,3,4} with {5,6}. Now, the centroid for the first cluster has

Dear R developers, I have a new class, which I called "Molecule", and have tried to define = a "+" operation for 2 objects of this class. This is what I have written so far, although the method is not complete = (I'm trying to look at it at intermediate stages): setMethod( f=3D"+", signature(x=3D"Molecule",y=3D"Molecule"),

Dear R gurus and users, I seem to have problem finding the right tool for plotting convex hulls over 2D plots, after a cluster analysis. In fact I would like to draw a convex hull in 2D for a generic group of points. I found a "convhulln", but this doesn't seem to give me a convex hull. Here is what I do: > library(mvtnorm) > Mean <- c(2,1) > Sigma <-

Dear RGL experts, I haven't been able to add greek letters to my rgl plot3d. I have tried "expression" with no success. Here is the interested bit: > library(rgl) > cb <- cube3d() > plot3d(cb,xlab=expression(alpha),ylab="",zlab="",box=FALSE,alpha=0.5) The expression(alpha) appears as "alpha", rather than as a greek symbol. I suspect greek

Dear R community (and Duncan more specifically), I can't work out how to make additional light sources work in rgl. Here is the example. First I create a cube and visualize it: > cubo <- cube3d(col="black") > shade3d(cubo) Next I position the viewpoint at theta=0 and phi=30: > view3d(theta=0,phi=30) Next, I want to create a 2nd light source which diffuses red light

Dear R community, I would like to source a file in my linux system to set a few environment variables. I have tried: > system("source /home/james/build//ccp4-6.1.13/include/ccp4.setup") and got: sh: source: not found In fact, using Sys.which("source") I get an empty string. How can I source that ccp4.setup file from within an R session? J Dr James Foadi PhD Membrane

Dear R community, I'm trying to develop a fast way of summing specific rows of a large data frame. Here is an example of the kind of data frames I'm dealing with: > refls H K L M/ISYM BATCH I SIGI 43247 1 0 5 21 79 61.44117 2.20553 1040 1 0 5 257 6 15.16316 0.54431 2324 1 0 5 257 5 46.76152 1.67858 31515 1 0 5 259 60 57.97305

Dear all, I have a very small script to plot a function. Here it is: ########################################## sinca <- function(N,th) { return(sin((N+0.5)*th)/sin(0.5*th)) } plot_sinca <- function(N) { x <- seq(-5*pi,5*pi,by=pi/100) y <- rep(0,length=length(x)) for (i in 1:length(x))y[i] <- sinca(N,x[i]) plot(x,y,type="l",ylim=c(0,2*N+4)) return(c(x,y)) }

Dear R gurus and users, I'm having problems with the use of write.table. I have a 28-variables data frame create at each cycle of a loop; it can contain between 2000 and 3000 rows for each cycle. After each cycle the data frame is written out to a file with the "append=TRUE" option and then removed from memory. These are the couple of lines involved: > data2 <-

Dear R users/gurus, I have recently installed Ubuntu 7.10 (Gutsy) on my Dell Latitude D520. First time for me, as I have always used Fedora. R installation using "apt" has proceeded fine, both for r-base and r-base-dev. Next, I need to install various other packages, "rgl" being first in the line. I have sorted some preliminaries by installing headers for X11 and Mesa: sudo

Hello All, I wanted to perform a fourier transform on high frequency financial data. I have searched and have not found much on this topic for R. I was wondering if anyone has used any libraries for it or have come across any papers I may read. Many Thanks, Neil Gupta [[alternative HTML version deleted]]

Hi, I am writing because I would like to know if centos will update the kernel-rt package from https://git.centos.org/rpms/kernel-rt/releases ? I ask because bz#1550584 https://bugzilla.redhat.com/show_bug.cgi?id=1550584 mentions "fixed in version kernel-rt-3.10.0-1063.rt56.1023.el7" And the latest tag in "releases" link above is for kernel-rt-3.10.0-1062.7.1.rt56.1030.el7

Hi I don't know where the announcements go. I cloned https://git.centos.org/git/rpms/kernel.git and grabbed the changelog section pasted below from the git diff in SPECS/kernel.spec (commit dated 4 Feb 2020) Regards +* Thu Dec 12 2019 Frantisek Hrbata <fhrbata at hrbata.com> [3.10.0-1062.12.1.el7] +- [powerpc] powerpc/pseries: Remove confusing warning message (Gustavo Duarte) [1780148

Dear all, While running some code requiring the "tcltk" package I have realised that my version of R was compiled with the Tcl/Tk libraries included in Fedora 6. It would be for me better to use the ActiveTcl libraries (which I have under /usr/local), and I'm aware that this probably means to recompile R with the proper configuration variables. But...is it by any chance

Dear R community, I'm in the process of writing a package. One of the functions reads a whole ascii file, full of data: > SYMINFO <- scan("/file/full/of/data.txt",what="character",sep="\n",quiet=TRUE) SYMINFO is now an array character with roughly 14000 elements. I thought the best way to use this file in the package was to list it, together with the

On our CentOS 7 servers I see there is a new kernel available when doing yum update: kernel 3.10.0-1062.12.1.el7 kernel-devel 3.10.0-1062.12.1.el7 However I am not able to find any announcement at redhat.com for this kernel. The only announcement I find is this from CentOS: https://www.mail-archive.com/centos-announce at centos.org/msg11573.html However that CentOS announcement links to

Hi, We run four-node Lustre 2.3, and I needed to both change hardware under MGS/MDS and reassign an OSS ip. Just the same, I added a brand new 10GE network to the system, which was the reason for MDS hardware change. I ran tunefs.lustre --writeconf as per chapter 14.4 in Lustre Manual, and everything mounts fine. Log regeneration apparently works, since it seems to do something, but

Some time ago I was allowed to discuss "Diamond Graphs", and whether they would be useful in R, in this mailing list. The August 2003 issue of The American Statistician has finally arrived here and I have been able to read the article. A number of points of interest arise. 1. The article is "A Diamond-Shaped Equiponderant Graphical Display of the Effects of Two

I apologise for starting a new thread, but we had a mail problem and I don't have the original message to refer to. Someone mentioned the new "Diamond Graphs" invented at Johns Hopkins. I haven't see the August 2003 issue of The American Statistician yet, but I _have_ read the press release. The press release is a bit of a stunner. I quote: "Who would have thought we