Re: [R] multiple separators in scan() functions

2008-07-30 Thread Prof Brian Ripley
On Thu, 31 Jul 2008, milton ruser wrote: Dear all, Is there a way of I deal with more than one kind of separators when reading a asc file using scan() function? The separators could be ":", "\t" and ";". Yes, pre-process the file (e.g. use a pipe connection and 'tr', or read the file with re

Re: [R] High resolution graphics from R Windows Vista

2008-07-30 Thread Prof Brian Ripley
On the TIFF point -- are you using 2.7.1 patched? It is likely that the crash is fixed there (it was intermittent and some of us never see it). Otherwise you have not described what is 'substandard' nor given any examples, and others are not reporting such differences. (I do have XP and Vist

Re: [R] predict rpart: new data has new level

2008-07-30 Thread Prof Brian Ripley
On Thu, 31 Jul 2008, Chua Siang Li wrote: Hi. I uses rpart to build a regression tree. Y is continuous. Now, I try to predict on a new set of data. In the new set of data, one of my x (call Incoterm, a factor) has a new level. I wonder why the error below appears as the guide say

[R] Generic plot function for GLM objects

2008-07-30 Thread Beat Huggler
Dear all, In R 2.7.1 on Windows it looks to me that the generic plot function for GLM objects uses standardized working residuals and not as labeled in the graph the standardized deviance residuals. It looks to me that from 2.7.0 to 2.7.1 there has been a bug introduced. Has anybody observ

[R] Rearranging Data

2008-07-30 Thread Jia Ying Mei
Hi everyone, A few days back Jim helped with out, helping me find an automated way of determining the mean and median of intervals between price changes. However, a new set of data I have is no longer in the same format. CDSmod is the file that was used with this code: x <- read.csv('/cdsmo

[R] stats question

2008-07-30 Thread Iasonas Lamprianou
Dear friends, I am not sure that this is the right place to ask,  but please feel free to suggest an alternative discussion group. My question is that I want to do a comparative study in order to compare the rate of incidence in two populations. I know that a pilot study was conducted a few wee

[R] High resolution graphics from R Windows Vista

2008-07-30 Thread Adam Langley
Hi I appreciate there has been some recent discussion regarding generating publication quality graphics in R. I am fairly familiar with these issues and have had no problems generating high quality images from R windows version. However, I have recently install Vista and what previously worked

[R] multiple separators in scan() functions

2008-07-30 Thread milton ruser
Dear all, Is there a way of I deal with more than one kind of separators when reading a asc file using scan() function? The separators could be ":", "\t" and ";". Kind regards miltinho [[alternative HTML version deleted]] __ R-help@r-project

[R] High resolution graphics from Windows Vista

2008-07-30 Thread Adam Langley
Hi I appreciate there has been some recent discussion regarding generating publication quality graphics in R. I am fairly familiar with these issues and have had no problems generating high quality images from R windows version. However, I have recently install Vista and what previously worked

[R] predict rpart: new data has new level

2008-07-30 Thread Chua Siang Li
Hi. I uses rpart to build a regression tree. Y is continuous. Now, I try to predict on a new set of data. In the new set of data, one of my x (call Incoterm, a factor) has a new level. I wonder why the error below appears as the guide says "For factor predictors, if an observ

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

2008-07-30 Thread Frank E Harrell Jr
Michal Figurski wrote: Tim, If I understand correctly, you are saying that one can't improve on estimating a mean by doing bootstrap and summarizing means of many such steps. As far as I understand (again), you're saying that this way one can only add bias without any improvement... Well, t

Re: [R] colnames from read.table

2008-07-30 Thread Erik Iverson
Chua Siang Li wrote: Hi. May I know why the colnames is NULL when reading only 1 column from a csv file? You are not "reading only 1 column from a csv file", you are subsetting one column of a data.frame. See ?Extract and the "drop" argument specifically. How do I get the colnam

Re: [R] S 3 generic method consistency warning please help

2008-07-30 Thread stephen sefick
changed data argument to x and it passes the check. FOR NOW. Thanks for all of the help On Wed, Jul 30, 2008 at 10:49 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > Mac OS X 10.5.4 > R 2.7.1 > the function works fine when I use it in a regular session. I can post > source code and the lik

Re: [R] S 3 generic method consistency warning please help

2008-07-30 Thread stephen sefick
Mac OS X 10.5.4 R 2.7.1 the function works fine when I use it in a regular session. I can post source code and the like if it is necessary to diagnose the problem. If ya'll need anything else just tell me. thanks Stephen On Wed, Jul 30, 2008 at 10:40 PM, stephen sefick <[EMAIL PROTECTED]> wrote

Re: [R] S 3 generic method consistency warning please help

2008-07-30 Thread Gabor Grothendieck
If you look at the generic formal arguments they look like this > args(window) function (x, ...) but your window method does not start with an argument called x so its inconsistent with its generic. On Wed, Jul 30, 2008 at 10:40 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > I would like to in

[R] S 3 generic method consistency warning please help

2008-07-30 Thread stephen sefick
I would like to include this in a package. The S3 methods on R CMD check says * checking S3 generic/method consistency ... WARNING window: function(x, ...) window.chron: function(data, day1, hour1, day2, hour2, ...) See section 'Generic functions and methods' of the 'Writing R Extensions' man

[R] colnames from read.table

2008-07-30 Thread Chua Siang Li
Hi. May I know why the colnames is NULL when reading only 1 column from a csv file? How do I get the colname then? Thanks. > xy = read.table("dataFile.csv",header=T, sep=",") > y <- xy[,1] > xDate <- xy[,2] > x <- xy[,3:8] > colnames(y) NULL > colnames(xDate) NULL

Re: [R] Re gression between adjacent columns - error with NAs

2008-07-30 Thread Krzysztof Sakrejda-Leavitt
A partial answer is... After doing debug(lm.fit) and debug(lm), and waiting to see where the messages come from, I can tell that when you do na.action=NULL, the R goes to call the Fortran routine dqrls (through .Fortran) and when R prepares the data to pass to the dqrls it triggers one of the

Re: [R] Sampling two exponentials

2008-07-30 Thread Moshe Olshansky
I am not sure that this is well defined. For a multivariate normal distribution (which is well defined), the covariance matrix (and the means vector) fully determine the distribution. In the exponential case, what is multivariate (bivariate) exponential distribution? I believe that knowing tha

Re: [R] Re gression between adjacent columns - error with NAs

2008-07-30 Thread Gabor Grothendieck
That's good. Try this: 1. put set.seed(1) at the top of the code to make it reproducible. 2. replace body of loop with: sel_col<-SourceMat[, i] out <- try(coef(lm(tt~sel_col, na.action=NULL))) if (!inherits(out, "try-error")) ResultMat[,i] <- out 3. email tends to wrap long

[R] curva PI

2008-07-30 Thread Bruno P Masi
Olá. Gostaria de uma ajuda sobre qual o pacote devo utilizar para realizar uma curva PI segundo Platt (1980). Achei um artigo que indica o pacote nlme, mas gostaria de saber se existem outros pacotes mais simples para modelos não lineares? Obrigado. Msc. Bruno Pereira Masi Doutorando em Ec

Re: [R] Sampling two exponentials

2008-07-30 Thread Ben Bolker
Zhang Yanwei - Princeton-MRAm munichreamerica.com> writes: > > Hi all, > I am going to sample two variables from two exponential distributions, but I want to specify a covariance > structure between these two variables. Is there any way to do it in R? Or is there a "Multivariate > Exponential"

[R] Is conditional evaluation of R code chunks possible in Sweave ?

2008-07-30 Thread Paul.Rustomji
Hi R mailing list Does anyone know of a good way to conditionally evaluate R code chunks when Sweaveing? What I'd kind of like to be able to do is set a number of variables in a list and if variable.a=TRUE then process code chunk labelled "methodA", if variable.b=TRUE then process code chunk labe

Re: [R] barplot for a categorial variable

2008-07-30 Thread Bill.Venables
Try with(yourData, barplot(table(Type))) or, if you really want percentages something like with(yourData, barplot(table(Type)/length(Type)*100, ylab = "Percentage")) (BTW, it would be nice to know who you are.) Bill Venables CSIRO Laboratories PO Box 120, Cleveland, 4163 AUSTRALIA Office Pho

Re: [R] Unexpected line type in lattice plot key on pdf device

2008-07-30 Thread Deepayan Sarkar
On Wed, Jul 30, 2008 at 6:41 AM, G. Draisma <[EMAIL PROTECTED]> wrote: > L.S., > With the code below, > on the Windows screen the line types in the key show > as solid and dashed as in the graph, > and in the pdf file they show > as solid in the key and solid and dashed in the graph. > I would no

[R] barplot for a categorial variable

2008-07-30 Thread kayj
Hi all, I have data with several columns and rows. one of the column is "Type" that can be either A, B, C or D. I need to plot the "Type" variable as a bar plot where the x axis is the type A, B ,C or D and the y axis shows the percentage. I was not able to get a bar plot. it seems that I need to

[R] Re gression between adjacent columns - error with NAs

2008-07-30 Thread rcoder
Hi Gabor, Thanks for your reply. I've written something that can be copied and pasted into your monitor to reproduce the error I am experiencing. Once the loop experiences a column full of NAs in SourceMat (column 3), it exits with errors, and ResultMat is only partially complete (up to column 2)

Re: [R] Lattice wireframe: How to avoid drawing lines around polygons when using shade=TRUE

2008-07-30 Thread Deepayan Sarkar
On Sun, Jul 27, 2008 at 1:41 AM, Oliver M. Haynold <[EMAIL PROTECTED]> wrote: > On Sun, 27 Jul 2008 05:00:48 +, Oliver M. Haynold wrote: >> I am using wireframe from the lattice package, with the shade option set >> to TRUE. When I output to PDF or Postscript, a line gets drawn around >> each

Re: [R] dotFirst option

2008-07-30 Thread Rolf Turner
On 31/07/2008, at 11:19 AM, Ferry wrote: Hi, I want to plot something the R opens. I created .First as follow: .First <- function() { plot(rnorm(100)) # do other stuff } But instead, I got error like below: Error in .First() : could not find function "plot" Any idea?

[R] dotFirst option

2008-07-30 Thread Ferry
Hi, I want to plot something the R opens. I created .First as follow: .First <- function() { plot(rnorm(100)) # do other stuff } But instead, I got error like below: Error in .First() : could not find function "plot" Any idea? Thanks. [[alternative HTML version deleted]]

[R] question about krige code in R

2008-07-30 Thread Alessandro
Hi All, I have this question. I wish to create a kriging map with R but I haven't a auxiliary map. I have only one txt file with X, Y and Z records. When I use this code: plot(variogram(Z~sqrt(X+Y), subground, cutoff=1800, width=10)) subground.or <- variogram(Z~sqrt(X+Y), subground,

Re: [R] john chambers R book

2008-07-30 Thread stephen sefick
amazon bonus bucks On Wed, Jul 30, 2008 at 4:05 PM, Marc Schwartz <[EMAIL PROTECTED]>wrote: > on 07/30/2008 01:27 PM [EMAIL PROTECTED] wrote: > >> Just to let everyone know: John Chambers R book is no longer at pre-order >> status on Amazon which is a good thing. The only problem >> is that now i

[R] Sampling two exponentials

2008-07-30 Thread Zhang Yanwei - Princeton-MRAm
Hi all, I am going to sample two variables from two exponential distributions, but I want to specify a covariance structure between these two variables. Is there any way to do it in R? Or is there a "Multivariate Exponential" thing corresponding to the multivariate normal? Thanks in advance.

Re: [R] setting editor environment variable EDITOR either when configuring R for installation or in .Rprofile

2008-07-30 Thread Mark Kimpel
Marc, That is exactly what I was looking for, although I am surprised that this variable does not appear to be configurable during installation from source as are some of the other environmental variables. Mark On Wed, Jul 30, 2008 at 2:04 PM, Marc Schwartz <[EMAIL PROTECTED]>wrote: > on 07/30/

Re: [R] Rprintf will not build in my C++ compiler

2008-07-30 Thread nmarti
Who cares that "An Introduction to the .C Interface to R" is not an "Official" manual. Does it matter? I was just trying to site where I got the code from and not claim it as my own. Would you prefer if I plagiarized stuff Professor? Sorry for not posting in devel, I thought this was a basic que

Re: [R] FFT - (STATS) - is this correct?

2008-07-30 Thread Rolf Turner
On 30/07/2008, at 11:14 PM, <[EMAIL PROTECTED]> wrote: Hello, I have calculated the fourier transform of the series enclosed at the end of this message, by doing: library(stats) x <- readLines("file1.txt") x.num <- as.numeric(x) ft.x.num <- fft(x.num) My question is: why is the first valu

Re: [R] john chambers R book

2008-07-30 Thread Marc Schwartz
on 07/30/2008 01:27 PM [EMAIL PROTECTED] wrote: Just to let everyone know: John Chambers R book is no longer at pre-order status on Amazon which is a good thing. The only problem is that now it is temporarily out of stock. There must have been many pre-orders so now I'm kicking myself for not ha

[R] Converting to subscripts and superscripts

2008-07-30 Thread naw3
Hi, I am reading in a CSV file of chemical reactions where the subscripts and superscripts are encoded in angle brackets, like below: 2 H+ + 2 O2- Is there a way to convert these to actual sub/superscripts and save them in another excel file? I greatly appreciate the help! Thanks, -Nina PS. I

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

2008-07-30 Thread Michal Figurski
Tim, If I understand correctly, you are saying that one can't improve on estimating a mean by doing bootstrap and summarizing means of many such steps. As far as I understand (again), you're saying that this way one can only add bias without any improvement... Well, this is in contradiction

Re: [R] read XML

2008-07-30 Thread Martin Morgan
"Paulo Cardoso" <[EMAIL PROTECTED]> writes: > Ok, > > I reach this with > setwd("E:/MBayRNP/output") > xml.mfd<-xmlTreeParse("myfile.xml") > str(xml.mfd$doc$children$layout) I usually parse but not edit XML, so the following might not be the best. But... > xml <- xmlRoot(xmlTreeParse("tmp.xml"))

Re: [R] john chambers R book

2008-07-30 Thread Christos Hatzis
I just received mine today (had pre-ordered it directly from Springer) -Christos > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of > [EMAIL PROTECTED] > Sent: Wednesday, July 30, 2008 2:27 PM > To: r-help@r-project.org > Subject: [R] john chambers

Re: [R] System exit codes

2008-07-30 Thread Ranjan Bagchi
Thank you -- I missed the bit about the invisible attribute. Can this happen with pipe, as well? x = readLines(pipe("perl -e 'print \"foo\\\n\";exit(-1);'")) x [1] "foo" I'd want to know that the command failed so I can process x differently. Thanks again, Ranjan Bagchi On Wed, 30 Jul 2008,

Re: [R] Re creating Procrustes Plot in Lattice

2008-07-30 Thread Gavin Simpson
On Wed, 2008-07-30 at 17:37 +0100, Gavin Simpson wrote: > On Wed, 2008-07-30 at 04:59 -0700, Andrewjohnclose wrote: > > Hi, I have been trying to create a function to generate a Procrustes plot, > > generated from package "vegan" > > in lattice. > > > > I can generate an empty plot, but not much

Re: [R] rolling regression between adjacent columns

2008-07-30 Thread Gabor Grothendieck
I was hoping for something reproducible. I can't run your code or make use of the data to save time in creating an example. Try this with builtin anscombe: > nc <- ncol(anscombe) > sapply(seq(1, nc, 2), function(i) coef(lm.fit(cbind(1, anscombe[,i]), > anscombe[,i+1]))) [,1] [,2] [,3]

Re: [R] History pruning

2008-07-30 Thread Marc Schwartz
on 07/30/2008 01:12 PM Ken Williams wrote: Hi, I find that a typical workflow for me looks something like this: 1) import some data from files 2) mess around with the data for a while 3) mess around with plotting for a while 4) get a plot or analysis that looks good 5) go back through my histor

