Re: [R] hist combines two lowest categories -- is there a workaround?

2008-01-31 Thread Dieter Menne
Ben Fairbank BEN at SSANET.COM writes: When preparing a series of histograms I found that hist was combining the two lowest categories or bins, 1 and 2. Specifying breaks, as illustrated below, resulted in the correct histogram: values - sample(10,500,replace=TRUE) hist(values)

Re: [R] re stricting points in a data frame

2008-01-31 Thread Petr PIKAL
Maybe what is seeked is rather between - function (x,y) x=range(y)[1] x=range(y)[2] A[between(A[,1], x1) between(A[,2], x2),] Regards Petr [EMAIL PROTECTED] [EMAIL PROTECTED] napsal dne 30.01.2008 18:14:55: Let's A = expand.grid(xk1=xk1,xk2=xk2) B = A[A[,1] %in% x1 A[,2] %in% x2, ]

[R] fastest way to compute the squared Euclidean distance between two vectors in R

2008-01-31 Thread Jason Liao
I have a program which needs to compute squared Euclidean distance between two vectors million of times, which the Rprof shows is the bottleneck. I wondered if there is any faster way than my own simple function distance2 = function(x1, x2) { temp = x1-x2 sum(temp*temp) } I have searched

[R] Bug with hist() ?

