[R] snow on OS X

2008-01-16 Thread Rajarshi Guha
Hi, I installed snow via R.app (from R 2.6.0) on OS X (10.4.11, Tiger), though I did not install Rmpi or Rsprng. My first question - has anybody installed Rsprng on OS X? The link from the CRAN package page seems to indicate that the sprng code was not installed rather than Rsprng not buildi

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-16 Thread Max Kuhn
Factors have huge benefits over character data in SAS. For a series regulatory filings, I had miles of SAS code to compute KxK tables where all the cells must show up. For example, if one of the levels of one of the variables was never observed, the corresponding row or column would not show up in

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-16 Thread Daniel Nordlund
> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf > Of Jeffrey J. Hallman > Sent: Wednesday, January 16, 2008 2:38 PM > To: [EMAIL PROTECTED] > Subject: Re: [R] things that are difficult/impossible to do in SAS or SPSS > but simple in > R > > SAS has no fa

[R] step-down bootstrap multiple comparisons in R?

2008-01-16 Thread Stewart, Gregory
I've been using SAS PROC MULTTEST to perform multiple comparisons on data that are not normally distributed by using the stepdown bootstrap procedures of Westfall and Young (1993). According to the SAS manual, "the bootstrap method creates pseudo-data sets by sampling observations with replacement

Re: [R] Problem with running installed package on Windows

2008-01-16 Thread James Ferguson
Thanks for everyone's help - I figured out the problem, or at least some way of making the package work properly. I noticed that everything was OK if I logged in as administrator. Not sure why that might be. Regards from Victoria, James Ferguson > On Wed, Jan 16, 2008 at 08:58:50AM -0800, James

Re: [R] Non linear regression with 2 explanatory variables

