Re: [R] Help with read.csv2

2005-11-17 Thread Marc Schwartz (via MN)
On Thu, 2005-11-17 at 17:09 +0100, Matthieu Cornec wrote: Hello, I am importing the following file ;aa;bb;cc 1988;12;12;12 1989;78;78;12 1990;78;78;12 1991;78;78;12 1992;78;78;12 1993;78;78;12 1994;78;78;12

Re: [R] invert y-axis in barplot

2005-11-16 Thread Marc Schwartz (via MN)
On Wed, 2005-11-16 at 16:46 +, Jörg Schlingemann wrote: Hi! This is probably a very trivial question. Is there an easy way to invert the y-axis (low values on top) when using the function barplot()? Thanks, Jrg You mean something like this?: barplot(1:10, ylim =

Re: [R] bug/feature with barplot?

2005-11-14 Thread Marc Schwartz (via MN)
On Mon, 2005-11-14 at 15:55 +0100, Karin Lagesen wrote: I have found a bug/feature with barplot that at least to me shows undesireable behaviour. When using barplot and plotting fewer groups/levels/factors(I am unsure what they are called) than the number of colors stated in a col statement,

Re: [R] Coercion of percentages by as.numeric

2005-11-14 Thread Marc Schwartz (via MN)
On Mon, 2005-11-14 at 19:07 +0200, Brandt, T. (Tobias) wrote: -Original Message- From: Gabor Grothendieck [mailto:[EMAIL PROTECTED] Sent: 14 November 2005 06:21 PM On 11/14/05, Brandt, T. (Tobias) [EMAIL PROTECTED] wrote: Hi Given that things like the following work a

Re: [R] as.integer with base other than ten.

2005-11-14 Thread Marc Schwartz (via MN)
On Mon, 2005-11-14 at 19:01 +, William Astle wrote: Is there an R function analogous to the C function strtol? I would like to convert a binary string into an integer. Cheers for advice Will There was some discussion in the past and you might want to search the archive for a more

Re: [R] Orientation of tickmarks labels in boxplot/plot

2005-11-02 Thread Marc Schwartz (via MN)
On Wed, 2005-11-02 at 16:06 -0600, Michal Lijowski wrote: Hi, I have been trying draw tickmark labels along the y - axis perpendicular to the y axis while labels along the x - axis parallel to x axis while making box plot. Here is my test dataset. TData ID Ratio 1 0 7.075 2

Re: [R] multiple boxplots

2005-10-28 Thread Marc Schwartz (via MN)
On Fri, 2005-10-28 at 10:26 -0700, J.M. Breiwick wrote: Hello, I want to plot 3 boxplots [ par(mfrow=c(3,1)) ] but the first one has 8 groups, the 2nd has 7 and the third has 6. But I the groups to line up: 1 2 3 4 5 6 7 8 2 3 4 5 6 7 8 3 4 5 6 7 8 where the numbers

Re: [R] encrypted RData file?

2005-10-27 Thread Marc Schwartz (via MN)
On Thu, 2005-10-27 at 16:15 -0500, Na Li wrote: On 27 Oct 2005, Duncan Temple Lang wrote: Yes, it is of interest and was sitting on my todo list at some time. If you want to go ahead and provide code to do it, that would be terrific. There are other areas where encryption would be

Re: [R] Graphics window always overlaps console window!

2005-10-25 Thread Marc Schwartz (via MN)
On Tue, 2005-10-25 at 11:55 -0600, [EMAIL PROTECTED] wrote: Does anyone know how I can set up R so that when I make a graphic, the graphics window remains behind the console window? It's annoying to have to reach for the mouse every time I want to type another line of code (e.g., to add

Re: [R] Graphics window always overlaps console window!

2005-10-25 Thread Marc Schwartz (via MN)
On Tue, 2005-10-25 at 13:07 -0500, Marc Schwartz (via MN) wrote: On Tue, 2005-10-25 at 11:55 -0600, [EMAIL PROTECTED] wrote: Does anyone know how I can set up R so that when I make a graphic, the graphics window remains behind the console window? It's annoying to have to reach for the mouse

Re: [R] Basic: setting resolution and size of an R graphic

2005-10-24 Thread Marc Schwartz (via MN)
On Mon, 2005-10-24 at 22:32 +0200, Dr. med. Peter Robinson wrote: Dear List, I am sorry if this perhaps a too basic question, but I have not found an answer in Google or in the R help system. I am trying to use R to do a very simple analysis of some data (RT-PCR and Western analysis) with a

Re: [R] make three plot to one plot

2005-10-21 Thread Marc Schwartz (via MN)
On Fri, 2005-10-21 at 21:13 +0200, Jan Sabee wrote: Dear all, I want to make three plot below to only one plot together with legend, how can I do that? I have tried with matplot function but I did not succeed. Thanks for your help. Sincerelly, Jan Sabee test.five.x -

Re: [R] peculiar matrices

2005-10-21 Thread Marc Schwartz (via MN)
On Fri, 2005-10-21 at 21:32 +, Ben Bolker wrote: As far as I can tell from reading The Fine Documentation (R Language Definition and Intro to R), matrices are supposed to be of homogeneous types. Yet giving matrix() an inhomogeneous list seems to work, although it produces a peculiar

Re: [R] Automatic rounding of values after factors , converted to numeric, are multipled by a real number

2005-10-20 Thread Marc Schwartz (via MN)
On Thu, 2005-10-20 at 16:20 +0200, Peter Dalgaard wrote: Nelson, Gary (FWE) [EMAIL PROTECTED] writes: Peter, Thank you for your response. I knew how close the values are to integers, but I still don't understand why I don't have control over how the numbers are displayed (rounded or

Re: [R] having scaling problems with a histogram

2005-10-20 Thread Marc Schwartz (via MN)
On Thu, 2005-10-20 at 17:31 +0200, Christian Jones wrote: Hello,?xml:namespace prefix = o ns = urn:schemas-microsoft-com:office:office /o:p/o:p I would like to create a histogram from a data collumn consisting of 4 classes (0; 0.05;0.5;25;75). Due to the difference in scale the classes

Re: [R] Boxplot labels

2005-10-20 Thread Marc Schwartz (via MN)
On Thu, 2005-10-20 at 14:46 -0400, Keith Sabol wrote: I am creating boxplots from a dataframe and would like to add to the standard output a marker representing the value from a particular row in the dataframe. .And I apologize if the solution is as trivial as it seems it should be.

Re: [R] Boxplot labels

2005-10-20 Thread Marc Schwartz (via MN)
On Thu, 2005-10-20 at 15:41 -0500, Marc Schwartz (via MN) wrote: On Thu, 2005-10-20 at 14:46 -0400, Keith Sabol wrote: I am creating boxplots from a dataframe and would like to add to the standard output a marker representing the value from a particular row in the dataframe

Re: [R] matching two plots

2005-10-19 Thread Marc Schwartz (via MN)
On Wed, 2005-10-19 at 18:27 +0100, Rui Cardoso wrote: Hi, I have a problem about graphics. I would like to plot two graphs: a barplot and curve. Here is the code: barplot(dpois(0:45,20),xlim=c(0,45),names=0:45) curve(dnorm(x,20,sqrt(20)),from=0,to=45,add=T) Both graphs are drawn

Re: [R] Efficient ways of finding functions and Breslow-Day test for homogeneity of the odds ratio

2005-10-18 Thread Marc Schwartz (via MN)
On Tue, 2005-10-18 at 14:27 +0100, MJ Price, Social Medicine wrote: Dear all, I have been trying to find a function to calculate the Breslow-Day test for homogeneity of the odds ratio in R. I know the test can be preformed in SAS but i was wondering if anyone could help me to perform this

Re: [R] p-value calculation

2005-10-18 Thread Marc Schwartz (via MN)
On Tue, 2005-10-18 at 17:13 +0200, richard mendes wrote: hello everybody i'm very new at using R so probably this is a very stupid question. I have a problem calculating a p-value. When i do this with excel i can use the method CHIDIST for 1.2654 with 1 freedom degree i get the answer 0.261

Re: [R] Efficient ways of finding functions and Breslow-Day test for homogeneity of the odds ratio

2005-10-18 Thread Marc Schwartz (via MN)
On Wed, 2005-10-19 at 06:47 +1000, Tim Churches wrote: Marc Schwartz (via MN) wrote: There is also code for the Woolf test in ?mantelhaen.test Is there? How is it obtained? The documentation on mantelhaen.test in R 2.2.0 contains a note: Currently, no inference on homogeneity of the odds

Re: [R] y axis in histograms

2005-10-17 Thread Marc Schwartz (via MN)
On Mon, 2005-10-17 at 16:51 +0200, Jorge Gaspar Sanz Salinas wrote: Hi all, This is my first post, I hope you will help me. I've some data to present with histograms. I have few values with almost 99% of the frequencies (thousands) and some other values with low frequencies (below one

Re: [R] apply and plot

2005-10-13 Thread Marc Schwartz (via MN)
On Thu, 2005-10-13 at 14:50 +0100, Luis Ridao Cruz wrote: R-help, I use the code below to plot some data by applying apply function. But I don't know how I can get the argument type or col on the plot function to distinguish the different lines in the graph: apply ( my.data, 2, function

Re: [R] Removing and restoring factor levels (TYPO CORRECTED)

2005-10-13 Thread Marc Schwartz (via MN)
On Thu, 2005-10-13 at 10:02 -0400, Duncan Murdoch wrote: Sorry, a typo in my previous message (parens in the wrong place in the conversion). Here it is corrected: I'm doing a big slow computation, and profiling shows that it is spending a lot of time in match(), apparently because I have

Re: [R] Removing and restoring factor levels (TYPO CORRECTED)

2005-10-13 Thread Marc Schwartz (via MN)
On Thu, 2005-10-13 at 14:31 -0400, Duncan Murdoch wrote: On 10/13/2005 1:07 PM, Marc Schwartz (via MN) wrote: On Thu, 2005-10-13 at 10:02 -0400, Duncan Murdoch wrote: Sorry, a typo in my previous message (parens in the wrong place in the conversion). Here it is corrected: I'm

Re: [R] How to generate for one vector matrix

2005-10-13 Thread Marc Schwartz (via MN)
On Thu, 2005-10-13 at 19:47 +0200, Jan Sabee wrote: Is there any routine to generate for one vector matrix. If I have X I want to generate start from zero to maximum value each vector. For example, I have a vector x = (4,2,3,1,4) I want to generate n=6 times, for 4, start 0 to 4, then 2

Re: [R] subsetting data frame using by() or tapply() or other

2005-10-13 Thread Marc Schwartz (via MN)
On Thu, 2005-10-13 at 14:28 -0600, Brian S Cade wrote: Ok so I see the problem that I'm having creating a new variable (LAG1DBC) in the example data transformation below is that tapply() is creating a list that is not dimensionally consistent with the data frame (data). So how do I go from

Re: [R] shell scripts in R

2005-10-13 Thread Marc Schwartz (via MN)
On Fri, 2005-10-14 at 00:04 +0200, Benedykt P. Barszcz wrote: Dnia czwartek, 13 października 2005 23:25, Andrew Robinson napisał: Marco, use the system command. ?system I hope that this helps, system(ls, intern = FALSE, ignore.stderr = TRUE) Error in as.character(args[[i]]) :

Re: [R] high resolution images for publication

2005-10-13 Thread Marc Schwartz (via MN)
On Thu, 2005-10-13 at 15:20 -0600, Chris Buddenhagen wrote: Dear all I am using R to produce ordinations library(vegan) and the plot function produced looks great on the screen but when I send it to jpg or pdf or eps the resolution is not so good. Can you tell me how to get high resolution

Re: [R] Hmisc latex function

2005-10-12 Thread Marc Schwartz (via MN)
On Wed, 2005-10-12 at 08:33 -0500, Charles Dupont wrote: Marc Schwartz (via MN) wrote: On Tue, 2005-10-11 at 10:01 -0400, Rick Bilonick wrote: I'm using R 2.2.0 on an up-to-date version of Fedora Core 4 with the latest version of Hmisc. When I run an example from the latex function I get

Re: [R] Correlation, by date, of two variables?

2005-10-12 Thread Marc Schwartz (via MN)
On Wed, 2005-10-12 at 12:01 -0700, t c wrote: I have a dataset with three variables: date, var1, var2 How can I calculate the correlation, by date, between var1 and var2? e.g. datevar1var2 1/1/200154 1/1/200185 1/1/200197 2/1/200172

Re: [R] Two factor (or more) non-parametric comparison of means

2005-10-11 Thread Marc Schwartz (via MN)
I suspect that John may be looking for ?friedman.test, which I believe will allow for a two-way non-parametric test. kruskal.test() will perform a non-parametric test on two or more samples (or factor levels) as a generalization of wilcox.test() for one or two samples (or factor levels). HTH,

Re: [R] Hmisc latex function

2005-10-11 Thread Marc Schwartz (via MN)
On Tue, 2005-10-11 at 10:01 -0400, Rick Bilonick wrote: I'm using R 2.2.0 on an up-to-date version of Fedora Core 4 with the latest version of Hmisc. When I run an example from the latex function I get the following: x - matrix(1:6, nrow=2, dimnames=list(c('a','b'),c('c','d','enLine 2')))

Re: [R] Two factor (or more) non-parametric comparison of means

2005-10-11 Thread Marc Schwartz (via MN)
Peter, Good points all around. Hopefully John will provide further clarification. Best regards, Marc On Tue, 2005-10-11 at 19:46 +0200, Peter Dalgaard wrote: Marc Schwartz (via MN) [EMAIL PROTECTED] writes: I suspect that John may be looking for ?friedman.test, which I believe

Re: [R] wildcards and removing variables

2005-10-10 Thread Marc Schwartz (via MN)
On Mon, 2005-10-10 at 10:37 -0400, Afshartous, David wrote: All, Is there are a wildcard in R for varible names as in unix? For example, rm(results*) to remove all variable or function names that begin w/ results? cheers, Dave ps - please

Re: [R] wildcards and removing variables

2005-10-10 Thread Marc Schwartz (via MN)
On Mon, 2005-10-10 at 16:01 +0100, Prof Brian Ripley wrote: On Mon, 10 Oct 2005, Marc Schwartz (via MN) wrote: On Mon, 2005-10-10 at 10:37 -0400, Afshartous, David wrote: All, Is there are a wildcard in R for varible names as in unix? For example, rm(results

Re: [R] spline.des

2005-10-05 Thread Marc Schwartz (via MN)
On Mon, 2005-10-03 at 15:28 -0500, lforzani wrote: Hello, I am using library fda and I can not run a lot of functions because I receive the error: Error in bsplineS(evalarg, breaks, norder, nderiv) : couldn't find function spline.des do you know how I can fix that? Thnaks.

Re: [R] Getting eps into Word documents.

2005-10-04 Thread Marc Schwartz (via MN)
On Tue, 2005-10-04 at 08:38 -0500, Robert Baer wrote: On 03-Oct-05 Marc Schwartz (via MN) wrote: On Mon, 2005-10-03 at 16:31 -0300, Rolf Turner wrote: A student in one of my courses has asked me about getting R graphics output (under Linux) into a Word document. I.e. she wants

Re: [R] gnomeGUI installation

2005-10-03 Thread Marc Schwartz (via MN)
On Mon, 2005-10-03 at 18:01 +0100, Prof Brian Ripley wrote: On Mon, 3 Oct 2005, Daniel Pick wrote: I have successfully downloaded the sources and built R as a shared library on a Red Hat Enterprise Level 3 box. I am now trying to build the GNOME GUI, but configure is barfing on glade.

Re: [R] Getting eps into Word documents.

2005-10-03 Thread Marc Schwartz (via MN)
On Mon, 2005-10-03 at 16:31 -0300, Rolf Turner wrote: A student in one of my courses has asked me about getting R graphics output (under Linux) into a Word document. I.e. she wants to do her R thing under Linux, but then do her word processing using Word. Scanning around the r-help archives

Re: [R] Getting eps into Word documents.

2005-10-03 Thread Marc Schwartz (via MN)
On Mon, 2005-10-03 at 22:00 +0100, Ted Harding wrote: Rolf ( Marc) On 03-Oct-05 Marc Schwartz (via MN) wrote: On Mon, 2005-10-03 at 16:31 -0300, Rolf Turner wrote: A student in one of my courses has asked me about getting R graphics output (under Linux) into a Word document. I.e. she

Re: [R] boxplot and xlim confusion?

2005-09-29 Thread Marc Schwartz (via MN)
On Thu, 2005-09-29 at 15:28 +0200, Karin Lagesen wrote: Petr Pikal [EMAIL PROTECTED] writes: Hi Karin I did not have seen any answer for your question yet so here is a try. I gues you want the horizotal layout or your boxplot. boxplot(split(rnorm(30), rep(1:3, each=10)),

Re: [R] Display values in piechart/barplot

2005-09-29 Thread Marc Schwartz (via MN)
On Thu, 2005-09-29 at 14:34 +0200, Volker Rehbock wrote: Is it possible to automatically display the underlying values of a piechart/barplot in the graphic? If so, which package/function/argument do I need for it? Thanks, Volker Using pie charts are not a particularly good way of

Re: [R] How to add frame (frame.plot=T) to Plot.Design?

2005-09-29 Thread Marc Schwartz (via MN)
On Thu, 2005-09-29 at 18:36 +0200, Jan Verbesselt wrote: Hi R-help, When using the package Design and the plot.Design function all the graphs of an lrm.fit (lrm()) are plotted without a frame (only with axes). How can the frame be added to these plots? In the plot.default

Re: [R] Binary Logit Regression with R

2005-09-29 Thread Marc Schwartz (via MN)
On Thu, 2005-09-29 at 18:08 -0400, Johann Park wrote: Hi to all, I am a PH.D Student doing statistical analysis. I am totally new to R. I previously use Stata and am changing into R. I ususally do with logit regression with binary dependent variable (war occurence:1 or 0). I just

Re: [R] correct syntax

2005-09-28 Thread Marc Schwartz (via MN)
On Thu, 2005-09-29 at 06:57 +1000, Stephen Choularton wrote: Hi I am trying to produce a little table like this: 0 1 0 601 408 1 290 2655 but I cannot get the syntax right. Can anyone help. Stephen Do you have the raw

Re: [R] Using unsplit - unsplit does not seem to reverse the effect of split

2005-09-27 Thread Marc Schwartz (via MN)
On Tue, 2005-09-27 at 19:12 +0200, Søren Højsgaard wrote: In data OME in MASS I would like to extract the first 5 observations per subject (=ID). So I do library(MASS) OMEsub - split(OME, OME$ID) OMEsub - lapply(OMEsub,function(x)x[1:5,]) unsplit(OMEsub, OME$ID) - which results in

Re: [R] Help: x11 position in the Unix environment

2005-09-26 Thread Marc Schwartz (via MN)
On Mon, 2005-09-26 at 17:45 +0200, Shengzhe Wu wrote: Hello, In the Unix environment, I open a window by x11(). May I specify the position of this window by specifying the position of the top left of the window as in Windows environment? Or some other parameters can be used to do that?

Re: [R] Fligner-Policello robust rank test

2005-09-22 Thread Marc Schwartz (via MN)
On Thu, 2005-09-22 at 10:44 -0400, Elisabetta Manduchi wrote: Can anybody tell me if there is an R implementation of the Fligner-Policello robust rank test? Thanks, Elisabetta I don't know about the Fligner-Policello test, but the Fligner-Killeen test has been implemented in fligner.test().

Re: [R] how to extract the column name or value from the numerical value of the matrix

2005-09-19 Thread Marc Schwartz (via MN)
On Mon, 2005-09-19 at 12:30 -0700, shanmuha boopathy wrote: Dear sir, i have a matrix like x-c(1:20) A-matrix(x,4,5) A [,1] [,2] [,3] [,4] [,5] [1,]159 13 17 [2,]26 10 14 18 [3,]37 11 15 19 [4,]48 12 16 20 I want to

Re: [R] How do I get the row indices?

2005-09-16 Thread Marc Schwartz (via MN)
On Fri, 2005-09-16 at 10:34 -0700, Martin Lam wrote: Hi, I was wondering if it's possible to get the row numbers from a filtering. Here's an example: # give me the rows with sepal.length == 6.2 iris[(iris[,1]==6.2),] # output Sepal.Length Sepal.Width Petal.Length Petal.Width

Re: [R] Graphics 'snapshots' in Linux?

2005-09-15 Thread Marc Schwartz (via MN)
On Thu, 2005-09-15 at 10:43 -0400, Tyler Smith wrote: Hi, I'm working on a MEPIS (Debian-based Linux) computer, using the emacs/ESS package to do my R work. I've got some plots that I label interactively using the locate function. With the Windows GUI there is an option to take a

Re: [R] Forcing hist()

2005-09-14 Thread Marc Schwartz (via MN)
On Wed, 2005-09-14 at 20:17 +0200, Par Leijonhufvud wrote: I'm trying to create histogram (using hist()) that fullfill the following criteria: * data is on a ordinal scale (1, 2, 3, 4, 5) * I want bars centered over the number on the x-axis * I want 5 bars of equal

Re: [R] axis of plot

2005-09-01 Thread Marc Schwartz (via MN)
On Thu, 2005-09-01 at 16:58 +0100, [EMAIL PROTECTED] wrote: hy, I need to have the 0 on the bottom left corner of the graph being joined , not with this little hole between x axis and y axis... I've saw this question with the answer one time but i'm unable to find it again.. thks.

Re: [R] png scaling problem

2005-09-01 Thread Marc Schwartz (via MN)
On Thu, 2005-09-01 at 21:51 +0200, Knut Krueger wrote: scaling-4 xywidth-480 resolution-150 png(filename = c:/r/anschluss/plots/4.png, width = xywidth*scaling, height = xywidth*scaling,pointsize = 12, bg = white, res = resolution*scaling) .. barplot(xrow,col =

Re: [R] crosstab for n-way contingency tables

2005-08-30 Thread Marc Schwartz (via MN)
On Tue, 2005-08-30 at 12:28 -0700, Isotta Felli wrote: Dear list. New to R, I'm looking for a way of using crosstab to output low-dimensional (higher than 2) contingency tables (frequencies, per-cents by rows, % by columns, mean, quantiles) I'm looking for something of the following

Re: [R] ylim for graphic

2005-08-29 Thread Marc Schwartz (via MN)
On Mon, 2005-08-29 at 12:58 -0400, Doran, Harold wrote: Dear list: I have some data for which I am generating a series of barplots for percentages. One issue that I am dealing with is that I am trying to get the legend to print in a fixed location for each chart generated by the data.

Re: [R] Matrix oriented computing

2005-08-26 Thread Marc Schwartz (via MN)
On Fri, 2005-08-26 at 15:25 +0200, Peter Dalgaard wrote: Marc Schwartz [EMAIL PROTECTED] writes: x - c(0.005, 0.010, 0.025, 0.05, 0.1, 0.5, 0.9, 0.95, 0.975, 0.99, 0.995) df - c(1:100) mat - sapply(x, qchisq, df) dim(mat) [1] 100 11 str(mat) num [1:100,

Re: [R] Matrix oriented computing

2005-08-26 Thread Marc Schwartz (via MN)
names-.default 0.29 0.02 0.29 0.02 names so almost all the time is in qchisq. On Fri, 26 Aug 2005, Marc Schwartz (via MN) wrote: On Fri, 2005-08-26 at 15:25 +0200, Peter Dalgaard wrote: Marc Schwartz [EMAIL PROTECTED] writes: x - c(0.005, 0.010, 0.025, 0.05

Re: [R] Unpaste Problem

2005-08-26 Thread Marc Schwartz (via MN)
On Fri, 2005-08-26 at 22:39 +0530, A Mani wrote: Hello, Easy ways to unpaste? xp - paste(x2, x3) # x2, x3 are two non-numeric columns. . . xfg - data.frame(xp,sc1, sc2, sc3) # sc1,sc2, sc3 are numeric cols. I want xp to be split up to form a

Re: [R] help on retrieving output from by( ) for regression

2005-08-25 Thread Marc Schwartz (via MN)
That works also (using the example in ?lmList) library(lme4) ?lmList fm1 - lmList(breaks ~ wool | tension, warpbreaks) However, one still would need to use either sapply() or lapply() as below to get the details that Krishna is looking for. 'fm1' above is a list of models (S4 class

Re: [R] help on retrieving output from by( ) for regression

2005-08-25 Thread Marc Schwartz (via MN)
, it would appear that the use of lmList would be advantageous over the use of by(). Marc On Thu, 2005-08-25 at 11:52 -0500, Marc Schwartz (via MN) wrote: That works also (using the example in ?lmList) library(lme4) ?lmList fm1 - lmList(breaks ~ wool | tension, warpbreaks) However

Re: [R] Remove NAs from Barplot

2005-08-24 Thread Marc Schwartz (via MN)
On Wed, 2005-08-24 at 13:06 -0400, Doran, Harold wrote: Dear List: I'm creating a series of barplots using Sweave that must assume a standard format. This is student achievement data and the x-axis must include all grades 3 to 8. In some cases, the data for a grade (or more than one grade)

Re: [R] Seeking help with an apparently simple recoding problem

2005-08-23 Thread Marc Schwartz (via MN)
On Tue, 2005-08-23 at 10:12 -0500, Greg Blevins wrote: Hello, I have struggled, for longer than I care to admit, with this seemingly simple problem, but I cannot find a solution other than the use of long drawn out ifelse statements. I know there has to be a better way. Here is stripped

Re: [R] plotting with same axes

2005-08-22 Thread Marc Schwartz (via MN)
On Mon, 2005-08-22 at 12:58 -0400, [EMAIL PROTECTED] wrote: I have used the 'par' command to overlay one plot on another. But how do I overlay it with the x-values plotted at the same points on the x-axis? Thank you, Steven The specific answer depends to an extent on the graphics

Re: [R] Copying rows from a matrix using a vector of indices

2005-08-17 Thread Marc Schwartz (via MN)
On Wed, 2005-08-17 at 12:02 -0700, Martin Lam wrote: Hi, I am trying to use a vector of indices to select some rows from a matrix. But before I can do that I somehow need to convert 'combinations' into a list, since 'mode(combinations)' says it's 'numerical'. Any idea how I can do that?

Re: [R] recategorizing a vector into discrete states

2005-08-16 Thread Marc Schwartz (via MN)
On Tue, 2005-08-16 at 13:18 -0400, Allan Strand wrote: Hi All, I'm trying to take a numerical vector and produce a new vector of the same length where each element in the first is placed into a category given by a 'breaks-like' vector. The values in the result should equal the lower

Re: [R] use different symbols for frequency in a plot

2005-08-08 Thread Marc Schwartz (via MN)
On Mon, 2005-08-08 at 11:57 -0700, Kerry Bush wrote: suppose I have the following data x-c(rep(.1,5),rep(.2,6),rep(.4,10),rep(.5,20)) y-c(rep(.5,3),rep(.6,8),rep(1.2,8),rep(2.5,18),rep(3,4)) If I plot(x,y) in R, I will only get seven distinct points. What I want to do is to use different

Re: [R] plotting horizontally

2005-07-26 Thread Marc Schwartz (via MN)
On Tue, 2005-07-26 at 08:58 +0700, [EMAIL PROTECTED] wrote: Hello, Is there any way to use plot() horizontally similar to boxplot(., horiz=TRUE)? I want to use to illustrate the distribution of y-values on an adjacent plot using layout(). Thanks in advance for any help. Regards,

Re: [R] elegant solution to transform vector into percentages?

2005-07-26 Thread Marc Schwartz (via MN)
On Tue, 2005-07-26 at 15:48 -0400, [EMAIL PROTECTED] wrote: Hi, I am looking for an elegant way to transform a vector into percentages of values that meet certain criteria. store-c(1,1.4,3,1.1,0.3,0.6,4,5) # now I want to get the precentages of values # that fall into the categories

Re: [R] problems with submitting an eps-file created in R

2005-07-22 Thread Marc Schwartz (via MN)
On Fri, 2005-07-22 at 16:26 +0200, Christian Bieli wrote: Dear all I've got some problems submitting a manuscript, because I can't manage creating the favourable eps-file of a graph created in R. The journal's graphic requirements are as followed: format: eps width: max. 6 inches

Re: [R] Question about 'text' (add lm summary to a plot)

2005-07-22 Thread Marc Schwartz (via MN)
Ok guys, So I played around with this a bit, going back to Dan's original requirements and using Thomas' do.call() approach with legend(). Gabor's approach using sapply() will also work here. I have the following: # Note the leading spaces here for alignment in the table # This could be

Re: [R] R:plot and dots

2005-07-21 Thread Marc Schwartz (via MN)
On Thu, 2005-07-21 at 16:18 +0200, Clark Allan wrote: hi all a very simple question. i have plot(x,y) but i would like to add in on the plot the observation number associated with each point. how can this be done? / allan If you mean the unique observation number associated with

Re: [R] Question about 'text' (add lm summary to a plot)

2005-07-21 Thread Marc Schwartz (via MN)
[Note: the initial posts have been re-arranged to attempt to maintain the flow from top to bottom] Dan Bolser writes: I would like to annotate my plot with a little box containing the slope, intercept and R^2 of a lm on the data. I would like it to look like...

Re: [R] using argument names (of indeterminate number) within a function

2005-07-19 Thread Marc Schwartz (via MN)
On Tue, 2005-07-19 at 17:29 +0200, Dirk Enzmann wrote: Although I tried to find an answer in the manuals and archives, I cannot solve this (please excuse that my English and/or R programming skills are not good enough to state my problem more clearly): I want to write a function with an

Re: [R] .gct file

2005-07-19 Thread Marc Schwartz (via MN)
On Tue, 2005-07-19 at 12:28 -0400, Duncan Murdoch wrote: On 7/19/2005 12:10 PM, mark salsburg wrote: I have two files to compare, one is a regular txt file that I can read in no prob. The other is a .gct file (How do I read in this one?) I tried a simple read.table(data.gct,

Re: [R] .gct file

2005-07-19 Thread Marc Schwartz (via MN)
On Tue, 2005-07-19 at 13:16 -0400, mark salsburg wrote: ok so the gct file looks like this: #1.2 (version number) 7283 19 (matrix size) Name Description Values ... .. How can I tell R to disregard the first two lines and start reading the 3rd line in this

Re: [R] .gct file

2005-07-19 Thread Marc Schwartz (via MN)
on the quotation marks.. thanks.. On 7/19/05, Marc Schwartz (via MN) [EMAIL PROTECTED] wrote: On Tue, 2005-07-19 at 13:16 -0400, mark salsburg wrote: ok so the gct file looks like this: #1.2 (version number) 7283 19 (matrix size) Name Description Values

Re: [R] .gct file

2005-07-19 Thread Marc Schwartz (via MN)
On Tue, 2005-07-19 at 13:08 -0500, Marc Schwartz (via MN) wrote: For the TAB delimited columns, adjust the 'sep' argument to: read.table(data.gct, skip = 2, header = TRUE, sep = \t) The 'quote' argument is by default: quote = \' which should take care of the quoted strings and bring

Re: [R] Can't get sample function from An Introduction to R to work

2005-07-15 Thread Marc Schwartz (via MN)
On Fri, 2005-07-15 at 15:40 -0500, David Groos wrote: I'm trying to figure out R, a piece at a time, hours at a time... I was trying to copy the sample function in, An Introduction to R (for version 2.1.0) by W. N. Venables, D. M. Smith, page 42. Section 10.1 Simple examples provides a

Re: [R] plot the number of replicates at the same point

2005-07-14 Thread Marc Schwartz (via MN)
On Thu, 2005-07-14 at 12:30 -0500, Deepayan Sarkar wrote: On 7/14/05, Kerry Bush [EMAIL PROTECTED] wrote: Thank you for thinking about the problem for me. However, I have found that your method doesn't work at all. You may test the following example:

Re: [R] plot the number of replicates at the same point

2005-07-14 Thread Marc Schwartz (via MN)
On Thu, 2005-07-14 at 15:08 -0500, Deepayan Sarkar wrote: On 7/14/05, Marc Schwartz (via MN) [EMAIL PROTECTED] wrote: On Thu, 2005-07-14 at 12:30 -0500, Deepayan Sarkar wrote: On 7/14/05, Kerry Bush [EMAIL PROTECTED] wrote: Thank you for thinking about the problem for me. However, I

Re: [R] question for IF ELSE usage

2005-07-12 Thread Marc Schwartz (via MN)
On Tue, 2005-07-12 at 16:22 +0200, ecoinfo wrote: Hi R users, Maybe the question is too simple. In a IF ... ELSE ... statement if(cond) cons.expr else alt.expr, IF and ELSE should be at the same line? For example, if (x1==12) { y1 - 5 }else { y1 - 3 } is right, while if

Re: [R] simple question: reading lower triangular matrix

2005-07-12 Thread Marc Schwartz (via MN)
On Tue, 2005-07-12 at 14:37 -0300, Rogério Rosa da Silva wrote: Dear list, I will like to learn how to read a lower triangular matrix in R. The input file *.txt have the following format: A B C D E A 0 B 10 C 250 D 3680 E 47

Re: [R] Plotting Character Variable

2005-07-07 Thread Marc Schwartz (via MN)
On Thu, 2005-07-07 at 11:47 -0400, [EMAIL PROTECTED] wrote: Any ideas about the following problem: I have a matrix (A) that looks like this: gene_names values hsa-mir-124 0.3 hsa-mir-234 0.1 hsa-mir-344 0.4 hsa-mir-333 0.7

Re: [R] xmat[1, 2:3] - NULL

2005-07-07 Thread Marc Schwartz (via MN)
On Thu, 2005-07-07 at 10:20 -0700, Mikkel Grum wrote: I have a situation where I'm filling out a dataframe from a database. Sometimes the database query doesn't get anything, so I end up trying to place NULL in the dataframe like below. temp - NULL xmat - as.data.frame(matrix(NA, 2, 3))

Re: [R] the format of the result

2005-07-01 Thread Marc Schwartz (via MN)
On Fri, 2005-07-01 at 19:40 +0800, ronggui wrote: I write a function to get the frequency and prop of a variable. freq-function(x,digits=3) {naa-is.na(x) nas-sum(naa) if (any(naa)) x-x[!naa] n-length(x) ta-table(x) prop-prop.table(ta)*100 res-rbind(ta,prop) rownames(res)-c(Freq,Prop)

Re: [R] How to rotate the axisnames in a BARPLOT

2005-06-30 Thread Marc Schwartz (via MN)
On Thu, 2005-06-30 at 23:05 +0200, [EMAIL PROTECTED] wrote: Hi all, - how can I do a barplot with rotated axis labels? I've seen the example for just a plot in the FAQ, but I'll missing the coordinates to plot my text at the right position beneath the bars. Is there any (easy?) solution?

Re: [R] faster algorithm for Kendall's tau

2005-06-28 Thread Marc Schwartz (via MN)
On Tue, 2005-06-28 at 13:03 -0400, ferdinand principia wrote: Hi, I need to calculate Kendall's tau for large data vectors (length 100'000). Is somebody aware of a faster algorithm or package function than cor(, method=kendall)? There are ties in the data to be considered (Kendall's

Re: [R] conversion

2005-06-28 Thread Marc Schwartz (via MN)
On Tue, 2005-06-28 at 17:23 -0400, Anna Oganyan wrote: Dear List, How can I convert a list with elements being character strings, like: c(1,2,3,4), “c(1,3,4,2) … to a list with elements as numerical vectors: c(1,2,3,4), c(1,3,4,2)…? Thanks! Anna l - list(c(1,2,3,4), c(1,3,4,2)) l [[1]]

<    1   2   3