Re: [R] System exit codes

2008-07-30 Thread Prof Brian Ripley
On Wed, 30 Jul 2008, Ranjan Bagchi wrote: Hi -- Is there a way of detecting the exit code while calling system or pipe? Yes, see the help page (you are ignoring the return value). print(system("perl -e 'print \"foo\\\n\";exit(-1);'")) foo [1] 65280 (= 256*255) eg (on Unix) system("per

Re: [R] Variable name in a String

2008-07-30 Thread Erik Iverson
?assign, ?get R_Learner wrote: Hi all, I have this string "year" and integer 2008 (both are inputs from the user), and I would like to make a variable called "year2008" that will store a vector of numbers. Does anyone know how to do this? Also, if the user later input "year" and "2008",

Re: [R] setting editor environment variable EDITOR either when configuring R for installation or in .Rprofile

2008-07-30 Thread Prof Brian Ripley
On Wed, 30 Jul 2008, Marc Schwartz wrote: on 07/30/2008 12:54 PM Mark Kimpel wrote: I'm running R on Linux and use emacs as my editor. When doing "edit(vignette(foo.vignette))" I would like to invoke emacs rather than the default vi. I am able to manually set this by editing $R_HOME/etc/Renviro

[R] System exit codes

2008-07-30 Thread Ranjan Bagchi
Hi -- Is there a way of detecting the exit code while calling system or pipe? eg (on Unix) system("perl -e 'print \"foo\\\n\";exit(-1);'") foo Any help appreciated, Ranjan __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/

[R] Setting fixed size for segement plot using stars() (axes size vs print size)

2008-07-30 Thread Jay Douillard
I have been making some segment plots with five variables. They work great, especially when I used a different scale function, which scaled them by area of the circle rather than radius scale <- function(x, Mr = 1 , Mx = 100) { ((x/Mx)^.5)*Mr} Where x is the the value, Mr is the Maximum radiu

Re: [R] Random subset

2008-07-30 Thread Marc Schwartz
on 07/30/2008 01:28 PM Marc Schwartz wrote: on 07/30/2008 01:18 PM Alessandro wrote: Hi all, I wish to do a random subset (i.e. 200 or 300 points) from a dataset, but I don't find the right code in R. See ?sample. Randomly select 10 elements from a 30 element vector: Vec <- 1:30 > sam

Re: [R] Random subset

2008-07-30 Thread Marc Schwartz
on 07/30/2008 01:18 PM Alessandro wrote: Hi all, I wish to do a random subset (i.e. 200 or 300 points) from a dataset, but I don't find the right code in R. See ?sample. Randomly select 10 elements from a 30 element vector: Vec <- 1:30 > sample(Vec, 10) [1] 16 5 10 29 27 30 21 1 12 28

[R] System exit codes

2008-07-30 Thread Ranjan Bagchi
Hi -- Is there a way of detecting the exit code while calling system or pipe? eg (on Unix) system("perl -e 'print \"foo\\\n\";exit(-1);'") foo Any help appreciated, Ranjan __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/

Re: [R] Random subset

2008-07-30 Thread Roland Rau
But probably you will need to set replace=FALSE if you want to sample from your original data set without replacement. Sorry for the possible confusion, Roland Roland Rau wrote: Hi, check ?sample n <- 200 mydata.set <- rnorm(10) my.random.subset <- sample(x=mydata.set, size=n, replace=T

Re: [R] rolling regression between adjacent columns

2008-07-30 Thread rcoder
Well, in this case I don't think my original code would have helped much... So, I've rewritten as below. I want to perform regression between one column in a matrix and all other columns in the same matrix. I have a for loop to achieve this, which succeeds in exporting intercept and slope coeffic

[R] Variable name in a String

2008-07-30 Thread R_Learner
Hi all, I have this string "year" and integer 2008 (both are inputs from the user), and I would like to make a variable called "year2008" that will store a vector of numbers. Does anyone know how to do this? Also, if the user later input "year" and "2008", how do I treat this as a variable;

[R] History pruning

2008-07-30 Thread Ken Williams
Hi, I find that a typical workflow for me looks something like this: 1) import some data from files 2) mess around with the data for a while 3) mess around with plotting for a while 4) get a plot or analysis that looks good 5) go back through my history to make a list of the shortest command sequ