2008-01-16 Thread Katharine Mullen
Dear Rolf, One thing that sometimes makes nls easier to apply is using the 'formula' argument like you would use the 'fn' argument of optim. That is, if you have a residual function that has arguments x, y, a, b and you need to optimize a and b, you would make a call like nls(~resid(x,y,a=astart

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-16 Thread Jeffrey J. Hallman
SAS has no facilities for date arithmetic and no easy way to build it yourself. In fact, that's the biggest problem with SAS: it stinks as a programming environment, so it's always much more difficult than it should be to do something new. As soon as you get away from the canned procs and have to

Re: [R] Rmpi on Linux x86_64 GNU/Linux

2008-01-16 Thread Prof Brian Ripley
On Wed, 16 Jan 2008, Brian O'Gorman wrote: > I'm having trouble with R CMD INSTALL Rmpi_0.5-5.tar.gz > --configure-args=~/lam > > lam is is installed locally. That's the problem. You need to build a shared library version, not a static version, to allow it to be compiled into shared objects. (

Re: [R] Can anyone help??

2008-01-16 Thread Cody Hamilton
Try replacing '\' with '/' Regards, -Cody Cody Hamilton Edwards Lifesciences -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of hoogeebear Sent: Wednesday, January 16, 2008 10:35 AM To: r-help@r-project.org Subject: [R] Can anyone help?? Hi, Im gettin

Re: [R] Rmpi on Linux x86_64 GNU/Linux

2008-01-16 Thread Martin Morgan
Hi Brian -- Prof. Ripley's reply to an earlier question http://tolstoy.newcastle.edu.au/R/e2/help/07/08/24306.html points to the need for mpi to be compiled as a shared library. Hope that helps, Martin "Brian O'Gorman" writes: > I'm having trouble with R CMD INSTALL Rmpi_0.5-5.tar.gz > --con

Re: [R] color ranges on a 2D plot

2008-01-16 Thread Marc Schwartz
dxc13 wrote: > useR's > > I am trying to color the points on a scatter plot (code below) with two > colors. Red for values 0.5 -1.0 and blue for 0.0 - .49. Does anyone know a > easy way to do this? > > x<- runif(100, 0, 1) > y<- runif(100, 0, 1) > plot(y ~ x, pch=16) > > Thanks, > dxc13 You did

Re: [R] color ranges on a 2D plot

2008-01-16 Thread Greg Snow
You did not say what variable you want to be in the ranges to determine the colors, but try any of these to see if they do what you want: > x <- runif(100) > y <- runif(100) > z <- runif(100) > plot(x,y, pch=16, col=ifelse( x < 0.5, 'blue','red')) > plot(x,y, pch=16, col=ifelse( y < 0.5, 'blue','r

[R] To Share And To Invest In

2008-01-16 Thread R. D. Hargaalga
THAT THEY MAY BE . . .. I am an insider who is an initiate into masking societies (masquerades) in West Africa. I am also an alumnus of the University of Alberta, Canada (1974). I have been engaged in research for over three decades on the causes, which prompted the remarks mentioned below

Re: [R] outer on matrices / arrays

2008-01-16 Thread Gabor Grothendieck
Check out: https://stat.ethz.ch/pipermail/r-help/2005-January/065085.html On Jan 16, 2008 3:59 PM, eariasca <[EMAIL PROTECTED]> wrote: > Hi, > > I would like to evaluate a function f(x,y) taking two vectors as > entries on all pairs of X rows and Y rows, where X and Y are n-by-p > and m-by-p matr

[R] outer on matrices / arrays

2008-01-16 Thread eariasca
Hi, I would like to evaluate a function f(x,y) taking two vectors as entries on all pairs of X rows and Y rows, where X and Y are n-by-p and m-by-p matrices respectively. Of course, I would like to avoid loops if possible, as the dataset I want to process is quite large. I was not able t

Re: [R] using table function

2008-01-16 Thread John Kane
Can you give a simple example. I am not clear on what you mean but table function in R console" when you talk about two vectors. --- Ricardo Perrone <[EMAIL PROTECTED]> wrote: > Hi, > > How to join two large vectors ordered, where one > has the variable's levels and another has the > frequenc

Re: [R] using an element of an array as a new object

2008-01-16 Thread Greg Snow
I would work with 2 different objects to be safe (you can rename the 2nd after everything works if you want). Try something like: > filesBox2 <- array( list(), dim(filesBox) ) > filesBox2[4,1] <- read.table(paste(dir2, filesBox[4,1], sep=""), header = FALSE) . . . Hope this helps, -- Gregory (

[R] color ranges on a 2D plot

2008-01-16 Thread dxc13
useR's I am trying to color the points on a scatter plot (code below) with two colors. Red for values 0.5 -1.0 and blue for 0.0 - .49. Does anyone know a easy way to do this? x <- runif(100, 0, 1) y <- runif(100, 0, 1) plot(y ~ x, pch=16) Thanks, dxc13 -- View this message in context: http:

Re: [R] Question concerning pdfs

2008-01-16 Thread Greg Snow
You may want to look at the logspline package. It is another way of estimating densities and allows you to specify boundries on the domain of x. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Message- > F

[R] Rmpi on Linux x86_64 GNU/Linux

2008-01-16 Thread Brian O'Gorman
I'm having trouble with R CMD INSTALL Rmpi_0.5-5.tar.gz --configure-args=~/lam lam is is installed locally. lamboot -d (or lamboot-d and also recon) works. make -k check from the lamtest suite passes all tests. Is this is problem with the -fPIC compiler as in the message? Should it be modifie

Re: [R] Grouping data

2008-01-16 Thread Marc Schwartz
K. Elo wrote: > Hi, > > I am quite new to R (but like it very much!), so please apologize if > this is a too simple question. > > I have a large data frame consisting of data from a survey. There is, > for example, information about age and education (a numeric value from > 1-9). Now I would like t

Re: [R] Grouping data

2008-01-16 Thread Andrew Robinson
Hi Kimmo, try cut() to create a factor with levels according to the range of values, and (among other options) table() to make the table. Cheers Andrew. On Wed, Jan 16, 2008 at 10:06:23PM +0200, K. Elo wrote: > Hi, > > I am quite new to R (but like it very much!), so please apologize if > thi

[R] Grouping data

2008-01-16 Thread K. Elo
Hi, I am quite new to R (but like it very much!), so please apologize if this is a too simple question. I have a large data frame consisting of data from a survey. There is, for example, information about age and education (a numeric value from 1-9). Now I would like to extract the total amoun

Re: [R] Pause on graphics

2008-01-16 Thread Vu Nguyen
Sorry, I meant par(ask=TRUE) before each plot. On MacOS, when I ran the my script, it did not wait for input at all. Only the last graph was remaining. Thanks, Vu - Original Message From: Rolf Turner <[EMAIL PROTECTED]> To: Vu Nguyen <[EMAIL PROTECTED]> Cc: R-help@r-project.org Sent:

Re: [R] R and Vista

2008-01-16 Thread Alberto Monteiro
Duncan Murdoch wrote: > >> Thank you for your answer. >> I get a very strange path on my R.home(): >> " C:\\PROGRA~2\\R\\R-25~1.0" >> And the command >> readLines(C:\\PROGRA~2\\R\\R-25~1.0\etc\Rprofile.site) returns an Error (a >> syntax one, I triple checked the syntax...) >> Would you know wha

Re: [R] Non linear regression with 2 explanatory variables

2008-01-16 Thread Rolf Turner
I have never had much success in using nls(). If you scan the archives you will find one or two postings from me on this topic. I have received no useful responses to these postings. I have found that anything that I tried (and failed) to do using nls() could be done quite easily using optim(

Re: [R] replace numbers in a column conditional on their value

2008-01-16 Thread Rolf Turner
On 17/01/2008, at 5:45 AM, Jabez Wilson wrote: > Dear R help, > > I have a data frame column in which I would like to replace some > of the numbers dependent on their value. > > data frame = zz > > AveExpr t P.Value FC > 7.481964 7.323950 1.778503e-04 2.218760 > 7.5857

Re: [R] Pause on graphics

2008-01-16 Thread Rolf Turner
On 17/01/2008, at 6:46 AM, Vu Nguyen wrote: > My R script needs to pause or wait for a key or mouse on each > graph. I used the following statement after each plot > > par(ask=TRUE) > > This works on Windows but not on MacOS. > Is there any way to pause on graph in MacOS? (1) Why do you set p

Re: [R] Can anyone help??

2008-01-16 Thread hoogeebear
Hi Guys, Thank you very much for the quick response. Adding the extra \ has fixed the errors. Thanks again, Regards, John. hoogeebear wrote: > > Hi, > > Im getting the following error when I attempt to read a csv file into R. I > have the file saved on my C drive. Id really appreciate it

Re: [R] Pause on graphics

2008-01-16 Thread mel
Vu Nguyen a écrit : > My R script needs to pause or wait for a key or mouse on each graph. I used > the following statement after each plot > par(ask=TRUE) > This works on Windows but not on MacOS. > Is there any way to pause on graph in MacOS? > Thanks, > Vu scan() should work __

Re: [R] Can anyone help??

2008-01-16 Thread Prof Brian Ripley
FAQ Q 7.8 rw-FAQ Q 2.16 should both help, as perhaps should the warning about > unrecognized escape removed from "C:\onlinesurvey.csv" You need to use forward slashes or double the backslashes. On Wed, 16 Jan 2008, hoogeebear wrote: > > Hi, > > Im getting the following error when I attempt to

Re: [R] using an element of an array as a new object

2008-01-16 Thread Julian Burgos
See ?assign. ppaarrkk wrote: > I have an array called filesBox. I want to take each element of the first > column and assign a dataset to it. > > For example : > > filesBox[4,1] returns > > [1] "fileR" > > Then I want to assign "fileR" which exists as a text file to the R object > "fileR" lik

Re: [R] Can anyone help??

2008-01-16 Thread Wittner, Ben, Ph.D.
The problem is that the backslash character (i.e., the "\" in "C:\online") is used to separate folder names on MS Windows, but it's used for something else in R. Try survey <- read.csv(file=file.path("C:", "onlinesurvey.csv"), head= TRUE, sep =",") -Ben p.s., As the posting guide instructs, plea

[R] Can anyone help??

2008-01-16 Thread hoogeebear
Hi, Im getting the following error when I attempt to read a csv file into R. I have the file saved on my C drive. Id really appreciate it if anyone can help My Command: > survey <- read.csv (file = "C:\onlinesurvey.csv",head= TRUE, sep =",") Error: Error in file(file, "r") : unable to ope

Re: [R] Probability weights with density estimation

2008-01-16 Thread Charles C. Berry
On Wed, 16 Jan 2008, David Winsemius wrote: > I am a physician examining an NHANES dataset available at the NCHS > website: > http://www.cdc.gov/nchs/about/major/nhanes/nhanes2005-2006/demo_d.xpt > http://www.cdc.gov/nchs/about/major/nhanes/nhanes2005-2006/hdl_d.xpt > http://www.cdc.gov/nchs/abo

Re: [R] odfWeave and xtable

2008-01-16 Thread Max Kuhn
2008/1/16 David Hajage <[EMAIL PROTECTED]>: > Ps : I would like to know if there is an R project to include all existing > format outputs (latex with Sweave, odf with odfWeave, html with rWeaveHTML) > and all the wonderful work of their author in a same package or in a same > project. All of these

Re: [R] Error in plot.new() : figure margins too large

2008-01-16 Thread Richard . Cotton
> I am getting error : > Error in plot.new() : figure margins too large > Can you please help > -code is > temp <- seq(550/2, 2200, 550/2) > sV2 <- c(1, 1+temp[-length(temp)]) > eV2 <- temp > i=1 > bitmap(paste("c:/vidh

Re: [R] degrees of freedom and random effects in lmer

2008-01-16 Thread Douglas Bates
I suggest this discussion be moved to the R-SIG-mixed-models mailing list which I am cc:ing on this reply. Please delete the R-help mailing list from replies to this message. On Jan 16, 2008 11:44 AM, Feldman, Tracy <[EMAIL PROTECTED]> wrote: > Dear All, > > > > I used lmer for data with non-norm

[R] Pause on graphics

2008-01-16 Thread Vu Nguyen
My R script needs to pause or wait for a key or mouse on each graph. I used the following statement after each plot par(ask=TRUE) This works on Windows but not on MacOS. Is there any way to pause on graph in MacOS? Thanks, Vu [[alternative HTML version deleted]]

[R] degrees of freedom and random effects in lmer

2008-01-16 Thread Feldman, Tracy
Dear All, I used lmer for data with non-normally distributed error and both fixed and random effects. I tried to calculate a "Type III" sums of squares result, by I conducting likelihood ratio tests of the full model against a model reduced by one variable at a time (for each variable separate

Re: [R] Own classes in "histogram"

2008-01-16 Thread David Hewitt
Denis Aydin wrote: > > I try to make a histogram from a variable that contains the number of > shoots from about 1000 individuals from a specific plant species (the > range is 1-110). > Those numbers are highly skewed to the right. > > My question is: how can I make my own classes with the latt

[R] "Halfing" Kernels in Kernel Density?

2008-01-16 Thread Nils Müller-Scheessel
hits=-2.5 tests=BAYES_00,FORGED_RCVD_HELO X-USF-Spam-Flag: NO Hello, I have got a question regarding the function "density": Is there a way in R to "half" kernels or to use kernels which are heavily skewed to one side? My problem is as follows: I have a set of archaeological sites which are d

Re: [R] Own classes in "histogram"

2008-01-16 Thread Deepayan Sarkar
On 1/16/08, Denis Aydin <[EMAIL PROTECTED]> wrote: > Hi, > > I try to make a histogram from a variable that contains the number of > shoots from about 1000 individuals from a specific plant species (the range > is 1-110). > Those numbers are highly skewed to the right. I would suggest you conside

Re: [R] [PS] Different results in calculating SD of 2 numbers

2008-01-16 Thread Ben Fairbank
And another problem, in addition to the points made by others, is that the formula for the SD gives a biased estimate (it underestimates it) of the population SD for small n when sampling from a normal distribution. When n is about twelve or so or more, the bias can usually be ignored (it is about

Re: [R] Problem with running installed package on Windows

2008-01-16 Thread Gabor Csardi
On Wed, Jan 16, 2008 at 08:58:50AM -0800, James Ferguson wrote: > Uwe, thanks for your response. The version of R was 2.6.1 and the version > of the package is the most recent 2.2-13. As mentioned, it works fine on > OS X but not on Windows XP. I will try using Rtools (from Duncan > Murdoch's si

Re: [R] replace numbers in a column conditional on their value

2008-01-16 Thread Richard . Cotton
> I have a data frame column in which I would like to replace some > of the numbers dependent on their value. > > data frame = zz > > AveExpr t P.Value FC > 7.481964 7.323950 1.778503e-04 2.218760 > 7.585783 12.233056 6.679776e-06 2.155867 > 6.953215 6.996525 2.353705e

Re: [R] Problem with running installed package on Windows

2008-01-16 Thread James Ferguson
Uwe, thanks for your response. The version of R was 2.6.1 and the version of the package is the most recent 2.2-13. As mentioned, it works fine on OS X but not on Windows XP. I will try using Rtools (from Duncan Murdoch's site) and install manually using the .zip file. Any other ideas? > Your

[R] Error in plot.new() : figure margins too large

2008-01-16 Thread Vidhu Choudhary
Hi I am getting error : Error in plot.new() : figure margins too large Can you please help -code is temp <- seq(550/2, 2200, 550/2) sV2 <- c(1, 1+temp[-length(temp)]) eV2 <- temp i=1 bitmap(paste("c:/vidhu/poster/poster_

Re: [R] replace numbers in a column conditional on their value

2008-01-16 Thread Erik Iverson
Try something like, ##Untested zz$FC <- ifelse(zz$FC > 2, zz$FC^2, zz$FC) and see ?ifelse -Erik Iverson Jabez Wilson wrote: > Dear R help, > > I have a data frame column in which I would like to replace some of the > numbers dependent on their value. > > data frame = zz > > A

[R] replace numbers in a column conditional on their value

2008-01-16 Thread Jabez Wilson
Dear R help, I have a data frame column in which I would like to replace some of the numbers dependent on their value. data frame = zz AveExpr t P.Value FC 7.481964 7.323950 1.778503e-04 2.218760 7.585783 12.233056 6.679776e-06 2.155867 6.953215 6.996525 2.35

[R] using an element of an array as a new object

2008-01-16 Thread ppaarrkk
I have an array called filesBox. I want to take each element of the first column and assign a dataset to it. For example : filesBox[4,1] returns [1] "fileR" Then I want to assign "fileR" which exists as a text file to the R object "fileR" like this : filesBox[4,1] <- read.table(paste(dir2, fi

[R] Supressing error messages on try function

2008-01-16 Thread Ignacio Perez
Hi! Why do I still get an error message on the try function?: > options(show.error.messages = FALSE) > > trab <-matrix(nrow=6,ncol=6) > trab[,] <-0.0001 > rownames(trab) <- c("A","B","C","D","E","I") > colnames(trab) <- c("A","B","C","D","E","I") > indice <- c(3,9,15,21,27,33) > for

Re: [R] Non linear regression with 2 explanatory variables

2008-01-16 Thread Roland Rau
Gavin Simpson wrote: > hits=-2.6 tests=BAYES_00 > X-USF-Spam-Flag: NO > > On Wed, 2008-01-16 at 11:02 +0100, Janice Kielbassa wrote: >> Hello! >> >> I want to do a non-linear regression with 2 explanatory variables >> (something like : length ~ a * time * exp( b* temperature)), having a >> data

Re: [R] R and Vista

2008-01-16 Thread Duncan Murdoch
On 1/16/2008 10:53 AM, Marcia Rocha wrote: > Hello Duncan, > Thank you for your answer. > I get a very strange path on my R.home(): > " C:\\PROGRA~2\\R\\R-25~1.0" > And the command > readLines(C:\\PROGRA~2\\R\\R-25~1.0\etc\Rprofile.site) returns an Error (a > syntax one, I triple checked the synta

Re: [R] Rprofile.site not executed?

2008-01-16 Thread Duncan Murdoch
On 1/16/2008 8:59 AM, Fernando Saldanha wrote: > It seems that my Rprofile.site file is not executed when I start R. To > test this I included the following code in that file: > > .First <- function(){ > cat("\nWelcome at", date(), "\n") > flush.console() > } > > When I start R the message abo

[R] exact method in coxph

2008-01-16 Thread Julian Wucherpfennig
I'm trying to estimate a cox proportional hazards regression for repeated events (in gap time) with time varying covariates. The dataset consists of just around 6000 observations (lines) (110 events). The (stylized) data look as follows: unit dur0 dur1 eventn event ongoing x 1 01

Re: [R] Rprofile.site not executed?

2008-01-16 Thread marciarr
I have exactly the same problem, just made the same question to the list 2 hours after you (your message was still not posted). Hope someone answers us! Fernando Saldanha wrote: > > It seems that my Rprofile.site file is not executed when I start R. To > test this I included the following code

Re: [R] R and Vista

2008-01-16 Thread Duncan Murdoch
On 1/16/2008 10:19 AM, marciarr wrote: > Hello R users, > I have been using R for a while on xp and have just updated my system to > vista. Have encountered some problems... First and main problem is that i > cannot make R recognize my own written functions by coping them into the > Rprofile.site o

[R] R and Vista

2008-01-16 Thread marciarr
Hello R users, I have been using R for a while on xp and have just updated my system to vista. Have encountered some problems... First and main problem is that i cannot make R recognize my own written functions by coping them into the Rprofile.site on the etc folder. For some reason, this very use

[R] how to calculate factor Scoring Coefficient s

2008-01-16 Thread Xingwang Ye
Dear R users, Recently, I am learning by myself to do factor analysis with R. Could some one guide me how to calculate 'factor scoring coefficients' or obtain 'factor scores' with the results from principal( ) or factor.pa( ) function in psych package? In SAS, it can be done with option "SCORE" b

Re: [R] affylmGUI

2008-01-16 Thread Prof Brian Ripley
See, from a few hours ago: https://stat.ethz.ch/pipermail/r-help/2008-January/150493.html On Wed, 16 Jan 2008, Marcos Pinho wrote: I have been trying to use the package affylmGUI but keep getting the error message that it cannot find the package "Bwidget". I have downloaded the activeTcl for

Re: [R] Using a data frame to create a legend

2008-01-16 Thread Henrique Dallazuanna
If I understand your question: x <- rnorm(100) plot(x) legend("topright", capture.output(t(t(summary(x) On 16/01/2008, Tom Snowden <[EMAIL PROTECTED]> wrote: > > Hi, > > I am currently in the process of translating a large script from S-Plus over > to R. It's going well so far and all the hi

Re: [R] Non linear regression with 2 explanatory variables

2008-01-16 Thread Gavin Simpson
hits=-2.6 tests=BAYES_00 X-USF-Spam-Flag: NO On Wed, 2008-01-16 at 11:02 +0100, Janice Kielbassa wrote: > Hello! > > I want to do a non-linear regression with 2 explanatory variables > (something like : length ~ a * time * exp( b* temperature)), having a > data set (length, time, temperature).

Re: [R] Rprofile.site not executed?

2008-01-16 Thread john seers (IFR)
Hi I am running under Windows XP using R2.6.1. I pasted your code in my Rprofile.site in the etc directory and it worked for me. All I can suugest is: 1. Check you really put it in the etc directory. (Do you have an old version of R somewhere ...?) 2. Is it really called Rprofile.site. (Not .

[R] Own classes in "histogram"

2008-01-16 Thread Denis Aydin
Hi, I try to make a histogram from a variable that contains the number of shoots from about 1000 individuals from a specific plant species (the range is 1-110). Those numbers are highly skewed to the right. My question is: how can I make my own classes with the lattice "histogram"? I tried it w

[R] affylmGUI

2008-01-16 Thread Marcos Pinho
I have been trying to use the package affylmGUI but keep getting the error message that it cannot find the package "Bwidget". I have downloaded the activeTcl for Windows and installed it in the default location C:\Tcl. The problem still persists. I have also tried to copy the files from C:\Tcl\l

[R] Rprofile.site not executed?

2008-01-16 Thread Fernando Saldanha
It seems that my Rprofile.site file is not executed when I start R. To test this I included the following code in that file: .First <- function(){ cat("\nWelcome at", date(), "\n") flush.console() } When I start R the message above is not displayed. I am running R 2.6.1 (rgui.exe) under Vista

Re: [R] Graphical output: dials and meters for a dashboard?

2008-01-16 Thread Parker Jones
> You can do a basic dial with something like this: > Thanks Jim. After installing plotrix I still had a bit of trouble getting this to work and I don't (yet) know enough R to fix it: > speedoplot(0.5,5,5,2,"Speedoplot") Error in polygon(xv, yv, border = border, col = col, lty = lty, lwd = lwd

Re: [R] Graphical output: dials and meters for a dashboard?

2008-01-16 Thread Parker Jones
> Here is a quick and dirty function to do so (things can easily be > changed to reflect preferences): > [...] Thanks Greg, that's working fine and is a fabulous starting point. Much appreciated! > Multiple dials can be created using par(mfrow= ...

[R] Using a data frame to create a legend

2008-01-16 Thread Tom Snowden
Hi, I am currently in the process of translating a large script from S-Plus over to R. It's going well so far and all the high level stuff is working fine, but as they say 'the devil is in the details' and my low level stuff is struggling. Specifically, I have enountered a problem translating fro

Re: [R] Graphical output: dials and meters for a dashboard?

2008-01-16 Thread Parker Jones
> > I'm new to R and am evaluating it to see whether it would be an appropriate > > tool to create a "dashboard" (a graphical statistical summary page). Could > > someone tell me if it is possible to display data as dials or meters. e.g. > > the four dials here: http://www.pentaho.com/images/

[R] between sum of squares for clusters

2008-01-16 Thread Serguei Kaniovski
Hello, how to compute the between sum of squares the for the clusters obtained by kmeans for the example below? I would like to compute the SSB-SSW based information measures for the clusters, as implemented in "cclust", but I am getting an error message. Serguei x1,x2,x3,x4,x5,cl 1.3,0.2,-1.2

Re: [R] Installing one channel GUI

2008-01-16 Thread Prof Brian Ripley
On Wed, 16 Jan 2008, Schraga Schwartz wrote: > Hello, > > I'm trying to install oneChannelGui on my computer. For this purpose, I've > tried to carefully follow the manual, and done the following: > > 1.I've installed the Affy power tools > 2.I've installed TCL from > http://www.activestat

Re: [R] Automatic traceback

2008-01-16 Thread Michael Hoffman
Prof Brian Ripley wrote: > On Wed, 16 Jan 2008, Michael Hoffman wrote: > >> I would like to have an automatic traceback printed on error. Is there >> a way to do this? > > There is a compact traceback printed automatically in non-interactive use, > and you can turn that on by the following options

Re: [R] help with reshaping data into long format (correct question)

2008-01-16 Thread Henrique Dallazuanna
try this: x[6, which(x[5,]=="y")] <- "y" levels(x$id) <- c(levels(x$id)[drop=T], "treat") x <- x[-5,] x[5, "id"] <- "treat" levels(x$id) <- gsub("^ques", "", levels(x$id)) x3 <- as.data.frame(t(x[,-1])) names(x3) <- x$id foo <- function(x, ...) { tmp <- as.numeric(as.character(unlist(x[,grep("_",

Re: [R] Automatic traceback

2008-01-16 Thread Prof Brian Ripley
On Wed, 16 Jan 2008, Michael Hoffman wrote: > I would like to have an automatic traceback printed on error. Is there > a way to do this? > > options(error=traceback) seems to print the previous error, not the > current error. I'm guessing this is because the traceback is not set > until the error

[R] Exact wilcoxon may differ in R and SPSS/StatXact (due to round off in the latter pair)

2008-01-16 Thread Nilsson Fredrik X
Dear R-users, If you use the exact Wilcoxon test in the coin package, I would like make you aware of that SPSS/StatXact MAY perform a round-off before doing their exact Wilcoxon-Mann-Whitney test (if you ever are unlucky enough not to use R). I have data from two treatments and was surprised to f

[R] Automatic traceback

2008-01-16 Thread Michael Hoffman
I would like to have an automatic traceback printed on error. Is there a way to do this? options(error=traceback) seems to print the previous error, not the current error. I'm guessing this is because the traceback is not set until the error handler is done running. > options(error=traceback) > s

Re: [R] Question concerning pdfs

2008-01-16 Thread Rainer M Krug
Thanks a lot - that gives the same result as the idea I had - only much easier. Rainer On Jan 16, 2008 12:18 PM, Dimitris Rizopoulos < [EMAIL PROTECTED]> wrote: > I think you need: > > dx2 <- exp(dx) > dp2 <- dp / dx2 > sum(c(0, diff(dx2)) * dp2) > > > I hope it helps. > > Best, > Dimitris > > -

Re: [R] Different results in calculating SD of 2 numbers

2008-01-16 Thread Heinz Tuechler
At 11:32 16.01.2008, Jim Lemon wrote: >(Ted Harding) wrote: > > On 16-Jan-08 08:45:04, Martin Maechler wrote: > > > >>>"RM" == Ron Michael <[EMAIL PROTECTED]> > >>>on Wed, 16 Jan 2008 00:14:56 -0800 (PST) writes: > >> > >>RM> Hi all, > >>RM> Can anyone tell me why I am getting d

[R] Installing one channel GUI

2008-01-16 Thread Schraga Schwartz
Hello, I'm trying to install oneChannelGui on my computer. For this purpose, I've tried to carefully follow the manual, and done the following: 1. I've installed the Affy power tools 2. I've installed TCL from http://www.activestate.com/Products/ActiveTcl/ (by following the link at t

Re: [R] plot labeling and expressions

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

[R] Analytical Optimization (Stat question)

2008-01-16 Thread francogrex
Dear Experts, this is more a general stat question, I tried to ask in other places but had no luck with answers (expect one that suggested numerical instead of analytical optimization): The likelihood below is a mixture of two negative binomial distributions: P*f(n;x1,x2,E) + (1-P)*f(n;x3,x4,E)

Re: [R] Quantile regression

2008-01-16 Thread Richard . Cotton
> I am trying to perform quantile regression analysis to analyse my work. I > could install the R package in windows xp. Now I am struggling > for the next work.I have *marks of students at the university > examinations*( say response variable Y) and their > *entrance examination marks* ( Independ

Re: [R] plot labeling and expressions

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

Re: [R] things that are difficult/impossible to do in SAS or SPSS but simple in R

2008-01-16 Thread Thomas Frööjd
As far as i know mixture modelling (sums of exponentials) cant be done in SAS or SPSS. For R there is the Rmix package that while not very user friendly at least works. On Jan 15, 2008 8:45 PM, Matthew Keller <[EMAIL PROTECTED]> wrote: > Hi all, > > I'm giving a talk in a few days to a group of ps

Re: [R] strange behaviour of is.factor()

2008-01-16 Thread Uwe Ligges
vito muggeo wrote: > Dear all, > It appears that the function is.factor() returns different results when > used inside the apply() function: that is, is.factor() fails to > recognize a factor.. > Where is the trick? > > many thanks, > vito > > > df1<-data.frame(y=1:10,x=rnorm(10),g=factor(c(

Re: [R] Problem with running installed package on Windows

2008-01-16 Thread Uwe Ligges
Your example works for me. Please read the posting guide that asks you to provide the following important information for us: Which version of R is this? Which version of the package? The recent version on CRAN master is: Package: clim.pact Version: 2.2-13 Date: January 10

[R] Derivation of data points

2008-01-16 Thread jessica . gervais
Dear All, I would like to make the derivation of a points. I tried to use D1tr (linear interpolation). I also tried to used D1ss (spline interpolation) but I do not realy know how to choose the smothing parameter the best. D1ss(x, y, xout = x, spar.offset = 0.1384, spl.spar=NULL) Here is an exam

Re: [R] Different results in calculating SD of 2 numbers

2008-01-16 Thread Jim Lemon
(Ted Harding) wrote: > On 16-Jan-08 08:45:04, Martin Maechler wrote: > >>>"RM" == Ron Michael <[EMAIL PROTECTED]> >>>on Wed, 16 Jan 2008 00:14:56 -0800 (PST) writes: >> >>RM> Hi all, >>RM> Can anyone tell me why I am getting different results in >>calculating SD of 2 numbers ?

[R] strange behaviour of is.factor()

2008-01-16 Thread vito muggeo
Dear all, It appears that the function is.factor() returns different results when used inside the apply() function: that is, is.factor() fails to recognize a factor.. Where is the trick? many thanks, vito > df1<-data.frame(y=1:10,x=rnorm(10),g=factor(c(rep("A",6),rep("B",4 > is.factor(df1

Re: [R] Question concerning pdfs

2008-01-16 Thread Dimitris Rizopoulos
I think you need: dx2 <- exp(dx) dp2 <- dp / dx2 sum(c(0, diff(dx2)) * dp2) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax:

[R] Non linear regression with 2 explanatory variables

2008-01-16 Thread Janice Kielbassa
Hello! I want to do a non-linear regression with 2 explanatory variables (something like : length ~ a * time * exp( b* temperature)), having a data set (length, time, temperature). Which function could I use (I tried nls but I think it doesn't work) Thanks a lot! Janice __

Re: [R] How to interrupt a loop by pressing a key?

2008-01-16 Thread mel
Bruno Jesus a écrit : > Does anyone know a way of interrupting a loop > by pressing a key (besides ctrl-c)? No hit key solution sorry, but may help : if (file.exists('zefile')) stop() allows to stop given zefile is created __ R-help@r-project.org mai

Re: [R] Different results in calculating SD of 2 numbers

2008-01-16 Thread Ted Harding
On 16-Jan-08 08:45:04, Martin Maechler wrote: >> "RM" == Ron Michael <[EMAIL PROTECTED]> >> on Wed, 16 Jan 2008 00:14:56 -0800 (PST) writes: > > RM> Hi all, > RM> Can anyone tell me why I am getting different results in > calculating SD of 2 numbers ? > > >> (1.25-0.95)/2

[R] Question concerning pdfs

2008-01-16 Thread Rainer M Krug
Hi I have a dataset of absolute growth rates g ranging from close to 0 to around whatever of which I want to calculate a pdf. If I use density(g) the pdf will extend to below 0 so I logtransform g and do d <- density(log(g)). Now I would like to transform this pdf back, i.e. d$x <- exp(d$x) but wh

Re: [R] Urgent : About R and 4D

2008-01-16 Thread Peter Danenberg
Quoth HAILE A. on Sweetmorn, Chaos 16, 3174: > Is it possible to interface R with 4D? > Thank you very much for a fast answer. Like the Law of Fives, it depends very much on the ingenuity of the user; how's that for a fast answer? __ R-help@r-project.or

Re: [R] Different results in calculating SD of 2 numbers

2008-01-16 Thread Martin Maechler
> "RM" == Ron Michael <[EMAIL PROTECTED]> > on Wed, 16 Jan 2008 00:14:56 -0800 (PST) writes: RM> Hi all, RM> Can anyone tell me why I am getting different results in calculating SD of 2 numbers ? >> (1.25-0.95)/2 RM> [1] 0.15 >> sd(c(1.25, 0.95)) RM> [1] 0.212

Re: [R] Different results in calculating SD of 2 numbers

2008-01-16 Thread Daniel Nordlund
> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf > Of Ron Michael > Sent: Wednesday, January 16, 2008 12:15 AM > To: [EMAIL PROTECTED] > Subject: [R] Different results in calculating SD of 2 numbers > > Hi all, > > Can anyone tell me why I am getting diff

[R] Urgent : About R and 4D

2008-01-16 Thread HAILE A.
Hello, Is it possible to interface R with 4D? Thank you very much for a fast answer. D. LEFLOCH eMail : [EMAIL PROTECTED] [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLE

  1   2   >