Re: [R] 1501-511 Compilation failed for file ch2inv.f - R on AIX 5.2.

2008-05-23 Thread Prof Brian Ripley
Look like you need to use a Fortran compiler set to accept Fortran 77 (which the R sources are). We don't know what 'f95' is, and this is not really an R issue so please ask your sysadmins for help. On Thu, 22 May 2008, Senthilkumaran wrote: Hi, While executing the make after successful

Re: [R] Where to download BRugs

2008-05-23 Thread Prof Brian Ripley
install.packages(BRugs) still works on Windows -- the binaries have been moved to the CRAN extras collections. If you want a direct URL http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.7/BRugs_0.4-1.zip On Thu, 22 May 2008, qqw wrote: Hi all, I tried to follow an online tutorial

[R] Est. Component Size with AIC/BIC under Gamma Distribution

2008-05-23 Thread Edward Wijaya
Dear all, I am trying to model number of samples from a given series. The series are modelled according Gamma function. In order to estimate the # samples, I use BIC/AIC with MLE (computed from dgamma function). Here is the code I have. __BEGIN__ mlogl - function( x_func, theta_func, samp) {

[R] Re : How to import package into R script

2008-05-23 Thread justin bem
Read documentation carefully ! see ?library (e.g library(tseries)) Justin BEM BP 1917 Yaoundé Tél (237) 99597295 (237) 22040246 - Message d'origine De : Edward Wijaya [EMAIL PROTECTED] À : r-help@r-project.org Envoyé le : Vendredi, 23 Mai 2008, 4h23mn 06s Objet : [R] How to

[R] New version of ecological modeling software using R

2008-05-23 Thread Bio7
Dear R users, Again with the hope that this application can also be useful for some R users, i like to announce a new Windows version of the ecological modeling software Bio7. In this release you can now easily transfer images from the well known image analysis tool ImageJ to R or create images

[R] Some problems with Sweave

2008-05-23 Thread martin . gaston
Dear R users, I'm working in a brief R-tutorial to a group of students. To make that I'm using Sweave but I've got two problems: First, I want show how R operates with the matrix type but, I write in the .rnw document the code echo=T,results=tex= matriz - matrix(vector,nrow=3,ncol=6) matriz @

Re: [R] Re : How to import package into R script

2008-05-23 Thread martin . gaston
Before call library(package) it is needed install the package. To make that, select Packages/InstallPkages... choose your nearest mirror and select the package from the list and accept. After that you can call library(tseries) for example Regards El Vie, 23 de Mayo de 2008, 9:00, justin bem

[R] Strange julian and/or strptime

2008-05-23 Thread Kåre Edvardsen
Hi r-helpers... Why do I get this strange huge jump of 36524 days when changing origin from 1969-01-01 to 1968-12-31. It should still be close to zero! This really messes up my calculations of follow-up times in my analyses. julian(strptime(010169, format = %d%m%y),origin = as.Date(1969-01-01))

Re: [R] plot 7 * 3 matrix on DIN A4 pdf

2008-05-23 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hello, I want to plot 21 scatter plots in a 7*3 matrix on a single A4 page (using mfcol). Below is an example, which -by and large- looks very close to the desired result. But I do not get it to fill out the whole page. If I specify the a different page size, the

[R] [slightly off topic] Sweave with markdown

2008-05-23 Thread baptiste Auguié
DeaR list, Has anyone tried to mix the Sweave paradigm with the Markdown[*] (and co.) syntax? Would this be hard to implement? My tiny understanding of Sweave is that one can define new drivers for the text part, while some functions that deal with the R code would not require any

Re: [R] 3D Area Diagram

2008-05-23 Thread Jim Lemon
AlGates wrote: Thank you all for your help. I will check the packages you suggested asap and maybe I find a good way to get a such a plot. The picture I posted was only an example taken from the microsoft-website, just to make clear what kind of diagram I am looking for. I know that one part is

Re: [R] Strange julian and/or strptime

2008-05-23 Thread Prof Brian Ripley
On my Linux box: strptime(010169, format = %d%m%y) [1] 1969-01-01 strptime(311268, format = %d%m%y) [1] 2068-12-31 From the help page: '%y' Year without century (00-99). If you use this on input, which century you get is system-specific. So don't! Often values

