search for: biostatmatt

...small patch (12 lines w/o documentation) that allows R's iconv to operate on raw objects, while not interfering or affecting the behavior of iconv on character vectors. To keep this message terse, I've put additional discussion, description of what the patch does, and examples here: http://biostatmatt.com/archives/456 Also, here is a link to the patch file: http://biostatmatt.com/R/R-devel-iconv-0.0.patch If this change is adopted, I'd be happy to submit a documentation patch also. -Matt Index: src/library/base/R/New-Internal.R ============================================================...

...l features could be reproduced using interpreted R code with little effort. So what do you think? Stay portable and expand to Windows? Exploit the full POSIX termios functionality? Expose connections API and write a package instead? Forget about the whole thing? -Matt Current patch URLs: http://biostatmatt.com/R/R-2.11.0-tty.patch http://biostatmatt.com/R/R-2.11.0-tty+conf.patch http://biostatmatt.com/archives/112 `diffstat -f0 R-2.11.0-tty.patch` : configure.ac | 2 1 + 1 - src/include/Internal.h | 2 2 + 0 - src/library/base/R/connections.R | 11 11...

Hello, I spent quite some time now searching for any hint that R can also be used to address the interfaces of a computer (i.e. RS232 or USB) to acquire data from measurement devices (like with the - I think it is the - devices or serial toolbox of Matlab). Is there any package available or a project going on that you know of? I would so much like to have never to work with Matlab again. The

Dear Sirs, I am a new user of the R package. When I try to use the curve function it confuses me. > curve(x^2) Works fine. > curve(x) Makes a complaint I don't understand. Why is x^2 valid and x is not? I check the documentation of curve, and it says the first argument must be an expression containing x. > expression(x) Is an expression containing x. > curve(expression(x))

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 found that the code only began to display the CPU time after quite a while and several iterations have finished. Is there a way to ask R to display time right after each iteration is finished? Thanks, -Jun [[alternative HTML version

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 *rank 4* but I want 5 rank 1 7 rank 2 7 rank 2 9 *rank 3* * * thanks SS [[alternative HTML version deleted]]

Dear R-users, Is there an R function to compute the multinomial beta function? That is, the normalizing constant that arises in a Dirichlet distribution. For example, with three parameters the beta function is Beta(n1,n2,n2) = Gamma(n1)*Gamma(n2)*Gamma(n3)/Gamma(n1+n2+n3) Thanks in advance for any assisstance. Regards, Greg [[alternative HTML version deleted]]

Hello list, Is there any way to bandpass filter in R thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version deleted]]

...lab('Drug') + ylab('Hazard ratio')+ > geom_hline(y=1, col='red', lty=2) > > > Abhijit > On Oct 30, 2010, at 5:31 PM, Matt Shotwell wrote: > >> 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&lt...

Dear R Developers: Is there a way to look at the underlying code from such items as R_setup_starma, please? Thanks, Erin Erin M. Hodgess, PhD Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: hodgesse@uhd.edu [[alternative HTML version deleted]]

...ml#x1-11000A.9 [2] https://stat.ethz.ch/pipermail/r-help/2008-June/164783.html [3] http://s000.tinyupload.com/index.php?file_id=22646520644600938098 [4] http://s000.tinyupload.com/index.php?file_id=73594116524789598605 [5] http://s000.tinyupload.com/index.php?file_id=00010585455841885948 [6] http://biostatmatt.com/archives/573/comment-page-1#comment-770 > > Hope that helps. If you have a specific application where you need > > this type of functionality, I may be able so suggest a workaround. > > > > -Cameron > >

Hi, I am trying to call a C function, that I wrote to parse a flat file, into R. The argument that will go into this function is an input file that I need to parse and write the desired output in an output file. I used some hit and trial approach but i keep on getting the "file not found" or "segmentation fault" error. I know that the error is in passing the argument but I

hi all - i'm working on an R package that makes use of my own shared library written in C. but i also am making use of another C-written library. (my package is for facilitating biological namespace translations via online (i.e. up-to-date) biological databases.) problem is, the library i'm using is not a standard library (i.e. i doubt it will be installed on most users' machines). i

Hello R users, I am trying to use R to do the low pass filter analysis for the tidal data. I am a novice in R and so far been doing only simple stuffs on R. I found a package called signal but couldn't find the proper tutorial for the low pass filter. Could anyone point me to the proper tutorial or starting point on how to do low pass filter analysis in R ? Thank you so much. Janesh

Hello UseRs, I need help on how to plot several raster images (such as those obtained from a kernel-smoothed intensity function) in a layout such as that obtained from the lattice package. I would like to obtain something such as obtained from using the "levelplot" or "xyplot" in lattice. I currently use: >par(mfrow=c(3,3) to set the workspace, but the resulting plots

I appreciate that this is OT, but I'd be grateful for pointers to examples of where Sweave has been used for web-based applications. In particular, examples of where reports/analyses are produced automatically through submission of data to a web-sever. I am mostly interested in situations where pdf reports have been produced rather than, say, a plot/table etc shown on a web page.

iterations <- 100 nvars <- 4 combined <- rbind(scaleMiceTrain, scaleMiceTest) reducedSample <- combined reducedSample <- subset(reducedSample, select = -pID50) reducedSample <- subset(reducedSample, select = -id) for (i in 1:iterations) { miceSample <- sample(combined[,-c(1,2)],nvars, replace=FALSE) miceSample$pID50 <- combined$pID50 miceTestSample <-

Hello, I would like to compress a long string (character vector), store the compressed string in the text field of a SQLite database (using RSQLite), and then load the text back into memory and decompress it back into the the original string. My character vector can be compressed considerably using standard gzip/bzip2 compression. In theory it should be much faster for me to compress/decompress

I have installed both 32-bit and 64-bit versions of R2.12.0 (2010-06-15 r52300) on my Ubuntu 10.04 64-bit system. I observe the following behavior when running the examples from base::connections. There appears to be a problem with seek() on a .gz file when using a 32-bit installation of R2.12.0, but the problem doesn't appear in the 64-bit installation. I realize that seek() has been

Hello, While assisting a fellow R-helper off list, I narrowed down an issue he was having to the following behavior of type.convert, called through read.table. This is using R 2.10.1, if newer versions don't exhibit this behavior, apologies. # generates numeric vector > type.convert(c("123.42", "NA")) [1] 123.42 NA # generates a numeric vector, notice the