Re: [R] Random subset

2008-07-30 Thread Roland Rau
Hi, check ?sample n <- 200 mydata.set <- rnorm(10) my.random.subset <- sample(x=mydata.set, size=n, replace=TRUE) my.random.subset I hope this helps, Roland Alessandro wrote: Hi all, I wish to do a random subset (i.e. 200 or 300 points) from a dataset, but I don't find the right

[R] Re gression between columns in matrix - with some/full NAs

2008-07-30 Thread rcoder
Hi Anna, Thanks for your reply and suggestions. I'm trying something different, based on regression this time, albeit giving similar problems. I want to regress a single column in a matrix with each other column in the same matrix, and o/p the intercept and slope coefficients to a results matrix.

Re: [R] Random subset

2008-07-30 Thread John Kane
?sample --- On Wed, 7/30/08, Alessandro <[EMAIL PROTECTED]> wrote: > From: Alessandro <[EMAIL PROTECTED]> > Subject: [R] Random subset > To: r-help@r-project.org > Received: Wednesday, July 30, 2008, 2:18 PM > Hi all, > > > > I wish to do a random subset (i.e. 200 or 300 points) from > a dat

[R] john chambers R book

2008-07-30 Thread markleeds
Just to let everyone know: John Chambers R book is no longer at pre-order status on Amazon which is a good thing. The only problem is that now it is temporarily out of stock. There must have been many pre-orders so now I'm kicking myself for not having done that. Mark

