Re: [R] Comparing rows of matrices with different dimensions

2005-11-21 Thread Marc Schwartz (via MN)
On Mon, 2005-11-21 at 16:57 +0100, Antje Döring wrote: > > Hi there, > > > > I have a question, which I thought is very easy to solve, but somehow > I can't find a solution. Probably someone could help me quickly? > > > > Here it is: > > > > I have two matrices: > > > > a > >

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 > > data<-read.csv2(

Re: [R] dev.copy legend problem

2005-11-17 Thread Marc Schwartz (via MN)
On Thu, 2005-11-17 at 17:03 +0100, Florence Combes wrote: > Dear all, > > We are facing this problem for long, and so ask for your help. > > We are plotting 2 graphs in a postscript device (left part -layout > function-), and the common legend for these graphs on the right part. > The legend in t

Re: [R] [Rd] Scan data from a .txt file

2005-11-17 Thread Marc Schwartz (via MN)
I have a feeling that Vasu wants (mistakenly) this: dat <- read.table("clipboard", header = FALSE) > dat V1 V2 V3 V4 1 Name Weight Height Gender 2 Anne150 65 F 3Rob160 68 M 4 George180 65 M 5 Greg205 69 M > str(dat)

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, y

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 mo

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 f

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] 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

Re: [R] R Graphs in Powerpoint

2005-11-01 Thread Marc Schwartz (via MN)
One other option, just to throw it out there, though it involves a few more steps. 1. Generate the R plots as EPS files. 2. Import them into Powerpoint onto the required slides. Resize and/or place as required. Recent versions of Powerpoint will auto-generate a bitmapped preview image upon import

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 number

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 > > wou

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 &

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 ano

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) wit

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 o

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 <- > c(0.02,0

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 > &g

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] 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, "urn:schemas-microsoft-com:office:office" /> > > 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 0;0.05 and 0.5 are displayed wit

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 display

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 a

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 i

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

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 t

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

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] 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.characte

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 f

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] 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 > >> convers

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

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.dat

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/2001

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

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 ?fr

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 >

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, Ma

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

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 >

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? Thn

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 > > >

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 > >> outp

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 arch

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

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 j

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.defa

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 disp

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

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 hav

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 resul

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 > >

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.Wi

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 sn

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 eq

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 = barcol

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.. > > >

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 followi

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] 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 t

Re: [R] Matrix oriented computing

2005-08-26 Thread Marc Schwartz (via MN)
.29 0.02 0.58 0.04 "names<-" >0.29 0.02 0.29 0.02 "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 (

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 > > >

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 ~ wo

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 'lmList'

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 gr

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 strip

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 grap

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 th

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 low

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 di

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 c

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. > >

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 automated

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 inch

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] 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 ass

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: >

Re: [R] .gct file

2005-07-19 Thread Marc Schwartz (via MN)
.. > > How do I get R to read the data as tab delemited, but read in the 2nd > coloumn as one value based on the quotation marks.. > > thanks.. > > On 7/19/05, Marc Schwartz (via MN) <[EMAIL PROTECTED]> wrote: > > On Tue, 2005-07-19 at

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 i

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.ta

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 a

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" provi

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

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: > > > > x1=c(0.6,0.4,.

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

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,

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(N

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

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

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?) soluti

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

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 (Ke

<    1   2   3