similar to: what.is(object)

Hi list Is there a possibility to filter certain values out of an R output? In my case: I want to create a vector of p-values in a for loop that invokes for every increment cor.test() on two vectors. I haven't found a way yet to tell cor.test() to only return the p- values instead of the whole text. Thanks, Tobi

cor.test and SuppDists give me different P-values for the same Spearman's rho. Which is correct, or am I doing something wrong? > x <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1) > y <- c( 2.6, 3.1, 2.5, 5.0, 3.6, 4.0, 5.2, 2.8, 3.8) > cor.test(x,y,method="spearman") Spearman's rank correlation rho data: x and y S = 48, p-value =

I tried 2 methods to estimate C.I. of correlation coefficient of variables x and y: > x <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1) > y <- c( 2.6, 3.1, 2.5, 5.0, 3.6, 4.0, 5.2, 2.8, 3.8) #METHOD 1: Pearson's ********************************************************** > cor.test(x, y, method = "pearson", conf.level = 0.95) Pearson's

hello help, i''ve searched through the manual pages and the only reference i can find to spearman rank correlation is cor.test, which only seems to give the significance value of the correlation. is there any way to get the actual value of rho? david. [[alternate HTML version deleted]]

To fix some problems with how condition registers are saved/restored for PowerPC, I need to override TargetRegisterInfo::hasReservedSpillSlot() in PPCRegisterInfo. I've had some difficulties because of the constness of the function, and I'm wondering what the best way to handle this would be. Essentially I need to add a field to PPCRegisterInfo, and modify that field in

Dear R experts--- after many years, I am planning to give in and write my first R package. I want to combine my collection of collected useful utility routines. as my guide, I am planning to use Friedrich Leisch's "Creating R Packages: A Tutorial" from Sep 2009. Is there a newer or better tutorial? this one is 4 years old. I also plan on one change---given that the

dear R experts--- I am experimenting with multicore processing, so far with pretty disappointing results. Here is my simple example: A <- 100000 randvalues <- abs(rnorm(A)) minfn <- function( x, i ) { log(abs(x))+x^3+i/A+randvalues[i] } ?## an arbitrary function ARGV <- commandArgs(trailingOnly=TRUE) if (ARGV[1] == "do-onecore") { ?library(foreach) ?discard <-

Dear Thomas, it is, unfortunately, not that simple. t.test() returns an object of class "htest" and not all such objects have standard errors. I'm not entirely sure what the point is since it's easy to compute the standard error of the difference from the information in the object (adapting an example from ?t.test): > (res <- t.test(1:10, y = c(7:20))) Welch Two Sample

Dear List, I am trying to get R's terminal output to a file and to the terminal at the same time, so that I can walk through some tests and keep a log concurrently. The function 'sink' with the option split=TRUE seems to do just that. It works fine for most output but for objects of class htest, the terminal output is incomplete (the lines are there but empty). Here is an

Dear R wizards: is there a clean way to assign to elements in a list? what I would like to do, in pseudo R+perl notation is f <- function(a,b) list(a+b,a-b) (c,d) <- f(1,2) and have c be assigned 1+2 and d be assigned 1-2. right now, I use the clunky x <- f(1,2) c <- x[[1]] d <- x[[2]] rm(x) which seems awful. is there a nicer syntax? regards, /iaw ---- Ivo Welch

dear r experts---Is there a multicore equivalent of by(), just like mclapply() is the multicore equivalent of lapply()? if not, is there a fast way to convert a data.table into a list based on a column that lapply and mclapply can consume? advice appreciated...as always. regards, /iaw ---- Ivo Welch (ivo.welch at gmail.com)

is it possible to throw a stop() that is so hard that it will escape even tryCatch? /iaw ---- Ivo Welch (ivo.welch at gmail.com)

Why does stop("we are done") print "Error in eval.with.vis(expr, envir, enclos) :" ? It would seem to me that a plain stop() is not an error, and that it would make more sense to have an error() function that is different from a stop(). Is there a rationale here that I am missing? sincerely, /iaw

hi: I searched the last 2 hours for a way to enter a data frame directly in my program. (I know how to read from a file.) that is, I would like to say something like d <- this.is.a.data.frame( c("obs1name", 0.2, 0.3), c("obs2name", 0.4, 1.0), c("obs3name", 0.6, 2.0) ,

dear R experts---I am trying to replicate a perl feature. I want to be able to embed R commands inside a character string, and have the string be printed with the command executed. my perl equivalent is my $a=10; my $teststring = "the expression, $a+1, is ::$a+1::, but add one more for ::$a+2::\n"; $teststring =~ s/::(.*?)::/$1/gee; print $teststring; of course, R does not use

dear R wizards: an operation I execute often is the deletion of all observations (in a matrix or data set) that have at least one NA. (I now need this operation for kde2d, because its internal quantile call complains; could this be considered a buglet?) usually, my data sets are small enough for speed not to matter, and there I do not care whether my method is pretty inefficient (ok, I

dear R experts---I am using the coef() function to pick off the coefficients from an lm() object. alas, I also need the standard errors and I need them fast. I know I can do a "summary()" on the object and pick them off this way, but this computes other stuff I do not need. Or, I can compute (X' X)^(-1) s^2 myself. Has someone written a fast se() function? incidentally, I think

every time I read the R release notes for the next release, I see many functions that I had forgotten about and many functions that I never knew existed to begin with. (who knew there were bibtex facilities in R? obviously, everyone except me.) I wonder whether there is a complete list of all R commands (incl the standard packages) somewhere, preferably each with its one-liner AND

this is not an important question, but I wonder why lm returns an error, and whether this can be shut off. it would seem to me that returning NA's would make more sense in some cases---after all, the problem is clearly that coefficients cannot be computed. I know that I can trap the lm.fit() error---although I have always found this to be quite inconvenient---and this is easy if I have only

Apologies, basic question on plot. y <- c(-4,3,-2,1); x <- c("time 1", "time 2", "time 3", "time 4"); plot(x,y, type="b"); of course fails. x <- 1:4 makes it succeed, but then I have too many ticks on my X axis. I want exactly 4 tickmarks. It would also be nicer if I could name the ticks. I looked at ?par and Venables&Ripley,