Re: [R] Where to download BRugs

2008-05-23 Thread qqw
Thanks, this works great! Thanks a lot! aldi-2 wrote: Hi, I just installed it from http://cran.wustl.edu/bin/windows/contrib/2.7/BRugs_0.4-1.zip Change 2.7 to 2.6 and you get the older version. HTH, Aldi Thanks Aldi, do you know anywhere we could download previous package of

[R] Create a data subset

2008-05-23 Thread MarySR
Hello, I want to select same data of a data frame that has information of fish catches: catch, vessel, trip, day, month, year where trip is the fishing operation number done in each day by each vessel, i.e. for each day, any vesel can have several register of catch corresponding to each

Re: [R] Some problems with Sweave

2008-05-23 Thread Charilaos Skiadas
Try results=verbatim instead of results=tex. Haris Skiadas Department of Mathematics and Computer Science Hanover College On May 23, 2008, at 4:16 AM, [EMAIL PROTECTED] wrote: Dear R users, I'm working in a brief R-tutorial to a group of students. To make that I'm using Sweave but I've

Re: [R] 3D Area Diagram

2008-05-23 Thread Patrick Connolly
On Thu, 22-May-2008 at 03:13PM -0700, AlGates wrote: | The one I rather need, is something like this (scan from a book): | http://algates.de/Scan.jpg | The only difference is the angle. That one is even stranger. There are two similar (but not parallel) lines for imports and another two

Re: [R] an unknown error message when using gamm function