[R] Random subset

2008-07-30 Thread Alessandro
Hi all, I wish to do a random subset (i.e. 200 or 300 points) from a dataset, but I don't find the right code in R. Thanks for help alessandro [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.et

Re: [R] bug in 'margins' behavior in reshape - cast

2008-07-30 Thread hadley wickham
On Wed, Jul 30, 2008 at 12:41 PM, mfrumin <[EMAIL PROTECTED]> wrote: > > I want the column that is never going to actually have the '(all)' level to > not become of type factor. continuing the example: > > chick_m$diet = as.integer(as.character(chick_m$diet)) > is.factor(chick_m$diet) #returns FAL

Re: [R] setting editor environment variable EDITOR either when configuring R for installation or in .Rprofile

2008-07-30 Thread Marc Schwartz
on 07/30/2008 12:54 PM Mark Kimpel wrote: I'm running R on Linux and use emacs as my editor. When doing "edit(vignette(foo.vignette))" I would like to invoke emacs rather than the default vi. I am able to manually set this by editing $R_HOME/etc/Renviron but would like to avoid doing this with ea

[R] setting editor environment variable EDITOR either when configuring R for installation or in .Rprofile

2008-07-30 Thread Mark Kimpel
I'm running R on Linux and use emacs as my editor. When doing "edit(vignette(foo.vignette))" I would like to invoke emacs rather than the default vi. I am able to manually set this by editing $R_HOME/etc/Renviron but would like to avoid doing this with each install. I assume this can be accomplishe

