[R] Function execution on package load?

2010-10-20 Thread Johannes Graumann
Hi, Can anyone give me a pointer on howto make a package execute a function at loading? Following an older post (http://bit.ly/cS1Go4), I'd like to do something along the lines of .localstuff - new.env() .localstuff$OftenUsedData - read.csv(...) upon loading the package ... Thanks, Joh

[R] Install Error

2010-10-28 Thread Johannes Graumann
Hi, I'm running into the error below when doing R CMD INSTALL MyPackage.tar.gz. This didn't use to be this way and I am at a loss as to where this might be coming from. Any pointers where to look? Joh ** building package indices ... Error in scan(file, what, nmax, sep, dec, quote, skip,

Re: [R] Install Error

2010-10-28 Thread Johannes Graumann
Duncan Murdoch wrote: On 28/10/2010 7:54 AM, Johannes Graumann wrote: Hi, I'm running into the error below when doing R CMD INSTALL MyPackage.tar.gz. This didn't use to be this way and I am at a loss as to where this might be coming from. Any pointers where to look? Joh ** building

[R] [SOLVED] Re: Install Error

2010-10-28 Thread Johannes Graumann
Johannes Graumann wrote: Hi, I'm running into the error below when doing R CMD INSTALL MyPackage.tar.gz. This didn't use to be this way and I am at a loss as to where this might be coming from. Any pointers where to look? Joh ** building package indices ... Error in scan(file, what

Re: [R] [SOLVED] Re: Install Error

2010-10-28 Thread Johannes Graumann
Johannes Graumann wrote: Johannes Graumann wrote: Hi, I'm running into the error below when doing R CMD INSTALL MyPackage.tar.gz. This didn't use to be this way and I am at a loss as to where this might be coming from. Any pointers where to look? Joh ** building package indices

[R] List to data frame

2010-07-26 Thread Johannes Graumann
Hi, Any ideas on how to efficiently convert list(c(1,2,3),c(4,5,6)) to data.frame(OriginalListIndex=c(1,1,1,2,2,2),Item=c(1,2,3,4,5,6)) Thanks for any hints, Joh __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] List to data frame

2010-07-26 Thread Johannes Graumann
(x = seq_along(mydata), times = unlist(lapply(mydata, length))), Item = unlist(mydata) ) HTH, Josh On Mon, Jul 26, 2010 at 7:46 AM, Johannes Graumann johannes_graum...@web.de wrote: Hi, Any ideas on how to efficiently convert list(c(1,2,3),c(4,5,6

[R] grid.table and expression in table body?

2010-08-03 Thread Johannes Graumann
Hi, Is there any way to get an expression into a data.frame, such that grid.table from gridExtra will plot it evaluated in the table body? The docu does it for the header, but is the body possible? Thanks, Joh __ R-help@r-project.org mailing list

Re: [R] grid.table and expression in table body?

2010-08-04 Thread Johannes Graumann
- data.frame(alpha, beta) grid.table(d, parse=T) you'll need revision 258 of gridExtra for this to work (googlecode now, r-forge in the following days, CRAN in the next stable version). HTH, baptiste On Aug 4, 2010, at 9:56 AM, Johannes Graumann wrote: Hi Baptiste, This is, I fear a bit

[R] Relation 1.5*IQR/Percentile in case of a normal Distribution

2010-08-14 Thread Johannes Graumann
Hi, can someone point me at material to understand how in http://upload.wikimedia.org/wikipedia/commons/8/89/Boxplot_vs_PDF.png the fivenum-corresponding percentages might be calculated? Thanks, Joh __ R-help@r-project.org mailing list

Re: [R] Relation 1.5*IQR/Percentile in case of a normal Distribution

2010-08-15 Thread Johannes Graumann
On Saturday 14 August 2010 23:08:31 Peter Dalgaard wrote: Johannes Graumann wrote: Hi, can someone point me at material to understand how in http://upload.wikimedia.org/wikipedia/commons/8/89/Boxplot_vs_PDF.png the fivenum-corresponding percentages might be calculated? Looks like

Re: [R] grid.table and expression in table body?

2010-08-17 Thread Johannes Graumann
. baptiste On 6 August 2010 17:11, Johannes Graumann johannes_graum...@web.de wrote: I updated the package from r-forge, but despite the fact that grid.table does not complain about the parse option if given, your example below is not being printed as a parsed expression. How can I check

[R] getNodeSet - what am I doing wrong?

2010-08-30 Thread Johannes Graumann
Hi, Why is the following retuning a nodset of length 0: library(XML) test - xmlTreeParse( http://www.unimod.org/xml/unimod_tables.xml,useInternalNodes=TRUE) getNodeSet(test,//modifications_row) Thanks for any hint. Joh __ R-help@r-project.org

Re: [R] getNodeSet - what am I doing wrong?

2010-08-31 Thread Johannes Graumann
() So getNodeSet(test, //x:modifications_row, x) gives you probably what you want. D. On 8/30/10 8:02 AM, Johannes Graumann wrote: library(XML) test - xmlTreeParse( http://www.unimod.org/xml/unimod_tables.xml,useInternalNodes=TRUE) getNodeSet(test,//modifications_row

Re: [R] getNodeSet - what am I doing wrong?

2010-08-31 Thread Johannes Graumann
Sorry about that - got dropped from my attempts yesterday (see the first example below, that has the useInternalNodes=TRUE) ... Thanks again, Joh Duncan Temple Lang wrote: Johannes Graumann wrote: Thanks! but: library(XML) xmlDoc - xmlTreeParse(http://www.unimod.org/xml

[R] choose.dir() gone?

2010-09-08 Thread Johannes Graumann
Hi, I fail to find choose.dir() in my current R install (see below)? Didn't that exist at some point? How to achieve file.choose() equivalent functionality for directories? Thanks for any hints, Joh sessionInfo() R version 2.11.1 (2010-05-31) x86_64-pc-linux-gnu locale: [1]

Re: [R] choose.dir() gone?

2010-09-08 Thread Johannes Graumann
OK. Just checked and choose.file/choose.dir exists in the windows version - apparently not in the linux one ... does anybody have a nice platform-agnostic solution for this? Thanks, Joh Johannes Graumann wrote: Hi, I fail to find choose.dir() in my current R install (see below)? Didn't

[R] Splitting a DF into rows according to a column

2010-10-04 Thread Johannes Graumann
Hi, I'm turning my wheels on this and keep coming around to the same wrong solution - please have a look and give a hand ... The premise is: a DF like so loremIpsum - Lorem ipsum dolor sit amet, consectetur adipiscing elit. Quisque leo ipsum, ultricies scelerisque volutpat non, volutpat et

Re: [R] Splitting a DF into rows according to a column

2010-10-05 Thread Johannes Graumann
Stupid Joh wants to give you a big hug! Thanks! Why rank works but order not, I have still to figure out, though ... Joh On Monday 04 October 2010 17:30:32 peter dalgaard wrote: On Oct 4, 2010, at 16:57 , Johannes Graumann wrote: Hi, I'm turning my wheels on this and keep coming around

Re: [R] loop and plot

2009-10-19 Thread Johannes Graumann
Rene wrote: Dear all, I am stuck at applying loop function for creating separated plots. I have coding like below: dataset.table - table(data.frame(var1=c(1,2,3,1,2,3,1),colour=c(a,b,c,c,a,b,b) )) kk = function(f) { ls=as.character(f)

[R] How to make XML support Expat?

2009-10-24 Thread Johannes Graumann
Hi, How can I make the result of the following lines TRUE? install.packages(XML) library(XML) supportsExpat() [1] FALSE I'm on linux, looked into the actual package, but don't seem to be able to wrap my head around how to compile this in ... Any pointers are welcome, Thanks Joh

Re: [R] How to make XML support Expat?

2009-10-25 Thread Johannes Graumann
the expat libraries, etc. HTH, D. Johannes Graumann wrote: Hi, How can I make the result of the following lines TRUE? install.packages(XML) library(XML) supportsExpat() [1] FALSE I'm on linux, looked into the actual package, but don't seem to be able to wrap my head

Re: [R] How to make XML support Expat?

2009-10-25 Thread Johannes Graumann
24 October 2009 23:31:46 Duncan Temple Lang wrote: Johannes Graumann wrote: Hi, I had heard that Expat is was faster. Your mail actually made me go check google for some comparisons and that does not seem the case ... do you have any insight into this? A couple of points.. i

Re: [R] How to make XML support Expat?

2009-10-25 Thread Johannes Graumann
On Sunday 25 October 2009 00:38:54 you wrote: xmlEventParse() is intended for handling files that we don't want to keep in memory. The branches parameter does make it easier to deal with sub-trees as the document is being parsed. And within these branches one can use XPath. Very

[R] Vector grouping challenge

2009-10-28 Thread Johannes Graumann
Dear all, Is there an efficient way to get this list testList - list(c(1),c(2,3,4),c(5,6),c(7),c(8),c(9,10,11,12),c(13)) from this vector testVector - c(12,32,NA,NA,56,NA,78,65,87,NA,NA,NA,90) ? Basically the vector should be grouped, such that non-NA and all following NAs end up in one

Re: [R] Vector grouping challenge

2009-10-28 Thread Johannes Graumann
` [1] 1 $`2` [1] 2 3 4 $`3` [1] 5 6 $`4` [1] 7 $`5` [1] 8 $`6` [1] 9 10 11 12 $`7` [1] 13 On Wed, Oct 28, 2009 at 7:57 AM, Johannes Graumann johannes_graum...@web.de wrote: Dear all, Is there an efficient way to get this list testList - list(c(1),c(2,3,4),c(5,6),c

[R] Regex matching that gives byte offset?

2009-10-28 Thread Johannes Graumann
Hi, Is there any way of doing 'grep' ore something like it on the content of a text file and extract the byte positioning of the match in the file? I'm facing the need to access rather largish (600MB) XML files and would like to be able to index them ... Thanks for any help or flogging, Joh

Re: [R] Regex matching that gives byte offset?

2009-11-02 Thread Johannes Graumann
the character offset: if you want to match in a MBCS and have byte offsets you will need to work a bit harder if useBytes=TRUE is not sufficient for you. On Wed, 28 Oct 2009, Johannes Graumann wrote: Hi, Is there any way of doing 'grep' ore something like it on the content of a text file

Re: [R] Regex matching that gives byte offset?

2009-11-02 Thread Johannes Graumann
On Monday 02 November 2009 13:41:45 Prof Brian Ripley wrote: On Mon, 2 Nov 2009, Johannes Graumann wrote: Hmmm ... that should do it, thanks. But how would one use this on a file without reading it into memory completely? ?file, ?readLines, ?readBin will tell you about connections

[R] Scanning grep through huge files

2009-11-03 Thread Johannes Graumann
Hi, I'm dealing which huge files I would like to index. On a linux system grep -buo PATTERN FILENAME hands me the byte offsets for PATTERN very quickly and I am looking to emulate that speed and ease with native R tools - for portability and elegance. gregexpr should be able to do that but I

[R] How to identify what is used as EOL in a given file?

2009-11-09 Thread Johannes Graumann
Hi, Is there any R-generic, OS-agnostic way to figure out what end-of-line character is being used in a file to be processed by readLines? Thanks, Joh __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Normal distribution test

2009-11-17 Thread Johannes Graumann
Markus Mehrwald wrote: Hi all, I am completely new to R and my knowledge of statistics is quite small so I hope you can help my. I have three dimensional point data which represents (and this is what I do not know for sure) a normal distribution. Now I want to test if this is true or not

Re: [R] symbol in the plot

2009-11-20 Thread Johannes Graumann
How about plot(sigma, delta1, ylim=range(-0.5, 2), xlab='sigma', ylab='delta', pch=22, type='o') points(sigma, delta2, col='red', axes=FALSE, type='o') legend(topleft,c(Delta1,Delta2),fill=TRUE,col=c(black,red)) Send runnable example next time. HTH, Joh gcheer3 wrote: TO be specific, here

Re: [R] symbol in the plot

2009-11-20 Thread Johannes Graumann
'. Also, it will be nice if there is a square symbol next to delta1 and a circle symbol next to delta 2, since sometimes I have to print the graph in a white and black paper. Thanks for any suggestions. Sorry for not asking question clearly. Johannes Graumann-2 wrote: How about plot(sigma

Re: [R] Method

2009-11-24 Thread Johannes Graumann
myVector - c(seq(10),23,35) length(myVector) myVector[length(myVector)] it's unclear to me which of the two you want ... HTH, Joh yonosoyelmejor wrote: Hello, i would like to ask you another question. Is exist anymethod to vectors that tells me the last element?That is to say,I have a

[R] Apply: Output matrix orientation

2010-05-27 Thread Johannes Graumann
Hi, Why is the result of below apply call rotated with respect to the input and how to remedy this? Thanks, Joh .ZScore - function(input){ #cat(input,\n) z - (input - mean(input))/sd(input) return(z) } apply(data.frame(x1=c(1,2,3,4,5),x2=c(2,3,4,5,6),x3=c(3,4,5,6,7)),1,.ZScore)

Re: [R] Apply: Output matrix orientation

2010-05-28 Thread Johannes Graumann
David Winsemius wrote: On May 27, 2010, at 7:24 AM, Johannes Graumann wrote: Hi, Why is the result of below apply call rotated with respect to the input and how to remedy this? Because the processing you requested is with respect to rows and the construction of matrices

[R] memDecompress and zlib compressed base64 encoded string

2010-01-14 Thread Johannes Graumann
Hi, I have zlib compressed strings (example is attached) and would like to decompress them using memDecompress ... I try this: connection - file(compressed.txt,r) compressed - readLines(connection) memDecompress(as.raw(compressed),type=g) Error in memDecompress(as.raw(compressed), type = g)

Re: [R] memDecompress and zlib compressed base64 encoded string

2010-01-15 Thread Johannes Graumann
was stripped, such files are 8-bit and yours is ASCII. Try x - 'Johannes Graumann' xx - charToRaw(x) xxx - memCompress(xx, g) rawToChar(xxx) [1] x\x9c\xf3\xca\xcfH\xcc\xcbK-Vp/J,\xcd\0052\001:\n\006\x90 to see what a real gzipped string looks like. and would like to decompress them using

[R] Lattice: How to implement varwidth analogous to graphics::boxplot in bwplot?

2010-02-12 Thread Johannes Graumann
Hello, Has anybody solved this? Thanks, Joh __ 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/posting-guide.html and provide commented, minimal, self-contained,

[R] lattice/ylim: how to fix ylim[1], but have ylim[2] dynamically calculated?

2010-02-14 Thread Johannes Graumann
Hello, When drawing barcharts, I find it not helpful if ylim[1] != 0 - bars for a quantity of 0, that do not show a length of 0 are quite non-intuitive. I have tried to study library(lattice) panel.barchart but am unable to figure out where ylim is taken care of and how one

Re: [R] lattice/ylim: how to fix ylim[1], but have ylim[2] dynamically calculated?

2010-02-14 Thread Johannes Graumann
David Winsemius wrote: On Feb 14, 2010, at 10:33 AM, Johannes Graumann wrote: Hello, When drawing barcharts, I find it not helpful if ylim[1] != 0 - bars for a quantity of 0, that do not show a length of 0 are quite non-intuitive. I have tried to study library(lattice

Re: [R] lattice/ylim: how to fix ylim[1], but have ylim[2] dynamically calculated?

2010-02-14 Thread Johannes Graumann
Deepayan Sarkar wrote: On Sun, Feb 14, 2010 at 7:33 AM, Johannes Graumann johannes_graum...@web.de wrote: Hello, When drawing barcharts, I find it not helpful if ylim[1] != 0 - bars for a quantity of 0, that do not show a length of 0 are quite non-intuitive. I have tried to study

Re: [R] lattice/ylim: how to fix ylim[1], but have ylim[2] dynamically calculated?

2010-02-14 Thread Johannes Graumann
Rolf Turner wrote: On 15/02/2010, at 9:40 AM, Johannes Graumann wrote SNIP (In response to some advice from David Winsemius): I am quite certain that this is the most elaborately worded version of RTFM I have ever come across. I nominate this as a fortune. (Despite Prof

Re: [R] Lattice: How to implement varwidth analogous to graphics::boxplot in bwplot?

2010-02-15 Thread Johannes Graumann
Johannes Graumann wrote: Has anybody solved this? For the benefit of others: after studying ?panel.bwplot I have to admit that bwplot(..., varwidth = TRUE) solves the issue. It's just not documented at ?bwplot Cheers, Joh

[R] Windows-GUI: Install Packages from local zip files and dependencies

2008-02-18 Thread Johannes Graumann
Hi, When using the Packages -- Install packages from local zip files menu item in the windows-gui: 1) is that supposed to automatically pull in dependencies (in that case I have to fix something in my package). 2) If that's not the default: is there a way to make it so? Thanks, Joh

[R] read.table and column class

2008-04-15 Thread Johannes Graumann
Hi, I read in some tabular data using this: read.table(file, quote = \, header = TRUE, sep = \t, stringsAsFactors = FALSE, comment.char = ) One slight problem I have now is that some columns in the data set contain either or +, which means FALSE or TRUE respectively. The command

Re: [R] read.table and column class

2008-04-15 Thread Johannes Graumann
, TRUE) y [1] FALSE FALSE TRUE FALSE TRUE Or define a custom class with an as() method and use that in colClasses. On Tue, Apr 15, 2008 at 5:10 AM, Johannes Graumann [EMAIL PROTECTED] wrote: Hi, I read in some tabular data using this: read.table(file, quote = \, header = TRUE, sep = \t

[R] R CMD check PACKAGE Error

2008-04-18 Thread Johannes Graumann
Hi, Can anyone give me a hint on what's wrong if R CMD check PACKAGE gives me this error: Rd files with unknown encoding: FILENAME.Rd I have checked that file over and over and can't find any weired characters or some such. How can I fix this? Thanks for your time, Joh

Re: [R] R CMD check PACKAGE Error

2008-04-22 Thread Johannes Graumann
Prof Brian Ripley wrote: One thing I should add for the archives -- R 2.7.0 and later have oNew function showNonASCII() in package tools to aid detection of non-ASCII characters in .R and .Rd files. Is there any chance of having the output of this show up in the test log when

[R] Feed list of vectors to vioplot()?

2008-04-22 Thread Johannes Graumann
Hi, I have a list of vectors and am trying to coerce them into something that vioplot will take as groups of data to be plotted independently. Can someone nudge me into the right direction? Thanks, Joh __ R-help@r-project.org mailing list

Re: [R] Feed list of vectors to vioplot()?

2008-04-23 Thread Johannes Graumann
Johannes Graumann wrote: Hi, I have a list of vectors and am trying to coerce them into something that vioplot will take as groups of data to be plotted independently. Can someone nudge me into the right direction? Thanks, Joh __ R-help@r

Re: [R] Feed list of vectors to vioplot()?

2008-04-23 Thread Johannes Graumann
a solution if we do not understand the problem to be solved. Sample data would be helpful along with an understanding of what you would expect for output. On Wed, Apr 23, 2008 at 10:11 AM, Johannes Graumann [EMAIL PROTECTED] wrote: Johannes Graumann wrote: Hi, I have a list of vectors

Re: [R] Feed list of vectors to vioplot()?

2008-04-23 Thread Johannes Graumann
Off-list it was pointed out to me that the trivial solution to this would look like this: myfunction2 - function(x, ...){ mylist - list(x, ...) # plenty of lapply stuff do.call(vioplot,mylist) } Thanks for everybodies patience, Joh Johannes Graumann wrote: posted mailed Thanks

Re: [R] Feed list of vectors to vioplot()?

2008-04-24 Thread Johannes Graumann
not # vioplot expects names as a plain vector, not in the # form names=c(One,Two,Three) ... names - c(One,Two,Three) do.call(vioplot,c(mylist,col=white,names)) Any help on how to do this? Thanks, Joh Johannes Graumann wrote: Off-list it was pointed out to me that the trivial solution to this would

Re: [R] Feed list of vectors to vioplot()?

2008-04-24 Thread Johannes Graumann
Dirty solution: switched off x-axis plotting via 'par' and added it the personalized way with 'axis'. Joh Johannes Graumann wrote: Hm, now I have trouble using additional vioplot parameters. mu-2 si-0.6 bimodal-c(rnorm(1000,-mu,si),rnorm(1000,mu,si)) uniform-runif(2000,-4,4) normal

Re: [R] vioplot kernel smooth via density kernel smooth

2008-04-29 Thread Johannes Graumann
posted mailed Johannes Graumann wrote: Vioplots have great appeal to me, as they manage to squeeze so much information into so little space ... Now some evaluation has made me suspicious about the implementation in the package vioplot and I would like to hear what you say about the appended

[R] heatmap on pre-established hclust output?

2008-05-16 Thread Johannes Graumann
Hi, Can someone please guide me towards how to produce heatmap output from the output of hclust run prior to the actual heatmap call? I have some rather lengthy clustering going on and tweeking the visual output with heatmap recalculating the clustering every time is not feasible. Thanks, Joh

[R] readBin: read from defined offset TO defined offset?

2009-05-11 Thread Johannes Graumann
Hello, With the help of seek I can start readBin from any byte offset within my file that I deem appropriate. What I would like to do is to be able to define the endpoint of that read as well. Is there any solution to that already out there? Thanks for any hints, Joh

Re: [R] readBin: read from defined offset TO defined offset?

2009-05-15 Thread Johannes Graumann
passing would be possible, since there's no way the thing will fit into memory). Thanks again, Joh Johannes Graumann wrote: Hello, With the help of seek I can start readBin from any byte offset within my file that I deem appropriate. What I would like to do is to be able to define

[R] Attempt at package documentation on debian: link to data.frame broken

2007-12-03 Thread Johannes Graumann
Hello, I'm trying to document this little package I'm working on and have the following issue: one of the man files (*.Rd) contains this bit \value{ Returns a \code{\link[base]{data.frame}} representing the tabular data. } The resulting link is broken on my debian system. It points to

Re: [R] Attempt at package documentation on debian: link to data.frame broken

2007-12-03 Thread Johannes Graumann
Prof Brian Ripley wrote: So you are talking about the HTML conversion of your help (.Rd) file? (Links appear in other versions too.) The HTML links are intended to be used via help.start(), not directly, and I think you find it actually points to ../../base/html/data.frame.html, which is

Re: [R] Attempt at package documentation on debian: link to data.frame broken

2007-12-03 Thread Johannes Graumann
posted mailed Prof Brian Ripley wrote: Using non-linked HTML file: style sheet and hyperlinks may be incorrect - which is what happens ... how to switch to linked HTML file? It all works via help.start() (as I did say), at least in a reliable R front end. Running help.start() is how you

Re: [R] Attempt at package documentation on debian: link to data.frame broken

2007-12-03 Thread Johannes Graumann
December 2007 14:08:47 Prof Brian Ripley wrote: On Mon, 3 Dec 2007, Johannes Graumann wrote: posted mailed Prof Brian Ripley wrote: Using non-linked HTML file: style sheet and hyperlinks may be incorrect - which is what happens ... how to switch to linked HTML file? It all works via

[R] R CMD Build feature searches or requests

2007-12-06 Thread Johannes Graumann
Hello, I'm missing two features in R CMD build: 1) Easy building of Windows/zip packaged package version alongside the *nix-style *.tar.gz. Right now I'm doing a scripted version of R CMD build PACKAGE R CMD INSTALL PACKAGE mkdir tmp cp -r

[R] Building package - tab delimited example data issue

2007-12-06 Thread Johannes Graumann
Hello, I'm trying to integrate example data in the shape of a tab delimited ASCII file into my package and therefore dropped it into the data subdirectory. The build works out just fine, but when I attempt to install I get: ** building package indices ... Error in scan(file, what, nmax, sep,

Re: [R] Building package - tab delimited example data issue

2007-12-06 Thread Johannes Graumann
Johannes Graumann wrote: On Thursday 06 December 2007 11:52:46 Peter Dalgaard wrote: Johannes Graumann wrote: Hello, I'm trying to integrate example data in the shape of a tab delimited ASCII file into my package and therefore dropped it into the data subdirectory. The build works out

Re: [R] Building package - tab delimited example data issue

2007-12-06 Thread Johannes Graumann
Peter Dalgaard wrote: Answering to myself and staying with the same example: system.file(data/stroke.csv,package=ISwR) allows direct access to the example file (name). Yes, but... This works right until you turn on LazyData for your package, then you end up with only

Re: [R] Extracting clusters from Data Frame

2007-12-11 Thread Johannes Graumann
Gustaf Rydevik wrote: On Dec 10, 2007 2:28 PM, Johannes Graumann [EMAIL PROTECTED] wrote: Hello, I have a large data frame (1006222 rows), which I subject to a crude clustering attempt that results in a vector stating whether the datapoint represented by a row belongs to a cluster

[R] All anchored series from a vector?

2007-12-18 Thread Johannes Graumann
Hi all, What may be a smart, efficient way to get the following result: myvector - c(A,B,C,D,E) myseries - miracle(myvector) myseries [1] [[1]] A [2] [[1]] A B [3] [[1]] A B [4] [[1]] A B C [5] [[1]] A B C D [6] [[1]] A B C D E Thanks for any hints, Joh

Re: [R] All anchored series from a vector?

2007-12-18 Thread Johannes Graumann
Debugged version: lapply(1:length(myvector), function(.length) { myvector[1:.length] }) Thanks for showing the direction! Joh [EMAIL PROTECTED] wrote: From: Johannes Graumann [EMAIL PROTECTED] Date: 2007/12/18 Tue PM 04:40:37 CST To: [EMAIL PROTECTED] Subject: [R] All anchored series from

Re: [R] All anchored series from a vector?

2007-12-18 Thread Johannes Graumann
Nothing to be sorry about. You suggested a viable solution untested ... my job to figure it out ;0) Joh [EMAIL PROTECTED] wrote: From: [EMAIL PROTECTED] Date: 2007/12/18 Tue PM 02:50:52 CST To: Johannes Graumann [EMAIL PROTECTED] Cc: r-help@r-project.org Subject: Re: [R] All anchored series

Re: [R] Factor Madness

2007-12-19 Thread Johannes Graumann
). This means that it will split matrix columns in data frame arguments, and convert character columns to factors unless stringsAsFactors = TRUE is passed. (I'm guessing 'spectrum' is a data.frame before the code fragment you've shown) hope this helps, Tony Plate Johannes Graumann wrote: Why

[R] Finding overlaps in vector

2007-12-21 Thread Johannes Graumann
posted mailed Dear all, I'm trying to solve the problem, of how to find clusters of values in a vector that are closer than a given value. Illustrated this might look as follows: vector - c(0,0.45,1,2,3,3.25,3.33,3.75,4.1,5,6,6.45,7,7.1,8) When using '0.5' as the proximity requirement, the

Re: [R] Finding overlaps in vector

2007-12-21 Thread Johannes Graumann
6.45 7.00 7.10 8.00 On Dec 21, 2007 4:56 AM, Johannes Graumann [EMAIL PROTECTED] wrote: posted mailed Dear all, I'm trying to solve the problem, of how to find clusters of values in a vector that are closer than a given value. Illustrated this might look as follows: vector - c

Re: [R] Finding overlaps in vector

2007-12-21 Thread Johannes Graumann
with more than 2 x.subsets - x.groups[sapply(x.groups, length) = 3] # print out the subsets invisible(lapply(x.subsets, function(a) print(vector[unique(a)]))) [1] 0.00 0.45 [1] 3.00 3.25 3.33 3.75 4.10 [1] 6.00 6.45 [1] 7.0 7.1 On Dec 21, 2007 4:56 AM, Johannes Graumann [EMAIL PROTECTED] wrote

Re: [R] Finding overlaps in vector

2007-12-22 Thread Johannes Graumann
. Split the output of cutree, instead. Continuing from the prior code: for(el in split(unname(vv), names(vv))) print(el) [1] 0.00 0.45 [1] 1 [1] 2 [1] 3.00 3.25 3.33 3.75 4.10 [1] 5 [1] 6.00 6.45 [1] 7.0 7.1 [1] 8 On Dec 21, 2007 3:24 PM, Johannes Graumann [EMAIL PROTECTED] wrote: Hm

Re: [R] Finding overlaps in vector

2007-12-22 Thread Johannes Graumann
of cutree, instead. Continuing from the prior code: for(el in split(unname(vv), names(vv))) print(el) [1] 0.00 0.45 [1] 1 [1] 2 [1] 3.00 3.25 3.33 3.75 4.10 [1] 5 [1] 6.00 6.45 [1] 7.0 7.1 [1] 8 On Dec 21, 2007 3:24 PM, Johannes Graumann [EMAIL PROTECTED] wrote: Hm, hm, rect.hclust

Re: [R] Finding overlaps in vector

2007-12-22 Thread Johannes Graumann
want to retain groups with 2+ elements then you can just Filter then out: twoplus - function(x) length(x) = 2 Filter(twoplus, split(seq_along(v), ct)) On Dec 22, 2007 5:12 AM, Johannes Graumann [EMAIL PROTECTED] wrote: But cutree does away with the indexes from the original input, which

[R] Plotting labeled impulses: label collision

2008-01-04 Thread Johannes Graumann
Dear all, As you can see from the attachment I'm using R to automatically annotate peptide fragmentation mass spectra, which are represented by impulse plots. I'd like to poll you on approaches of how to deal as generally as possible with the two biggest annotation issues I run into: 1) very

[R] Efficient way to substract all entries in two vectors from each other

2008-01-07 Thread Johannes Graumann
Hi all, I'm to inexperienced to come up with the matrix solution elusively appearing on the horizon for the following problem and would appreciate if you could give me a nudge ... I have two vectors a, and b and need to find the closest match for each value of a in b. How to do that efficiently?

Re: [R] Efficient way to substract all entries in two vectors from eachother

2008-01-07 Thread Johannes Graumann
: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Johannes Graumann Sent: Monday, January 07, 2008 11:19 AM To: [EMAIL PROTECTED] Subject: [R] Efficient way to substract all entries in two vectors from eachother Hi all, I'm to inexperienced to come up with the matrix solution

Re: [R] Efficient way to substract all entries in two vectors from each other

2008-01-07 Thread Johannes Graumann
Ido M. Tamir wrote: matchpt Thanks for this hint. It is exactly what I'm looking for. Cheers, Joh __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Plotting labeled impulses: label collision

2008-01-09 Thread Johannes Graumann
Jim, I finally got back to this implementation of mine and dude, this function is amazing! Thank you so much! Joh On Saturday 05 January 2008 11:42:30 Jim Lemon wrote: Johannes Graumann wrote: Dear all, As you can see from the attachment I'm using R to automatically annotate peptide

[R] plot labeling and expressions

2008-01-15 Thread Johannes Graumann
Hi all, I'm trying to use an iteration that appends something to vector 'peaklabels' like so: peaklabelNames - append( peaklabelNames, substitute(i^{z+phantom()}* (*a*AMU),list(i=y2,z=2,a=0)) ) The vector is supposed to be used with 'text' to put labels on signals in an impulse

Re: [R] plot labeling and expressions

2008-01-16 Thread Johannes Graumann
Uwe Ligges wrote: Johannes Graumann wrote: Hi all, I'm trying to use an iteration that appends something to vector 'peaklabels' like so: peaklabelNames - append( peaklabelNames, substitute(i^{z+phantom()}* (*a*AMU),list(i=y2,z=2,a=0)) ) The vector is supposed

[R] multiple text placements and expressions revisited

2008-01-17 Thread Johannes Graumann
Hi all, I asked something like this earlier but decided that a proper minimal example might be helpfull ;0) Why does this work with regards to the expression (substitution): require(stats) plot(cars) text(5,120,labels=substitute(i^{z+phantom()}*(*a* AMU),list(i=yx,z=2,a=0)))

Re: [R] multiple text placements and expressions revisited

2008-01-17 Thread Johannes Graumann
Thanks for your help! Works like a charm now - I can even append to an expression abject as if it was plain 'c()' ... Joh Prof Brian Ripley wrote: On Thu, 17 Jan 2008, Johannes Graumann wrote: Hi all, I asked something like this earlier but decided that a proper minimal example might

[R] non-plot plotting

2008-01-17 Thread Johannes Graumann
I really do not know ho to else title this ... I want to draw something like the attached png with R and would like to poll you on how to start ... make an empty plot first and then start positioning the characterstring by 'text' and then drawing the lines ... Joh attachment:

Re: [R] non-plot plotting

2008-01-18 Thread Johannes Graumann
Marc, you are kind of crazy ... I asked for a discussion of the path to take to get there, not to have it done ... but hey, what am I complaining about? ;0) Thank you very much! Joh Marc Schwartz wrote: Johannes Graumann wrote: I really do not know ho to else title this ... I want to draw

[R] Regex magic anyone?

2008-01-18 Thread Johannes Graumann
Hi again, how to elegantly split s - ABCDEFGT(P)HIJK into A B C D E F G T(P) H I J K (independently of which letters 'T' or 'P' actually represent ...). Please jumstart my regexing, Joh __ R-help@r-project.org mailing list

Re: [R] Regex magic anyone?

2008-01-18 Thread Johannes Graumann
at: http://gsubfn.googlecode.com and vignette: library(gsubfn); vignette(gsubfn) On Jan 18, 2008 9:09 AM, Johannes Graumann [EMAIL PROTECTED] wrote: Hi again, how to elegantly split s - ABCDEFGT(P)HIJK into A B C D E F G T(P) H I J K (independently of which letters 'T

[R] Package Building and Name Space

2008-01-23 Thread Johannes Graumann
Hello, I just don't get this and would appreciate if someone could write a line or two: I'm trying to build this package and it stops installing after I add the following to the NAMESPACES file: importFrom(gsubfn,strapply) The error during the package test is: Error in MyPackage::MyFunction :

Re: [R] Package Building and Name Space

2008-01-23 Thread Johannes Graumann
need to rely on the old mechanism using require(). Apart from that, it also seems that 'strapply' package is not even in the R search path. Do you see it in the listing created by a call to library()? HTH, Michal Johannes Graumann-2 wrote: Hello, I just don't get this and would

[R] layout: how to check whether a pending plot is going to happen within layout or not?

2008-01-25 Thread Johannes Graumann
How can I test whether a plot I'm about to call will be within a layout or not? I need to fix some par options based on that ... Thanks, Joh __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] xmlEventParse returning trimmed content?

2009-04-09 Thread Johannes Graumann
Hello, I wrote the function below and have the problem, that the text bit returns only a trimmed version (686 chars as far as I can see) of the content under the fetchPeaks condition. Any hunches why that might be? Thanks for pointer, Joh xmlEventParse(fileName, list(

Re: [R] xmlEventParse returning trimmed content?

2009-04-10 Thread Johannes Graumann
= FALSE in the call to xmlEventParse(). If you specify this, you might well get the results you expect. If not, can you post the actual file you are reading so we can reproduce your results. D. Johannes Graumann wrote: Hello, I wrote the function below and have the problem

[R] gsubfn, strapply, REGEX Problem

2008-10-28 Thread Johannes Graumann
Hi all, I swear this used to work: library(gsubfn) strapply(S(AC,P)TVDK(8)EELVQK(8), .[(].{1,2}[)]|.)[[1]] But somewhere along the update path it stopped ... now giving me this Error in base::gsub(pattern, rs, x, ...) : invalid backreference 2 in regular expression Can't figure it out.

Re: [R] gsubfn, strapply, REGEX Problem

2008-10-28 Thread Johannes Graumann
in no error message. On Tue, Oct 28, 2008 at 7:39 AM, Johannes Graumann [EMAIL PROTECTED] wrote: Hi all, I swear this used to work: library(gsubfn) strapply(S(AC,P)TVDK(8)EELVQK(8), .[(].{1,2}[)]|.)[[1]] But somewhere along the update path it stopped ... now giving me this Error in base

Re: [R] gsubfn, strapply, REGEX Problem

2008-10-28 Thread Johannes Graumann
. The NEWS file in the gsubfn distribution does mention the change. On Tue, Oct 28, 2008 at 8:32 AM, Johannes Graumann [EMAIL PROTECTED] wrote: -BEGIN PGP SIGNED MESSAGE- Hash: SHA512 Thanks for looking at this. The \ was an oversight for the example, but the backref bit solves my

  1   2   >