Re: [R] Re peated Measure ANOVA-Old question

2008-07-30 Thread Mark Difford
Hi Chunhao, I google the website and I found that there are three ways to perform repeated measure ANOVA: aov, lme and lmer. It's also a good idea to search through the archives. I use the example that is provided in the above link and I try

Re: [R] Urgent

2008-07-30 Thread Moshe Olshansky
Hi Yunlei, Is your problem constrained or not? If it is unconstrained and your matrix is not positive definite, the minimum is unbounded (unless you are extremely lucky and the matrix is positive semi-definite and the vector which multiplies the unknowns is exactly perpendicular to all the

Re: [R] Re peated Measure ANOVA-Old question

2008-07-30 Thread ctu
Tank you (Anna Mark) very much for this super information Let me confirm that IF I want to perform the RM-ANOVA I should use lmer and perform the post-hoc test by using glht, right? Because I am not so familiar with lmer so I have two more questions. Let me use that following example again, 1.

Re: [R] R command history -- can it be like Matlab's?

2008-07-30 Thread Martin Maechler
l == losemind [EMAIL PROTECTED] on Tue, 29 Jul 2008 12:23:36 -0700 (PDT) writes: l Hi all, l In R GUI window, if you use up and down key, you will be able to recall l the previous and next command that has been used and stored in the command l history cache. l But

Re: [R] What's the best way to operate R from within Excel?

2008-07-30 Thread Patrizio Frederic
have a look at http://sunsite.univie.ac.at/rcom/server/doc/RExcel.html regards, PF 2008/7/30 losemind [EMAIL PROTECTED]: Hi all, How do I operate R from within Excel? I mean, highlight a bunch of cells, and send to R, and do some statistics in R, and return back the numbers and do some

[R] problem with nested loops

2008-07-30 Thread Jens Oldeland
Dear all, I have a problem with constructing a nested loop. I have two matrices: pnr: 800 rows 14 columns where rows are 40x20 meaning that 40 rows belong to one of twenty objects in the matrix pnr mvp: 20 rows and 14 columns I want to: calculate a distance value with the first 40 rows of

[R] ~ . help

2008-07-30 Thread renata.camargo05
Hello, What the tilde and dot means? I believe is something to do with formula, but how do they work? An example: (fk5 ~ ., size = 10) Do you know what the fk5 ~ . means? Thank you! Renata - Email sent from www.virginmedia.com/email Virus-checked

Re: [R] Re peated Measure ANOVA-Old question

2008-07-30 Thread Mark Difford
Hi Chunhao, If you carefully read the posting that was referred to you will see that lme() and not lmer() was used as an example (for using with the multcomp package). lmer() was only mentioned as an aside... lmer() is S4 and doesn't work with multcomp, which is S3. Apropos of specifying random

Re: [R] Coarsening the Resolution of a Dataset

2008-07-30 Thread Steve Murray
Hi - thanks for the advice - I am however applying this to the whole data frame. And the code that I'm using is just to read in the data (using read.table) and then the code that you supplied. I could send you the actual dataset if you don't mind a file ~50MB?! Thanks again, Steve Date:

Re: [R] Coarsening the Resolution of a Dataset

2008-07-30 Thread Steve Murray
Hi - thanks for the advice - I am however applying this to the whole data frame. And the code that I'm using is just to read in the data (using read.table) and then the code that you supplied. I could send you the actual dataset if you don't mind a file ~50MB?! Thanks again, Steve Date:

[R] Eclipse/Statet: How to set breakpoints

2008-07-30 Thread [EMAIL PROTECTED]
Hi, I'm using the IDE Eclipse with the Statet-Plugin to develop R-code. Is there a possibility to set breakpoints for debugging in R-script files (just like in Java-code)? Cheers, Eli __ R-help@r-project.org mailing list

Re: [R] R command history -- can it be like Matlab's?

2008-07-30 Thread Gad Abraham
Prof Brian Ripley wrote: A patch to do this was posted on 2007-09-29 by Glenn Davis. Some people not addicted to Matlab find the behaviour very inconvenient and prefer the getline/readline behaviour (triggered by ^R/^S) of Rterm and R on Unixen. Now that you've mentioned readline, it

[R] need some help

2008-07-30 Thread Kurapati, Ravichandra (Ravichandra)
Hi Can any body explain the meaning of the following R code if(interactive()) stop(msg, call. = FALSE) else quit() Thanks K.Ravichandra [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] Howto Draw Bimodal Gamma Curve with User Supplied Parameters

2008-07-30 Thread Patrizio Frederic
can you write the bimodal gamma pdf? if so create your own pdf: dbigamma=function(x,alpha1,alpha2,beta1,beta2) { # ... the bimodal gamma pdf (I can't find it) } hist(x,probability=T) curve(dbigamma(x,alpha1,alpha2,beta1,beta2)) and look at the result regards, PF 2008/7/29 Gundala Viswanath

[R] Need help

2008-07-30 Thread renata.camargo05
Hello, Can someone help me to understand the meaning of the following R line? list(fk5 ~ .) Thanks, Renata - Email sent from www.virginmedia.com/email Virus-checked using McAfee(R) Software and scanned for spam

Re: [R] R command history -- can it be like Matlab's?

2008-07-30 Thread bartjoosen
Here is the post where Prof Ripley is reffering to: http://tolstoy.newcastle.edu.au/R/e2/help/07/09/26668.html Bart losemind wrote: Prof Brian Ripley wrote: A patch to do this was posted on 2007-09-29 by Glenn Davis. Some people not addicted to Matlab find the behaviour very

Re: [R] Need help

2008-07-30 Thread Richard . Cotton
Can someone help me to understand the meaning of the following R line? list(fk5 ~ .) fk5 ~. is a formula, most likely used by a regression (or similar model). fk5 is the response variable, and the dot is a placeholder for 'whatever was on that side of the formula before'. (This is usually

Re: [R] Re peated Measure ANOVA-Old question

2008-07-30 Thread Gavin Simpson
On Wed, 2008-07-30 at 00:49 -0700, Mark Difford wrote: Hi Chunhao, If you carefully read the posting that was referred to you will see that lme() and not lmer() was used as an example (for using with the multcomp package). lmer() was only mentioned as an aside... lmer() is S4 and doesn't

Re: [R] Coarsening the Resolution of a Dataset

2008-07-30 Thread jim holtman
Can you at least send the output from the script so that we can see exactly how you are reading in the data and then processing it. Doing 'str(yourData)' would help to understand what the data looks like. The exactly error messages would also help. Is the data entirely numeric, or does it also

Re: [R] need some help

2008-07-30 Thread jim holtman
I assume that it is testing to see if you are in the interactive mode (typically using RGUI or Rterm) and not a batch file. If it is interactive, it issues a 'stop' (this lets the person take action on the problem). If in the batch, it just terminates the execution of the script. On Wed, Jul

Re: [R] problem with nested loops

2008-07-30 Thread Kenn Konstabel
What you need is chapter 9 in http://cran.r-project.org/doc/manuals/R-intro.pdf . And you could also use chapter 2, especially 2.2. specdist -matrix(NA,nrow=40,ncol=20) for (j in 1:40){ for(i in 1:20){

[R] read XML

2008-07-30 Thread Paulo Cardoso
I have a xml exported by Manifold GIS but I'm not being able to import it into R using XLM package. The file have this structure: ?xml version=1.0 encoding=UTF-8 ? - layout nameLayout 2/name pagesByX1/pagesByX pagesByY1/pagesByY - elements legend refParent=r3

[R] function to transform response of a formula

2008-07-30 Thread Paul Emberson
Hi, I am trying to write a function which takes a formula as input and outputs a new formula with a different response while keeping the rest of the formula the same. I.e. respapply : ( y ~ a+b ) - ( f(y) ~ a + b ) I have tried the following but it doesn't work. The terms become invalid

Re: [R] problem with read.table()

2008-07-30 Thread Birgitle
I am sorry I just found the stupid mistake. I did not specify dec=, because I usually use . Anyway thanks for having the opportunity to get help in this list. B. Birgitle wrote: Hello R-User I have a table as tab-delimited textfile (291 rows, 83 columns). The first row are labels and

Re: [R] problem with read.table()

2008-07-30 Thread Philipp Pagel
Table-read.table(Table.txt,header=TRUE,row.names=1, na.strings=NA ,colClasses = Classe82) Fehler in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : scan() erwartete 'a real', bekam '8,5' Error in scan(file, what, nmax, sep, dec, quote, skip, nlines,

[R] Bizarre - R crashes on merge

2008-07-30 Thread Williams, Robin
Hi all, I have a large data.frame, 1530 observation with 6 columns. I want to merge a 7th column, a transformation of the response variable (hospital admissions), namely trans-sqrt(copd$admissions+0.25) trans-data.frame(trans) And now when I do copd2-merge(copd,trans) (copd being my original

Re: [R] Need help

2008-07-30 Thread Prof Brian Ripley
See ?formula ?terms.formula or ask the author of the line (it doesn't make much sense to use '.' unless you have a 'data' argument to pass to terms()). On Wed, 30 Jul 2008, [EMAIL PROTECTED] wrote: Hello, Can someone help me to understand the meaning of the following R line? list(fk5 ~

Re: [R] Need help

2008-07-30 Thread Prof Brian Ripley
On Wed, 30 Jul 2008, [EMAIL PROTECTED] wrote: Can someone help me to understand the meaning of the following R line? list(fk5 ~ .) fk5 ~. is a formula, most likely used by a regression (or similar model). fk5 is the response variable, and the dot is a placeholder for 'whatever was on that

Re: [R] Eclipse/Statet: How to set breakpoints

2008-07-30 Thread Duncan Murdoch
On 30/07/2008 4:28 AM, [EMAIL PROTECTED] wrote: Hi, I'm using the IDE Eclipse with the Statet-Plugin to develop R-code. Is there a possibility to set breakpoints for debugging in R-script files (just like in Java-code)? I believe the answer is no, but R does contain (some of?) the necessary

Re: [R] Bizarre - R crashes on merge

2008-07-30 Thread Philipp Pagel
I have a large data.frame, 1530 observation with 6 columns. I want to merge a 7th column, a transformation of the response variable (hospital admissions), namely trans-sqrt(copd$admissions+0.25) trans-data.frame(trans) And now when I do copd2-merge(copd,trans) (copd being my original

Re: [R] table questions

2008-07-30 Thread Jim Lemon
On Tue, 2008-07-29 at 11:28 -0500, Edna Bell wrote: Hi again! Suppose I have the following: xy - round(rexp(20),1) xy [1] 0.1 3.4 1.6 0.4 1.0 1.4 0.2 0.3 1.6 0.2 0.0 0.1 0.1 1.0 2.0 0.9 2.5 0.1 1.5 0.4 table(xy) xy 0 0.1 0.2 0.3 0.4 0.9 1 1.4 1.5 1.6 2 2.5 3.4 1 4 2

[R] glm error: factor levels

2008-07-30 Thread Birgitle
Hello R-Users, again me! I have a data.frame with 291 rows, 82 columns. Tha variables in the columns are factors, numerics and ordered factors. The response variable is a factor with two levels. I would like to find the best model by trying every possible variable combination using a logistic

Re: [R] function to transform response of a formula

2008-07-30 Thread email
Its ok. I've just read about update.formula in another message. Paul Paul Emberson wrote: Hi, I am trying to write a function which takes a formula as input and outputs a new formula with a different response while keeping the rest of the formula the same. I.e. respapply : ( y ~ a+b ) -

Re: [R] Bizarre - R crashes on merge

2008-07-30 Thread Philipp Pagel
On Wed, Jul 30, 2008 at 12:24:22PM +0200, Philipp Pagel wrote: If I understand correctly, all you want to do is add another column to your data.frame. Unless I have overlooked something, a simple assignement should do: copd$trans - copd$admissions+0.25 Sorry, that should have been:

[R] png() does not generate a 24 bit PNG file

2008-07-30 Thread Fán Lóng
Hi there, When I am using the png() function as the following to generate a png file on Mac OSX Tiger, I found that the bit depth is 32. But I want a 24 bit depth png. What should I do? png(filename=/Users/usg/Desktop/new.png) plot(1:20) dev.off() To avoid using cairo, and R-2.7.0 does not

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

2008-07-30 Thread renata.camargo05
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 value (Real) of ft.x.num that big? (954.833870) all the other

Re: [R] placing a text in the corner of a plot

2008-07-30 Thread Jim Lemon
On Wed, 2008-07-30 at 00:04 +0200, Jörg Groß wrote: Is there an easy way to add a text into an R-plot and place it in the upper left corner? Like that: | (a) | | | |

Re: [R] Eclipse/Statet: How to set breakpoints

2008-07-30 Thread Duncan Murdoch
On 7/30/2008 7:15 AM, [EMAIL PROTECTED] wrote: Hi Duncan, thank you for your fast answering. Is there any function implemented in R to stop it? (something like q() but i only want to stop a script not shutdown the whole R-console at a certain point) That would be a good alternative for

[R] dens function

2008-07-30 Thread Shubha Vishwanath Karanth
Hi R, Is there a function called dens in R? Is it similar to density? I found this function in an RExcel demo... Or is it in a specific other package? The function is used as (=rapply(dens,A11:A1010,100,D8,-5,5)) Thanks, Shubha This e-mail may contain confidential and/or privileged

[R] Connecting R to PostgreSQL on Windows.

2008-07-30 Thread Andre Cesta
Hi All, It took me 30minutes to connect GNU R on Windows XP to a PostgreSQL database running on the same machine and then generate the first graphics. I've taken screenshots of all the steps and documented them here: http://code.google.com/p/jwanalytics/wiki/GnuRe Thanks for making it easy

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

2008-07-30 Thread M Ensbey
Hi, I have searched the archives and can't quite confirm the answer to this. I appreciate your time... I have 4 treatments (fixed) and I would like to know if there is a significant difference in metal volume (metal) between the treatments. The experiment has 5 blocks (random) in each

[R] Successfully connected GNU R to PostgreSQL on windows XP.

2008-07-30 Thread Andre Cesta
Hi All, I guess it is also important to write when you did not have any problems doing something. This morning it took me 30minutes to connect GNU R on Windows XP to a PostgreSQL database running on the same machine and then generate the first graphics. I've taken screenshots of all the

[R] R GUI help

2008-07-30 Thread Tom.O
Hi R gurus Im currently doing some GUI coding with the rpanel package. I have come across a very interesting problem, how do I close a panel that I have created. This sounds silly for I can always click the X up in the right corner but how do I do it by code? for example; when I create a control

[R] adding lines to multiple plots

2008-07-30 Thread Gary W
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 Matlab there is a hold function, is there

[R] Re creating Procrustes Plot in Lattice

2008-07-30 Thread Andrewjohnclose
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) Now, here is the code for the function that I have failed to get to

Re: [R] read XML

2008-07-30 Thread Martin Morgan
Paulo Cardoso [EMAIL PROTECTED] writes: I have a xml exported by Manifold GIS but I'm not being able to import it into R using XLM package. I can library(XML) xml - xmlTreeParse(tmp.xml, useInternal=TRUE) xpathApply(xml, //component/@ymin, xmlValue) [[1]] [1] 0.32371064328183174 [[2]] [1]

Re: [R] function to transform response of a formula

2008-07-30 Thread Paul Emberson
Hi, I am now trying to use update.formula from within a function but I am not able to get the result I want. Here is an example f - function(fm, n) { update(fm, n*. ~ .); } If I run fm - formula( y ~ a + b ) f(fm,3) n * y ~ a + b Now when I try to use the formula, n does not exist.

[R] adding lines to multiple plot

2008-07-30 Thread Gary W
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 Matlab there is a hold function, is there

Re: [R] What's the best way to operate R from within Excel?

2008-07-30 Thread losemind
HI all, Is that the best R link from within Excel currently? Thanks! Patrizio Frederic wrote: have a look at http://sunsite.univie.ac.at/rcom/server/doc/RExcel.html regards, PF 2008/7/30 losemind [EMAIL PROTECTED]: Hi all, How do I operate R from within Excel? I mean,

Re: [R] R command history -- can it be like Matlab's?

2008-07-30 Thread losemind
Prof Brian Ripley wrote: A patch to do this was posted on 2007-09-29 by Glenn Davis. Some people not addicted to Matlab find the behaviour very inconvenient and prefer the getline/readline behaviour (triggered by ^R/^S) of Rterm and R on Unixen. On Tue, 29 Jul 2008, losemind wrote:

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

2008-07-30 Thread Ajay ohri
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 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

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

2008-07-30 Thread G. Draisma
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 not expect that, but may be I get something wrong. Gerrit. ---

Re: [R] read XML

2008-07-30 Thread Paulo Cardoso
Ok, I reach this with setwd(E:/MBayRNP/output) xml.mfd-xmlTreeParse(myfile.xml) str(xml.mfd$doc$children$layout) and get this: layout nameLayout 2/name pagesByX1/pagesByX pagesByY1/pagesByY elements legend refParent=r3 xmax=0.54156171284634758 xmin=0.020151133501259435

[R] odds ratios in multiway tables (stratified)

2008-07-30 Thread Wim Bertels
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 calculated for every 2*2 subtable in it.

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

2008-07-30 Thread stephen sefick
#ts with a freq=1 d - ts(x) #it seems that this is an oscillating signal of ~ a day plot(d) #spectrogram which shows that the main frequency component is ~ 0.7. this translates to about a day (1.3 or so) sec.pgram(d) On Wed, Jul 30, 2008 at 7:14 AM, [EMAIL PROTECTED] wrote: Hello, I have

Re: [R] adding lines to multiple plot

2008-07-30 Thread Ben Bolker
Gary W g.watmough at 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

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

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

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

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

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

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

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

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

[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

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

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

[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] ~ . 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

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

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) Now,

[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

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

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 GPL

[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

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

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

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

[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

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) pro=protest(P1,

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

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,

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

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

[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

[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

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 dataset, but I

[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

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

[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; how can

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

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,

[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

  1   2   >