[R] using a variable for a column name in a formula

2013-10-13 Thread David Epstein
lm(height ~ ., data=X) works fine. However nnn <- "height" ; lm(nnn ~ . ,data=X) fails How do I write such a formula, which depends on the value of a string variable like nnn above? A typical application might be a program that takes a data frame containing only numerical data, and figures ou

[R] knitr: Was previously " Sweave: printing an underscore in the output from an R command"

2013-09-02 Thread David Epstein
} > \Sexpr{sub("_", "\\\\_", version$platform)} > \end{document} > > Regards, > Yihui > -- > Yihui Xie > Web: http://yihui.name > Department of Statistics, Iowa State University > 2215 Snedecor Hall, Ames, IA > > > On Mon, Sep

Re: [R] Sweave: printing an underscore in the output from an R command

2013-09-02 Thread David Epstein
ot;, version$platform))} When Sweave is run, this expression is evaluated to illegal LaTeX David. On 2 Sep 2013, at 16:47, ONKELINX, Thierry wrote: > You have to escape the underscore > > \Sexpr{gsub("_", "\_", print(version$platform))} > > Best regards,

[R] Meaning of "Integer,19"

2013-09-02 Thread David Epstein
I tried example('apply'). Among the various examples, there was the following: apply> z <- array(1:24, dim = 2:4) apply> zseq <- apply(z, 1:2, function(x) seq_len(max(x))) apply> zseq ## a 2 x 3 matrix [,1] [,2] [,3] [1,] Integer,19 Integer,21 Integer,23 [2,] Integer,20 In

[R] Sweave: printing an underscore in the output from an R command