2008-05-23 Thread Simon Wood
Could you try upgrading to the current version of mgcv, please? looking in the `changeLog' there were a couple of fixed problems that could be causing this. On Wednesday 21 May 2008 16:14, Lung-Chang Chien wrote: The version is mgcv 1.3-27. Last night, I dropped out many arguments from my

[R] Cleaning a time series

2008-05-23 Thread tolga . i . uzuner
Dear R Users, Was wondering if anyone can give me pointers to functionality in R that can help clean a time series ? For example, some kind of package/functionality which identifies potential errors and takes some action, such as replacement by some suitable value (carry-forward, average of

Re: [R] colorspace package does not compile on ubuntu 7.04 32 bit

2008-05-23 Thread Mikhail Spivakov
Hi again, I have managed to solve this problem myself - indeed, the trouble was with the package being incompatible with gcc v4. Here what I've done - perhaps it's going to be useful for someone else: 1) installed gcc v3.4: sudo apt-get install gcc-3.4 2) opened .../R/etc/Makeconf (changing

Re: [R] Some problems with Sweave

2008-05-23 Thread Andrew Robinson
Martin, try omitting the results=tex argument. Andrew On Fri, May 23, 2008 at 10:16:33AM +0200, [EMAIL PROTECTED] wrote: Dear R users, I'm working in a brief R-tutorial to a group of students. To make that I'm using Sweave but I've got two problems: First, I want show how R operates with

Re: [R] Pros and Cons of R

2008-05-23 Thread Neil Shephard
Installation under Gentoo is straightforward too (emerge dev-lang/R). Updating has never really been a problem. CRAN packages are rebuilt if needed when updating R, and periodically all you need to do is fire up R and use update.packages() to update any packages you've installed. Another pro

Re: [R] Strange julian and/or strptime

2008-05-23 Thread Kåre Edvardsen
I'm on an Ubuntu Linux PC, and get the same wrong result as you. I could not work out what the description of '%y' really ment, so I did not realize this was operating system specific in this sense. Anyway, I'll find a way to work around this bug. Have a nice weekend, Kare On Fri, 2008-05-23 at

Re: [R] Create a data subset

2008-05-23 Thread jim holtman
PLEASE do read the posting guide http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. You should at least provide a sample of the data that you are working with and what you might think the

Re: [R] Strange julian and/or strptime

2008-05-23 Thread jim holtman
It is not a bug. Y2K strikes again, a couple of years late. On Fri, May 23, 2008 at 7:17 AM, Kåre Edvardsen [EMAIL PROTECTED] wrote: I'm on an Ubuntu Linux PC, and get the same wrong result as you. I could not work out what the description of '%y' really ment, so I did not realize this was

[R] Evaluation conflict in combination of replicate() and rexp()as variable inside a function

2008-05-23 Thread Gerrit Eichner
Dear userRs, playing around with combinations of replicate() and random number generating functions inside a self-defined wrapper function I encounterd a puzzling behaviour. The following are intentionally simple (and rather nonsense-) examples to isolate the relevant aspects. Please, note

Re: [R] Pros and Cons of R

2008-05-23 Thread seanpor
Neil Shephard wrote: Another pro to consider is the cost, you can obtain R for free, SAS/S-Plus/Stata all have licenses of some sort that require purchasing. Neil Which has the side effect of *not* restricting how many machines are available for use or where; e.g. I was running big

Re: [R] Cleaning a time series

2008-05-23 Thread Gabor Grothendieck
The zoo package has six na.* routines for carrying values forward, etc. library(zoo) ?zoo describes them. Also see the vignettes. On Fri, May 23, 2008 at 6:55 AM, [EMAIL PROTECTED] wrote: Dear R Users, Was wondering if anyone can give me pointers to functionality in R that can help clean

Re: [R] [slightly off topic] Sweave with markdown

2008-05-23 Thread Gabor Grothendieck
You could check out the brew package: https://stat.ethz.ch/pipermail/r-packages/2007/000327.html On Fri, May 23, 2008 at 4:37 AM, baptiste Auguié [EMAIL PROTECTED] wrote: DeaR list, Has anyone tried to mix the Sweave paradigm with the Markdown[*] (and co.) syntax? Would this be hard to

[R] Advise in R- plotting graphs

2008-05-23 Thread Jason Lee
Hi I have a data table matrix,data which looks like below:- V1V2 V3 1 -6382.719 -1096.554 6998994 2 -some values- 3 -some values- 4 -some values- 5 -some values- Querying dim of data gives me 3 columns and 5 rows. And currently I want to plot V1 against V2 with the

Re: [R] Advise in R- plotting graphs

2008-05-23 Thread jim holtman
try this: x - matrix(runif(15, 100, 300), ncol=3) plot(x[,1], x[,2], type='l') points(x[5,1], x[5,2], pch=17, cex=3, col='red') # plot point #5 On Fri, May 23, 2008 at 8:50 AM, Jason Lee [EMAIL PROTECTED] wrote: Hi I have a data table matrix,data which looks like below:- V1

Re: [R] [slightly off topic] Sweave with markdown

2008-05-23 Thread Max Kuhn
On Fri, May 23, 2008 at 4:37 AM, baptiste Auguié [EMAIL PROTECTED] wrote: (although the zero compilation time is a plus), while for the latter I do not have a decent compatible editor (on a Mac, I tried Openoffice and Abiword but the fonts look like my handwriting for some obscure reason).

[R] van der Corput sequences

2008-05-23 Thread Alberto Monteiro
In package fOptions, there are functions that generate Halton sequences. The van der Corput sequence for base 2 is a particular case of the Halton sequence generated by: n - 8 # anything here... x - runif.halton(n, 1) In fact, x - runif.halton(n, dim) will generate the van der Corput sequences

Re: [R] extracting columns from a list

2008-05-23 Thread N. Lapidus
Hi Mohamed Try: lapply (NameOfYourList, function (dat, NumCol) dat[,NumCol], c(12,13)) But there must be a shorter way to write this. Nael On Fri, May 23, 2008 at 3:37 PM, mohamed nur anisah [EMAIL PROTECTED] wrote: Dear all, i have 2 lists of data with each of the list contain 14

Re: [R] extracting columns from a list

2008-05-23 Thread N. Lapidus
I forgot to answer to the last part of your question. I think what you call a list is actually an element of a list, right? If so, the command you want depends on the way you want to combine these elements. For example, the following lines will extract columns 12 and 13 of any array-like element

Re: [R] Fit a sine to data

2008-05-23 Thread Juan Manuel Barreneche
i guess nls or gnls should work JM El Viernes, 23 de Mayo de 2008 10:47, Zroutik Zroutik escribió: Dear R-users, I'd like to fit a sine function to my data. The result should have a format (and thus the formula, too) y ~ a + sin(x+b) where y and x are vectors, and a and b are (yet)

Re: [R] Pros and Cons of R

2008-05-23 Thread Roland Rau
Monica Pisica wrote: - You can save scripts, but not *.exe. If you want to contrast R with statistical packages like SPSS or Stata (and if your audience has rather a background in those than in general purpose languages), I think this is not really a problem unless I missed something

Re: [R] extracting columns from a list

2008-05-23 Thread Vincent Goulet
Le ven. 23 mai à 09:37, mohamed nur anisah a écrit : Dear all, i have 2 lists of data with each of the list contain 14 columns. No, you have one list with two elements; each is a 14-column data frame. How am i going to extract column 12 and 13 from each of the list ?? Let's call your

[R] Preparing high quality figures with tiff as end result

2008-05-23 Thread Gustaf Rydevik
Hi all, I'm currently preparing some figures that will be submitted to PloS One. In their guidelines they state that they will only accept figures in tiff or eps format, with the warning that eps figures will be converted to tiff format ( see http://www.plosone.org/static/figureGuidelines.action

Re: [R] Preparing high quality figures with tiff as end result

2008-05-23 Thread Gustaf Rydevik
On Fri, May 23, 2008 at 4:40 PM, Gustaf Rydevik [EMAIL PROTECTED] wrote: Hi all, I'm currently preparing some figures that will be submitted to PloS One. In their guidelines they state that they will only accept figures in tiff or eps format, with the warning that eps figures will be

[R] Percentages for categorical data by group

2008-05-23 Thread Economics Guy
I can think of several ways to blunt force hard code what I want but I imagine there is a command or two that can be easily combined to do this: I have a data frame with about 23000 observations. There first variable is the group to which the observation belongs (about 500 different groups). The

[R] Line Breaks and Axis breaks..

2008-05-23 Thread Megan J Bellamy
Hello all, I have two questions. One probably has a very simple answer but I have checked the FAQ, other websites and still have not found an answer. I am new to using R. My very simple question is how to do line breaks when creating an axis title (xlab) ? I have a few that are too long and

Re: [R] Fit a sine to data

2008-05-23 Thread Peter Dalgaard
Juan Manuel Barreneche wrote: i guess nls or gnls should work Yes. Notice though that this is really a linear problem because of the trigonometric identity sin(x+b) = cos(b) * sin(x) + sin(b) * cos(x) which means that you can reparametrize c*sin(x+b) (you forgot the multiplier, did you

Re: [R] Line Breaks and Axis breaks..

2008-05-23 Thread Doran, Harold
plot(1:10, xlab='Hi Megan, \n Is this what you want?') -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Megan J Bellamy Sent: Friday, May 23, 2008 10:56 AM To: r-help@r-project.org Subject: [R] Line Breaks and Axis breaks.. Hello all, I have

[R] Question about Runuran (Triangular Distribution)

2008-05-23 Thread Alberto Monteiro
The Runuran package includes a simulation for the Triangular Distribution, urtriang. But I didn't find an analogue of the q* functions (rnorm - qnorm, runif - qunif), that would invert the CDF (?uqtriang). Are there any such functions? Alberto Monteiro

Re: [R] Percentages for categorical data by group

2008-05-23 Thread Marc Schwartz
on 05/23/2008 09:51 AM Economics Guy wrote: I can think of several ways to blunt force hard code what I want but I imagine there is a command or two that can be easily combined to do this: I have a data frame with about 23000 observations. There first variable is the group to which the

Re: [R] Question about Runuran (Triangular Distribution)

2008-05-23 Thread Yasir Kaheil
Please refer to: qtriangle function in VGAM library thanks y Alberto Monteiro wrote: The Runuran package includes a simulation for the Triangular Distribution, urtriang. But I didn't find an analogue of the q* functions (rnorm - qnorm, runif - qunif), that would invert the CDF

[R] problem using embedded R

2008-05-23 Thread Lobanov, Victor S. [PRDUS]
I have a peculiar problem with an application that uses embedded R. The application comes with its own copy of R installation and it sets up correct environment path for the process during the startup. However, when I execute R code that loads an additional package (e.g. library(XYZ)), the

Re: [R] Fit a sine to data

2008-05-23 Thread apjaworski
Milan, This is a fairly standard trick. Let us generalize your equation slightly: y ~ a + c*sin(x+b) so the amplitude of the sine wave is adjustable (otherwise, you assume (or know) that the amplitude is 1). Then y ~ a + c*sin(b)*cos(x) + c*cos(b)*sin(x) or y ~ b0 + b1*x1 + b2*x2 which is

Re: [R] automation of R? running an R script at a certain time each night?

2008-05-23 Thread Felipe Carrillo
Hi: Edward, were you able to automate the process? if so, do you mind giving me a hint on how you did it? I am facing the same problem. I created a batch file which it runs fine using the task scheduler but only opens Tinn-R and R but it doesn't execute my script. My task scheduler executes

Re: [R] Percentages for categorical data by group

2008-05-23 Thread Michael Conklin
tapply(example.data$responseVar,example.data$groupVar,function(x){prop.t able(table(x))}) Michael Conklin Chief Methodologist - Advanced Analytics MarketTools, Inc. 6465 Wayzata Blvd. Suite 170 Minneapolis, MN 55426 Tel: 952.417.4719 | Mobile:612.201.8978 [EMAIL PROTECTED]

Re: [R] Preparing high quality figures with tiff as end result

2008-05-23 Thread Prof Brian Ripley
Try bitmap(): this produces PostScript and converts with ghostscript. BTW, everything (including text) in TIFF files *is* 'low-res bitmaps' by definition, so you cannot avoid this. Anti-aliasing *may* help, and the bitmap() device in R-devel supports it. On Fri, 23 May 2008, Gustaf Rydevik

[R] Boxplot width

2008-05-23 Thread Antje
Hi there, assume that you have data with different sampling like d1 - rnorm(100) d2 - rnorm(150) now, I'd like to create two boxplots in one graph but each plot located at the sampling number at the x-axis. This, I can do with at l - list(d1,d2) boxplot(l, at=c(length(d1), length(d2)),

Re: [R] automation of R? running an R script at a certain time each night?

2008-05-23 Thread Felipe Carrillo
Thank you Prof Ripley --- Prof Brian Ripley [EMAIL PROTECTED] wrote: On Fri, 23 May 2008, Felipe Carrillo wrote: Hi: Edward, were you able to automate the process? if so, do you mind giving me a hint on how you did it? I am facing the same problem. I created a batch file which

Re: [R] Preparing high quality figures with tiff as end result

2008-05-23 Thread Prof Brian Ripley
One bit you didn't tell us was that the journal asks for 600dpi TIFF files, but you also forgot to tell R that. Try bitmap(type=tifflzw, res=600) -- you probably will not need anti-aliasing. On Fri, 23 May 2008, Prof Brian Ripley wrote: Try bitmap(): this produces PostScript and converts

Re: [R] problem using embedded R

2008-05-23 Thread Prof Brian Ripley
All the usual questions: What OS? What version of R? What is meant by 'environment path'? (The environment of a process is not a path.) I have a hunch that if this is Windows = XP it was solved in R 2.7.0. Otherwise the footer of this message applies. On Fri, 23 May 2008, Lobanov, Victor

Re: [R] automation of R? running an R script at a certain time each night?

2008-05-23 Thread Tony Plate
Felipe Carrillo wrote: Hi: Edward, were you able to automate the process? if so, do you mind giving me a hint on how you did it? I am facing the same problem. I created a batch file which it runs fine using the task scheduler but only opens Tinn-R and R but it doesn't execute my script. My task

Re: [R] Alternatives to rJava and JRI

2008-05-23 Thread Tobias Verbeke
Roy Mendelssohn wrote: I do not know if this is better as we are just starting to look into it (and better is in the eye of the beholder) but look at: http://tolstoy.newcastle.edu.au/R/e3/help/07/11/4253.html in the archives. The link to the software is

Re: [R] problem using embedded R

2008-05-23 Thread Lobanov, Victor S. [PRDUS]
It is Windows 2000 still... I use latest version of R 2.7.0. There is a system environment variable PATH that lists directories where system looks for dlls when needed. Each running process can modify this environment variable, but changes apply to this process only. Instead of modifying system

Re: [R] 1501-511 Compilation failed for file ch2inv.f - R on AIX 5.2.

2008-05-23 Thread Senthilkumaran Ramaiyan
Hi, Thanks for the response.  We do have a IBM XL Fortran Enterprise Edition for AIX, V11.1 coupled with IBM XL C/C++ Enterprise Edition for AIX, V9.0 specifically on our IBM-AIX 5.2 system for the 'R' software. The 'f95' represents Fortran 95 compiler and I presume the 'R' sources can be

Re: [R] Alternatives to rJava and JRI

2008-05-23 Thread Munir, Danish
Thanks a lot everyone. These seem like very useful resources. I will check all of these out, and hopefully provide feedback so others can benefit too. -Original Message- From: Tobias Verbeke [mailto:[EMAIL PROTECTED] Sent: Friday, May 23, 2008 1:03 PM To: Roy Mendelssohn Cc: Munir,

[R] IDE

2008-05-23 Thread Alexandra Almeida
People, I'm a ubunto user and I used to write my scipts in Java Gui for R, but it is a very slow tool to run my scripts... Do you know some efficient IDE for R? Thank!!! Alexandra Almeida -- Alexandra R M de Almeida [[alternative HTML version deleted]]

[R] Colspan rowspan

2008-05-23 Thread Arshavir
Hi, What's the best way to build tables with multiple sets of column/row headings, each spanning 1 or more cols/rows, similar to the colspan/ rowspan options in HTML. I am using Sweave but couldn't find an option to do this in xtable. I can obviously build it by hand, just wanted to make

[R] GARCH-like

2008-05-23 Thread Renato Costa
I need to change the code of Garch to the FCGARCH (a non-linear multi-regime GARCH). I don't know nothing about R. I'd like to know how can I get the code of the garch in order to change it and make the fit for the FC-GARCH. Any non-linear code will be helpfull because if doesn't help in the

Re: [R] Percentages for categorical data by group

2008-05-23 Thread Economics Guy
Thanks Michael that works! Now I am having a problem getting the results into a format I can use. prop.table generates a contingency table that I tried to turn into a data frame I could graph from with as.data.frame() however the results are not a true data.frame. The rows have commas between

Re: [R] Boxplot width

2008-05-23 Thread Jorge Ivan Velez
Hi Antje, Try this: d1 - rnorm(100) d2 - rnorm(150) boxplot(c(d1,d2) ~ rep(c(1,2),c(100,150)),col=c(4:5)) HTH, Jorge On Fri, May 23, 2008 at 11:52 AM, Antje [EMAIL PROTECTED] wrote: Hi there, assume that you have data with different sampling like d1 - rnorm(100) d2 - rnorm(150) now,

Re: [R] Percentages for categorical data by group

2008-05-23 Thread Jorge Ivan Velez
Hi there, Try this: do.call(rbind,data.table) HTH, Jorge On Fri, May 23, 2008 at 1:31 PM, Economics Guy [EMAIL PROTECTED] wrote: Thanks Michael that works! Now I am having a problem getting the results into a format I can use. prop.table generates a contingency table that I tried to

Re: [R] Pros and Cons of R

2008-05-23 Thread GREGOR Brian J
Monica, here are some other Pros to consider about R: 1) IMHO, the most important reason for using R is that expressed by John Chambers as the aim of the S language: to turn ideas into software, quickly and faithfully. The broad capabilities of R facilitate the integration of data maintenance and

Re: [R] Percentages for categorical data by group

2008-05-23 Thread Economics Guy
I appreciate all the help. The trouble is that in my real data set each group does not always have an observation that choose each response. This results in some of the rows returned from prop.table() to be shorter than others so I get: Warning message: In function (..., deparse.level = 1) :

Re: [R] Pros and Cons of R

2008-05-23 Thread Antonio_Paredes
Yes, but R is not a spectator sport and that is the beauty of it. GREGOR Brian J [EMAIL PROTECTED] Sent by: [EMAIL PROTECTED] 05/23/2008 01:00 PM To r-help@r-project.org cc Subject Re: [R] Pros and Cons of R Monica, here are some other Pros to consider about R: 1) IMHO, the most

Re: [R] 3D Area Diagram

2008-05-23 Thread John Kane
--- Patrick Connolly [EMAIL PROTECTED] wrote: On Thu, 22-May-2008 at 03:13PM -0700, AlGates wrote: | The one I rather need, is something like this (scan from a book): | http://algates.de/Scan.jpg | The only difference is the angle. That one is even stranger. There are two similar

Re: [R] Pros and Cons of R

2008-05-23 Thread John Kane
I even carry a copy of R on a USB so that I can do a bit of work on a machine that does not have it installed. It's a bit slow and one would not want to do anything major with it but it's handy to show someone a quick graph or check something when far from the office. --- seanpor [EMAIL

Re: [R] Pros and Cons of R

2008-05-23 Thread Spencer Graves
GREGOR Brian J wrote: Monica, here are some other Pros to consider about R: snip 2) R facilitates documentation and replication. snip When I finally made the effort to learn how to make R packages, I experienced a substantial increase in my productivity. Before I start coding a

[R] Seeking help with trellis: log scales on xyplot

2008-05-23 Thread hobie perry
Good afternoon. The basic plot function can automatically generate log scales as follows: plot(calcium ~ soil_ph, log=y) Here is my basic model in xyplot... xyplot(calcium+magnesium ~ soil_ph|depth*region) I would like the calcium and magnesium values to be reported on a log scale.

Re: [R] Percentages for categorical data by group

2008-05-23 Thread Michael Conklin
prop.table(table(factor(x,levels=1:5))) Michael Conklin Chief Methodologist - Advanced Analytics MarketTools, Inc. 6465 Wayzata Blvd. Suite 170 Minneapolis, MN 55426 Tel: 952.417.4719 | Mobile:612.201.8978 [EMAIL PROTECTED] MarketTools(r)http://www.markettools.com This

[R] Installing R on Linux is easy. (was Re: Pros and Cons of R)

2008-05-23 Thread Patrick Connolly
On Thu, 22-May-2008 at 09:28PM +0200, Hans Ekbrand wrote: | On Thu, May 22, 2008 at 02:07:01PM -0400, R P Herrold wrote: | On Thu, 22 May 2008, Monica Pisica wrote: | | [...] | | When a new R version is in place you | cannot up-grade your old R version, you have to do a new | installation

[R] maximizing the gamma likelihood

2008-05-23 Thread markleeds
for learning purposes and also to help someone, i used roger peng's document to get the mle's of the gamma where the gamma is defined as f(y_i) = (1/gammafunction(shape)) * (scale^shape) * (y_i^(shape-1)) * exp(-scale*y_i) ( i'm defining the scale as lambda rather than 1/lambda. various

[R] exrpression(), double subscript

2008-05-23 Thread Afshartous, David
All, The first use of expression produces a double subscript while the second does not. I searched in ?plotmath and the archives and didn't find anything. I'm sure I'm missing something obvious but can't seem to find. any suggestions much appreciated. Cheers, David plot(1:10, 1:10)

[R] van der Corput sequences

2008-05-23 Thread Carlos Ungil
Alberto, I think the functions below do what you want: vanDerCorput(12,6) [1] 0.1667 0. 0.5000 0.6667 0.8333 0.0278 [7] 0.1944 0.3611 0.5278 0.6944 0.8611 0.0556 Regards, Carlos number2digits=function(n,base){ #first digit in output is

[R] dividing the data

2008-05-23 Thread kayj
hello, I have data of individuals; with their age and other infromation. Each row corresponds to a single individual. I would like to divide the data and store them into two files. one file with people less than 50 years of age and the other with people more than 50 years of age. I am new into R

[R] SVD on a matix

2008-05-23 Thread kayj
Hi All, I performed an svd on a matrix X and saved the first three column of the left singular matrix U. ( I assume that they correspond to the projection of the matrix on the first three eigen vectors that corresponds to the first three largest eigenvalues). I would like to know how much

Re: [R] dividing the data

2008-05-23 Thread Marianne Promberger
On 05/23/08 14:34, kayj wrote: hello, I have data of individuals; with their age and other infromation. Each row corresponds to a single individual. I would like to divide the data and store them into two files. one file with people less than 50 years of age and the other with people more

[R] Shared Library Error

2008-05-23 Thread what's in the name?
Can anyone help me resolve this? A part of the R function looks like this: print(is.loaded('merge_xtabs_patterns_file')) print(is.loaded('merge_xtabs_patterns_file_')) .Fortran('merge_xtabs_patterns_file_',ydim[1],ydim[2],x=as.integer(as.matrix(y)),na=as.integer(c),

Re: [R] Seeking help with trellis: log scales on xyplot

2008-05-23 Thread Deepayan Sarkar
On 5/23/08, hobie perry [EMAIL PROTECTED] wrote: Good afternoon. The basic plot function can automatically generate log scales as follows: plot(calcium ~ soil_ph, log=y) Here is my basic model in xyplot... xyplot(calcium+magnesium ~ soil_ph|depth*region) I would like the

Re: [R] exrpression(), double subscript

2008-05-23 Thread Bert Gunter
Not sure what you mean by double subscript. 11 is the number eleven and that is what is shown, as expected. 0i is the complex number 0+0i (type it at your console), and that is what is shown. Does that solve the mystery? Maybe try it with 11i as the subscript to help you see it. -- Bert Gunter

[R] Spatial Overlay - SQL - Geoprocessing

2008-05-23 Thread Paulo Cardoso
Can we by any change perform in R something similar geoprocessing or, more advanced, Spatial SQL over vectorial data? Paulo [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] How to pass variable of for loop on read table text

2008-05-23 Thread Charilaos Skiadas
?paste paste('text',y,'.txt', sep=) and you likely need data[[y]] instead of data$y. On May 23, 2008, at 9:05 PM, Jason Lee wrote: Hi, I have a couple of text and would like to automate of reading these multiple files using (namely; text1.txt, text2.txt) for(y in 3:10){

Re: [R] How to pass variable of for loop on read table text

2008-05-23 Thread markleeds
maybe this should also be a FAQ because i've seen it a lot since I've been on this list ? assign each data.frame to a component of a list, rather than to a component of a dataframe. you also need paste to get the number. listDFs-list() for ( y in 3:10 ) {

Re: [R] Colspan rowspan

2008-05-23 Thread Gabor Grothendieck
latex in Hmisc can do that with the n.cgroup= and n.rgroup= args although if your tables are getting very complicated you might be best off just generating the latex yourself since then you get all the flexibility of latex. On Fri, May 23, 2008 at 1:26 PM, Arshavir [EMAIL PROTECTED] wrote: Hi,

[R] Importing data in text file into R

2008-05-23 Thread amarjit singh sethi
Dear all, I am quite new to R; facing certain problems: Say, I have a text file( named as try): YearC1 C2 C3 C4 C5 C6 Y1 3.5 13.89.5 6.8 0.4 24.2 Y2 3.8 13.99.9 7.6 0.7 12.8 Y3 4.5 14.514.29.2 0.6

Re: [R] Importing data in text file into R

2008-05-23 Thread John Fox
Dear AJSS, The problem is that the way you've read the data, the first column of the data frame is a factor, not a numerical variable, and thus is not suitable for computing correlations. You could use the command cor(idt[,-1]) to compute correlations on all but the first column, but

Re: [R] Importing data in text file into R

2008-05-23 Thread Patrick Connolly
On Fri, 23-May-2008 at 08:40PM -0700, amarjit singh sethi wrote: | Dear all, | I am quite new to R; facing certain problems: | Say, I have a text file( named as try): | idt=read.table(df,header=T, sep=\t) | idt | Year C1 C2 C3 C4 C5 C6 | 1 Y1 3.5 13.8 9.5 6.8 0.4 24.2 | 2