Re: [R] Opening SAS file using read.sas7bdat() function in sas7bdat library.

Thanks for the helpful comments from others. The KNOWNHOST variable lists the types of file that are known to work with the read.sas7bdat function. It's likely that most files written on Windows platforms will work, even if not listed in KNOWNHOST. If you're feeling experimental, you might just co

Re: [R] Any good R server-with connection examples

> I want to connect R with HTML/PHP pages to take input from user,do > some > statistical processing on it & show results to HTML page again. > I search on net,i got Rserve package,but examples are mainly for java > langaure & not for PHP > i am wondering how to connect it to PHP-Apache-MySQL > I

Re: [R] Nested brew call yields Error in .brew.cat(26, 28) : unused argument(s) (26, 28)

On Wed, 2012-03-28 at 11:40 +0100, Chris Beeley wrote: > I am writing several webpages using the brew package and R2HTML. I would > like to work off one script so I am using nested brew calls. The > documentation for brew states that: > > "NOTE: brew calls can be nested and rely on placing a fun

[R] Help with Matrix code optimization

The chol and solve methods for dpoMatrix (Matrix package) are much faster than the default methods. But, the time required to coerce a regular matrix to dpoMatrix swamps the advantage. Hence, I have the following problem, where use of dpoMatrix is worse than a regular matrix. library(Matrix) x <

Re: [R] function restrictedparts

That's because the number of partitions of 281 items of order 10 is quite large: R> library('partitions') R> R(10,281) [1] 1218681472 Without thinking about this too hard, the result of restrictedparts(281,10) should require around R> 1218681472 * 10 * 4 / 10^9 [1] 48.74726 gigabytes of storage

Re: [R] Bayesian data analysis recommendations

On Thu, 2012-01-19 at 19:23 -0500, C W wrote: > Thanks, Rich, I will look at the book. > > I agree, there are many nice packages, but what if the package changes in a > few years? I would have no idea what is going on! I've heard > from predecessor in the industry who emphasize the learning, not

Re: [R] R logo in eps formt

See this earlier post for SVG logos: http://tolstoy.newcastle.edu.au/R/e12/devel/10/10/0112.html Using Image Magick, do something like convert logo.svg logo.eps On Thu, 2011-12-01 at 10:56 +0700, Ben Madin wrote: > G'day all, > > Sorry if this message has been posted before, but searching fo

Re: [R] contact person for UseR 2012, please?

The contact person is: Stephania McNeal-Goddard email: stephania.mcneal-godd...@vanderbilt.edu phone: (615)322-2768 Vanderbilt University School of Medicine Department of Biostatistics S-2323 Medical Center North Nashville, TN 37232-2158 On Tue, 2011-10-18 at 12:41 -0400, David Winsemius wrote:

Re: [R] [Related Topic] need help on read.spss

Would it be worthwhile to update the read.spss implementation using the more recent discoveries from the PSPP group? I don't mean to copy their code; but to use the ideas in their code. Is anyone working on this? I wouldn't want the effort to be duplicated. On Thu, 2011-10-13 at 16:22 +0200, Uwe L

Re: [R] Rweb and setting up R on a server

Erin, I haven't used Rweb recently. The URL is http://www.math.montana.edu/Rweb/ . If you have a server, you could set up the server version of RStudio: http://rstudio.org/download/server . It worked well when I tried it. Best, Matt On Tue, 2011-09-06 at 17:07 -0500, Erin Hodgess wrote: > Dea

Re: [R] readBin fails to read large files

On Thu, 2011-09-01 at 17:36 +0100, Prof Brian Ripley wrote: > readBin is intended to read a few items at a time, not 10^9. You are > probably getting 32-bit integer overflow inside your OS, since the > number of bytes you are trying to read in one go exceeds 2GB. > > Don't do that: read say a m

Re: [R] Bootstrap