Re: [R] About clustering techniques

2008-07-30 Thread pacomet
Hi Christian I've been reading about daisy and think I need to do something like.. > mydaisydata <- daisy(mydata,metric=c("euclidean"),stand=FALSE) Error en vector("double", length) : tamaño del vector especificado es muy grande(which means, specified vector size is too big) mydata is an

Re: [R] best quad core configuration?

2008-07-30 Thread Vincent Goulet
Le mer. 30 juil. à 13:04, Alan Swanson a écrit : I am purchasing a quad-core 64-bit system. My proposed configuration is an intel core duo quad with 8Gb memory and an ubuntu OS. I'm choosing ubunto because I have good local support, but it doesn't seem as well supported as debian for R.

Re: [R] bug in 'margins' behavior in reshape - cast

2008-07-30 Thread mfrumin
I want the column that is never going to actually have the '(all)' level to not become of type factor. continuing the example: chick_m$diet = as.integer(as.character(chick_m$diet)) is.factor(chick_m$diet) #returns FALSE is.factor(cast(subset(chick_m, time == 0), diet + chick ~ time, mean, margin

Re: [R] Mixed effects model where nested factor is not the repeated across treatments lme???

2008-07-30 Thread Mark Difford
Hi Miki and Chunhao, << Rusers (Anna, and Mark {thank you guys}) provide me a vary valuable << information. Also see Gavin Simpson's posting earlier today: apparently multcomp does now work with lmer objects (it's gone through phases of not working, then working: it's still being developed). B