2013-09-02 Thread David Epstein
I am working with Sweave and would like to print out into my latex document the result of the R command version$platform So what I first tried in my .Rnw document was \Sexpr{print(version$platform)}. However, the output from this command is the string "x86_64-apple-darwin10.8.0" (without the quo

[R] Unsuccessful beginner's struggle with lm

2013-08-29 Thread David Epstein
I have two data frames, "train" and "response". Here is my attempt to do a linear regression. All entries of both data frames are numeric. I am expecting the intercept value to lie between 2 and 3 (in particular, non-zero). Here is a record of my interaction with R: > class(response) [1] "data.fr

[R] Mac Os X 10.6.8, R and edit/vi

2013-08-08 Thread David Epstein
I tried using various versions of the 'edit' command. Here is an account of how this failed. I hope I have included all relevant information. I haven't used R for a couple of years. Before restarting with R, I downloaded the latest version I could find in its binary version, and installed it witho

[R] lmer (LME4) and survey standard error

2013-03-09 Thread David Epstein
Hello, I am using lmer (LME4) to build a model from data for 19 different neighborhoods drawn, in part, from the American Communities Survey (ACS). The ACS data is static while other variables change over the five years under investigation. I am new to mixed effects models and was hoping that some

Re: [R] Recoding multiple TRUE/FALSE columns into a single list of TRUE columns

2011-12-25 Thread David Epstein
Jim, Wow, that does it! I think I can use strsplit and unlist to convert the string of row names into a R list. thank you! -david __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http:

[R] Recoding multiple TRUE/FALSE columns into a single list of TRUE columns

2011-12-25 Thread David Epstein
Hi everyone, I need to recode multiple columns in a dataframe into a single column in a variety of different ways. Often the values will be TRUE/FALSE and I want a list of the columns that are true as in the Result column below: P1 P2 P3 P4 Result 1 0011P3,P4 2

[R] side-by-side map with different geographies using spplot

2011-12-02 Thread David Epstein
Hello, I want to create side-by-side maps of similar attribute data in two different cities using a single legend. To simply display side-by-side census block group boundary (non-thematic) maps for Minneapolis & Cleveland I do the following: library(rgdal) library(sp) Minneapolis=readOGR(".

Re: [R] replacing values in one vector with corresponding index from another vector

2011-10-19 Thread David Epstein
thank you! that is straight forward. On Wed, 2011-10-19 at 22:37 +, William Dunlap wrote: > > match(b, a) > [1] 4 5 2 3 1 5 5 1 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http:/

[R] replacing values in one vector with corresponding index from another vector

2011-10-19 Thread David Epstein
A quick question for the gurus... Given: a=c(58,73,100,40,70) b=c(40,70,73,100,58,70,70,58) How can I replace the elements of "b" with the corresponding index numbers from "a" that start at 1? All values in "a" are unique. So, I end up with: b=c(4,5,2,3,1,5,5,1) I believe I need to use one of

Re: [R] Using SQL "IN" with RJDBC or RMySQL

2011-10-19 Thread David Epstein
Gabor, thank you. I just got the script working using paste but your solution looks good as well. I was not familiar with gsubfn. Appears very useful. -david > With gsubfn if you preface your function with fn$ as shown below then > it turns on a quasi-perl style string interpolation: > > libr

Re: [R] Using SQL "IN" with RJDBC or RMySQL

2011-10-19 Thread David Epstein
thank you! I was able to get it to work with collapse="," On Wed, 2011-10-19 at 01:52 -0500, Jeff Newmiller wrote: > paste("SELECT * FROM this_table WHERE this_column IN (", > paste(org_table$id, collapse=TRUE),")",sep="") __ R-help@r-project.org maili

[R] Using SQL "IN" with RJDBC or RMySQL

2011-10-18 Thread David Epstein
Hello, The code below works fine up until I try to use the "IN" statement in the last line. The proper SQL format is: SELECT * FROM this_table WHERE this_column IN (1,2,3,4,5) But, I think I may be getting something like: SELECT * FROM this_table WHERE this_column IN c(1,2,3,4,5) Which makes

[R] Problems with Treemaps in VCD & PORTFOLIO

2010-08-08 Thread David Epstein
Hi, I am trying to create a treemap for a non-symmetric tree where some branches have more sub-branches than others. I have tried tools in two packages, but have encountered some problems that I hope more experienced users on this list can help address. The PORTFOLIO package offers MAP.MARKET, b

Re: [R] trying to install plotrix

2009-09-22 Thread David Epstein
David Epstein wrote: > > Here are the error messages: >> install.packages('plotrix') > Warning in install.packages("plotrix") : > argument 'lib' is missing: using '/Users/dbae/Library/R/library' > --- Please select a CRAN mirror f

Re: [R] trying to install plotrix

2009-09-22 Thread David Epstein
David Epstein wrote: > > Here are the error messages: ...snip ... > Sorry. I rebooted my machine and then plotrix loaded without complaint. D -- View this message in context: http://www.nabble.com/trying-to-install-plotrix-tp25530666p25530766.html Sent from the R help mailing lis

[R] trying to install plotrix

2009-09-22 Thread David Epstein
Here are the error messages: > install.packages('plotrix') Warning in install.packages("plotrix") : argument 'lib' is missing: using '/Users/dbae/Library/R/library' --- Please select a CRAN mirror for use in this session --- Error in .install.macbinary(pkgs = pkgs, lib = lib, contriburl = contri

Re: [R] eliminating control characters from formatted data files

2009-02-05 Thread David Epstein
Murray Cooper wrote: > > This may be a case of "If all you have is a hammer, everything looks like > a > nail". > If all you want to do is remove the last line if it contains a CONTROL-Z, > why > not use something like perl to process the files? > My first thought was to use perl, and this w

[R] eliminating control characters from formatted data files

2009-02-05 Thread David Epstein
I have a few hundred files of formatted data. Unfortunately most of them end with a spurious CONTROL-Z. I want to rewrite the files without the spurious character. Here's what I've come up with so far, but my code is unsafe because it assumes without justification that the last row of df contains

[R] Difference between a[[i]] and a[i]

2009-02-02 Thread David Epstein
I'm sure I've read about the difference between a[[i]] and a[i] in R, but I cannot recall what I read. Even more disturbing is the fact that I don't know how to search the newsgroup for this. All the different combinations I tried were declared not to be valid search syntax. 1. What sort of objec

[R] Environmental variables

2009-02-01 Thread David Epstein
I use a Mac (10.4.11 Mac Os X). In my .tcshrc I define an environmental variable MY. Is it possible to find out its value from inside R? When one loads R for Mac OS X Cocoa GUI written by: Simon Urbanek Stefano M. Iacus are files like .tcshrc read by R? Can I make the value of th

Re: [R] simultaneous plots

2008-12-03 Thread David Epstein
Prof Brian Ripley wrote: > > How about dev.copy2pdf ? > I put the following into a script: quartz(...) #many graphics commands dev.copy2pdf(pdf,file="ratio1.pdf") and then sourced the script. I got the error message Error in pdfFonts(family) : invalid arguments in 'pdfFonts' (must be

[R] Making a picture that is wide and small height

2008-12-03 Thread David Epstein
How do I make a picture that is a horizontal strip? I tried > plot(x=c(1,2,3,4),y=c(1,1,1,1)) #works but screen image is square. > pdf("ratio.pdf",height=1,width=6) > plot(x=c(1,2,3,4),y=c(1,1,1,1)) I got the following error message: Error in plot.new() : figure margins too large Is it possible

[R] simultaneous plots

2008-12-02 Thread David Epstein
Is there a good and concise way of making simultaneous plots that are identical, but directed to different devices? I'm writing an R-script that produces a pdf file. I would really like to check visually whether the pdf file shows what I expect. So I would like the same commands to produce a plot

[R] library(gdata) fails

2008-11-05 Thread David Epstein
I wanted to install the package "gdata". Here are the commands I gave and the responses: install.packages("gdata") Warning in install.packages("gdata") : argument 'lib' is missing: using '/Users/dbae/Library/R/library' trying URL 'http://cran.uk.r-project.org/bin/macosx/universal/contrib/2.8/gd

[R] re flecting a line

2008-09-10 Thread David Epstein
Suppose x and y are numeric vectors of the same length. plot(x,y) #scatterplot lmObj1 <- lm(y~x) # best fit line abline(lmObj1) # good lmObj2 <- lm(x~y) #get best fit but with axes interchanged abline(lmObj2) # not what I want. I want the correct line, drawn on the same graph, but with

[R] options("contrasts")

2008-08-26 Thread David Epstein
ingle item in a list?" My answer to the meta-meta-question is to post to this list. I hope that at least that part is correct. Thanks for any help. David Epstein -- View this message in context: http://www.nabble.com/options%28%22contrasts%22%29-tp

Re: [R] asp and ylim

2008-07-21 Thread David Epstein
A really great answer to my concerns! I'll get hold of the Paul Murrell book, and see how far I can get. On 21 Jul, 2008, at 10:48, Martin Maechler wrote: Play around resizing your graphics window.. This is very instructive, with an 'asp = .' using "traditional graphics" plot(). OK, but I don

Re: [R] asp and ylim

2008-07-21 Thread David Epstein
elp. David On 20 Jul, 2008, at 23:14, Prof Brian Ripley wrote: Take a look at eqscplot() in package MASS for a different approach. You last para forgets that once you have set the device region and the margins the physical plot region and hence its aspect ratio is determined -- see the figur

Re: [R] asp and ylim

2008-07-21 Thread David Epstein
ion and the margins the physical plot region and hence its aspect ratio is determined -- see the figures in 'An Introduction to R'. On Sun, 20 Jul 2008, David Epstein wrote: #See David Williams' book "Weighing the odds", p286 y <- c(1.21, 0.51, 0.14, 1.62, -

Re: [R] asp and ylim

2008-07-20 Thread David Epstein
On 20 Jul, 2008, at 22:00, Rolf Turner wrote: I just tried > set.seed(42) > x <- runif(10) > y <- runif(10) > plot(x,y,ylim=c(-1,2),asp=1) and it seemed to give results as expected/desired. cheers, Ro

[R] asp and ylim

2008-07-20 Thread David Epstein
#See David Williams' book "Weighing the odds", p286 y <- c(1.21, 0.51, 0.14, 1.62, -0.8, 0.72, -1.71, 0.84, 0.02, -0.12) ybar <- mean(y) ylength <- length(y) ybarv <- rep(ybar, ylength) x <- 1:ylength plot(x,y,asp=1,xlab="position",ylab="ybar",type="n",ylim=c(-1,1)) segments(x[1], ybar, x

Re: [R] drawing segments through points with pch=1

2008-07-20 Thread David Epstein
t;n") segments(x[1], ybar, x[ll], ybar) points(x, ybarv, pch=21, bg="white") Uwe Ligges David Epstein wrote: Please excuse me for asking such basic questions: Here is my code y=c(1.21,0.51,0.14,1.62,-0.8,0.72,-1.71,0.84,0.02,-0.12) ybar=mean(y) ll=length(y); ybarv=rep(ybar,ll) x=1:

[R] drawing segments through points with pch=1

2008-07-20 Thread David Epstein
hing inside them? I tried putting the segments command first, but then no segment appeared at all. In general, is there a method of laying a drawing "on top" of another. I tried inserting add="T" as an argument to plot, and R objected s

Re: [R] Error in edit(name,file,title,editor)

2008-07-20 Thread David Epstein
list.files()' to see what is addressable in whatever directory you are in. On Sun, Jul 20, 2008 at 7:19 AM, David Epstein <[EMAIL PROTECTED]> wrote: Can anyone help me with the following attempt to use an external editor from within R vi(file="p286.R") Error in edit(name, file

[R] Error in edit(name,file,title,editor)

2008-07-20 Thread David Epstein
pose I could use the function source(), but I haven't tried this. Thanks for any help David Epstein __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/p