2008-01-31 Thread Ng Stanley
Hi, Is this a bug from hist() ? The total density is greater than 1. test [1] 0.05077802 -0.50585520 -0.98053648 -0.35513059 0.86767129 0.61736097 [7] 0.14264062 0.26243841 -0.41477782 0.64172618 0.87230953 -0.62570629 [13] 0.43596467 -0.93214106 -0.26781599 h - hist(test, freq =

Re: [R] R-help going to become subscriber-only

2008-01-31 Thread Dieter Menne
Martin Maechler maechler at stat.math.ethz.ch writes: The consequence of this will be that .. before the weekend .. only __e-mail addresses__ subscribed to R-help will be allowed to post to R-help without any moderation, and we are even considering outright rejection all non-subscribers,

Re: [R] fastest way to compute the squared Euclidean distance betweentwo vectors in R

2008-01-31 Thread Dimitris Rizopoulos
hits=0.5 tests=BAYES_00,ETHZ_GEOCITIES,FORGED_RCVD_HELO X-USF-Spam-Flag: NO try this: distance2 - function (x1, x2) { temp - x1 - x2 sum(temp * temp) } x1 - rnorm(1e06) x2 - rnorm(1e06) system.time(for (i in 1:100) distance2(x1, x2)) system.time(for (i in 1:100) crossprod(x1 - x2)) I

[R] Error handling in nlme call

2008-01-31 Thread Rob Forsyth
In some trial simulation work I need to create batch files that will repeatedly generate pseudoreplicate datasets and then create non- linear mixed effects models using nlme. Inevitably these models sometimes fail to converge but I need the batch file to simply move on to another simulation

[R] PDF Different Page Sizes

2008-01-31 Thread Samuel Kemp
Hi, I am using the pdf() function to write out a pdf file that contains multiple plots (a plot per page). Is there a way to specify different page sizes for each plot in the file e.g. on page 1 in the pdf file I want (width=10, height=10) on page 2 (width=20, height=20), etc. Any help would be

Re: [R] PDF Different Page Sizes

2008-01-31 Thread Prof Brian Ripley
On Thu, 31 Jan 2008, Samuel Kemp wrote: Hi, I am using the pdf() function to write out a pdf file that contains multiple plots (a plot per page). Is there a way to specify different page sizes for each plot in the file e.g. on page 1 in the pdf file I want (width=10, height=10) on page 2

Re: [R] Error handling in nlme call

2008-01-31 Thread Prof Brian Ripley
On Thu, 31 Jan 2008, Rob Forsyth wrote: In some trial simulation work I need to create batch files that will repeatedly generate pseudoreplicate datasets and then create non- linear mixed effects models using nlme. Inevitably these models sometimes fail to converge but I need the batch file

Re: [R] R-help going to become subscriber-only

2008-01-31 Thread Martin Maechler
DM == Dieter Menne [EMAIL PROTECTED] on Thu, 31 Jan 2008 08:58:54 + (UTC) writes: DM Martin Maechler maechler at stat.math.ethz.ch writes: The consequence of this will be that .. before the weekend .. only __e-mail addresses__ subscribed to R-help will be allowed to

[R] permutates and/or samples a matrix

2008-01-31 Thread Ng Stanley
Hi, sample(x) only permutates and/or samples from a vector, which I can't use for a matrix. Please help. Thanks Stanley [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] Odp: Delete/insert rows in data matrix

2008-01-31 Thread Petr PIKAL
Hi data[, -5] without fifth column data[-5,] without fifth row rbind(data[1:2,], NA, data[3:4,]) insert row of NA Regards Petr [EMAIL PROTECTED] [EMAIL PROTECTED] napsal dne 31.01.2008 13:52:02: Hello all, I have the following problem. I want to insert blank rows in an data frame so as

Re: [R] R-help going to become

2008-01-31 Thread Dieter Menne
Duncan Murdoch murdoch at stats.uwo.ca writes: It's worth mentioning that you can be subscribed and have no mail delivered: just check Disabled in the mail delivery options on the subscription page. That's what anyone using gmane should do, so that they can still post without requiring

Re: [R] Delete/insert rows in data matrix

2008-01-31 Thread ronggui
How about this: DF-data.frame(x=1:3,y=2:4) DF[4,]-c(4,5) DF x y 1 1 2 2 2 3 3 3 4 4 4 5 DF-DF[,-4] DF x y 1 1 2 2 2 3 3 3 4 4 4 5 On 1/31/08, Eleni Christodoulou [EMAIL PROTECTED] wrote: Hello all, I have the following problem. I want to insert blank rows in an data frame so as to

Re: [R] Memory problem?

2008-01-31 Thread Jay Emerson
Elena, Page 23 of the R Installation Guide provides some memory guidelines that you might find helpful. There are a few things you could try using R, at least to get up and running: - Look at fewer tumors at a time using standard R as you have been. - Look at the ff package, which leaves the

[R] Newbie: Using R to analyse Apache logs

2008-01-31 Thread Raj Mathur
hits=-2.5 tests=BAYES_00,FORGED_RCVD_HELO X-USF-Spam-Flag: NO Hi, I have a requirement to scan Apache logs and discover ``exceptions''. Exceptions can be of two types: 1. A single IP generating a large amount of traffic within a given time frame (for definable values of ``large'' and ``time

Re: [R] permutates and/or samples a matrix

2008-01-31 Thread David Winsemius
On 31 Jan 2008, you wrote in gmane.comp.lang.r.general: Hi, sample(x) only permutates and/or samples from a vector, which I can't use for a matrix. Please help. You just feed sample(.) a redimensioned object. If you need to do it programmatically just multiply the dimensions of the

Re: [R] R-help going to become subscriber-only

2008-01-31 Thread hadley wickham
If you are *subscribed* to R-help, there definitely won't be a problem. Otherwise, such posting may have to wait for moderator approval, but I think that maybe could be changed to work automatically {at least as long as the spammers don't misuse it}. Unless nabble/gmane have hefty anti-spam

Re: [R] permutates and/or samples a matrix

2008-01-31 Thread Gabor Csardi
What exactly is the question? Selecting/permuting rows? M[sample(length=nrow(M), count), ] Selecting/permuting columns? M[ , sample(length=ncol(M), count) ] Permuting elements? structure(sample(M), dim=dim(M)) Selecting elements? sample(M, count) Gabor On Thu, Jan 31, 2008 at 09:10:11PM

Re: [R] R-help going to become subscriber-only

2008-01-31 Thread Dieter Menne
hadley wickham h.wickham at gmail.com writes: Unless nabble/gmane have hefty anti-spam protection in place, spammers will misuse those services. A centralised target is just too tempting. Works well for gmane. You must be subscribed, and you have to enter a word from a blurbed

Re: [R] R-help going to become subscriber-only

2008-01-31 Thread Martin Maechler
DM == Duncan Murdoch [EMAIL PROTECTED] on Thu, 31 Jan 2008 08:07:06 -0500 writes: DM On 1/31/2008 6:48 AM, Martin Maechler wrote: DM == Dieter Menne [EMAIL PROTECTED] on Thu, 31 Jan 2008 08:58:54 + (UTC) writes: DM Martin Maechler maechler at stat.math.ethz.ch

Re: [R] Newbie: Using R to analyse Apache logs

2008-01-31 Thread Zembower, Kevin
Raj, I've been experimenting with R to compute simple statistics from my web logs somewhat similar to what you're describing. For instance, I'm working on trying to classify a unique IP or domain name requestor as 'human' or 'robot' based on the number of seconds between requests for pages. I've

Re: [R] permutates and/or samples a matrix

2008-01-31 Thread Gabor Csardi
I mean On Thu, Jan 31, 2008 at 02:32:20PM +0100, Gabor Csardi wrote: What exactly is the question? Selecting/permuting rows? M[sample(length=nrow(M), count), ] M[sample(seq(length=nrow(M)), count), ] Selecting/permuting columns? M[ , sample(length=ncol(M), count) ] M[ ,

Re: [R] R-help going to become subscriber-only

2008-01-31 Thread François Pinard
[Dieter Menne] hadley wickham h.wickham at gmail.com writes: Unless nabble/gmane have hefty anti-spam protection in place, spammers will misuse those services. A centralised target is just too tempting. Works well for gmane. You must be subscribed, and you have to enter a word from a

Re: [R] dates in French format

2008-01-31 Thread Denis Chabot
(I've put the R Mac list in cc because of the crashes I have experienced trying some of the suggestions below) Hi Gabor and Prof Ripley, Le 31 janv. 08 à 02:11, Prof Brian Ripley a écrit : The output from sessionInfo() the posting guide asked for would have been very helpful here. You

[R] ESS 5.3.7

2008-01-31 Thread 宋时歌
The ESS manual is about version 5.3.7, but I can only find version 5.3.6 on the official web site. Am I missing anything here? Shige __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] R-help going to become subscriber-only

2008-01-31 Thread Marc Schwartz
François Pinard wrote: [Dieter Menne] hadley wickham h.wickham at gmail.com writes: Unless nabble/gmane have hefty anti-spam protection in place, spammers will misuse those services. A centralised target is just too tempting. Works well for gmane. You must be subscribed, and you have to

[R] Delete Columns with zeros

2008-01-31 Thread Jens Oldeland
Hi, I know there have been some discussions on that topic. all their solutions failed in my case... My problem is that I have a dataframe with many zeros. but while plotting they are not useful. so I want to get rid of column 1,3,4,..,n i have 628 columns 1000 rows, so I can´t look

[R] How to calculate Intraclass-coefficient in 2-level Linear Mixed-Effects models?

2008-01-31 Thread Christine A.
Dear R-users, consider a 2-level linear mixed effects model (LME) with random intercept AND random slope for level 1 AND 2. Does anybody know how to calculate Intraclass-coefficient (ICC) for highest (innermost) level 2 ??? In the literature, I did not find an example for these kind of komplex

[R] Attributes of lm, gee and lme functions

2008-01-31 Thread Lucy Radford
Hi I am doing a simulation study as part of my PhD which inovles fitting an lm model, a gee model (via geeglm from geepack package) and a lme model (via the nlme package). After fitting the models i ideally would like to save the p-value of the beta coefficient, so that i can see how many were

Re: [R] Delete Columns with zeros

2008-01-31 Thread Uwe Ligges
Jens Oldeland wrote: Hi, I know there have been some discussions on that topic. all their solutions failed in my case... My problem is that I have a dataframe with many zeros. but while plotting they are not useful. so I want to get rid of column 1,3,4,..,n i have 628 columns 1000

[R] how to customize boxplot

2008-01-31 Thread Cornelis de Gier
Dear List, I'd like to make boxplots of a large number of observations (+/- 20.000), which are distributed log-normal and right skewed. The problem is that with standard boxplots a too large number of observations are displayed as outliers. I also tried to display the log of the observations, but

Re: [R] permutates and/or samples a matrix

2008-01-31 Thread David Winsemius
David Winsemius [EMAIL PROTECTED] wrote : On 31 Jan 2008, Stanley wrote: sample(x) only permutates and/or samples from a vector, which I can't use for a matrix. Please help. You just feed sample(.) a redimensioned object. If you need to do it programmatically just multiply the dimensions

Re: [R] difficulties computing a simple anova

2008-01-31 Thread Will Holcomb
Excellent. I was able to do the analysis with no problem. I just had to do some twiddling with the inputs to get things from the CSV to the format you described. When the data is in that format, I would like to be able to say, create a column which is each square of the difference of each

Re: [R] ESS 5.3.7

2008-01-31 Thread Marc Schwartz
宋时歌 wrote: The ESS manual is about version 5.3.7, but I can only find version 5.3.6 on the official web site. Am I missing anything here? Shige Presumably 5.3.7 is imminent, but not yet rolled up for release. BTW, there is a dedicated ESS e-mail list. More info here:

[R] useR! 2008: submission registration started!

2008-01-31 Thread Uwe Ligges
We are happy to inform you that the online abstract submission and registration for `useR! 2008' is now available online from http://www.R-project.org/useR-2008/ This meeting of the R user community will take place at the Technische Universität Dortmund, Dortmund, Germany, August 12 to 14,

Re: [R] how to customize boxplot

2008-01-31 Thread Marc Schwartz
Cornelis de Gier wrote: Dear List, I'd like to make boxplots of a large number of observations (+/- 20.000), which are distributed log-normal and right skewed. The problem is that with standard boxplots a too large number of observations are displayed as outliers. I also tried to display

Re: [R] Direct adjusted survival?

2008-01-31 Thread Terry Therneau
The lines that I hoped to be the survival probabilities for each edtrt-group adjusted for confounding by log(bili) are nearly identical to the KM-lines, and they certainly don't appear adjusted for the very strong confounding by log(bili). I'm not quite sure what they are, though. Yes,

Re: [R] Delete Columns with zeros

2008-01-31 Thread Gavin Simpson
hits=-2.6 tests�YES_00 X-USF-Spam-Flag: NO On Thu, 2008-01-31 at 15:33 +0100, Jens Oldeland wrote: Hi, I know there have been some discussions on that topic. all their solutions failed in my case... My problem is that I have a dataframe with many zeros. but while plotting they are not

[R] Log rank test power calculations

2008-01-31 Thread Daniel Brewer
Does anyone have any ideas how I could do a power calculation for a log rank test. I would like to know what the suggested sample sizes would be to pick a difference when the control to active are in a ratio of 80% to 20%. Thanks Dan --

[R] Date time, MSSQL and RODBC

2008-01-31 Thread Daniele Amberti
Hi, I have to connect to MSSQL server to get some TimeSeries value (datetime, float) on output. Datetime data type in R is a POSIXct date with information about timezone, daylight or solar time. In the DB I have not this kind of information and this cause some problem at the change of time for

Re: [R] Log rank test power calculations

2008-01-31 Thread Wittner, Ben, Ph.D.
I suspect that someone will write saying that what you want is in one or a few packages, but you might also find the code below useful. (Note, the first two functions require the third.) -Ben ## ## Do power calculation of Freedman for Cox PH as set forth on page 733 ## of 5th edition of Rosner,

Re: [R] Log rank test power calculations

2008-01-31 Thread Marc Schwartz
Daniel Brewer wrote: Does anyone have any ideas how I could do a power calculation for a log rank test. I would like to know what the suggested sample sizes would be to pick a difference when the control to active are in a ratio of 80% to 20%. Thanks Dan See the cpower() and spower()

[R] Is there a formal test for long vs short memory processes?

2008-01-31 Thread tom soyer
Hi, Does anyone know if there are formal tests for long vs short memory processes? i.e., quantitative tests instead of visual examination of corellograms produced by acf. Thanks! -- Tom [[alternative HTML version deleted]] __

[R] Ryacas

2008-01-31 Thread Doran, Harold
I'm tinkering around in Ryacas trying to find an easy way to get the first and second partial derivatives of mu and sigma from a normal distribution (actually a bitterly ugly likelihood but this example works for now). I've done all of the work in Mathematica, but I then need to manually write R

[R] R for iPhone and iPod touch

2008-01-31 Thread Hisaji ONO
Hi. Any plan for developing of R for iPhone and iPod touch underway? Regards. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] Confidence intervals for PCA scores/eigenvalues

2008-01-31 Thread Markus Loecher
Dear all, I have read various descriptions of employing resampling techniques, such as the bootstrap, to estimate the uncertainties of the eigenvectors computed by PCA. When I try __ R-help@r-project.org mailing list

[R] decrease amount of ticks on y axis in lattice levelplot

2008-01-31 Thread Bram Kuijper
Hi all, How can I decrease the number of ticks on the y-axis in a lattice levelplot()? I have as many ticks displayed on the y-axis as I have columns of data (1000 columns), how can I decrease this amount of ticks, while still properly displaying all the data? Note that I get my data from a

[R] R2WinBUGS is broken

2008-01-31 Thread John Smith
Dear R-users, I am trying to use the following code to reproduce results from Prof. Gelman's book, but have the listed error for R2WinBUGS version (the openbugs version is good). I am using R-2.6.1 on windows XP, and all the R packages are most current ones. schools.bug can be found at

Re: [R] decrease amount of ticks on y axis in lattice levelplot

2008-01-31 Thread Weidong Gu
levelplot(, scale=list(y=list(tick.number=))) ?xyplot() Weidong Gu, Department of Medicine University of Alabama, Birmingham 1900 University Blvd., Birmingham, Alabama 35294 Email: [EMAIL PROTECTED] PH: (205)-975-9053 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]

Re: [R] difficulties computing a simple anova

2008-01-31 Thread Greg Snow
Here are a couple of ways: tmp.df - data.frame( group=factor(sample(LETTERS[1:3], 25, replace=TRUE)), val = rnorm(25, 50, 5) ) tmp.df - transform(tmp.df, ssq = ave(val, group, FUN=function(x) x-mean(x))^2 ) # another way tmp.means - tapply(tmp.df$val, tmp.df$group, FUN=mean) tmp.df$ssq2

Re: [R] Direct adjusted survival?

2008-01-31 Thread Peter Jepsen
Dear Terry Therneau: Thank you for replying. Please forgive me for carrying on. Here is the example I gave, now with some output shown: require(survival) data(pbc) coxph(Surv(time,status)~edtrt,data=pbc) --- OUTPUT Call: coxph(formula = Surv(time, status) ~ edtrt, data = pbc) coef

[R] Box Plot With Groups being numbers

2008-01-31 Thread stephen sefick
I would like to Summarize values that are repeated measures at a certain river mile with box plot i.e. The data matrix looks like this 123 124 125 #fiver mile 0.5 0.6 0.7 0.4 0.5 0.6 ...... ...#values I would like to make a boxplot with the

Re: [R] Ryacas

2008-01-31 Thread Gabor Grothendieck
You can only use those functions that Ryacas knows about or that Yacas 1.0.63 knows about: library(Ryacas) View(transtab) and ls(package:Ryacas) The Ryacas home page http://ryacas.googlecode.com has links to the yacas home page. Sometimes using the first few terms of the Taylor representation

Re: [R] Box Plot With Groups being numbers

2008-01-31 Thread Uwe Ligges
stephen sefick wrote: I would like to Summarize values that are repeated measures at a certain river mile with box plot i.e. The data matrix looks like this 123 124 125 #fiver mile 0.5 0.6 0.7 0.4 0.5 0.6 ...... ...#values I

[R] [R-pkgs] Update of X2R (with FishGraph) sent to CRAN, 30 Jan 2008

2008-01-31 Thread Michael H. Prager
X2R is a bundle of three software libraries allowing the user to pass structured data easily from Fortran, C/C++, or AD Model Builder to R. An update to X2R has been sent to CRAN and should be available at mirrors shortly. We also have uploaded FishGraph, a set of R functions to generate

Re: [R] Log rank test power calculations

2008-01-31 Thread Greg Snow
Here is one approach using simulation: library(survival) lrsim1 - function(n, diff=0) { n1 - round( .8 * n ) n2 - n - n1 tmp.df - data.frame( group=rep(c('a','b'), c(n1,n2) ) ) tmp.time - rexp( n, 1/rep(c(3,3+diff),c(n1,n2)) ) tmp.cens - rexp( n, 1/4 ) tmp.df$time - pmin(tmp.time,tmp.cens)

[R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Wensui Liu
Dear All, Is there a emacs/xemacs for unix without compiling it? If it is possible, I prefer something that I can use immediately after putting it on the unix server. Is it possible or am I too stupid to ask this kind of question? Thanks for your insight! -- ===

[R] box plot

2008-01-31 Thread stephen sefick
the data is attached I want to boxplot this data y-read.table(afdmmgs.txt, header=T) boxplot(y) Error in oldClass(stats) - cl : adding class factor to an invalid object I get this error when I try and plot something with na in it can boxplot just overlook the na and boxplot the values?

[R] xreg in ARIMA function

2008-01-31 Thread E C
Hi everyone, I'm trying to include an external regressor in my ARIMA model but am having some problems with the data format in R. I've named my dependent variable of interest count and the external regressor abc. The external regressor is a binary variable. Here are the contents of

Re: [R] box plot

2008-01-31 Thread Rolf Turner
Note that R is case sensitive. (Hint: ``na'' is not the same thing as ``NA'' !!!) cheers, Rolf Turner On 1/02/2008, at 9:21 AM, stephen sefick wrote: the data is attached I want to boxplot this data y-read.table(afdmmgs.txt, header=T) boxplot(y) Error in

Re: [R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Wittner, Ben, Ph.D.
I don't know what unix you're using, but if it's fedora 8, the command yum install emacs did the trick for me. You probably need to have root (a.k.a. superuser) privileges to do that. If you like to use ESS with emacs to edit R scripts, yum install emacs-ess My guess is this would work for any

Re: [R] box plot

2008-01-31 Thread Douglas Bates
2008/1/31 stephen sefick [EMAIL PROTECTED]: the data is attached Thanks. I want to boxplot this data y-read.table(afdmmgs.txt, header=T) boxplot(y) At this point you should use str(y) to examine the structure of the data. str(y) 'data.frame': 4 obs. of 12 variables: $ X215: num

Re: [R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Wensui Liu
I think the unix is SunOS. the secret is I don't have root priviledge. ^_^. So is there a possibility? Thanks. On Jan 31, 2008 3:35 PM, Wittner, Ben, Ph.D. [EMAIL PROTECTED] wrote: I don't know what unix you're using, but if it's fedora 8, the command yum install emacs did the trick for me.

Re: [R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Roland Rau
Hi, Wensui Liu wrote: Dear All, Is there a emacs/xemacs for unix without compiling it? If it is possible, I prefer something that I can use immediately after putting it on the unix server. Is it possible or am I too stupid to ask this kind of question? Thanks for your insight! I guess

[R] need help with my time series analysis problem

2008-01-31 Thread vibhava
hi, I've got two time series sorted datewise. Both the time series have lots of missing data (sometime on same day and sometime on different dates). I am trying to find cross correlation betwen the two time series but don't know exactly how to go about it in R (with missing data values). I can't

Re: [R] Box Plot With Groups being numbers

2008-01-31 Thread Longinus
Try: plot(as.factor(rivermile), measurements) On Jan 31, 2:02 pm, stephen sefick [EMAIL PROTECTED] wrote: I would like to Summarize values that are repeated measures at a certain river mile with box plot i.e. The data matrix looks like this 123 124 125 #fiver mile 0.5

Re: [R] numeric coercion when one or more elements is non numerice

2008-01-31 Thread Arthur Steinmetz
That did it. Thanks! What I was getting was temp GDP CPIYOY 23832 3108.2 garbage 23923 garbage 1.8 24015 3214.1 1.8 24107 3291.8 2 str(temp) `data.frame': 4 obs. of 2 variables: $ GDP :Error in importIntoEnv(impenv, impnames, ns, impvars) : objects

Re: [R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Gabor Csardi
You don't need root access to compile a program (in mose cases). Just a compiler and enough space. Gabor On Thu, Jan 31, 2008 at 03:45:51PM -0500, Wensui Liu wrote: I think the unix is SunOS. the secret is I don't have root priviledge. ^_^. So is there a possibility? Thanks. On Jan 31,

Re: [R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Wensui Liu
SunOS 5.10 Generic_118833-24 sun4u sparc SUNW,Sun-Fire-V890 On Jan 31, 2008 4:10 PM, Horan, Brian [EMAIL PROTECTED] wrote: What does 'uname -a' say? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Wensui Liu Sent: Thursday, January 31, 2008 3:46

Re: [R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Prof Brian Ripley
On Thu, 31 Jan 2008, Wensui Liu wrote: I think the unix is SunOS. the secret is I don't have root priviledge. ^_^. So is there a possibility? At least two projects have Emacs for Solaris aka SunOS using the Solaris pkgadd format. Request your sysadmins to add emacs for you. Thanks. On

Re: [R] dates in French format

2008-01-31 Thread Denis Chabot
Hi all, The crashes I reported earlier were cause by R 2.6.1 for Mac not liking the OS date setting french canada, an issue that has been solved (by Simon Urbanek). The crashes did not occur when the OS was set to use normal french formats for dates. With that setting, the suggestions by

Re: [R] [OT] emacs / xemacs for unix without compile

2008-01-31 Thread Horan, Brian
What does 'uname -a' say? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Wensui Liu Sent: Thursday, January 31, 2008 3:46 PM To: Wittner, Ben, Ph.D. Cc: r-help@r-project.org Subject: Re: [R] [OT] emacs / xemacs for unix without compile I think the unix

Re: [R] dates in French format

2008-01-31 Thread Gabor Grothendieck
On Jan 31, 2008 4:25 PM, Denis Chabot [EMAIL PROTECTED] wrote: Hi all, The crashes I reported earlier were cause by R 2.6.1 for Mac not liking the OS date setting french canada, an issue that has been solved (by Simon Urbanek). The crashes did not occur when the OS was set to use normal

[R] importance measures Random Forest

2008-01-31 Thread Gonçalo Carrera
Hi, I'm using Random Forest e R enviorment (2.4.0) and i want to interprete %IncMSE and IncNodePurity, what are the definitions and the formulas used to obtain this to importance parameters? What is the more reliable importance measure? Thanks in advance Gonçalo [[alternative HTML

Re: [R] difficulties computing a simple anova

2008-01-31 Thread Will Holcomb
That's excellent. I somehow missed that variable names were case sensitive. I was using all lowercase names and it was throwing things off. The error I would get using fun instead of FUN was along the lines of: function (x) x - mean(x) environment: 0x84be554 Error in unique.default(x) : unique()

[R] display pictures

2008-01-31 Thread Roland Rau
Dear all, is there a possibility to display pictures (jpegs, pngs, ...) in R? Maybe the following (hypothetical) code makes it clearer what I have in mind: plot(x=0:100, y=0:100, type=n) rect(xleft=25, xright=75, ybottom=30, ytop=95, fill=apicture.png) Thank you, Roland [[alternative

Re: [R] numeric coercion when one or more elements is non numerice

2008-01-31 Thread Gavin Simpson
hits=-2.6 tests=BAYES_00 X-USF-Spam-Flag: NO On Thu, 2008-01-31 at 12:57 -0800, Arthur Steinmetz wrote: That did it. Thanks! What I was getting was temp GDP CPIYOY 23832 3108.2 garbage 23923 garbage 1.8 24015 3214.1 1.8 24107 3291.8 2 str(temp)

Re: [R] display pictures

2008-01-31 Thread Roland Rau
Dear all, thank you very much for your fast help (via list and also off-list to me). Once again: the R community rocks! Fast and absolutely helpful advice. Thank you all, Roland On Jan 31, 2008 7:30 PM, Gabor Grothendieck [EMAIL PROTECTED] wrote: See:

Re: [R] display pictures

2008-01-31 Thread Gabor Grothendieck
See: http://wiki.r-project.org/rwiki/doku.php?id=tips:graphics-misc:translucency On Jan 31, 2008 7:17 PM, Roland Rau [EMAIL PROTECTED] wrote: Dear all, is there a possibility to display pictures (jpegs, pngs, ...) in R? Maybe the following (hypothetical) code makes it clearer what I have

[R] package could not be loaded

2008-01-31 Thread Suhaila Zainudin
Hi, I am an R newbie. I am running R2.5.1 on WinXP. I am trying to run a pacakge (BNArray1.0) that depends on deal 1.2-26 and dynamicGraph 0.2.0.1 . These are available on http://www.cls.zju.edu.cn/binfo/BNArray/#OLE9 . I have installed the 3 packages (using local files). When I try to run

Re: [R] Newbie: Using R to analyse Apache logs

2008-01-31 Thread Raj Mathur
hits=-2.5 tests=BAYES_00,FORGED_RCVD_HELO X-USF-Spam-Flag: NO Hi Kevin, On Thursday 31 Jan 2008, Zembower, Kevin wrote: Raj, I've been experimenting with R to compute simple statistics from my web logs somewhat similar to what you're describing. For instance, I'm working on trying to

[R] Converting a character string into an object variable

2008-01-31 Thread Brant Inman
R-helpers: Assume that I want to create a series of sequentially named R objects. For example, I might want to call these objects V1, V2, V3 ... V50. To do this, I thought of some sort of looping function like: input.value - seq(1:50) * 3 for(i in 1:50){ paste(V, i, sep=) - input.value[i]

Re: [R] Converting a character string into an object variable

2008-01-31 Thread Benilton Carvalho
assign(paste(V, i, sep=), input.value[i]) b On Jan 31, 2008, at 11:28 PM, Brant Inman wrote: R-helpers: Assume that I want to create a series of sequentially named R objects. For example, I might want to call these objects V1, V2, V3 ... V50. To do this, I thought of some sort of

Re: [R] Converting a character string into an object variable

2008-01-31 Thread Benilton Carvalho
You mean the list has elements called V1 ... V50? assign(paste(output, i, sep=), mylist[[paste(V, i, sep=)]]) b On Jan 31, 2008, at 11:44 PM, Inman, Brant A., M.D. [EMAIL PROTECTED] wrote: Thanks B. Now what if V1, V2, V3 ...V50 are elements of a list and I wanted to extract them? For

Re: [R] Converting a character string into an object variable

2008-01-31 Thread Bill.Venables
This is pretty standard. Here are two solutions, there are many more. input.value - 3*(1:50) for(i in 1:50) assign(paste(V, i, sep = ), input.value[i]) --- input.value - 3*(1:50) for(i in 1:50) eval(substitute( V - input.value[i], list(V =

Re: [R] package could not be loaded

2008-01-31 Thread Prof Brian Ripley
I suggest you update your R, install deal and dynamicGraph from CRAN and then try BNArray. If that then fails, consult its provider for help. Please note the posting guide makes clear that we don't provide support for obsolete versions of R, nor binary package distributions. On Fri, 1 Feb

[R] Cannot save histogram picture when run from script

2008-01-31 Thread Ng Stanley
Hi, I am faced with a strange problem. My picture file is empty when the following statements were run from R script. But the picture file shows up correctly, when the commands are individually run from Rgui. png(file=histogram_correlation.png, bg = white) require(stats) histogram( ~ height |