Re: [R] vectorized sub, gsub, grep, etc.

2008-10-08 Thread Christos Hatzis
Hi John, As I mentioned in our private exchange, this is well known in R, i.e. vectorized versions are not always faster or more efficient than straight loops. It is a misconception that loops should be avoided at any cost. See John Fox's illuminating article on Rnews (p. 46) on this subject. ht

Re: [R] Creating GUIs for R

2008-10-08 Thread Greg Snow
http://rcom.univie.ac.at/ Or just install the RExcelInstaller package and run the installRExcel function. A book on the interface will be coming out sometime in Winter. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] 801.408.8111 > -Origi

Re: [R] RApache

2008-10-08 Thread Jeffrey Horner
lamack lamack wrote: Dear all, I have a question about RApache. Is it possible install it under XAMPP server? I would like of some directions. Is it needed install the R before? As long as your XAMPP server runs some variant of UNIX, rapache should install. and Yes, you should install R fi

Re: [R] vectorized sub, gsub, grep, etc.

2008-10-08 Thread john
Hello Christos, To my surprise, vectorization actually hurt processing speed! #Example X <- c("ab", "cd", "ef") patt <- c("b", "cd", "a") repl <- c("B", "CD", "A") sub2 <- function(pattern, replacement, x) { len <- length(x) if (length(pattern) == 1) pattern <- rep(pattern, len

[R] Adding bars to barplots

2008-10-08 Thread Kay, Winston
I'm trying to add stacked bars to a third column of a barplot and using the following code: year <- c(2004, 2005, 2006, 2007, 2008) nlic <- c(5,8,6,4,2) allow <- c(5,8,6,4,2) shot <- c(1,6,0,1,0) farm <- c(0,0,0,0,0) dfs <-matrix(data=c(nlic, allow), nrow=5, ncol=2) barplot(t(dfs), beside=T, nam

Re: [R] error installing lattice package

2008-10-08 Thread Waichler, Scott R
> > I just updated my Redhat EL systems to R-2.7.2, and tried > to update my > > packages as well. Lattice is one that failed. What do I > need to do? > > Install the R-devel RPM? (Assuming you installed R from an RPM.) Yes, I installed R from an RPM. Installing the R-devel RPM worked on o

Re: [R] NA's in multiple choice summary table in Hmisc

2008-10-08 Thread Frank E Harrell Jr
Erich Studerus, Psychiatrische Uni-Klinik wrote: Hi, I have a set of 30 binary variables measuring side effects after drug treatment. Since each subject can have multiple side effects, I want to display these side effects in a multiple choice table. I'm using the summary and mChoice functions

Re: [R] Creating a matrix