Re: [R] Re creating Procrustes Plot in Lattice

2008-07-30 Thread Dylan Beaudette
On Wednesday 30 July 2008, Gavin Simpson wrote: > On Wed, 2008-07-30 at 04:59 -0700, Andrewjohnclose wrote: > > Hi, I have been trying to create a function to generate a Procrustes > > plot, generated from package "vegan" > > in lattice. > > > > standard vegan code as follows > > library(vegan) > >

[R] best quad core configuration?

2008-07-30 Thread Alan Swanson
I am purchasing a quad-core 64-bit system. My proposed configuration is an intel core duo quad with 8Gb memory and an ubuntu OS. I'm choosing ubunto because I have good local support, but it doesn't seem as well supported as debian for R. Will ubuntu limit me in any significant way? I have not

Re: [R] odds ratios in multiway tables (stratified)

2008-07-30 Thread Marc Schwartz
on 07/30/2008 08:48 AM Wim Bertels wrote: Hi, does anyone know of a function to calculate odds ratios in multiway tables (stratified) (+ the other usual statistics involved) i mean: say we have a table r*c*d, For every d (depth) we have a r*c table, and in this table the odds ratio's are calcul

Re: [R] bug in 'margins' behavior in reshape - cast

2008-07-30 Thread hadley wickham
On Wed, Jul 30, 2008 at 11:47 AM, mfrumin <[EMAIL PROTECTED]> wrote: > > according to the documentation of the cast function in the reshape function, > I would expect this bit of code from the examples to calculate marginal > means over only the 'diet' variable. > > #Chick weight example > names(Ch

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread Gavin Simpson
On Wed, 2008-07-30 at 21:05 +0530, [EMAIL PROTECTED] wrote: > now which license to use? Use whatever you or your company want to use. As long as it is in compliance with GPL2 (as that is the licence R currently uses) then you are fine. As to whether what you intend to do is in compliance with the

