similar to: sourcing a linux file with "system"

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 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 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'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 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]]

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 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

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

Hi, I have a data frame consisting of coordinates on a 10*10 grid, i.e. > example x y 1 4 5 2 6 7 3 6 6 4 7 5 5 5 7 6 6 7 7 4 5 8 6 7 9 7 6 10 5 6 What I would like to do is return an 10*10 matrix consisting of counts at each position, so in the above example I would have a matrix where, for example, cell [4,5] contains 2 and [6,7] contains 3. At the

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

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 all, I have looked for an answer for a couple of days, but can't come with any solution. I have a set of functions, say: > t0 <- function(x) {1} > t1 <- function(x) {x} > t2 <- function(x) {x^2} > t3 <- function(x) {x^3} I would like to find a way to add up the previous 4 functions and obtain a new function: > rrr <- function(x) {1+x+x^2+x^3} without,

Dear all, I have tried installing "rgl" with the usual command: R CMD INSTALL rgl_0.67-2.tar.gz Differently from what happened last time I have succesfully installed this package, this time there was a failure: ... ...g++ -I/usr/lib/R/include -I/usr/lib/R/include -I -DHAVE_PNG_H -I/usr/include/libpng12 -I/usr/local/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2

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

Hello. I wonder if anyone can help me with this. I'm performing cluster analysis by using hclust in stats package. My data are contained in a data frame with 10 columns, named "drops". Firs I create a distance matrix using dist: distanxe <- dist(drops) Then I perform cluster analysis via hclust: clusters <- hclust(distanze) At this point I want to view the tree

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