2008-10-08 Thread jim holtman
You can read it in as a data.frame and then use 'aggregate' to get the mean: > x <- read.table(textConnection(' Seller Art. Unit_Price + 1 v1 p1 9.148352 + 2 v1 p1 2.858073 + 3 v1 p2 3.775315 + 4 v1 p2 9.821429 + 5 v1 p3 3.286827 +

Re: [R] question of R to do perl matching and matched string extraction

2008-10-08 Thread jim holtman
Is this what you want: > x <- "/pages-cell.net/deepan/sony/" > z <- gsub("/(.*)/", "\\1", x) > > z [1] "pages-cell.net/deepan/sony" On Wed, Oct 8, 2008 at 8:20 PM, Waverley <[EMAIL PROTECTED]> wrote: > Hi, > > I want to extract some of the substring via pattern recognition. But > I don't know

Re: [R] Choose subset for plot use (bwplot)

2008-10-08 Thread Michael Just
Deepayan, It is empty. > subset(dat11, Pri_No %in% c(1,2)) [1] tallyscapenumidbbContag bbED [6] bbENN_MN bbLPIbbLSIbbPAFRAC bbPD [11] bbPROX_MNpfor year extent bi_ca [16] rex_bin PriNopri1234

Re: [R] Strange horns on notched box plots

2008-10-08 Thread Ben Bolker
Frostygoat gmail.com> writes: > > Hi I'm getting a weird result when I try to switch from a normal box > plot to a notched one. The ends of the box fold down toward the > median giving a horned appearance. Is just the sample itself? It is > small, but the un-notched plot looks okay. Anyway t

[R] interpreting Shapiro-Wilks test result

2008-10-08 Thread Halizah Basiron
Hi all, I am newbie in using R software and also doing statistical test. I want to know if my data in in normal distribution. I have 2 groups of data and I did calculate Shapiro Wilks using R software. Here is the results: Group 1: W = 0.9206, p-value = 0.01683 Group 2: W = 0.9626, p-value =

Re: [R] Choose subset for plot use (bwplot)

2008-10-08 Thread Deepayan Sarkar
On Wed, Oct 8, 2008 at 1:56 PM, Michael Just <[EMAIL PROTECTED]> wrote: > Hello, this code below was from a helpful R-help user. > > dat <- read.csv("Resid_fix2.csv", sep="," , header=T) > dat11 <- dat[1:413,] > # convert ambiguous columns to factors: > dat11$Pri_No <- factor(dat11$Pri_No) > dat11

Re: [R] Histogram colours in lattice.

2008-10-08 Thread Deepayan Sarkar
On Wed, Oct 8, 2008 at 12:15 PM, Rolf Turner <[EMAIL PROTECTED]> wrote: > > On 8/10/2008, at 5:16 PM, Deepayan Sarkar wrote: > >> On 10/7/08, Rolf Turner <[EMAIL PROTECTED]> wrote: >>> >>> I am trying to do a histogram lattice plot and I would like the >>> histogram to be filled with a different

Re: [R] How to join the two tables based on one overlapped column

2008-10-08 Thread ss
Hi Dan, Thanks a lot for this! It works but I have some minor issues. Let's take the working example for instance, which x is of 100 rows and y of 120 rows. After merge(), z is of 120 rows and is sorted by the ascending order of "SEARCH_KEY1". The following is what I got out of the working examp

[R] RApache

2008-10-08 Thread lamack lamack
Dear all, I have a question about RApache. Is it possible install it under XAMPP server? I would like of some directions. Is it needed install the R before? Best regards. _ Confira vídeos com notícias do NY Times, gols direto

[R] question of R to do perl matching and matched string extraction

2008-10-08 Thread Waverley
Hi, I want to extract some of the substring via pattern recognition. But I don't know how to do it in R. In perl: my $url = "/pages-cell.net/deepan/sony/"; if($url =~ m/\/(.*)\//g) { my @result = $1; return @result; } How does the same work in R? Thanks much in advance -- Waverley @ Palo Al

Re: [R] Lattice question: plotting two sets of data, defining groups for the second set

2008-10-08 Thread Deepayan Sarkar
On Wed, Oct 8, 2008 at 3:41 PM, Alex Karner <[EMAIL PROTECTED]> wrote: > R friends, > > I'm running R 2.7.2 on Windows XP SP2. > > I have some data that's amenable to smoothing, and some that's not. I'm > trying to plot smoothed lines for the former along with just points for the > latter in a sing

Re: [R] Creating GUIs for R

2008-10-08 Thread Michael Bibo
Wade Wall gmail.com> writes: > What I am wanting to do is learn to build some simple GUIs for a limited > number of functions. Basically, I am envisioning a screen with check boxes, > a drop down menu etc. that users could select to run analyses on imported > data. The R Gui Generator (http:/

[R] NA's in multiple choice summary table in Hmisc

2008-10-08 Thread Erich Studerus, Psychiatrische Uni-Klinik
Hi, I have a set of 30 binary variables measuring side effects after drug treatment. Since each subject can have multiple side effects, I want to display these side effects in a multiple choice table. I'm using the summary and mChoice functions of the Hmisc package, because it produces nicely

[R] Lattice question: plotting two sets of data, defining groups for the second set

2008-10-08 Thread Alex Karner
R friends, I'm running R 2.7.2 on Windows XP SP2. I have some data that's amenable to smoothing, and some that's not. I'm trying to plot smoothed lines for the former along with just points for the latter in a single panel. The problem comes when trying to break out the points by group. My sample

Re: [R] detect repeated number in a vector

2008-10-08 Thread liujb
Duncan, Thank you so much. It is exactly what I was looking for. However, for some reason, indices() does not work ("Error: could not find function "indices"). I used which(), it worked. Thanks again for the help, Julia Duncan Murdoch-2 wrote: > > On 08/10/2008 2:36 PM, liujb wrote: >> Dear

[R] Creating a matrix

2008-10-08 Thread
Good evening. I have this following table and I would like to turn it into a matrix in which my rows would be filled with de "Sellers", my columns with my "Articles" and my data would be the mean unitary price used by each seller in each produt. Seller Art. Unit Price 1 v1 p1

[R] [randomForest]: display decision trees

2008-10-08 Thread Christian Sturz
Hi, I'm using the package randomForest to generate a classifier for the exemplary iris data set: data(iris) iris.rf<-randomForest(Species~.,iris) Is it possible to print all decision trees in the generated forest? If so, can the trees be also written to disk? What I actually need is to translat

Re: [R] How to join the two tables based on one overlapped column

2008-10-08 Thread Daniel Malter
Hi Allen, look at the following working example: x=data.frame(rnorm(100),1:100) names(x)=c("PROBE_ID1","SEARCH_KEY1") y=data.frame(rnorm(120),1:120) names(y)=c("PROBE_ID2","SEARCH_KEY2") z=merge(x,y,by.x="SEARCH_KEY1",by.y="SEARCH_KEY2",all.x=T,all.y=T) z # check the resulting data frame it sho

Re: [R] How to join the two tables based on one overlapped column

2008-10-08 Thread Daniel Malter
Hi Allen, look at the following working example: x=data.frame(rnorm(100),1:100) names(x)=c("PROBE_ID1","SEARCH_KEY1") y=data.frame(rnorm(120),1:120) names(y)=c("PROBE_ID2","SEARCH_KEY2") z=merge(x,y,by.x="SEARCH_KEY1",by.y="SEARCH_KEY2",all.x=T,all.y=T) z # check the resulting data frame it sho

[R] Applying an R script to data within MySQL? How to?

2008-10-08 Thread Ted Byers
I am trying something I haven't attempted before and the available documentation doesn't quite answer my questions (at least in a way I can understand). My usual course of action would be to extract my data from my DB, do whatever manipulation is necessary, either manually or using a C++ program,

Re: [R] trouble with character \u00e2

2008-10-08 Thread Charles Annis, P.E.
PS - After installing rcom from CRAN as a zipped binary, my home-brew package runs on R2.8.0 beta as it did on r releases before R2.7.x Thank you, Professor Ripley. Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Origina

Re: [R] print results from R

2008-10-08 Thread Liviu Andronic
Hello, On 10/8/08, gallon li <[EMAIL PROTECTED]> wrote: > I want to print the following multiple boxes of output from R. > [..] > Instead of manually typing all these numbers, can I have an easy way to > output these in the right format from R? > It would help to know what commands you issued, a

Re: [R] trouble with character \u00e2

2008-10-08 Thread Charles Annis, P.E.
Professor Ripley: I'm sorry. I misunderstood. R version 2.8.0 beta (2008-10-07 r46631) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type

Re: [R] error installing lattice package

2008-10-08 Thread Prof Brian Ripley
On Wed, 8 Oct 2008, Waichler, Scott R wrote: I just updated my Redhat EL systems to R-2.7.2, and tried to update my packages as well. Lattice is one that failed. What do I need to do? Install the R-devel RPM? (Assuming you installed R from an RPM.) R version 2.7.2 (2008-08-25) install.

Re: [R] ParallelR

2008-10-08 Thread Markus Schmidberger
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi Chi, ParallelR is a commercial software to run R in parallel. It is working very well. We tested it at a small linux cluster. You also can use R and the parallel packages (snow, Rmpi, nws, ...) packages. All Open Source and for free. You probably

Re: [R] trouble with character \u00e2

2008-10-08 Thread Prof Brian Ripley
On Wed, 8 Oct 2008, Charles Annis, P.E. wrote: Thank you Professor: Here is an example using R2.8.0 beta. It shows the coding to be "latin1" But you did not use file.choose or basename here. I installed my package which requires rcom, RODBC, RColorBrewer, survival I was unable to find rcom

[R] error installing lattice package

2008-10-08 Thread Waichler, Scott R
I just updated my Redhat EL systems to R-2.7.2, and tried to update my packages as well. Lattice is one that failed. What do I need to do? R version 2.7.2 (2008-08-25) > install.packages("lattice", repos = "http://cran.fhcrc.org/";) Warning in install.packages("lattice", repos = "http://cran.fh

Re: [R] How to join the two tables based on one overlapped column

2008-10-08 Thread ss
Dear Daniel, Thank you very much for the help! I tried the code and got the following: > John<-read.table(file="John_probe.txt",header=TRUE,row.names=NULL,fill=TRUE) > Susan<-read.table(file="Susan_probe.txt",header=TRUE,row.names=NULL,fill=TRUE) > dim(John) [1] 48701 2 > dim(Susan) [1] 4671

Re: [R] detect repeated number in a vector

2008-10-08 Thread N. Lapidus
Can this be an answer ? which(v %in% names(table(v)[table(v)>1])) [1] 2 5 Nael On Wed, Oct 8, 2008 at 8:36 PM, liujb <[EMAIL PROTECTED]> wrote: > > Dear R users, > > I have this vector that consists numeric numbers. Is there a command that > detects the repeated numbers in a vector and returns

[R] Choose subset for plot use (bwplot)

2008-10-08 Thread Michael Just
Hello, this code below was from a helpful R-help user. dat <- read.csv("Resid_fix2.csv", sep="," , header=T) dat11 <- dat[1:413,] # convert ambiguous columns to factors: dat11$Pri_No <- factor(dat11$Pri_No) dat11$RecovUnit <- factor(dat11$RecovUnit) # plot: require(lattice) bwplot(bbED~ Pri_No |

Re: [R] Re move repeated values

2008-10-08 Thread Erik Iverson
kathie wrote: Dear R users, I'd like to make this data rem.y = c(-1,0,2,4,5) from y = c(-1,-1,0,2,2,2,2,4,4,5,5,5,5,5). That is, I need to remove repeated values. Here is my code, but I don't think it is efficient. How could I improve this? By using the 'unique' fu

Re: [R] Re move repeated values

2008-10-08 Thread Jorge Ivan Velez
Dear Kathie, See ?unique. y = c(-1,-1,0,2,2,2,2,4,4,5,5,5,5,5) unique(y) [1] -1 0 2 4 5 HTH, Jorge On Wed, Oct 8, 2008 at 3:12 PM, kathie <[EMAIL PROTECTED]> wrote: > > Dear R users, > > I'd like to make this data > >rem.y = c(-1,0,2,4,5) > > from > >y = c(-1,-1,0,2,2,2,2,

[R] Announcing R-PLUS 3.3 Experience-Rplus

2008-10-08 Thread Sue Turner
Interactive, Comprehensive and Highly Visual ! R-PLUS 3.3 Rocks :) At a click, import data from different formats, use Excel-like spreadsheet for manipulating your data, create publication-quality reports, generate editable graphics ... and click click to run your favorite models through dialogs

Re: [R] detect repeated number in a vector

2008-10-08 Thread Duncan Murdoch
On 08/10/2008 2:36 PM, liujb wrote: Dear R users, I have this vector that consists numeric numbers. Is there a command that detects the repeated numbers in a vector and returns the index of the repeated numbers (or the actual numbers)? For example, v <- c(3,4,5,7,4). The command would return me

Re: [R] Re move repeated values

2008-10-08 Thread Peter Alspach
Kathryn ?unique and see also duplicated. > unique(c(-1,-1,0,2,2,2,2,4,4,5,5,5,5,5)) [1] -1 0 2 4 5 HTH Peter Alspach > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of kathie > Sent: Thursday, 9 October 2008 8:12 a.m. > To: r-help@r-proje

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Michael Just
Manuel, Thanks, this worked well. I was also toying around with other options in bargraph.CI per your suggestion. Thanks, Michael bargraph.CI(RecovUnit, bbED, group = year, data =scape234, + xlab = "Recovery Unit", ylab = "Edge Density", cex.lab = 1.5, x.leg = 1, + density

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Manuel Morales
On Wed, 2008-10-08 at 09:49 -0700, Dylan Beaudette wrote: > On Wednesday 08 October 2008, Manuel Morales wrote: > > Another option is bargraph.CI or lineplot.CI from the package sciplot. > > > > See http://mutualism.williams.edu/sciplot for examples. > > > > On Tue, 2008-10-07 at 23:31 -0500, Micha

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Manuel Morales
On Wed, 2008-10-08 at 12:01 -0500, Michael Just wrote: > Thank you all for you suggestions. They are all helpful. However, I have > come to a more fundamental problem. Preparing my data to even make such a > graph. I thought I was ready. I will obviously need to find the n, mean, > and confidence

[R] detect repeated number in a vector

2008-10-08 Thread liujb
Dear R users, I have this vector that consists numeric numbers. Is there a command that detects the repeated numbers in a vector and returns the index of the repeated numbers (or the actual numbers)? For example, v <- c(3,4,5,7,4). The command would return me index 2 and 5 (or the repeated number

[R] Re move repeated values

2008-10-08 Thread kathie
Dear R users, I'd like to make this data rem.y = c(-1,0,2,4,5) from y = c(-1,-1,0,2,2,2,2,4,4,5,5,5,5,5). That is, I need to remove repeated values. Here is my code, but I don't think it is efficient. How could I improve this? #---

[R] Thank you very much for all your possible solutions!

2008-10-08 Thread mentor_
I also managed to get the right result but within a for loop ;) So I really appreciate your solutions! Thanks a lot! -- View this message in context: http://www.nabble.com/Using-grep-tp19881017p19882769.html Sent from the R help mailing list archive at Nabble.com. _

[R] Strange horns on notched box plots

2008-10-08 Thread Frostygoat
Hi I'm getting a weird result when I try to switch from a normal box plot to a notched one. The ends of the box fold down toward the median giving a horned appearance. Is just the sample itself? It is small, but the un-notched plot looks okay. Anyway to fix this? e7=as.vector(c(234,37,98,116,4

Re: [R] par(new = TRUE) - overplotting

2008-10-08 Thread Antje
Hi Dieter, thanks a lot for looking inside my code though it was not executable...(sorry for that). Finally, I found a rather stupid mistake. My original code did not use the variable i for the second boxplot. So the second round actually plotted two different data at the two calls... so it h

Re: [R] Ecological Niche Modelling on R

2008-10-08 Thread Kingsford Jones
On Wed, Oct 8, 2008 at 6:17 AM, milton ruser <[EMAIL PROTECTED]> wrote: [snip] > > Finally, in fact GRASP do what I am looking for, and I am starting to > compare the results of this packages with other very wel- know softwares > (DescktopGarp, Maxent, OpenModeller). If someone of you have suggesti

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Dylan Beaudette
On Wednesday 08 October 2008, Manuel Morales wrote: > On Wed, 2008-10-08 at 09:49 -0700, Dylan Beaudette wrote: > > On Wednesday 08 October 2008, Manuel Morales wrote: > > > Another option is bargraph.CI or lineplot.CI from the package sciplot. > > > > > > See http://mutualism.williams.edu/sciplot

Re: [R] Quantiles of weighted sample

2008-10-08 Thread Thomas Lumley
On Wed, 8 Oct 2008, Giovanni Petris wrote: I am wondering if there is a function in R that returns quantiles of a weighted sample, i.e., a set of values, each coming with a different weight. The function quantile() does that only for the case when the weights are all equal. In other words, I am

Re: [R] lme and lmer df's and F-statistics again

2008-10-08 Thread Rolf Turner
On 9/10/2008, at 12:34 AM, Julia S. wrote: Hm. Bert Gunter wrote: that even the most technical aspects of the discipline can be made manifest to anyone with half a brain and a stat 101 course under their belt. I don't think this is something I can use in a rebuttal. The reviewer may

Re: [R] question from Braun/Murdoch book

2008-10-08 Thread Rolf Turner
On 8/10/2008, at 6:48 PM, Erin Hodgess wrote: Hi R People: I am looking at the Braun/Murdoch book, " A First Course in Statistical Programming in R", and I have a question about a function there. It's on page 52, Example 4.5; the sieve of Erastosthenes. There is a line: primes <- c() Is the

Re: [R] Histogram colours in lattice.

2008-10-08 Thread Rolf Turner
On 8/10/2008, at 5:16 PM, Deepayan Sarkar wrote: On 10/7/08, Rolf Turner <[EMAIL PROTECTED]> wrote: I am trying to do a histogram lattice plot and I would like the histogram to be filled with a different colour in each panel. Note: I want every bar in each histogram to be the same colour

Re: [R] Axes limits in rgl.surface (again)

2008-10-08 Thread Duncan Murdoch
On 10/8/2008 2:28 PM, Kevin Wright wrote: An earlier post asked about limiting the z-axis range. I have the opposite problem. I need to expand the z-axis range in order to reduce the vertical relief on the surface plot (i.e. make it flatter). Any suggestions on how to do that? persp3d suppor

Re: [R] trouble with character \u00e2

2008-10-08 Thread Charles Annis, P.E.
Thank you Professor: Here is an example using R2.8.0 beta. It shows the coding to be "latin1" I installed my package which requires rcom, RODBC, RColorBrewer, survival I was unable to find rcom in the packages drop-down menu. I tried mirrors USA(PA) and USA(PA2). rcom does appear in the menu

Re: [R] R seven years ago

2008-10-08 Thread Peter Dalgaard
(Ted Harding) wrote: > On 08-Oct-08 18:00:27, Liviu Andronic wrote: > >> Hello everyone, >> >> As some may know, today Google unveiled its 2001 search index [1]. I >> was curious to see how was R like at that time, and was not >> disappointed. Compared to today's main page [2], seven years ago t

Re: [R] R seven years ago

2008-10-08 Thread Prof Brian Ripley
On Wed, 8 Oct 2008, Liviu Andronic wrote: Hello everyone, As some may know, today Google unveiled its 2001 search index [1]. I was curious to see how was R like at that time, and was not disappointed. Compared to today's main page [2], seven years ago the page looked [3] a bit rudimentary, espe

Re: [R] R seven years ago

2008-10-08 Thread Ted Harding
On 08-Oct-08 18:00:27, Liviu Andronic wrote: > Hello everyone, > > As some may know, today Google unveiled its 2001 search index [1]. I > was curious to see how was R like at that time, and was not > disappointed. Compared to today's main page [2], seven years ago the > page looked [3] a bit rudim

[R] Axes limits in rgl.surface (again)

2008-10-08 Thread Kevin Wright
An earlier post asked about limiting the z-axis range. I have the opposite problem. I need to expand the z-axis range in order to reduce the vertical relief on the surface plot (i.e. make it flatter). Any suggestions on how to do that? Kevin __ R-hel

Re: [R] extracting a subset of sublists form a list; vectorized form

2008-10-08 Thread Claudia Beleites
Dear Javier, > sublists, to be passed to a function, sintaxis like this won't work: > >> sublist <- main.lst[[1:4]] are you looking for: sublist <- main.lst[1:4] HTH Claudia __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-

Re: [R] Fitting weibull, exponential and lognormal distributions to left-truncated data.

2008-10-08 Thread Göran Broström
The package 'eha' fits these distributions (and more) with general left truncation and right censoring, and also regression models a la survreg. Look at 'phreg' for parametric proportional hazards models and 'aftreg' for accelerated failure time models. In your case of no covariates, the two functi

Re: [R] Space between bars in barplot

2008-10-08 Thread Dieter Menne
mentor_ gmx.net> writes: > with the "space" parameter it is possible to change the gap / distance > between the bars, but is it also > possible to change the "space" after each 6th bar? > So for example you have bars from 1 to 6 then a large gap and then the next > six bars from 7 to 12 Space c

[R] R seven years ago

2008-10-08 Thread Liviu Andronic
Hello everyone, As some may know, today Google unveiled its 2001 search index [1]. I was curious to see how was R like at that time, and was not disappointed. Compared to today's main page [2], seven years ago the page looked [3] a bit rudimentary, especially the graphic. (It is wort noting that s

[R] Trouble with word-wrapping and other Tinn-R options

2008-10-08 Thread Glen A Sargeant
I'm running R 2.7.1 with Tinn-R 2.0.0.7 and Windows XP Professional. I have write privileges to the folders containing my R installation, my Tinn-R installation, and my Tinn-R initialization settings. I installed R first, followed by Tinn-R. I have modified my RProfile.site file (in Tinn-R,

Re: [R] Quantiles of weighted sample

2008-10-08 Thread roger koenker
see http://www.nabble.com/weighted-quantiles-to19864562.html#a19865869 url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Champ

[R] Quantiles of weighted sample

2008-10-08 Thread Giovanni Petris
Hello! I am wondering if there is a function in R that returns quantiles of a weighted sample, i.e., a set of values, each coming with a different weight. The function quantile() does that only for the case when the weights are all equal. In other words, I am looking for a quantile function tha

[R] Observed responses in 'augPred' data frame - Wrong order ?

2008-10-08 Thread Marco Barbàra
Dea-R community. I'd like to draw your attention to an issue I have recently encountered while doing my current data analysis. I've got an unexpected (to me) result from the command: > augPred(lmList(my.object)), 'my.object' being a grouped data frame of class: > class(my.object) [1] "nfnGroup

Re: [R] rgl_081.708: rgl.snapshot fails, causing persistent problems

2008-10-08 Thread Duncan Murdoch
On 10/8/2008 12:42 PM, Michael Friendly wrote: Thanks, John What you say about mixing rgl.* and *3d calls may be true, but that is not my problem. After this occurred, I tried many things, but reduced it to the smallest, most basic example (in my original post), using only plot3d(), that worke

Re: [R] trouble with character \u00e2

2008-10-08 Thread Prof Brian Ripley
Can you please try a 2.8.0 beta build? I have a suspicion as to what might be going on, and it cannot happen there. If my guess is correct, nfile <- paste("diagnostic â vs a ", file.label, ".jpg", sep = "") savePlot(path.expand(nfile), type="jpg") may work for you in 2.7.2 (but as I said, I w

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Michael Just
Thank you all for you suggestions. They are all helpful. However, I have come to a more fundamental problem. Preparing my data to even make such a graph. I thought I was ready. I will obviously need to find the n, mean, and confidence interval for my data before I can plot them. for some of these

Re: [R] help for compile sjava

2008-10-08 Thread Martin Morgan
Jacky -- Duncan has kindly allowed me to maintain SJava. A current version is not readily available, but I will provide you with some help off-list (and arrange for more public dissemination soon). Martin Jacky Huang wrote: Dear duncan, I'm writing to ask you for some help about compiling

Re: [R] lme and lmer df's and F-statistics again

2008-10-08 Thread Peter Dalgaard
Julia S. wrote: > Hi there, > > thanks for your help. I did read Bates statement several times, and I am > very glad and thankful that many statisticians spend much time on this. The > problem is, as Dieter pointed it out, that many "end users" often have to > use statistics without being able to f

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Dylan Beaudette
On Wednesday 08 October 2008, Manuel Morales wrote: > Another option is bargraph.CI or lineplot.CI from the package sciplot. > > See http://mutualism.williams.edu/sciplot for examples. > > On Tue, 2008-10-07 at 23:31 -0500, Michael Just wrote: > > Hello, > > I'd appreciate a suggestion on how to co

Re: [R] Using grep

2008-10-08 Thread Ted Harding
On 08-Oct-08 15:19:02, mentor_ wrote: > Hi, > I have a vector A with (200, 201, 202, 203, 204, ... 210) and > a vector B with (201, 204, 209). > Now I would like to get the position in vector A matches with > the entries in vector B > So what I want to have is the following result: > [1] 2 5 10 Fi

Re: [R] rgl_081.708: rgl.snapshot fails, causing persistent problems

2008-10-08 Thread Michael Friendly
Thanks, John What you say about mixing rgl.* and *3d calls may be true, but that is not my problem. After this occurred, I tried many things, but reduced it to the smallest, most basic example (in my original post), using only plot3d(), that worked perfectly in both the CRAN and R-Forge version

[R] extracting a subset of sublists form a list; vectorized form

2008-10-08 Thread jgarcia
Hello; I'll put my real problem through a simple example: I've got a main list: > main.lst <- lst() With a number of sublists: > for(i in 1:1000){ main.lst[[i]] <- list() main.lst[[i]]$first <- runif(1,0,1) main.lst[[i]]$second <- runif(2,3,4) } If later on I need to split this list

Re: [R] Using grep

2008-10-08 Thread Doran, Harold
A <- seq(200,210,1) B <- c(201,204,209) which(A %in% B) > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of mentor_ > Sent: Wednesday, October 08, 2008 11:19 AM > To: r-help@r-project.org > Subject: [R] Using grep > > > Hi, > > I have a vector A wi

Re: [R] Using grep

2008-10-08 Thread Erin Hodgess
Here is a possible solution: > A [1] 200 201 202 203 204 205 206 207 208 209 210 > B [1] 201 204 209 > which(!is.na(match(A,B))) [1] 2 5 10 > Hope this helps, Sincerely, Erin On Wed, Oct 8, 2008 at 10:19 AM, mentor_ <[EMAIL PROTECTED]> wrote: > > Hi, > > I have a vector A with (200, 201, 202

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread hadley wickham
On Tue, Oct 7, 2008 at 11:31 PM, Michael Just <[EMAIL PROTECTED]> wrote: > Hello, > I'd appreciate a suggestion on how to construct plots (barplots?) that use > means on the Y axis instead of density/count. I'd also like to use groups > and plot error or confidence interval bars on these graphs. I

Re: [R] Using grep

2008-10-08 Thread Christos Hatzis
which(A %in% B) -Christos > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of mentor_ > Sent: Wednesday, October 08, 2008 11:19 AM > To: r-help@r-project.org > Subject: [R] Using grep > > > Hi, > > I have a vector A with (200, 201, 202, 203, 204, .

[R] help for compile sjava

2008-10-08 Thread Jacky Huang
Dear duncan, I'm writing to ask you for some help about compiling SJava. I encounter some errors when I compiled SJava on Window XP: Building JNI header files... Extracting the classes from Environment.jar /jdk1.3/bin/jar: not found RForeignReference After executing command "sh ./co

Re: [R] Plot means with error bars - A novice needs help

2008-10-08 Thread Manuel Morales
Another option is bargraph.CI or lineplot.CI from the package sciplot. See http://mutualism.williams.edu/sciplot for examples. On Tue, 2008-10-07 at 23:31 -0500, Michael Just wrote: > Hello, > I'd appreciate a suggestion on how to construct plots (barplots?) that use > means on the Y axis instead

[R] Using grep

2008-10-08 Thread mentor_
Hi, I have a vector A with (200, 201, 202, 203, 204, ... 210) and a vector B with (201, 204, 209). Now I would like to get the position in vector A matches with the entries in vector B So what I want to have is the following result: [1] 2 5 10 I tried the following: grep(B, A) grep(c(B), A) A

[R] Space between bars in barplot

2008-10-08 Thread mentor_
Hi, with the "space" parameter it is possible to change the gap / distance between the bars, but is it also possible to change the "space" after each 6th bar? So for example you have bars from 1 to 6 then a large gap and then the next six bars from 7 to 12 Thanks a lot! -- View this message in

Re: [R] lme and lmer df's and F-statistics again

2008-10-08 Thread Julia S.
Hi there, thanks for your help. I did read Bates statement several times, and I am very glad and thankful that many statisticians spend much time on this. The problem is, as Dieter pointed it out, that many "end users" often have to use statistics without being able to fully understand the math

Re: [R] rgl_081.708: rgl.snapshot fails, causing persistent problems

2008-10-08 Thread John Fox
Dear Michael, I haven't tried rgl.snapshot() in the development version of rgl, so I can't comment on that, but I believe that these are two unrelated problems. I think that the version of rgl on CRAN fails to display text when rgl.* and *3d function calls are mixed, while this works in the develo

Re: [R] lme and lmer df's and F-statistics again

2008-10-08 Thread Julia S.
Hm. Bert Gunter wrote: > > that even the most technical > aspects of the discipline can be made manifest to anyone with half a brain > and a stat 101 course under their belt. > > I don't think this is something I can use in a rebuttal. The reviewer may be offended and reviewers are people one

Re: [R] par(new = TRUE) - overplotting

2008-10-08 Thread Dieter Menne
Antje yahoo.de> writes: > I want to create some boxplots (as png) within an lapply method. To get > nice gridlines behind the boxplot, I plotted it twice and therefore I > set par(new=TRUE). > This works nicely for the first plot but the second does plot on the > first plot > too and creates

Re: [R] ParallelR

2008-10-08 Thread Erin Hodgess
Hi! Have you looked at snow, snowfall, Rmpi, or rparallel, please? Hope this helps, Sincerely, Erin On Wed, Oct 8, 2008 at 10:14 AM, Chi Chan <[EMAIL PROTECTED]> wrote: > Anyone using or has access to ParallelR? I was looking at the page and > found nothing really useful! > > http://www.revoluti

Re: [R] How to join the two tables based on one overlapped column

2008-10-08 Thread Daniel Malter
dataToMerge=data.frame(yourtablename2$PROBE_ID2, yourtablename2$SEARCH_KEY2) ##Puts the two columns of interest in dataset 2 in a separate data frame. mergedData=merge(yourtablename1,dataToMerge,by.x=SEARCH_KEY1,by.y=SEARCH_KEY 2,all.x=T,all.y=F) ##merges the first table with the data frame just c

[R] How to join the two tables based on one overlapped column

2008-10-08 Thread ss
Dear list, I need some clues on this. I have two excel files and I basically want to map one to the other one. Can you give me some hints how to do it? The first excel file, named as "Susan_probe.xls", there are two columns, "PROBE_ID1" and "SEARCH_KEY1" PROBE_ID1 SEARCH_KEY1 ILMN_30212 ILMN

[R] ParallelR

2008-10-08 Thread Chi Chan
Anyone using or has access to ParallelR? I was looking at the page and found nothing really useful! http://www.revolution-computing.com/sitegenius/topic.php?id=195 I want to see if I can run R on a cluster of workstation, and use batch systems like Grid Engine or Xgrid: http://gridengine.sunsour

Re: [R] trouble with character \u00e2

2008-10-08 Thread Charles Annis, P.E.
Thank you Professor: After reading in the file this is what I see: > file.label [1] "EXAMPLE 1 â vs a.xls" charToRaw(file.label) [1] 45 58 41 4d 50 4c 45 20 31 20 c3 a2 20 76 73 20 61 2e 78 6c 73 > Encoding(file.label) [1] "UTF-8" > Encoding(paste("diagnostic â vs a ", file.label, ".jpg", sep

[R] rgl_081.708: rgl.snapshot fails, causing persistent problems

2008-10-08 Thread Michael Friendly
*Summary*: The latest Windows binary version of rgl_081.708 from R-Forge has some problem that causes rgl.snapshot() to fail, at least on my system. Thereafter, *all* rgl 3D graphics are rendered without any text labels. The last problem remains even after (a) removing rgl and re-installing fro

[R] Error in spdep: system is computationally singular

2008-10-08 Thread Alberto Jiménez Valverde
Hi all, I am trying to run an autologistic model using the function errorsarlm from spdep package. **I built an XY matrix extracting the two colums from matriz** coords1<-matriz[matriz$casos1==1, c(4,5)] coords1<-as.matrix(coords1) **I identify neighbours of region points** nb20<-dnearneigh(

[R] par(new = TRUE) - overplotting

2008-10-08 Thread Antje
Hi everybody, I want to create some boxplots (as png) within an lapply method. To get nice gridlines behind the boxplot, I plotted it twice and therefore I set par(new=TRUE). This works nicely for the first plot but the second does plot on the first plot too and creates a mess... How can I for

Re: [R] trouble with character \u00e2

2008-10-08 Thread Prof Brian Ripley
That also works without a hitch on my box, even in vanilla 2.7.2. What exactly is in file.label as given by charToRaw(file.label) Encoding(file.label) ? It should be in UTF-8, and so should paste("diagnostic â vs a ", file.label, ".jpg", sep = "") It looks like the latter is not being treat

  1   2   >