[R] Hello,

2008-07-30 Thread Kiran Annaiah
Hello, A newbie to R. I am trying to use the exonmap package in R. According to the docs, the xmapDatabase() command should read the config file with all the connection parameters and connect to the DB. But i get the error shown below > xmapDatabase("Human") Error in mysqlNewConnection(dbDri

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread Gavin Simpson
On Wed, 2008-07-30 at 09:50 -0400, Stevens, Martin Henry H. Dr. wrote: > The limit of use is that you cannot profit from it. It is Open Source. I think you mean you are unlikely to profit from it if you base your business model on "selling" the "product" and the product is the code; because it is

[R] bug in 'margins' behavior in reshape - cast

2008-07-30 Thread mfrumin
according to the documentation of the cast function in the reshape function, I would expect this bit of code from the examples to calculate marginal means over only the 'diet' variable. #Chick weight example names(ChickWeight) <- tolower(names(ChickWeight)) chick_m <- melt(ChickWeight, id=2:4,

[R] OT: Data mining on demand .Using cloud computers.

2008-07-30 Thread Ajay ohri
Dear Lists, There are many softwares (with many passionate adherents). Some of them are very good upto say 2 million rows but lose steam after that despite having very good focus in innovating algorithms. Some of them are very good on same PC for much more data but become comfortable with the copy

Re: [R] Re creating Procrustes Plot in Lattice

2008-07-30 Thread Gavin Simpson
On Wed, 2008-07-30 at 04:59 -0700, Andrewjohnclose wrote: > Hi, I have been trying to create a function to generate a Procrustes plot, > generated from package "vegan" > in lattice. > > standard vegan code as follows > library(vegan) > pro=protest(P1, P8, permutations=4999,choices=1:4) > plot(pro)

Re: [R] R -Legality Question about R's Open Source GNU GPL Licen

2008-07-30 Thread Csardi Gabor
Ted, yes, there were a couple of cases, both in Europe and in the US. See more here: http://www.softwarefreedom.org/ As far as I can remember, all lawsuits ended with an agreement, and the publication of the source code, of course. Best, Gabor On Wed, Jul 30, 2008 at 05:19:31PM +0100, Ted Har

Re: [R] ~ . help

2008-07-30 Thread Michael Dewey
At 08:44 30/07/2008, [EMAIL PROTECTED] wrote: Hello, What the tilde and dot means? I believe is something to do with formula, but how do they work? You are already half way there ?formula should help you find out what ~ means ?update.formula should help you find out what . means in the contex

[R] FOURIER TRANSFORM HELP

2008-07-30 Thread renata.camargo05
Hello, I have a series and I need to run a Fourier Transform for that series. I have done that using the function fft from stats package. However I am not sure whether the result I am getting is correct or not. Seems that the first value of the Fourier Transform list is the sum of all elements

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread ohri2007
well, i am planning to test rattle (r gui) or even R on a cloud computer. IMO SAS (software ) scores over R because it handles bigger files more efficiently on SAME local hardware. but prices a huge premium on licenses. with R and a cloud, you pay for time spent (1+profit P). cloud ramps up proce

Re: [R] R -Legality Question about R's Open Source GNU GPL Licen

2008-07-30 Thread Ted Harding
In this midst of all this debate, can I try a simple question? Have the enforceability and/or interpretation of the GPL ever been tested in court? If so, then are there any references to where one can read about the case (or cases)? With thanks, Ted.

Re: [R] Rprintf will not build in my C++ compiler

2008-07-30 Thread Prof Brian Ripley
That is a linker (not compiler) message, and you apparently forgot to link against R.dll. See README.packages in the R distribution for how to use VC++ with R under Windows (which is neither recommended nor supported). Three other comments. 1) See the posting guide. R-devel is the list for n

[R] Axes in perspective plot - persp()

2008-07-30 Thread David Afshartous
All, I'm making some plots with the persp() function and am having trouble sorting out the following: - How does one avoid the axes label overlapping the tick values? - Is doesn't seem possible to independently control the number of ticks of the x,y, and z-axes, e.g., I'd like say only 4 ticks

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread Marc Schwartz
on 07/30/2008 10:35 AM [EMAIL PROTECTED] wrote: now which license to use? it is bad enough explaining SaaS (software as a service), and open source and data mining to lawyers, without the experts fighting it as well. and who pays for the cloud computing expenses for me then if there is no mone

[R] Rprintf will not build in my C++ compiler

