search for: iverson

Hi, How to rbind these vectors from a list?: > l <- list(a = c(1, 2), b = c(1, 2, 3)) > l $a [1] 1 2 $b [1] 1 2 3 > do.call(rbind, l) [,1] [,2] [,3] a 1 2 1 b 1 2 3 Warning message: In function (..., deparse.level = 1) : number of columns of result is not a multiple of vector length (arg 1) > -J

...nd not computed from an algorithm) is an integer in the sense of FAQ 7.31. This seems to be the case on my machine. > identical(2.0, 2) [1] TRUE Apologies that this is such a trivial question, it seems so obvious on the surface, I just want to be sure I am understanding it correctly. Erik Iverson iverson at biostat.wisc.edu sessionInfo() R version 2.7.0 (2008-04-22) i686-pc-linux-gnu locale: LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8...

...e are other misalignments). To resolve this, I am currently going into the /tmp/Rtmp****/ directory where the .tex files are created by the dvi.latex function, and rerunning latex on the appropriate file****.tex file. Is there any thing I'm missing that would make this easier? Thanks, Erik Iverson iverson at biostat.wisc.edu

...pr("HH", ht))+1 #mean of HH position, 6 FYI, this is not homework, I have not been in school in years. I saw a similar problem posed in a blog post on the Revolutions R blog, and although I believe the answer, I'm having a hard time figuring out why this should be? Thanks, Erik Iverson

...erms.default(formula, special, data = data) : no terms component It appears from debugging that the formula argument passed into coxph is simply `~` in this case, with class "name" instead of "formula". Any ideas on how I can get this to work using this approach? Best, Erik Iverson iverson at biostat.wisc.edu > sessionInfo() R version 2.5.1 (2007-06-27) i686-pc-linux-gnu locale: LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=en_US.UTF-8;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASURE...

...ireless network lists around my laptop ,but , i can't connect to the wireless network even under the right authority my centos 5.2 dectected my wirless network card as eth1 , it's really strange , i think any help ? thx and wish you a happy new year :) -- Keep It Simple Stupid Thomas X. Iverson

...result without defining my own dummy function. Can anyone think of a way? As I said, I already have a way that works, I'm just curious if there is a more 'elegant' solution that does not rely on my having to define another function. Seems like it should be possible. Thanks, Erik Iverson

...on$People pop2=as.numeric(pop) popdich = ifelse(pop2 < 100000, "red", "blue") popnames.lower = tolower(popnames) cols=popdich[match(mapnames,popnames.lower)] plot(nz,fill=TRUE,col=cols,proj="GCS_NZGD_2000") thanks in advance for any assistance Kind regards Iverson -- View this message in context: http://r.789695.n4.nabble.com/R-help-using-R-to-build-choropleth-tp4634686.html Sent from the R help mailing list archive at Nabble.com.

...I'd really like is to name the components, and get the results of the following line without having to specify the right-hand side individually for each component: names(list.example) <- c("df.example$a > 7", "df.example$b < 4") Any ideas? Best Regards, Erik Iverson

Hello, Can anyone think of a non-iterative way to generate a decreasing geometric sequence in R? For example, for a hypothetical function dg, I would like: > dg(20) [1] 20 10 5 2 1 where I am using integer division by 2 to get each subsequent value in the sequence. There is of course: dg <- function(x) { res <- integer() while(x >= 1) { res <- c(res, x) x

...} subwrap3(testdf, sex == 'F') #does not work # Using a string, this works, however. subwrap4 <- function(x, sb) { coxph(Surv(times,event)~trt, data = x, subset = eval(parse(text=sb))) } subwrap4(testdf, "sex == 'F'") ### END R SAMPLE CODE Thanks so much, Erik Iverson iverson at biostat.wisc.edu > sessionInfo() R version 2.5.1 (2007-06-27) i686-pc-linux-gnu locale: LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=en_US.UTF-8;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASURE...

Dear all I was willing to use argument 'exclude' in function xtabs to remove some levels of factors (xtabs help page says '"exclude: a vector of values to be excluded when forming the set of levels of the classifying factors"). I tried: > mydata <- data.frame( treatment = c("B", "A", "C", "C", "B",

...C=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base Thanks, Erik Iverson

Perhaps I'm being even thicker than usual, but I can't find the tcltk package on CRAN. There is a tcltk2 package, which says that it is a collection of supplements to tcltk, but I cannot see a just-plain tcltk anywhere. If I try to install tcltk2 (from the Linux command line, or using install.packages() in R) it complains that it needs tcltk. If I try to install tcltk using

...stdf' does not exist in the global environment, but an object called 'mydata' does. I understand why, but I'd like 'wrapsave' to save the object passed to it with the name of the object it was called with, in this case 'testdf'. Is this possible? Thank you, Erik Iverson

...e groups". By my reading, this makes it impossible to summarize a continuous variable with, for example, its correlation with the response variable. Is there some sort of functionality I'm missing here, or is this just not possible with how summary.formula is written now? Thanks, Erik Iverson

Hi, Is there a function taking a function as an argument, which returns all the functions it calls, and all the the functions those functions call, and so on? I could use Rprof, but that would involve executing the function, which may miss some branches of code. I'd really like a function which looks at the source code to work out all the functions that could possibly be called. When I

Hi list. I have a numerical dataset 22,000 rows deep and 43 columns wide. I would like to remove those rows which contain only values less than 100 (ie if any value in the row is greater than 100 the row stays in the dataset). I am unsure how to test each individual value across the rows and then identify the rows which meet my criteria. Can anyone help? Thank you. Iain [[alternative HTML

...etting the geometry in .Xresources, but it does not seem to be the same "class" as WM_CLASS referenced above. Unfortunately, this is about as far as my knowledge of X goes, so hopefully this makes sense. If anyone has any pointers, they would be welcome. Thank you! Best Regards, Erik Iverson

I am attempting to generate a series of plots with cubic B-Splines and the code I am using is generating the following error message that I have been unable to resolve. Error in grid.Call.graphics("L_xspline", x$x, x$y, x$shape, x$open, x$arrow, : add_point - reached MAXNUMPTS (25200) This is the code that generates this error. plot.results <- function(dfn) { dt <-