[R] substituting elements in vector according to sample(unique(vector))

2007-08-29 Thread Arthur Wuster
Hello! Assuming I have a vector, such as v - c(1,2,1,2,3,3,1) This vector has three unique elements: 1, 2, and 3. unique(v) [1] 1 2 3 If I shuffle this vector of unique elements, I get something like this: sample(unique(v)) [1] 3 2 1 In the vector v I started with, I would now like to

[R] sql query over local tables

2007-08-29 Thread Jorge Cornejo Donoso
Hi i have to table with IDs in each one. I want to make a join (as in sql) by the ID. Is any way to use the RODBC package (or other) in local tables (not a access, mysql, sql, etc. ) and made the join? Thanks in advance __ R-help@stat.math.ethz.ch

[R] adding deviances as meta-analysis?

2007-08-29 Thread Peter Claessens
Dear all, Not a problem that is very specific to R, but I think that it is not only of interest to me... So I hope that someone finds the time to provide me with some clues on this probably rather basic issue. I'm performing an analysis of experimental data with a categorical response

[R] How to let the program know the end of a file is reached?

2007-08-29 Thread Yuchen Luo
Dear collegues. I am using scan( ) to read from a table (a csv file). I am wondering how to let the program know that the end of the file is reached? Your help will be highly appreciated! Best Wishes Yuchen Luo [[alternative HTML version deleted]]

[R] Age-Length key with kimura algorith

2007-08-29 Thread Jorge Cornejo Donoso
Hi, I’m looking for information to implement the kimura method (kimura Chikuni. 1987. Mixtures of empirical distributions: an interative application of the age length key) for calculation of fisheries age-length key. If someone have a manuals, methodology or examples about it I will be really

Re: [R] sql query over local tables