2008-07-30 Thread nmarti
I have searched through the threads and "Rprintf" causing a build error dosen't seem to be a problem for anyone else. And I've read through "R Extensions" and "An Introduction to the .C Interface to R" and there doesn't seem to be any troubleshooting for my problem. My code is straight from "An

Re: [R] Mixed effects model where nested factor is not the repeated across treatments lme???

2008-07-30 Thread ctu
Hi Miki, I just got the same problem with you couple hours ago. Rusers (Anna, and Mark {thank you guys}) provide me a vary valuable information. link to following address. http://www.nabble.com/Tukey-HSD-(or-other-post-hoc-tests)-following-repeated-measures-ANOVA-td17508294.html#a17559307 for

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread ohri2007
now which license to use? it is bad enough explaining SaaS (software as a service), and open source and data mining to lawyers, without the experts fighting it as well. and who pays for the cloud computing expenses for me then if there is no money in it. cloud computing doesnt exist in India rig

[R] Percent Change After Each Interval

2008-07-30 Thread Jia Ying Mei
Hi, I have a data set (attached, and which I've asked about before) containing countries, prices and dates. I am trying to find the mean and median (by country) percentage change of each price change. For example: Date Australia 6/30/07 2 3/31/07

Re: [R] adding lines to multiple plot

2008-07-30 Thread Ben Bolker
Gary W gmail.com> writes: > > > Hi ben > > Thanks for that, the density lines have now worked, however i was hoping > that the labels would be the names of the columns in my original data set > not numeric. any advice how to do this? > > gary w > mapply(function(v,n) { hist(v,prob=TRUE

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread Roland Rau
Hi, Stevens, Martin Henry H. Dr. wrote: The limit of use is that you cannot profit from it. It is Open Source. I am not sure what you mean exactly. If you mean that you can not earn any money selling free software (in terms of the GPL2, which is R's licence), this seems to be not the case:

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread Marc Schwartz
Absolutely FALSE. Read this FAQ and the following 3: http://www.fsf.org/licensing/licenses/gpl-faq.html#DoesTheGPLAllowMoney Let's not spread FUD. Marc Schwartz on 07/30/2008 08:50 AM Stevens, Martin Henry H. Dr. wrote: The limit of use is that you cannot profit from it. It is Open Source. _

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread Stevens, Martin Henry H. Dr.
The limit of use is that you cannot profit from it. It is Open Source. From: [EMAIL PROTECTED] [EMAIL PROTECTED] On Behalf Of Ajay ohri [EMAIL PROTECTED] Sent: Wednesday, July 30, 2008 9:39 AM To: r-help@r-project.org Subject: [R] R -Legality Question about

Re: [R] Most often pairs of chars across grouping variable

2008-07-30 Thread svga
Hi Marc, many thanks, that is exactly what I was looking for. Best, Sven - Original Nachricht Von: Marc Schwartz <[EMAIL PROTECTED]> An: [EMAIL PROTECTED] Datum: 29.07.2008 17:15 Betreff: Re: [R] Most often pairs of chars across grouping variable > on 07/29/2008 09:51 AM [

Re: [R] adding lines to multiple plot

2008-07-30 Thread Gary W
Hi ben Thanks for that, the density lines have now worked, however i was hoping that the labels would be the names of the columns in my original data set not numeric. any advice how to do this? gary w Ben Bolker wrote: > >> >> ##find unique column names in table >> variablenames<-unique(n

[R] model mix problem. FALSE CONVERGENCE

2008-07-30 Thread mergullo
Hi R users I try to use the lme but I can´t! My script is (some words in french, sorry!!): rm(list=ls(all=TRUE)) #Efface tous les objets en mémoire pour éviter des erreurs library(MASS) #Chargement des Librairies library(car) library(Hmisc) library(tkWidgets) library(svDialogs) librar

Re: [R] R -Legality Question about R's Open Source GNU GPL License

2008-07-30 Thread Marc Schwartz
on 07/30/2008 08:39 AM Ajay ohri wrote: whats the license policy of using r as a SaaS with other apps, priced for a fee. I am not sure which license to invoke. regards, ajay See this post: https://stat.ethz.ch/pipermail/r-help/2008-June/165906.html HTH, Marc Schwartz _

Re: [R] adding lines to multiple plot

2008-07-30 Thread Ben Bolker
Gary W gmail.com> writes: > > > Hi > > I am trying to add density lines to multiple histograms created using R. I > have managed to get all 18 lines on the last histogram in the series, > however, i cannot figure out how to get each line appear on the correct > histogram plot. I know in Mat

  1   2   >