similar to: Add row to data frame

Hi all, I tried to save a complete log of a R session we had in a seminar today... but I didn't succeed. 1) R | tee session.log This saves both input and output, but I do get the cursor key escape sequences from editing (cursor-up to get last command etc) instead of the actual command line executed. 2) savehistory Gets commands only, not the output 3) sink Gets output, without the

Hi all, Given a number of points (x,y) in a plane, I'd like to plot a map of polygons, so that 1) each polygon contains exactly one point 2) the polygon defines the area for which this specific point is closer than any other point. It's a bit like a map of areas "influenced" by that point, and it's obviously a matter of intersecting the perpendicular bisectors

Dear all, I hope this is the right mailing list for my question -- I felt that this was too technical for R-help. I am trying to compile R-2.4.0 on a HP-UX system: ./configure MAKE=gmake --prefix=$HOME --without-x --without-tcltk --disable-R-profiling --without-readline --disable-multibyte R is now configured for ia64-hp-hpux11.23 Source directory: . Installation directory:

Hi all, Is there a function to check if a particular value is contained in a vector? I've looked at grep in the hope that I could use a Perl-like syntax, but obviously it's different... I'd like to do something like: y <- c("a","b","c") if("a" in y) { # "a" is not in y } Also, is there a way to

I am trying to plot a figure within a figure (an inset that shows a closeup of part of the data set). I have searched R-help and other sources but not found a solution. What I would like to do is (1) produce a plot (2) specify a window that will be used for the next plot (in inches or using the coordinate system of the plot produced in (1) (3) overlay a new plot in the window specified under (2)

Hi all, I could not find a way to add the dashed horizontal separator line to the top panel of a coplot (the one indicating which panel belongs to which range of the given variable a in y~x | a). This line should separate the ranges plotted in a row of panels in the panel matrix below (hmmm, hope this is reasonably clear). I like it because it enhances readability. I'll try a sketch: >

Hi all, I wonder how to correctly write the following expression (it's the axis label in a plot command): ylab=expression(y' == y - bar(y) ) Somehow the single quote in y' is causing the problems, I guess because it is interpreted as a quote... Does it have to be escaped? But how? Thanks for your help Pascal

Hi all, I am playing around with latin squares, and wrote a recursive function that searches for valid combinations. Apart from the fact that there are very many, I run into troubles beginning with size 10x10 because the recursion depth becomes too large (max of 10x9-1=89 in this case). Why is this a problem? Isn't there enough space allocated to the stack? Can this be increased? The

Hi all, I repeatedly encounter the following problem: After importing a data set into a data frame, I wish to set a column with numeric values to be a factor, but can't figure out how to do this. Also, I do not wish to write as.factor(x) all the time. I can create a new vector with x <- factor(x), but the new vector resides outside the attached data frame. Pascal > attach(ngrad)

Hi all, I'd like to extract a value from an ANOVA table, but experience the following problem: ### This works: > s.pseudo <- summary(aov(m ~ block + mix*graz,data=split1)) > s.pseudo Df Sum Sq Mean Sq F value Pr(>F) block 2 1114.66 557.33 4.4296 0.04192 * mix 1 6.14 6.14 0.0488 0.82956 graz 2 1.45 0.72 0.0057 0.99427 mix:graz

Given a set of data: > names(data) [1] "city" "house" "visit" "value" I am looking for a way to compute the variance components of the nested model (ie, visit 1 at house 2 at city 3 isn't related to visit 1 and house 2 at city 4), but different houses in the same city may be related, and different visits to the same house are probably

Dear list, My apologies if a solution / explanation to this already exists on the list, but it is difficult to assign it to a certain keyword. test<-c(1:3) testfct <- function(x) {test[1]<-100} test [1] 1 2 3 testfct(1) [1] 1 2 3 Basically, I would like to write data into an external table that the function does not know. Why is this not working / what alternatives exist? Thanks,

hi, this may be trivial, but we can't seem to find anything adequate, (although there is a work around with match() ). We are looking for something along the lines of plot(table1[table1$var2==or("a","b","c","d"),"var1"]) would be handy, with the potential or() function leading to what plot(table1[table1$var2=="a" |

Dear R-list I am (forced) to change from Linux to Mac and am now looking for a new editor for R. I would like one that features a split window (console + editor) as well as syntax highlighting. Can anyone help? Especially the split-window feature does not seem to be easily available in the editors desribed on the R-help site, except Emacs, which I am reluctant to start using. I am on a MacBook Air

Hi all, I posted this question several days ago, but did not get any answer until now. Since I still have no clue about the source of this error message, I repost a description of the problem including some code: A student at our institute fitted an aov model, and got the following error message: > plot(p.aov) Hit <Return> to see next plot: Hit <Return> to see next plot:

Full_Name: Pascal Version: 1.9.0 OS: GNU/Linux Submission from: (NULL) (212.152.21.2) In order to compile under SuSE 9.1, I needed to change #define NeedFunctionPrototypes 0 to #define NeedFunctionPrototypes 1 in line 29 of src/modules/X11/dataentry.c; /* don't use X11 function prototypes (which tend to ...): */ #define NeedFunctionPrototypes 1 #include <X11/X.h> HTH Pascal

As many of you have probably realized, the spam filtering at @stat.math.ethz.ch has been dead for since yesterday (09-25) ~16:50 till today ~08:30. The sudden death may have been caused by unrelated installation of some perl modules (spamassassin *is* running on perl) by our IT staff. We are very sorry for this event. On the bright side: You have been able to get a glimpse of what you are

Hi All: I am clustering 500 genes using hclust of R. Visualizing cluster membership becomes difficult with so many genes in each cluster...Is there a way of printing the dendrogram in multiple pages so that I can clearly see what is in each cluster? Thanks in advance. Karthi.

Dear all I imported a Stata .dta file with the read.dta-function from the foreign-package. The dataframe's dimensions are > dim(d.apc) [1] 15806 1300 Importing needs up to 15 min and calculations with these data are rather slow (although I subset the data before starting analyses). My questions are: 1. Has someone experiences importing Stata files (alternatives to read.dta) ? 2.

Dear R-list member, I'm using the function pmnorm() (-->library(mnormt)) within a for-loop. Certain parameter values leads to an error message: "(In sqrt(diag(S)) : NaNs produced, In sqrt(1/diag(V)) : NaNs produced, In cov2cor(S) : diag(.) had 0 or NA entries; non-finite result is doubtful)" obviously because "NaNs" were produced. Is it possible to tell R that it