2007-08-29 Thread Prof Brian Ripley
On Tue, 28 Aug 2007, Jorge Cornejo Donoso wrote: Hi i have to table with IDs in each one. And what is a 'table'? If these are data frames, see ?merge. If they are tables (which are arrays in R), then still use merge() if they can be converted to data frames. I want to make a join (as in

[R] xeon processor and ATLAS

2007-08-29 Thread hui xie
hi everyone: I have a Dell Server that has a Xeon processor, and I would like to use the best ATLAS posted in the R website. I find that R has ATLAS for core2duo and P4. I am not sure which one of these two is best suited for Xeon processor, or is that neither of these two is good and I should

Re: [R] Excel

2007-08-29 Thread christian.ritter
Hmm, Excel bashing always brings joy ... but then again, it's a user's community of more than 100 million people, and if one is careful one can do quite a few interesting things with excel, in particular if it can be extended by R (using R(D)COM by T Baier and RExcel by E Neuwirth).

[R] Stratified weighted statistics

2007-08-29 Thread Vikas Rawal
I need to compute weighted descriptive statistics (mean, median, variance) for a vector of data disaggregated by a set of index variables. The weights to be used are vector in the same data frame and have to be disaggregated by the same index variables. Package Hmisc has functions for computing

Re: [R] sql query over local tables

2007-08-29 Thread Ido M. Tamir
Hi i have to table with IDs in each one. I want to make a join (as in sql) by the ID. Is any way to use the RODBC package (or other) in local tables (not a access, mysql, sql, etc. )  and made the join? With RODBC it would be normal SQL, so I guess with local tables you mean data frames? then

Re: [R] substituting elements in vector according tosample(unique(vector))

2007-08-29 Thread Dimitris Rizopoulos
try the following: v - c(1, 2, 1, 2, 3, 3, 1) x - c(3, 2, 1) fv - factor(v, levels = x) as.vector(unclass(fv)) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35,

Re: [R] sql query over local tables

2007-08-29 Thread ONKELINX, Thierry
Jorge, I'm assuming that you mean dataframes of matrices. Then you could use merge to join both dataframes into a new one. See ?merge for more detail. HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en

Re: [R] sql query over local tables

2007-08-29 Thread Uwe Ligges
Jorge Cornejo Donoso wrote: Hi i have to table with IDs in each one. I want to make a join (as in sql) by the ID. Is any way to use the RODBC package (or other) in local tables (not a access, mysql, sql, etc. ) and made the join? I guess you are looking for ?merge Uwe Ligges

Re: [R] Recoding multiple columns consistently

2007-08-29 Thread Uwe Ligges
Ron Crump wrote: Hi, I have a dataframe that contains pedigree information; that is individual, sire and dam identities as separate columns. It also has date of birth. These identifiers are not numeric, or not sequential. Obviously, an identifier can appear in one or two columns,

Re: [R] attempt at making a polygon class failed

2007-08-29 Thread Uwe Ligges
Leeds, Mark (IED) wrote: I was reading a presentation of Professor Peng's and typed the presentation code into R but I changed it to make plot.polygon a separate function instread of defining the function in SetMethod itself as he did. Is that the problem with the code below because

Re: [R] Nodes edges with similarity matrix

2007-08-29 Thread Birgit Lemcke
Hello, sorry that I interfer, but I tried to do this The SymTest matrix: A B C D E F G A 0.3 0.1 0.5 0.7 0.2 0.8 0.1 B 0.1 0.5 0.4 0.6 0.7 0.3 0.7 C 0.5 0.4 0.3 0.5 0.1 0.1 0.5 D 0.7 0.6 0.5 0.4 0.3 0.5 0.7 E 0.2 0.7 0.1 0.3 0.5 0.2 0.1 F 0.8 0.3 0.1 0.5 0.2 0.3 0.5 G 0.1 0.7

Re: [R] sql query over local tables

2007-08-29 Thread Gabor Grothendieck
I assume that by local tables you mean data frames in R. You can use the merge function in the base of R, as others have already mentioned, or if you want to use SQL syntax you can use the sqldf package. See example 4 on the sqldf home page: http://sqldf.googlecode.com On 8/28/07, Jorge

Re: [R] Age-Length key with kimura algorith

2007-08-29 Thread Alberto Murta
Hi Jorge In fact I have been working with that method. I made a simple R function to apply it, and I'm sending you the code below. If you have any questions about it just send me a message. Cheers Alberto On Tuesday 28 August 2007 8:44 pm, Jorge Cornejo Donoso wrote: Hi, I�m looking for

Re: [R] FW: How to fit an linear model withou intercept

2007-08-29 Thread Eik Vettorazzi
Hi Mark, as last comment you may also take a look at ?summary.lm where you will notice, that R reports two different R squares depending on the presence or absence of an intercept term. For comparison issues you should ensure that you use the same mathematical object. There was a thread about

Re: [R] Recoding multiple columns consistently

2007-08-29 Thread Jim Lemon
Ron Crump wrote: Hi, I have a dataframe that contains pedigree information; that is individual, sire and dam identities as separate columns. It also has date of birth. These identifiers are not numeric, or not sequential. Obviously, an identifier can appear in one or two columns,

Re: [R] Excel (off-topic, sort of)

2007-08-29 Thread Alberto Monteiro
Chris wrote: Typically, people in the R community are not used to the spreadsheet paradigm and need some time to be able to take advantage of automatic recalculation, (...) Do you know what's in my wish list? I wish spreadsheets and computer languages had gone one step further. I mean,

[R] Month end calculations

2007-08-29 Thread Shubha Vishwanath Karanth
Hi R users, Is there a function in R, which does some calculation only for the month end in a daily data?... In other words, is there a command in R, equivalent to last. function in SAS? BR, Shubha [[alternative HTML version deleted]]

Re: [R] Month end calculations

2007-08-29 Thread Henrique Dallazuanna
Hi, Perhaps if object is of the type 'ts', the command 'end' can usefully. -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 29/08/2007, Shubha Vishwanath Karanth [EMAIL PROTECTED] wrote: Hi R users, Is there a function in R, which does some calculation only for

Re: [R] Excel

2007-08-29 Thread John Kane
--- Rolf Turner [EMAIL PROTECTED] wrote: On 28/08/2007, at 7:16 PM, J Dougherty wrote: snip PS, I quit using Excel for most important work after it returned a negative variance on some data I was collecting descriptive statistics on. Those of you who have not seen it

Re: [R] Month end calculations

2007-08-29 Thread Shubha Vishwanath Karanth
But my dataset is a data frame BR, Shubha From: Henrique Dallazuanna [mailto:[EMAIL PROTECTED] Sent: Wednesday, August 29, 2007 5:12 PM To: Shubha Vishwanath Karanth Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Month end calculations Hi, Perhaps if

[R] Plot command drops part of the plot for large plots in multiple figure environment

2007-08-29 Thread Lindveld, Charles
I have run into some surprising behaviour when plotting data in a 3x2 grid: random parts of the data on my graphs aren't rendered. Not on screen, not on file, and not when I send the plot to a printer. I don't get any error messages either. When I do the individual plots (not in a grid) they

Re: [R] Month end calculations

2007-08-29 Thread Henrique Dallazuanna
Then, you can try: For example: head(df) datevalue 1 1990-01-01 4.716572 2 1990-02-01 3.418355 3 1990-03-01 1.982209 4 1990-04-01 3.974942 5 1990-05-01 6.624222 6 1990-06-01 5.729092 max(as.Date(df$date)) or if you interest is in the month only: max(format.Date(as.Date(df$date),

[R] Histogram Title SizeFont

2007-08-29 Thread livia
Hello, I would like to plot a histogram with title Return, and I would like the font for the title to be Bold and the size to be 8( as in Excel). I tried the following code, but it does not make any change. Could anyone give me some advice? hist (preturn, seq(-0.05,0.05,0.005),freq = FALSE,

Re: [R] Stratified weighted statistics

2007-08-29 Thread Frank E Harrell Jr
Vikas Rawal wrote: I need to compute weighted descriptive statistics (mean, median, variance) for a vector of data disaggregated by a set of index variables. The weights to be used are vector in the same data frame and have to be disaggregated by the same index variables. Package Hmisc has

Re: [R] Excel (off-topic, sort of)

2007-08-29 Thread Philippe Grosjean
Take a look at Mathematica or Maple. This is the kind of thing you do with these languages. Best, Philippe Grosjean Alberto Monteiro wrote: Chris wrote: Typically, people in the R community are not used to the spreadsheet paradigm and need some time to be able to take advantage of

[R] Strage result with an append/strptime combination

2007-08-29 Thread Ptit_Bleu
Hi, I keep on trying to write some small scripts in order to learn R but even with basic scripts I have problems ... I start with the name of a file which is in fact the time the file has been generated (I cannot change the format). Then I convert namefile with strptime. The problem occurs when

Re: [R] How to let the program know the end of a file is reached?

2007-08-29 Thread Ptit_Bleu
Hi Yuchen, I'm not sure my message will help you but who knows (I'm a newbie : I discovered R one month ago). To load data file, I use donnees-read.table(datafile.dat, quote=\, sep=;, dec=,, skip=18) and I don't need to tell when to stop. R knows it. Good luck, Ptit Bleu. Yuchen Luo wrote:

Re: [R] Strage result with an append/strptime combination

2007-08-29 Thread Gabor Grothendieck
Try chron: library(chron) namefile - 070707050642.dat#day-month-year-hour-minute-second.dat x - chron(substr(namefile, 1, 6), substr(namefile, 7, 12), + format = c(dmy, hms), out.format = c(m/d/y, h:m:s)) c(x, x) [1] (07/07/07 05:06:42) (07/07/07 05:06:42) See R News 4/1 Help Desk

Re: [R] xeon processor and ATLAS

2007-08-29 Thread Prof Brian Ripley
On Tue, 28 Aug 2007, hui xie wrote: hi everyone: I have a Dell Server that has a Xeon processor, and I would like to use the best ATLAS posted in the R website. I find that R has ATLAS for core2duo and P4. I am not sure which one of these two is best suited for Xeon processor, or is that

Re: [R] Strage result with an append/strptime combination

2007-08-29 Thread Ptit_Bleu
Thanks Gabor ! It works. Just one more thing : is there a possibility to remove ( and ) before I copy the data to a MySQL database. Again thank you for the tip. Ptit Bleu. Gabor Grothendieck wrote: Try chron: library(chron) namefile - 070707050642.dat

Re: [R] Interpreting the eigen value of a population matrix (2nd try)

2007-08-29 Thread Simon Blomberg
To get a confidence interval on lambda, you need to have measures of variability in the elements of the transition matrix. If you have that, you can use a parametric bootstrap to get approximate confidence intervals. I have done this, and it seems to work. Alternatively, you could calculate a

Re: [R] Interpreting the eigen value of a population matrix (2nd try)

2007-08-29 Thread Ben Bolker
Simon Blomberg s.blomberg1 at uq.edu.au writes: To get a confidence interval on lambda, you need to have measures of variability in the elements of the transition matrix. If you have that, you can use a parametric bootstrap to get approximate confidence intervals. I have done this, and it

Re: [R] Nodes edges with similarity matrix

2007-08-29 Thread Seth Falcon
Birgit Lemcke [EMAIL PROTECTED] writes: Hello, sorry that I interfer, but I tried to do this The SymTest matrix: A B C D E F G A 0.3 0.1 0.5 0.7 0.2 0.8 0.1 B 0.1 0.5 0.4 0.6 0.7 0.3 0.7 C 0.5 0.4 0.3 0.5 0.1 0.1 0.5 D 0.7 0.6 0.5 0.4 0.3 0.5 0.7 E 0.2 0.7 0.1 0.3 0.5

Re: [R] Strage result with an append/strptime combination

2007-08-29 Thread Ptit_Bleu
In french, I wouls say Chapeau bas. In english, no idea. So I just say Thanks again. Ptit Bleu (reading chron reference manual :-). Try fmt - function(x) with(month.day.year(x), sprintf(%02d/%02d/%02d %02d:%02d:%02d, month, day, year, hours(x), minutes(x),

Re: [R] Excel

2007-08-29 Thread Duncan Murdoch
[EMAIL PROTECTED] wrote: ... Here are a few additional comments related to the representation issue in .csv files: What is said about the .csv files with respect to rounding also holds for the windows clipboard but not for the office clipboard. If you format data in an excel range,

Re: [R] Strage result with an append/strptime combination

2007-08-29 Thread Gabor Grothendieck
Try fmt - function(x) with(month.day.year(x), sprintf(%02d/%02d/%02d %02d:%02d:%02d, month, day, year, hours(x), minutes(x), seconds(x))) fmt(x) On 8/29/07, Ptit_Bleu [EMAIL PROTECTED] wrote: Thanks Gabor ! It works. Just one more thing : is there a possibility to

[R] a faster and shorter way to perform calculations?

2007-08-29 Thread dxc13
This is a continuation from a previous posting of mine: The following algorithm below is what I want to accomplish: Z(xk) = Average(Yi, i belongs to Ik), where Ik contains all i such that for each j, |Xi,j - xkj│≤ 2. Here, j = 1, 2 and i corresponds to the elements in each X and/or xk data x1

Re: [R] xeon processor and ATLAS

2007-08-29 Thread Uwe Ligges
Prof Brian Ripley wrote: On Tue, 28 Aug 2007, hui xie wrote: hi everyone: I have a Dell Server that has a Xeon processor, and I would like to use the best ATLAS posted in the R website. I find that R has ATLAS for core2duo and P4. I am not sure which one of these two is best suited for

Re: [R] Histogram Title SizeFont

2007-08-29 Thread Uwe Ligges
livia wrote: Hello, I would like to plot a histogram with title Return, and I would like the font for the title to be Bold and the size to be 8( as in Excel). I tried the following code, but it does not make any change. Could anyone give me some advice? See ?par. Uwe Ligges hist

Re: [R] Plot command drops part of the plot for large plots in multiple figure environment

2007-08-29 Thread Uwe Ligges
Your examples are not reproducible. Hence hard to tell what goes wrong. There never was an R version 1.51. Please use a current version of R and read the pposoting guide. Uwe Ligges Lindveld, Charles wrote: I have run into some surprising behaviour when plotting data in a 3x2 grid: random

Re: [R] The l1ce function in lasso2: The bound and absolute.tparameters.

2007-08-29 Thread Greg Snow
This is my understanding of what is happening. 1. Standardize all the x variables to have mean 0 and variance 1 (possibly y as well). 2. Compute the unconstrained least squares regression. 3. Sum the abs values of the b's. That sum is the scaling factor. A bound of 1 means the sum above (and

Re: [R] xeon processor and ATLAS

2007-08-29 Thread hui xie
Thanks very much for all your advice. To be clear, my OS is window XP. I bought this server last year. It's Dell Precision PWS690. THe processor is Xeon(TM) CPU 3GHZ, 2G RAM. I am not sure how to check more details of processor on my computer. But I went to Dell website and from what I can

Re: [R] xeon processor and ATLAS

2007-08-29 Thread Uwe Ligges
hui xie wrote: Thanks very much for all your advice. To be clear, my OS is window XP. I bought this server last year. It's Dell Precision PWS690. THe processor is Xeon(TM) CPU 3GHZ, 2G RAM. I am not sure how to check more details of processor on my computer. But I went to Dell website and

Re: [R] Excel (off-topic, sort of)

2007-08-29 Thread François Pinard
[Alberto Monteiro] Maybe I'll write a letter to Santa Claus [there are people who write to congressman; they must have more faith than me]. :-) :-) I wish a language where I can write a = b + 10 and then when I write a = 20 the language automatically assigns b = 10. METAFONT does

[R] How to cite the the PDF user's guide for LIMMA

2007-08-29 Thread Vallejo, Roger
Dear All, How to cite the PDF user's guide for the LIMMA package? This is not about how to cite the LIMMA package. Roger Roger L. Vallejo, Ph.D. Computational Biologist Geneticist U.S. Department of Agriculture, ARS National Center for Cool Cold Water Aquaculture 11861 Leetown Road

[R] retrieve p-value from a cox.obj

2007-08-29 Thread clearsky
I have a cox.obj named obj, obj - coxph( Surv(time, status) ~ group, surv.data) now I want to retrieve the p-value from obj, so that I can run this hundreds of times and plot out the distribution of the p-value. could anyone tell me how to get p-value from obj? thanks, -- View this message in

Re: [R] How to cite the the PDF user's guide for LIMMA

2007-08-29 Thread Vallejo, Roger
Dear John, In the link that you provide, I do not see an example on referencing the LIMMA package. Perhaps I am not being clear on my question. Again, I am publishing a paper (Genetics journal) and I would like to reference the User's Guide for the Limma Package which is this PDF document: limma:

Re: [R] How to cite the the PDF user's guide for LIMMA

2007-08-29 Thread John Kane
This depends on what style manual you , your publisher, or your institution uses. One common syle (APA ) recommends this: http://library.osu.edu/sites/guides/apagd.php#sid --- Vallejo, Roger [EMAIL PROTECTED] wrote: Dear All, How to cite the PDF user's guide for the LIMMA package? This

Re: [R] retrieve p-value from a cox.obj

2007-08-29 Thread Bert Gunter
str() is your friend. It tells you about the structure of any R object, from which you can usually glean what you need to know to get what you want. It is often useful to use it on summary(object) rather than on the object, as the summary method for an (S3) classed object often contains what

Re: [R] retrieve p-value from a cox.obj

2007-08-29 Thread clearsky
I tried obj$p, obj$pvalue, but they both returned null~~ clearsky wrote: I have a cox.obj named obj, obj - coxph( Surv(time, status) ~ group, surv.data) now I want to retrieve the p-value from obj, so that I can run this hundreds of times and plot out the distribution of the p-value.

Re: [R] retrieve p-value from a cox.obj

2007-08-29 Thread Henrique Dallazuanna
See the source code of function: survival:::print.coxph -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 29/08/2007, clearsky [EMAIL PROTECTED] wrote: I tried obj$p, obj$pvalue, but they both returned null~~ clearsky wrote: I have a cox.obj named obj,

Re: [R] Excel (off-topic, sort of)

2007-08-29 Thread Albicelli, Nicholas \(Exchange\)
Except for the ability to perform circular recalculation, I believe that the closest programming analogy to a spreadsheet is a functional programming language. Check out Haskell (or LISP or Erlang) to do what you describe. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

[R] Alternating Logistic Regression

2007-08-29 Thread Rusell Barbour Ph.D.
Dear All, I am trying to see if there is a package that will allow ALR within GEE for R. I have found V J Carey of Harvard's page, but it seems out of data. Does anyone know of such a package? Sincerely Russell Barbour Ph.D. Co-Director for Statistical Analysis Yale Center for

Re: [R] retrieve p-value from a cox.obj

2007-08-29 Thread clearsky
Thanks, Bert . It is a lot help! gunter.berton wrote: str() is your friend. It tells you about the structure of any R object, from which you can usually glean what you need to know to get what you want. It is often useful to use it on summary(object) rather than on the object, as the

[R] OT: distribution of a pathological random variate

2007-08-29 Thread Horace Tso
Folks, I wonder if anything could be said about the distribution of a random variate x, where x = N(0,1)/N(0,1) Obviously x is pathological because it could be 0/0. If we exclude this point, so the set is {x/(0/0)}, does x have a well defined distribution? or does it exist a distribution

[R] Q: how to interrupt long calculation?

2007-08-29 Thread D. R. Evans
The subject says it all really: I've tried hitting control-C (multiple times), but that doesn't seem to be a reliable way to interrupt a long calculation. What is the right way to interrupt a calculation that has been proceeding for several minutes and shows no sign of finishing soon?

[R] sum(data.frame(),na.rm=TRUE) still throws an Error

2007-08-29 Thread Jacob Colvin
Hello, sum(data.frame(),na.rm=TRUE) still throws an Error, is this really the intended behavior? sum(data.frame()) Error in FUN(X[[1L]], ...) : only defined on a data frame with all numeric or complex variables sum(data.frame(),na.rm=TRUE) Error in FUN(X[[1L]], ...) : only defined on a data

[R] breaking the x-axis and having two different x-axis labels

2007-08-29 Thread Georg Ehret
Dear R community, I have two questions concerning barplots that I struggle to resolve: 1) How can I break (interrupt) the x-axis (e.g.: have it display values from -100 to -90 and 90 to 100 only)? 2) I overlay two horizontal barplots: one with negative values only and one with positive

Re: [R] OT: distribution of a pathological random variate

2007-08-29 Thread Daniel Lakeland
On Wed, Aug 29, 2007 at 10:39:17AM -0700, Horace Tso wrote: Folks, I wonder if anything could be said about the distribution of a random variate x, where x = N(0,1)/N(0,1) Obviously x is pathological because it could be 0/0. If we exclude this point, so the set is {x/(0/0)}, does x

Re: [R] OT: distribution of a pathological random variate

2007-08-29 Thread Ted Harding
On 29-Aug-07 17:39:17, Horace Tso wrote: Folks, I wonder if anything could be said about the distribution of a random variate x, where x = N(0,1)/N(0,1) Obviously x is pathological because it could be 0/0. If we exclude this point, so the set is {x/(0/0)}, does x have a well defined

Re: [R] OT: distribution of a pathological random variate

2007-08-29 Thread Charles C. Berry
On Wed, 29 Aug 2007, Horace Tso wrote: Folks, I wonder if anything could be said about the distribution of a random variate x, where x = N(0,1)/N(0,1) Instead of asking this off topic question here try googling 'gaussian ratio' Obviously x is pathological because it could be

Re: [R] Q: how to interrupt long calculation?

2007-08-29 Thread D. R. Evans
On 29/08/2007, D. R. Evans [EMAIL PROTECTED] wrote: The subject says it all really: I've tried hitting control-C (multiple times), but that doesn't seem to be a reliable way to interrupt a long calculation. What is the right way to interrupt a calculation that has been proceeding for several

Re: [R] Excel

2007-08-29 Thread Erich Neuwirth
Excel bashing can be fun but also can be dangerous because you are makeing your life harder than necessary. Statisticians meanwhile know that the numerics of statistical computation can be quite bad, therefore one should not use them. But using our (we = Thomas Baier + Erich Neuwirth) RExcel addin

Re: [R] OT: distribution of a pathological random variate

2007-08-29 Thread Horace Tso
Thank you Ted, Daniel, and Charles. I thought Cauchy distribution has to do with resonance. Didn't know this nice extension. H. Ted Harding [EMAIL PROTECTED] 8/29/2007 11:02 AM On 29-Aug-07 17:39:17, Horace Tso wrote: Folks, I wonder if anything could be said about the distribution of a

Re: [R] Excel

2007-08-29 Thread Bert Gunter
Erich: This is not a comment either for or against the use of Excel. I only wish to point out that AFAICS, Hadley Wickham's reshape package offers all the pivot table functionality and more. If I am wrong about this, please let me and everyone else know. Bert Gunter Genentech Nonclinical

Re: [R] Excel

2007-08-29 Thread Gabor Grothendieck
You would still need the interactive GUI to get to the point where its at all comparable to Excel. Using rpad you could construct such an interface although its a bit of work. Here is an example using rpad and reshape: http://www.rpad.org/Rpad/DataExplorer.Rpad On 8/29/07, Bert Gunter [EMAIL

Re: [R] Excel

2007-08-29 Thread Greg Snow
Erich Neuwirth said: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Erich Neuwirth Sent: Wednesday, August 29, 2007 12:43 PM To: r-help Subject: Re: [R] Excel Excel bashing can be fun but also can be dangerous because you are makeing your

Re: [R] Excel

2007-08-29 Thread Rolf Turner
On 30/08/2007, at 8:49 AM, Greg Snow wrote: Erich Neuwirth said: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Erich Neuwirth Sent: Wednesday, August 29, 2007 12:43 PM To: r-help Subject: Re: [R] Excel Excel bashing can be fun but also can be

[R] a new-bie question about obtaining certain value from the print out

2007-08-29 Thread kurt Zhao
Hi everyone, I am quite new to R. my command is t.test(c(1:5,7:11), y=c(1:10),alternative = c(two.sided), paired = TRUE) The output is Paired t-test data: c(1:5, 7:11) and c(1:10) t = 3, df = 9, p-value = 0.01496 alternative hypothesis: true difference in means is not equal to 0 95

Re: [R] a new-bie question about obtaining certain value from the printout

2007-08-29 Thread Greg Snow
You need to save the output in order to do something with it (the default if you don't save the output is to call the print method, so what you are seeing is the results of calling print.htest on the return value from t.test). Save the output by doing something like: out - t.test(c(1:5,7:11),

Re: [R] a new-bie question about obtaining certain value from the print out

2007-08-29 Thread Rolf Turner
On 30/08/2007, at 8:53 AM, kurt Zhao wrote: Hi everyone, I am quite new to R. my command is t.test(c(1:5,7:11), y=c(1:10),alternative = c(two.sided), paired = TRUE) The output is Paired t-test data: c(1:5, 7:11) and c(1:10) t = 3, df = 9, p-value = 0.01496 alternative

Re: [R] combining datasets by row

2007-08-29 Thread tkobayas
Rusers, I am trying to append multiple .csv files of different dimensions (but # of columns are the same for all .csv files). I do have .csv files whose names are CA1.csv ~ CA100.csv. CA means california and 1 means the first file. So what I have been doing (after googling how to append by

Re: [R] combining datasets by row

2007-08-29 Thread Rolf Turner
On 30/08/2007, at 9:41 AM, [EMAIL PROTECTED] wrote: Rusers, I am trying to append multiple .csv files of different dimensions (but # of columns are the same for all .csv files). I do have .csv files whose names are CA1.csv ~ CA100.csv. CA means california and 1 means the first file. So

Re: [R] Q: how to interrupt long calculation?

2007-08-29 Thread Paul Smith
On 8/29/07, D. R. Evans [EMAIL PROTECTED] wrote: The subject says it all really: I've tried hitting control-C (multiple times), but that doesn't seem to be a reliable way to interrupt a long calculation. What is the right way to interrupt a calculation that has been proceeding for several

Re: [R] Excel

2007-08-29 Thread Erich Neuwirth
Greg Snow wrote: Or do you trust all of your clients to know to use R(D)COM as well as how to install and use it? Do you trust your clients to be fluent enough in R to use it? For most of my clients, that is not true. For this kind of users, the following strategy works. They have their

Re: [R] Stratified weighted statistics

2007-08-29 Thread Vikas Rawal
Please look again at the help files for summarize and wtd.mean which show how to do this. This is from the summarize help file: FUN: a function of a single vector argument, used to create the statistical summaries for 'summarize'. 'FUN' may compute any number of

Re: [R] Stratified weighted statistics

2007-08-29 Thread Erich Neuwirth
In the help file for summarize in the examples section search for the following line # Get stratified weighted means Vikas Rawal wrote: Please look again at the help files for summarize and wtd.mean which show how to do this. -- Erich Neuwirth, University of Vienna Faculty of Computer

Re: [R] Excel

2007-08-29 Thread Frank E Harrell Jr
Rolf Turner wrote: On 30/08/2007, at 8:49 AM, Greg Snow wrote: Erich Neuwirth said: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Erich Neuwirth Sent: Wednesday, August 29, 2007 12:43 PM To: r-help Subject: Re: [R] Excel Excel bashing can be

[R] Linear modelling confusion.

2007-08-29 Thread Rolf Turner
I'm puzzled by a phenomenon that I have just encountered, and was hoping that someone out there might be able to give me some insight. It's not an R problem as such, but . well, there are R aspects. I have a rather messy data set, with a response, a fixed effect, and 3 random effects:

[R] piecewise linear approximation

2007-08-29 Thread Naxerova, Kamila
Dear list, I have a series of data points which I want to approximate with exactly two linear functions. I would like to choose the intervals so that the total deviation from my fitted lines is minimal. How do I best do this? Thanks! Kamila The information transmitted in this electronic

Re: [R] Excel

2007-08-29 Thread Philippe Grosjean
I think that one big problem of Excel regarding readability of the calculation done is its non linear nature: mixing data and calculation on a 2D grid could look nice at first sight, but indeed, it turns into a nightmare to understand all the calculations actually done. On the counterpart,

Re: [R] piecewise linear approximation

2007-08-29 Thread Achim Zeileis
On Wed, 29 Aug 2007, Naxerova, Kamila wrote: Dear list, I have a series of data points which I want to approximate with exactly two linear functions. I would like to choose the intervals so that the total deviation from my fitted lines is minimal. How do I best do this? From the information

Re: [R] piecewise linear approximation

2007-08-29 Thread Rolf Turner
On 30/08/2007, at 12:45 PM, Naxerova, Kamila wrote: Dear list, I have a series of data points which I want to approximate with exactly two linear functions. I would like to choose the intervals so that the total deviation from my fitted lines is minimal. How do I best do this?

[R] sneding my email

2007-08-29 Thread Luke Neraas
hello, I would like ot post questions to this list so i am sending my email Thanks Luke Neraas [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] How to multiply all dataframe rows by another dataframe's columns

2007-08-29 Thread Luke Neraas
Hello, I have two data frames, X and Y, with two columns each and different numbers of rows. # creation of data frame X Loc1.alleles - c(1,5,6,7,8) Loc1.Freq- c(0.35, 0.15, 0.05, 0.10, 0.35) Loc1 - cbind( Loc1.alleles,Loc1.Freq) X- data.frame(Loc1) #creation

Re: [R] How to multiply all dataframe rows by another dataframe's columns

2007-08-29 Thread Anders Nielsen
matrix(Y[,2]%o%X[,2],ncol=1) On Wednesday 29 August 2007 03:21 pm, Luke Neraas wrote: Hello, I have two data frames, X and Y, with two columns each and different numbers of rows. # creation of data frame X Loc1.alleles - c(1,5,6,7,8) Loc1.Freq- c(0.35, 0.15, 0.05,

Re: [R] How to multiply all dataframe rows by another dataframe's columns

2007-08-29 Thread Moshe Olshansky
Below is one way to do this: Loc1.alleles - c(1,5,6,7,8) Loc1.Freq- c(0.35, 0.15, 0.05, 0.10, 0.35) Loc1 - cbind( Loc1.alleles,Loc1.Freq) X- data.frame(Loc1) Loc2.alleles - c(1,4,6,8) Loc2.Freq - c(0.35, 0.35, 0.10, 0.20) Loc2 -

[R] filehash package: error in unserialize(con)

2007-08-29 Thread Adrian Dragulescu
This question is related to the filehash package only. I would like to use the filehash package because of it's ability to deal with large amounts of data. I've set up a filehash database on a network drive and populated it with several GB of data. Somehow, I must have corrupted the files

Re: [R] Interpreting the eigen value of a population matrix (2nd try)

2007-08-29 Thread Simon Blomberg
[I've cc'ed this to the list because I think that it may be of value to other useRs, at least for archiving purposes.] Sorry, I did that work a long time ago in XLispStat (before I switched to R). The delta method and the bootstrap method for standard errors for eigenvalues are described in

Re: [R] Month end calculations

2007-08-29 Thread Jim Porzak
Hi Shubha, By using the tautology that the end of a month is immediately followed by the first of a month, the following returns a TRUE when the date is the last day of a month IsMonthEnd - format(MyDates + 1, %d) == 01 where MyDates is a vector, or column in a data frame, typed as Date (eg

[R] customizing the color and point shape for each line drawn using lattice's xyplot

2007-08-29 Thread Gen
Description of what I am trying to do: I am using the xyplot code below to plot the variable “MeanBxg” against the variable “PercentVarExplained” for all 9 possible combinations of variables “bdg” and “bdx”. Within each of these 9 scenarios I am plotting a separate line for each of up to 9