In order to apply the bootstrap, you must resample, uniformly at random from the independent units of measurement in your data. Assuming that these represent the rows of 'data', consider the following: est <- function(y, x, obeta = c(1,1), verbose=FALSE) { n <- length(x) X <- cbind(rep(1,

Re: [R] How to capture console output in a numeric format

On Fri, 2011-06-24 at 12:09 -0400, David Winsemius wrote: > On Jun 24, 2011, at 11:27 AM, Matt Shotwell wrote: > > > Ravi, > > > > Consider using an environment (i.e. a 'reference' object) to store the > > results, avoiding string manipulation, and

Re: [R] How to capture console output in a numeric format

Ravi, Consider using an environment (i.e. a 'reference' object) to store the results, avoiding string manipulation, and the potential for loss of precision: fr <- function(x, env) { ## Rosenbrock Banana function x1 <- x[1] x2 <- x[2] f <- 100 * (x2 - x1 * x1)^2 + (1 - x1)^2 if(

Re: [R] Elbow criterion

On Mon, 2011-06-20 at 13:38 +0200, Dominik P.H. Kalisch wrote: > Hi, > > I would like to cluster a dataset with the ward algorithm. I'm assuming that this refers to the agglomerative partitioning method [1]. That is, the number of clusters is selected according to the data partition that is seque

Re: [R] Can we prepare a questionaire in R

As Mike had written, there are frameworks for web-development with R. RApache http://www.rapache.net is one. Also, see the R package Rook: http://cran.r-project.org/web/packages/Rook/index.html . On Wed, 2011-06-08 at 17:26 +0530, amrita gs wrote: > How can we create HTML forms in R Wouldn't

Re: [R] Question about curve function

On Tue, 2011-06-07 at 16:17 +0200, Uwe Ligges wrote: > > On 07.06.2011 11:57, peter dalgaard wrote: > > > > On Jun 6, 2011, at 11:22 , Prof Brian Ripley wrote: > > > >> As a further example of the trickiness, the "function" method of plot() > >> relies on curve(x, ...) being a request to plot the

Re: [R] Question about curve function

I think there is trouble because expr in curve(expr) may be the name of a function, and it's ambiguous whether 'x' should be interpreted as a mathematical expression involving x, or the name of a function. Here are some examples that work: curve(I(x)) curve(1*x) On Sun, 2011-06-05 at 12:07 -0500,

Re: [R] creating a vector from a file

On Tue, 2011-05-31 at 16:19 +0200, heimat los wrote: > On Tue, May 31, 2011 at 4:12 PM, Matt Shotwell > wrote: > On Tue, 2011-05-31 at 15:36 +0200, heimat los wrote: > > Hello all, > > I am new to R and my question should be trivial. I need to &g

Re: [R] creating a vector from a file

On Tue, 2011-05-31 at 15:36 +0200, heimat los wrote: > Hello all, > I am new to R and my question should be trivial. I need to create a word > cloud from a txt file containing the words and their occurrence number. For > that purposes I am using the snippets package [1]. > As it can be seen at the

Re: [R] blank space escape sequence in R?

adLines), then the six characters would be read into memory, and printed like this: "a\\x20b" That is, not with a space character substituted for \x20. So, now I'm not sure this is a solution. On Mon, 2011-04-25 at 12:24 -0500, Matt Shotwell wrote: > You can embed hex escapes in s

Re: [R] blank space escape sequence in R?

You can embed hex escapes in strings (except \x00). The value(s) that you embed will depend on the character encoding used on you platform. If this is UTF-8, or some other ASCII compatible encoding, \x20 will work: > "foo\x20bar" [1] "foo bar" > For other locales, you might try charToRaw(" ") to

Re: [R] Converting 16-bit to 8-bit encoding?

s (I think) is to be read in and resaved in the less memory-efficient 16-bit format. Thanks! Sorry if the explanation isn't clear. On Thu, Apr 21, 2011 at 9:50 AM, Matt Shotwell mailto:matt.shotw...@vanderbilt.edu>> wrote: On 04/21/2011 10:36 AM, Brian Buma wrote: Hel

Re: [R] Converting 16-bit to 8-bit encoding?

On 04/21/2011 10:36 AM, Brian Buma wrote: Hello all- I have a question related to encoding. I'm using a seperate program which takes either 16 bit or 8 bit (flat binary files) as inputs (they are raster satellite imagery and the associated quality files), but can't handle both at the same time.

Re: [R] print.raw - but convert ASCII?

On Tue, 2011-04-19 at 03:14 -0400, Duncan Murdoch wrote: > On 11-04-18 9:51 PM, Matt Shotwell wrote: > > Does anyone know if there is a simple way to print raw vectors, such > > that ASCII characters are printed for bytes in the ASCII range, and > > their hex representation

[R] print.raw - but convert ASCII?

Does anyone know if there is a simple way to print raw vectors, such that ASCII characters are printed for bytes in the ASCII range, and their hex representation otherwise? rawToChar doesn't work when we have something like c(0x00, 0x00, 0x44, 0x00). -Matt

Re: [R] integer and floating-point storage

Hi Mike, There are some facilities for storing and manipulating small (2 bit) integers. See here: http://cran.r-project.org/web/packages/ff/index.html -Matt On 04/14/2011 01:20 PM, Mike Miller wrote: I note that "current implementations of R use 32-bit integers for integer vectors," but I a

[R] understanding dump.frames; typo;

When a function I have stop()s, I'd like it to return its evaluation frame, but not halt execution of the script. In experimenting with this, I became confused with dump.frames. From ?dump.frames: If ‘dump.frames’ is installed as the error handler, execution will continue even in non-

Re: [R] Examples of web-based Sweave use?

That's an interesting idea. I had written a long email describing a proof-of-concept, but decided to post is to the website below instead. http://biostatmatt.com/archives/1184 Matt On 04/04/2011 07:31 AM, carslaw wrote: I appreciate that this is OT, but I'd be grateful for pointers to example

Re: [R] library(foreign) read.spss warning

There is some information about this subtype in the PSPP source code, and for other subtypes not yet implemented by read.spss. The PSPP source code indicates that this subtype consists of "Value labels for long strings", which isn't very illuminating to me (probably because I don't use PSPP, or SPS

Re: [R] Venn Diagram corresponding to size in R

Try here: https://stat.ethz.ch/pipermail/r-help/2003-February/029393.html On Tue, 2011-03-08 at 20:25 -0500, Shira Rockowitz wrote: > I was wondering if anyone could help me figure out how to make a Venn > diagram in R where the circles are scaled to the size of each dataset. I > have looked at

Re: [R] assignment by value or reference

On 03/08/2011 07:20 AM, Xiaobo Gu wrote: > On Wed, Sep 15, 2010 at 5:05 PM, Uwe Ligges > wrote: >> See the R Language Definition manual. Since R knows about lazy evaluation, >> it is sometimes neither by reference nor by value. >> If you want to think binary, then "by value" fits better than "by

Re: [R] Rapache ( was Developing a web crawler )

On Sun, 2011-03-06 at 08:06 -0500, Mike Marchywka wrote: > > > > > > > > Date: Thu, 3 Mar 2011 13:04:11 -0600 > > From: matt.shotw...@vanderbilt.edu > > To: r-help@r-project.org > > Subject: Re: [R] Developing a web crawler / R "webkit" or something >

Re: [R] Developing a web crawler / R "webkit" or something similar? [off topic]

On 03/03/2011 08:07 AM, Mike Marchywka wrote: Date: Thu, 3 Mar 2011 01:22:44 -0800 From: antuj...@gmail.com To: r-help@r-project.org Subject: [R] Developing a web crawler Hi, I wish to develop a web crawler in R. I have been using the functionalities available under the RCurl package. I

Re: [R] Robust variance estimation with rq (failure of the bootstrap?)

hich uses a local estimate of the sparsity, to fail? > > -- > Jim > > > > On Mon, Feb 28, 2011 at 8:59 PM, Matt Shotwell wrote: > > Jim, > > > > If repeated measurements on patients are correlated, then resampling all > > measurements independently ind

Re: [R] Robust variance estimation with rq (failure of the bootstrap?)

Jim, If repeated measurements on patients are correlated, then resampling all measurements independently induces an incorrect sampling distribution (=> incorrect variance) on a statistic of these data. One solution, as you mention, is the block or cluster bootstrap, which preserves the correlatio

Re: [R] Visualizing Points on a Sphere

That's interesting. You might also like: http://en.wikipedia.org/wiki/Von_Mises%E2%80%93Fisher_distribution I'm not sure how to plot the wireframe sphere, but you can visualize the points by transforming to Cartesian coordinates like so: u <- runif(1000,0,1) v <- runif(1000,0,1) theta <- 2 * pi *

Re: [R] Problem with writing a file in UTF-8

Thomas, I wasn't able to reproduce your finding. The last two characters in my 'out.txt' file were just as expected. But, I'm in an UTF-8 locale. Your locale affects the encoding of characters on your platform. If you're not in a UTF-8 locale, then characters are converted from your native encodi

Re: [R] non-ascii characters in R output

On Fri, 2011-02-18 at 19:50 -0500, Duncan Murdoch wrote: > On 18/02/2011 5:58 PM, Matt Shotwell wrote: > > OK, looks like my web browser does render non-ascii characters output by > > R when it's given the encoding explicitly. This works for me: > http-equiv="Cont

Re: [R] non-ascii characters in R output

OK, looks like my web browser does render non-ascii characters output by R when it's given the encoding explicitly. This works for me: . So that's another solution, but not a general one. -Matt On Fri, 2011-02-18 at 12:47 -0600, Matt Shotwell wrote: > All, > > I'd like

[R] non-ascii characters in R output

All, I'd like to automatically output text from R to HTML. In doing this I've run into trouble with non-ascii characters, as my browser (and presumably others) does not render such characters correctly. For example, the 'fancy' single quotes associated with summary.lm are multi-byte character

Re: [R] Revolution Analytics reading SAS datasets

On Thu, 2011-02-10 at 10:44 -0800, David Smith wrote: > The SAS import/export feature of Revolution R Enterprise 4.2 isn't > open-source, so we can't release it in open-source Revolution R > Community, or to CRAN as we do with the ParallelR packages (foreach, > doMC, etc.). Judging by the language

Re: [R] clustering with finite mixture model

There are quite a few packages that work with finite mixtures, as evidenced by the descriptions here: http://cran.r-project.org/web/packages/index.html These might be useful: http://cran.r-project.org/web/packages/flexmix/index.html http://cran.r-project.org/web/packages/mclust/index.html -Ma

Re: [R] User input in R program

Martyn Plummer's 'coda' package has some nice interactive menus. The package appears to be written entirely in R. You could start with the codamenu() function in the package source: http://cran.r-project.org/web/packages/coda/index.html -Matt On Fri, 2011-01-21 at 14:26 +0200, christiaan pauw wr

Re: [R] Encoding problem - I fails to read Hebrew text from online

- > Contact me: tal.gal...@gmail.com | 972-52-7275845 > Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) > | www.r-statistics.com (English) > ------ &g

Re: [R] Encoding problem - I fails to read Hebrew text from online

Tal, It looks like the data you received has HTML special hex characters. That is, 'ש' is just an ASCII HTML representation of a hex character. It's not encoded in a special manner. The trick is to substitute the HTML encoded hex character for its binary representation, or "decode" the character

Re: [R] statistical test for comparison of two classifications (nominal)

will also have a look at profdpm, that did not pop-up in my searches. > Indeed, the interpretation is going to be critical... Could you please > elaborate on what you mean by the bootstrap process? > > Thanks a lot for your helps, > Martin > > On 11/17/2010 3:50 PM, Matt Shotwel

Re: [R] statistical test for comparison of two classifications (nominal)

There are several statistics used to compare nominal classifications, or _partitions_ of a data set. A partition isn't quite the same in this context because partitioned data are not restricted to a fixed number of classes. However, the statistics used to compare partitions should also work for the

Re: [R] Fatal Error R

Please see below. On Wed, 2010-11-17 at 04:41 -0500, Ted Harding wrote: > On 17-Nov-10 00:02:39, José Fernando Zea Castro wrote: > > Hello. > > First, I'm thankful about your wonderful project. > > > > However, I have serious worries about the reliability of R. > > I found the next bug which I co

Re: [R] log-transformed linear regression

Servet, These data do look linear in log space. Fortunately, the model log(y) = a + b * log(x) does have intercept zero in linear space. To see this, consider log(y) = a + b * log(x) y = 10^(a + b * log(x)) y = 10^a * 10^(b * log(x)) y = 10^a * 10^(log(x^b)) y = 10^a * x^b

Re: [R] data acquisition with R?

R implements (almost) all IO through its 'connections'. Unfortunately, there is no API (public or private) for adding connections, and therefore no packages that implement connections. You will find more discussion of connections and hardware (serial, USB) interface in the R-devel list archives. T

Re: [R] Lattice plots for images

Have you tried using the 'mai' argument to par()? Something like: par(mfrow=c(3,3), mai=c(0,0,0,0)) I've used this in conjunction with image() to plot raster data in a tight grid. -Matt On Wed, 2010-11-03 at 11:13 -0400, Neba Funwi-Gabga wrote: > Hello UseRs

Re: [R] ForestPlot or similar

Here is a small function for forest plots in R, with an example: http://biostatmatt.com/wiki/r-credplot -Matt On Sat, 2010-10-30 at 11:40 -0400, Mestat wrote: > Here is one example: > I have three vectors (mean,lower interval, upper interval) > mean<-c(2,4,6,8) > l<-c(1,2,3,4) > u<-c(4,8,12,16)

Re: [R] Regular expression to find value between brackets

Here's a shorter (but more cryptic) one: > gsub("^([^\\(]+)(\\((.+)\\))?", "\\2", tests) [1] """(%)" "(%)" "(mg/ml)" > gsub("^([^\\(]+)(\\((.+)\\))?", "\\3", tests) [1] "" "%" "%" "mg/ml" -Matt On Wed, 2010-10-13 at 14:34 -0400, Henrique Dallazuanna wrote: > Try th

Re: [R] puzzle with integrate over infinite range

You could try pnorm also: shiftedGaussR <- function(x0 = 500) { sd <- 100/sqrt(2) int <- pnorm(0, x0, sd, lower.tail=FALSE, log.p=TRUE) exp(int + log(sd) + 0.5 * log(2*pi)) } > shiftedGaussR(500) [1] 177.2454 > shiftedGauss(500) [1] 177.2454 -Matt On Tue, 2010-09-21 at 09:38 -0400

Re: [R] Is there a bisection method in R?

I was just reading about the merge sort algorithm last night (BTW, here is a fun link http://www.youtube.com/watch?v=t8g-iYGHpEA). There are some interesting similarities in this context. Here's a recursive method for bisection: bisectMatt <- function(fn, lo, hi, tol = 1e-7, ...) { flo <- fn(

Re: [R] R Founding

On Thu, 2010-09-16 at 17:30 -0400, Tal Galili wrote: > Hello dear Jaroslaw, > I strongly agree with you that the R foundation should have an easier method > of enabling people to give donations. > At the same time, I feel there is a (friendly) disagreement between us on > how such money should be u

Re: [R] Reproducible research

Well, the attachment was a dud. Try this: http://biostatmatt.com/R/markup_0.0.tar.gz -Matt On Thu, 2010-09-09 at 10:54 -0400, Matt Shotwell wrote: > I have a little package I've been using to write template blog posts (in > HTML) with embedded R code. It's quite small but v

Re: [R] Reproducible research

I have a little package I've been using to write template blog posts (in HTML) with embedded R code. It's quite small but very flexible and extensible, and aims to do something similar to Sweave and brew. In fact, the package is heavily influenced by the brew package, though implemented quite diffe

Re: [R] Which language is faster for numerical computation?

For the compiled languages, it depends heavily on the compiler. This sort of comparison is rendered moot by the huge variety of compiler and hardware specific optimizations. My suggestion is to use C, or possibly C++ in conjunction with Rcpp, as these are most compatible with R. Also, C and C++ are

Re: [R] Uncompressing data from read.socket

Have a look at gzcon, for decompressing data as they arrive. From the help file: ‘gzcon’ provides a modified connection that wraps an existing connection, and decompresses reads or compresses writes through that connection. Standard ‘gzip’ headers are assumed. There is no indicati

Re: [R] remove accents in strings

quot;E" "I" "O" "U" "a" "e" "i" "o" "u" "y" "A" "E" "I" [20] "O" "U" "Y" By the way, I'll take this moment to remind anyone interested that

Re: [R] remove accents in strings

If you know the encoding of the string, or if its encoding is the current locale encoding, then you can use the iconv function to convert the string to ASCII. Something like: iconv(accented.string, to="ASCII//TRANSLIT") While 7-bit ASCII does not permit accented characters, extended (8-bit) ASCII

Re: [R] [semi-OT] Using fortune() in an email signature

Or using R & GNU tools: m...@max:~$ R -e "fortunes::fortune()" | gawk '/^[^>]/ {print}' It's not a question of trying variations, rather of following instructions. -- Brian D. Ripley (about using 'Writing R Extensions') R-help (January 2006) -Matt On Wed, 2010-09-01 at 16:49 -0400, Stu

Re: [R] documentation to upgrade R-package from 32 to 64bit

Try this: http://www.stats.ox.ac.uk/~ripley/Win64/W64porting.html -Matt On Wed, 2010-09-01 at 07:40 -0400, Hayes, Daniel wrote: > Dear all, > > I am working with the an R-package named GAMLSS > (www.gamlss.com) it is currently only functional under > the 32-bit version o

Re: [R] Why is vector assignment in R recreates the entire vector ?

Tal, For your first example, x is not duplicated in memory. If you compile R with --enable-memory-profiling, you have access to the tracemem() function, which will report whether x is duplicate()d: > x <- rep(1,100) > tracemem(x) [1] "<0x8f71c38>" > x[10] <- NA This does not result in duplicati

Re: [R] Sweave.sty

Here is one: http://svn.r-project.org/R/trunk/share/texmf/tex/latex/Sweave.sty -Matt On Tue, 2010-08-24 at 15:40 -0400, r.ookie wrote: > Does anyone know where I can download the latest version of Sweave.sty? I > have looked all over the site http://www.stat.umn.edu/~charlie/Sweave/ with > no

Re: [R] on abort error, always show call stack?

to STDERR > automatically on a program abort? > > /iaw > > On Sun, Aug 22, 2010 at 11:50 AM, Matt Shotwell wrote: > > On Sun, 2010-08-22 at 11:41 -0400, ivo welch wrote: > >> Dear R Wizards---is it possible to get R to show its current call > >> stack (sys.c

Re: [R] on abort error, always show call stack?

On Sun, 2010-08-22 at 11:41 -0400, ivo welch wrote: > Dear R Wizards---is it possible to get R to show its current call > stack (sys.calls()) upon an error abort? I don't use ESS for > execution, and it is often not obvious how to locate how I triggered > an error in an R internal function. Seein

Re: [R] Pass By Value Questions

On Thu, 2010-08-19 at 14:27 -0400, Duncan Murdoch wrote: > On 19/08/2010 12:57 PM, li...@jdadesign.net wrote: > > I understand R is a "Pass-By-Value" language. I have a few practical > > questions, however. > > > > I'm dealing with a "large" dataset (~1GB) and so my understanding of the > > nuances

Re: [R] Does R always insist on sending plot output to a file?

d terminal $ R > postscript(file="Rfifo") > plot(0) > dev.off() -Matt On Wed, 2010-08-18 at 23:21 -0400, Matt Shotwell wrote: > Donald, > > At least for the PDF device (I know you asked about png, but I believe > they are similar), the answer no. Ultimately, th

Re: [R] Does R always insist on sending plot output to a file?

Donald, At least for the PDF device (I know you asked about png, but I believe they are similar), the answer no. Ultimately, this device calls the standard C function fopen, and writes its data to the resulting file stream. If you're using GNU Linux, you might trick R into writing to a fifo (a na

Re: [R] Delete rpart/mvpart cross-validation output

Or, if using GNU Linux or other UNIX-like system: sink("/dev/null") # Issue commands sink() -Matt On Wed, 2010-08-18 at 09:14 -0400, Gabor Grothendieck wrote: > On Fri, Aug 13, 2010 at 1:52 PM, Marie-Hélène Ouellette > wrote: > > Dear all, > > > > I was wondering if there is a simple way to avo

Re: [R] when to use textConnection ??

Also, many R functions are designed to operate on R connections, to input and output text. Alternatively, we may wish to provide the input text as an R character vector, or output text to a character vector. The textConnection makes a character vector look like a connection, so R routines that oper

Re: [R] ASCI characters

How about: > rawToChar(as.raw(82)) [1] "R" -Matt On Sun, 2010-08-15 at 19:50 -0400, Orvalho Augusto wrote: > Hello guys! > > Is there any function that permits me to get an ASCI character from its > code? Eg. ascifunction(34) would give me ' > or ascifunction(92) gives \ > > Thanks > Caveman >

Re: [R] reading a text file, one line at a time

Walt, Something like: con <- file("your-large-file.txt", "rt") readLines(con, 1) # Read one line -Matt On Sun, 2010-08-15 at 10:58 -0400, Data Analytics Corp. wrote: > Hi, > > I have an upcoming project that will involve a large text file. I want to > >1. read the file into R one line

Re: [R] band pass filter

nuncio, If you already have a filter kernel, you can use the filter function. Of course, convolution filters can be applied directly using the discrete Fourier transform via the fft function. For an example of filtering (lowpass) with R, see http://biostatmatt.com/archives/78 , and the associated

Re: [R] Numerical Methods Course

TGS, Given that you have to pay an outrageous $155.86 for that book, it seems reasonable to look for a free environment for numerical computing (like R!). If your instructor says that such a variety of programming languages would work, you could probably make a good argument to use R. But why not

Re: [R] Changing downloaded source code into a package

See comments below. On Mon, 2010-08-09 at 10:22 -0400, JH wrote: > I am wanting to change some lines of code in the R package named "nlme" > http://cran.r-project.org/web/packages/nlme/index.html > To do this I have downloaded the Package source named nlme_3.1-96.tar.gz, > opened up the file and c

Re: [R] Good Book To Work Through This Summer

There are some book-length documents (downloadable for free) at the contributed documentation section of the R project website here: http://cran.r-project.org/other-docs.html In particular, the book “Practical Regression and Anova using R” by Julian Faraway looks to have the content you want, tho

Re: [R] grep with search terms defined by a variable

Daniel, If you want to search for each term at the beginning of a sting, using the regular expression construct '^', you might use the following > search.terms <- c("Emil", "Meryl") > names <- c("Emil Jannings", +"Charles Chaplin", +"Katherine Hepburn", +"M

Re: [R] Meaning of following function

Ron, In arithmetic, '-' and '+' are binary _and_ unary operators. That is, both -1 and 1-1 are valid arithmetic expressions, the former negates its argument, and the latter subtracts the second from the first. Since much of R is designed do arithmetic, R honors the unary _and_ binary versions of '

Re: [R] how to code it??

If I take your meaning correctly, you want something like this. > x <- c(0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, + 1) > easy <- function(x) { + state <- 0 + for (i in 1:length(x)) { + if (x[i] == 0) + x[i] <- state + state <- 0 + if (x

Re: [R] what is a vignette?

Alex, Vignettes are optional supplemental documentation. That is, they are in addition to the required boilerplate documentation for R functions and datasets. Vignettes are written in the spirit of sharing knowledge, and assisting new users in learning the purpose and use of a package. Maybe the

Re: [R] sink function

I had addressed a problem similar to this only a few days ago. Please see the following URL: http://tolstoy.newcastle.edu.au/R/e11/help/10/07/1677.html On Fri, 2010-07-23 at 08:45 -0400, nuncio m wrote: > I have the following code to write the output from auto.arima function. The > issue is not

Re: [R] Help with Sink Function

Your code between calls to sink() does not generate any output. Hence, nothing will be diverted to the file. To illustrate this point, consider for(i in 1:10) i This produces no output. However, for(i in 1:10) print(i) produces output as expected. -Matt On Fri, 2010-07-16 at 13:34 -0400, Ad

Re: [R] calling a c function from R

Fahim, Please see the Writing R Extensions manual http://cran.r-project.org/doc/manuals/R-exts.pdf There are simple instructions in this document under the heading "System and foreign language interfaces." -Matt On Wed, 2010-07-14 at 01:21 -0400, Fahim Md wrote: > Hi, > I am trying to call a

Re: [R] Fast string comparison

acter" ) > strcmp2 <- cxxfunction(sig, code2, settings = settings) > system.time(strcmp2( lhs, rhs) ) > > > > I get: > > $ Rscript strings.R > Le chargement a nécessité le package : methods > utilisateur système écoulé >0.002 0.000

Re: [R] Fast string comparison

On Tue, 2010-07-13 at 01:42 -0400, Hadley Wickham wrote: > strings <- replicate(1e5, paste(sample(letters, 100, rep = T), collapse = > "")) > system.time(strings[-1] == strings[-1e5]) > # user system elapsed > # 0.016 0.000 0.017 > > So it takes ~1/100 of a second to do ~100,000 string c

Re: [R] Compress string memCompress/Decompress

lieve it is stored as a BLOB regardless. -Matt > Thanks again!, > Erik > > > > On Jul 9, 2010, at 3:21 PM, Matt Shotwell wrote: > > > Erik, > > > > Can you store the data as a blob? For example: > > > >> #create string, compre

Re: [R] Compress string memCompress/Decompress

Erik, Can you store the data as a blob? For example: > #create string, compress with gzip, convert to SQLite blob string > string <- "gzip this string, store as blob in SQLite database" > string.gz <- memCompress(string, type="gzip") > string.sqlite <- paste("x'",paste(string.gz,collapse=""),"'"

Re: [R] Calling Gnuplot from R

I recently wrote a small R function to draw simple ASCII scatterplots. http://biostatmatt.com/archives/491 Bill Harris commented that the plots reminded him of the "dumb" terminal of Gnuplot. I think it would be really neat to have an R graphics driver to Gnuplot in order to generate more complete

Re: [R] Function to compute the multinomial beta function?

How about this? mbeta <- function(...) { exp(sum(lgamma(c(...)))-lgamma(sum(c(... } > gamma(5)*gamma(6)*gamma(7)/gamma(18) [1] 5.829838e-09 > mbeta(5,6,7) [1] 5.829838e-09 On Mon, 2010-07-05 at 17:10 -0400, Gregory Gentlemen wrote: > Dear R-users, > > Is there an R function to compu

Re: [R] Plot with whispers

It looks like read.table is reading the first line as a data value, which is the default for read.table. Try using read.table with the argument header=TRUE. Also, consider using a box and whiskers plot for these data (?boxplot, ?lattice::bwplot). -Matt On Mon, 2010-07-05 at 12:08 -0400, Ian Bentl

Re: [R] is there a way to do dense rank in R

> x <- c(5,7,7,9) > rank(unique(x))[match(x, unique(x))] [1] 1 2 2 3 On Thu, 2010-07-01 at 21:30 -0400, Suresh Singh wrote: > I have not been able to find a way to do dense rank in R > > Here is an example of what I need > > rank() gives the following > > 5 rank 1 > 7 rank 2 > 7 rank 2 > 9 *ran

Re: [R] how to display the clock time in the loop

Try to flush output after printing: cat(paste(Sys.time()),"\n"); flush(stdout()) On Thu, 2010-07-01 at 16:17 -0400, Jack Luo wrote: > Hi, > > I am doing some computation which is pretty time consuming, I want R to > display CPU time after each iteration using the command Sys.time(). However, > I

Re: [R] left end or right end

Suku, Just to clarify, in your table and each of your images, it appears that the start position of P (start1) is _after_ or at the start position of Q (start2), and the end position of P (end1) is _before_ or at the end position of Q (end2). If these positions represent increasing integers, then

Re: [R] left end or right end

Suku, It looks like you might want to consult with a [bio]statistician, but I'm interested in what these distances represent. Can you give some additional context for your problem? How were these distances collected? Is it a collection of pairs of intervals, like this: P

Re: [R] advice on package devel with external libs

Some ideas, 1. Wrap the library as an R package, as you said, and check for the library at configure time (i.e. with autoconf or custom script). But if you do, it would be great to provide an R-level API so that we can all use it. This is the strategy of the 'cairo', 'RGtk', 'rgl', and 'gsl' packa

Re: [R] integration of two normal density

Isn't it equally trivial to demonstrate that the product of two pdfs _may_ be a normalized pdf? For example, the uniform (0,1) pdf: f(x) = 1 for x in (0, 1), and 0 otherwise Hence, g(x) = f(x)*f(x) = 1 for x in (0, 1), and 0 otherwise _is_ a normalized pdf. But this is a little silly. Rather th

  1   2   >