Re: [R] Two Noobie questions

2009-01-08 Thread Prof Brian Ripley
On Thu, 8 Jan 2009, Neil Beddoe wrote: You can also use subscripts to get at things with a bit of playing around. You can, but it is easier not to fight R to do so. Much more transsparent is: cf - coef(summary(lm fit) cf[2,2] # index as a matrix. You are a) indexing a matrix as a vector

[R] library error on start up, error on update packages

2009-01-08 Thread Greg Blevins
Hello, I am unable to update packages. I am running R version 2.8.1 (2008-12-22) under VISTA. When I start R I get the following error: Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. Error in

[R] VCOV Source Code

2009-01-08 Thread Yang Wan
Dear R Help, I wonder the way to show the source code of [vcov] command. Usually, it can show the source code after input the command and enter. But for [vcov], it shows function (object, ...) UseMethod(vcov) environment: namespace:stats I appreciate for your help. Best wishes.

Re: [R] Regarding Books on R

2009-01-08 Thread Gavin Simpson
On Thu, 2009-01-08 at 10:39 +0530, Kishore wrote: Hi, I have good understanding on Econometrics and statistical techniques. However, I am new to R. What would be the best way to learn R as I would be one of the few in my team started exploring R in your team. I have got a few downloads on

Re: [R] VCOV Source Code

2009-01-08 Thread Carlos J. Gil Bellosta
Hello, You can do stats:::vcov.lm to see the source code for that particular method. In order to see which are the methods supported by vcov, write methods(vcov) Best regards, Carlos J. Gil Bellosta http://www.datanalytics.com On Wed, 2009-01-07 at 21:37 -0600, Yang Wan wrote: Dear R

Re: [R] Dataframe with unequal rows

2009-01-08 Thread Carlos J. Gil Bellosta
Hello, You are not very precise there. Do you mean that the rows in your text file do not all have the same number of separators (commas, in your case)? Best regards, Carlos J. Gil Bellosta http://www.datanalytics.com On Thu, 2009-01-08 at 04:38 -0500, rahul-a.agar...@ubs.com wrote: I have a

Re: [R] New York Times Article: Data Analysts Captivated by R's Power

2009-01-08 Thread Patrick Burns
Simon Blomberg wrote: Some people familiar with R describe it as a supercharged version of Microsoft's Excel spreadsheet software... Maybe it's my dry Australian humour, but I think this should go into the fortunes package. It is somewhat humourous even in damper London. But I think

Re: [R] Dataframe with unequal rows

2009-01-08 Thread anna freni sterrantino
Hi Rahul, try to see ?read.table you may substitute commas with a NA and then sum the NA. I hope that it can help. Cheers Anna Anna Freni Sterrantino Ph.D Student Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO.

Re: [R] library error on start up, error on update packages

2009-01-08 Thread Prof Brian Ripley
Please do read the rw-FAQ (as requested in the posting guide): it will enlighten you. On Wed, 7 Jan 2009, Greg Blevins wrote: Hello, I am unable to update packages. I am running R version 2.8.1 (2008-12-22) under VISTA. When I start R I get the following error: I presume that is

[R] VaR-Monte carlo Simulation, Historic simulation, Variance-Covariance Simulation

2009-01-08 Thread Maithili Shiva
Dear R helpers Suppose I have a portfolio of securities with exposure to Equity, Bonds and Forex (say $ 100 each). Is there any fucntion in R that will help me calculate Value at Risk (VaR) using Monte carlo Simulation , Historic simulation and Variance - Covariance Simulation. With

Re: [R] VaR-Monte carlo Simulation, Historic simulation, Variance-Covariance Simulation

2009-01-08 Thread Carlos J. Gil Bellosta
Yes, there are: replicate and quantile are your friends. You will find better support in the R-Finance list, though. Best regards, Carlos J. Gil Bellosta http://www.datanalytics.com On Thu, 2009-01-08 at 01:36 -0800, Maithili Shiva wrote: Dear R helpers Suppose I have a portfolio of

Re: [R] resolution (dpi) problem

2009-01-08 Thread Gabor Grothendieck
On Tue, Apr 29, 2008 at 7:56 AM, Ted Harding ted.hard...@manchester.ac.uk wrote: Certainly for Windows users, the suggested change could be wise, since Windows (for reasons best known to someone who should have known better) conceals the extension from the I think the reason for hiding

[R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Gundala Viswanath
Dear all, I found that printing with 'cat' is very slow. For example in my machine this snippet __BEGIN__ # I need to resolve to use this type of loop. # because using write(), I need to create a matrix which # consumes so much memory. Note that foo, bar, qux object # is already very large

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Henrik Bengtsson
On Thu, Jan 8, 2009 at 3:12 AM, Gundala Viswanath gunda...@gmail.com wrote: Dear all, I found that printing with 'cat' is very slow. For example in my machine this snippet __BEGIN__ # I need to resolve to use this type of loop. # because using write(), I need to create a matrix which #

Re: [R] troubles performing Moran.I test

2009-01-08 Thread Roger Bivand
Barbara.Spillmann at agrar.uni-giessen.de writes: dear R users, I have troubles performing Moran.I test as suggested on http://www.ats.ucla.edu/stat/r/faq/morans_i.htm Using the ape package which is written for a different application area is not necessarily a good idea - the web page

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread jim holtman
What exactly is the problem you are trying to solve. What is going to be done with the data? Is it going to be read by some other program? How much physical memory do you have on your machine? Is there paging occuring due to the size of the objects? Have you consider creating a structure with

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Gundala Viswanath
Dear Jim and Henrik, What exactly is the problem you are trying to solve. Is it going to be read by some other program? I simply want to print the data out. Surely, this data will be manipulated (with Excel or other programming languages) by other people suit to their purpose. Typically the

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 08.01.2009 14:26:32: Dear Jim and Henrik, What exactly is the problem you are trying to solve. Is it going to be read by some other program? I simply want to print the data out. Surely, this data will be manipulated (with Excel or other

Re: [R] R in the NY Times

2009-01-08 Thread Marc Schwartz
on 01/07/2009 09:47 PM Gabor Grothendieck wrote: On Wed, Jan 7, 2009 at 10:26 PM, Dirk Eddelbuettel e...@debian.org wrote: On 7 January 2009 at 18:24, Gabor Grothendieck wrote: | By running the code below we see that the: | - sum of the three seems to be rising at a constant rate | - S is

Re: [R] Dataframe with unequal rows

2009-01-08 Thread David Winsemius
I doubt that you have a dataframe with those features, since R would not allow such an event; more likely you have data in a file. If your goal is to determine the number of items, then you should definitely look at: ?count.fields count.fields(filename, sep=,) # would give you 1 + the

Re: [R] Faster Printing Alternatives to 'cat'

2009-01-08 Thread jim holtman
Here is one way of doing it. To write out 1 million rows on my system took 21 seconds. # create some data dataSize - 1e6 foo - runif(dataSize) bar - runif(dataSize) n - 1000 # number of items to write out each time output - file('/output.txt', 'w') # now split the indices into groups of

[R] Reading UTF-8 from MySQL in Windows

2009-01-08 Thread Tomas Krilavičius
Dear R-help, I am reading strings (varchar) from MySQL database and instead of some characters get garbled chars. Strings are encoded in UTF-8. I have tried to change encoding in MySQL to latin7, but it does not help. I was not able to find solution, how to get strings in proper encoding/form.

[R] NAs are not allowed in subscripted assignments

2009-01-08 Thread Rune Schjellerup Philosof
Browse[1] j - c(1,2,NA) Browse[1] j[j==1][-1] [1] NA Browse[1] j[j==1][-2] [1] 1 Browse[1] j[j==1][-2] - 2 Error during wrapup: NAs are not allowed in subscripted assignments As far as I can see, I have no NA in the lhs (not after the second subscript anyway). Besides, I have a single value on

Re: [R] R in the NY Times

2009-01-08 Thread Max Kuhn
More commentary on Slashdot: http://developers.slashdot.org/article.pl?sid=09/01/07/2316227 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

[R] ftp connections for uploading files

2009-01-08 Thread Thomas Loridan
Hi all, I would like to upload some plots I create wth R via ftp or something similar but I don t really understand which command/syntax I should use: should I go for make.socket + write.socket or try and create environment variables like frp_proxy_user and then ftp my files? how? many thanks

Re: [R] xtable-longtable question

2009-01-08 Thread Felipe Carrillo
Thanks Duncan, I'll try that with my data. --- On Wed, 1/7/09, Duncan Mackay mac...@northnet.com.au wrote: From: Duncan Mackay mac...@northnet.com.au Subject: Re: [R] xtable-longtable question To: mazatlanmex...@yahoo.com Cc: r-help@r-project.org Date: Wednesday, January 7, 2009, 1:15 PM I

Re: [R] R in the NY Times

2009-01-08 Thread Doran, Harold
The open-source mentality is invaluable, as most on this list know. That is what keeps the R evolution progressing at a pace that SAS cannot keep up with. On a side note (a very side note), I am a zealot for an exercise program called Crossfit. Crossfit has adopted the same open-source mentality

Re: [R] ftp connections for uploading files

2009-01-08 Thread Prof Brian Ripley
Try system() with curl or a decent ftp client (I don't see that package RCurl covers this, but it might despite its description only mentioning HTTP). From 'man curl' curl offers a busload of useful tricks like proxy support, user authen- tication, ftp upload, HTTP post, SSL

[R] interpolation to abscissa

2009-01-08 Thread e-letter
Readers, I have looked at various documents hosted on the web site; I couldn't find anything on interpolation. So I started r and accessed the help (help.start()). (by the way is it possible to configure r to open help in opera instead of firefox?) Initially I read the help for the akima package

Re: [R] Generating GUI for r-scripts

2009-01-08 Thread Mike Prager
Daren Tan dare...@hotmail.com wrote: I have developed some scripts that basically ask for input tab-limited format files, do some processing, and output several pictures or csv. Now I need to have some gui to wrap on top of the scripts, [...] Please advice me if there is any tools or project

Re: [R] R in the NY Times

2009-01-08 Thread Stas Kolenikov
On 1/7/09, Gabor Grothendieck ggrothendi...@gmail.com wrote: Here is the same number of messages/posts data for each of S, SAS, R: - reworked into a 3 column ts class time series - with Jan 2009 removed since its not complete - leading and trailing NA rows removed My software of choice

Re: [R] ftp connections for uploading files

2009-01-08 Thread Thomas Loridan
Many thanks for that. I am running R under windows and I am not sure what you mean by this curl? Is'nt there a way to ftp via an intrinsic R function? Thanks again Thomas 2009/1/8 Prof Brian Ripley rip...@stats.ox.ac.uk: Try system() with curl or a decent ftp client (I don't see that package

[R] how to plot histogram plot and fitted distributions on the same graph

2009-01-08 Thread Xin Shi
Dear: I am trying to plot the histogram graph for my observed data. Then plot fitted distribution on the same graph of histogram plot in R. 1. histogram plot y. 2. based on 1, plotting y1 v. x; 3. based on 1, plotting y2 v. x; 4. based on 1, plotting y3 v. x; All of

[R] Letter-based representation of pairwise comparisons

2009-01-08 Thread retama
Hi! I have been working several years with R but it's my first public question. I hope I'll be clear :) . This question is related to obtaining letter-based representation of non-parametric pairwise comparisons. I have a dataframe with this structure (but with quite more rows and cols): A B C

Re: [R] Converting data frame to symmetric matrix

2009-01-08 Thread poastpd
Thank you. It worked! Quoting Henrique Dallazuanna www...@gmail.com: Try this: reshape(d, timevar = x, idvar = y, direction = wide) On Tue, Jan 6, 2009 at 6:29 AM, poas...@umich.edu wrote: Dear Sir or Madam, I have the following data frame (which is just a toy example of my larger

[R] cosinor analysis

2009-01-08 Thread Anne Berger
Hallo, I didn´t found any facilities for Halbergs cosinor analysis in R. This analysis is well known in the Chronobiology as the least square approximation of time series using cosine function of known period (in my case of 24hours-period). I tried to write a script but crashed... Can you give

Re: [R] Numerical Integration Problems

2009-01-08 Thread Christos Argyropoulos
Hi, You may want to try the double exponential transformation on the numerator and the denominator on this one. The method is described in detail here: http://projecteuclid.org/DPubS?service=UIversion=1.0verb=Displayhandle=euclid.prims/1145474600 If you want to give it a shot outside R

[R] ftp connections for uploading files

2009-01-08 Thread Thomas Loridan
Hi all, I would like to upload some plots I create wth R via ftp or something similar but I don t really understand which command/syntax I should use: should I go for make.socket + write.socket or try and create environment variables like frp_proxy_user and then ftp my files? how? many thanks

Re: [R] Generating GUI for r-scripts

2009-01-08 Thread Bert Gunter
If you are on Windows (you didn't say!) see also ?winMenuAdd, ?winDialog ?choose.files etc. -- Bert Gunter Genentech -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Mike Prager Sent: Thursday, January 08, 2009 8:29 AM To:

[R] Why Excel is our friend

2009-01-08 Thread Stavros Macrakis
Some people familiar with R describe it as a supercharged version of Microsoft's Excel spreadsheet software... It is easy to ridicule this line from the NYT article. But this is not only a very sensible comment by a smart reporter, but also one that is good for R: It is good for R because

Re: [R] ftp connections for uploading files

2009-01-08 Thread Gábor Csárdi
On Thu, Jan 8, 2009 at 5:42 PM, Thomas Loridan thomas.lori...@gmail.com wrote: Many thanks for that. I am running R under windows and I am not sure what you mean by this curl? http://curl.haxx.se/ Is'nt there a way to ftp via an intrinsic R function? That I don't know, perhaps some else

Re: [R] Excluding data with apply

2009-01-08 Thread David Winsemius
Using indexing and putting a minus sign in front of a vector of column names that you want to exclude would be a typical approach: df - data.frame(a=LETTERS[1:4], b= rnorm(4), c=rnorm(4), d= letters[5:9]) apply(df[ , -c(a,d)], 2, sum) (Pretty sure this will run properly but don't have R

Re: [R] ftp connections for uploading files

2009-01-08 Thread Prof Brian Ripley
On Thu, 8 Jan 2009, Gábor Csárdi wrote: On Thu, Jan 8, 2009 at 5:42 PM, Thomas Loridan thomas.lori...@gmail.com wrote: Many thanks for that. I am running R under windows and I am not sure what you mean by this curl? http://curl.haxx.se/ And it comes as standard on Mac OS X and many Linux

Re: [R] Letter-based representation of pairwise comparisons

2009-01-08 Thread hadley wickham
On Thu, Jan 8, 2009 at 5:40 AM, retama retama1...@gmail.com wrote: Hi! I have been working several years with R but it's my first public question. I hope I'll be clear :) . This question is related to obtaining letter-based representation of non-parametric pairwise comparisons. I have a

[R] longtable example

2009-01-08 Thread Felipe Carrillo
Hi: Could someone try to run this example? I can't make it work. ## Not run: \begin{small} \setlongtables \begin{longtable}{ results=tex,fig=FALSE= cat(paste(c('c', rep('cc', 34/2-1), 'c'), collapse='@{\hspace{2pt}}')) @ } \hline \endhead \hline \endfoot results=tex,fig=FALSE=

Re: [R] NAs are not allowed in subscripted assignments

2009-01-08 Thread Prof Brian Ripley
On Thu, 8 Jan 2009, Rune Schjellerup Philosof wrote: Browse[1] j - c(1,2,NA) Browse[1] j[j==1][-1] [1] NA Browse[1] j[j==1][-2] [1] 1 Browse[1] j[j==1][-2] - 2 Error during wrapup: NAs are not allowed in subscripted assignments As far as I can see, I have no NA in the lhs (not after the second

Re: [R] how to plot histogram plot and fitted distributions on the same graph

2009-01-08 Thread Albyn Jones
You are plotting the histogram in the frequency scale. A quick look at the doc page for hist() would reveal the freq option: hist(x,freq=FALSE) then you can add the densities with lines() albyn Quoting Xin Shi jasonshi...@hotmail.com: Dear: I am trying to plot the histogram

Re: [R] Regarding Books on R

2009-01-08 Thread David M Smith
In addition to the other excellent examples already given, I maintain a list of resources for newcomers to R here: http://www.revolution-computing.com/community/resources.php#TOC-Resources-for-Newcomers-to-R I can also suggest this blog category (which will grow further with time):

[R] lattice question: independent per-row or per-column scaling?

2009-01-08 Thread René J.V. Bertin
Hello - and happy newyear to all of you! I've got some data that I'm plotting with bwplot, a 3x2x3 design where the observable decreases with the principle independent factor, but at different rates. I'd like to get lattice to impose not a single set of axes ranges identical for all panels, but

[R] R and Excel

2009-01-08 Thread ohri2007
Even using the VBA back of Excel to create interfaces with R would make a lot of sense. Suppose I could have access to VBA macros that import and export data into R , it would be great. The R GUI series like Rattle come even closer to Excel...so a VBA _R_ExCel package might be useful to ordinary

Re: [R] VaR-Monte carlo Simulation, Historic simulation, Variance-Covariance Simulation

2009-01-08 Thread Liviu Andronic
On 1/8/09, Maithili Shiva maithili_sh...@yahoo.com wrote: Is there any fucntion in R that will help me calculate Value at Risk (VaR) using Monte carlo Simulation , Historic simulation and Variance - Covariance Simulation. There are some un-published Crystal Ball functions for R [1], which

Re: [R] ftp connections for uploading files

2009-01-08 Thread Duncan Temple Lang
Prof Brian Ripley wrote: Try system() with curl or a decent ftp client (I don't see that package RCurl covers this, but it might despite its description only mentioning HTTP). It does support FTP, and all of the protocols that are supported in the installed libcurl, so it depends the

Re: [R] R and Excel

2009-01-08 Thread Erich Neuwirth
There is RExcel (available by downloading the CRAN package RExcelInstaller. It allows to transfer data between R and Excel, and run R code from within Excel. So you can start with your data in Excel, let R do an analysis, and transfer the results back to Excel. You can write VBA macros which do

Re: [R] interpolation to abscissa

2009-01-08 Thread Greg Snow
If you want to just linearly interpolate, then use the functions approx or approxfun from the stats package (one of those that is loaded by default). See the function TkApprox in the TeachingDemos package for an interactive way to plot the approximations with the interpolations plotted.

Re: [R] R Stacked Histogram

2009-01-08 Thread hadley wickham
# Unfortunately due to a small bug in the current version of ggplot, you need to use # jpeg with an e to get jpegs ggsave(file = ggsave_qplot_output.jpeg) # But I'd recommend using png instead - much better quality ggsave(file = ggsave_qplot_output.png) Here's a cute illustration of why I

Re: [R] Excluding data with apply

2009-01-08 Thread Stavros Macrakis
On Thu, Jan 8, 2009 at 5:52 AM, Christian Kamenik christian.kame...@giub.unibe.ch wrote: 'Apply' is a great thing for running functions on rows or columns of a matrix: X - rnorm(20, mean = 0, sd = 1) dim(X) - c(5,4) apply(X,2,sum) Is there a way to use apply for excluding rows or

[R] Convert to as.Date

2009-01-08 Thread Shruthi Jayaram
Hi, I have an vector object that looks like DA - c(1991q1, 1993q2) (first quarter of 1991 etc) and I want to convert it into a date object using as.Date(). I did this for montly data but am stumped when it comes to dealing with quarterly data and as.Date. Would anyone be able to help?

Re: [R] R and Excel

2009-01-08 Thread Ajay ohri
Hi Erich, I saw that it uses a remote server ( which can be the same machine ) to compute. Here is the question- What is the remote server is Amazon EC2 which has upscalaing and downscaling facillity for RAM and CPU... Will it work ? is there a SaaS version of this? Regards, Ajay On Thu,

Re: [R] R and Excel

2009-01-08 Thread Ajay ohri
Hi Erich, I would like to share and embed the RExcel Training video (just like youtube allows me to) . How can I do that ? Regards, Ajay www.decisionstats.com On Fri, Jan 9, 2009 at 12:20 AM, Ajay ohri ohri2...@gmail.com wrote: Hi Erich, I saw that it uses a remote server ( which can be

Re: [R] R and Excel

2009-01-08 Thread Erich Neuwirth
Can't you just set a link to our server? I do not know how to embed it. I used Camtasia to produce this video, and it does all the embedding. Ajay ohri wrote: Hi Erich, I would like to share and embed the RExcel Training video (just like youtube allows me to) . How can I do that ?

Re: [R] how to plot histogram plot and fitted distributions on the same graph

2009-01-08 Thread Greg Snow
Here is a direct quote from the help page for hist: if 'FALSE', probability densities, component 'density', are plotted (so that the histogram has a total area of one). That looks like it may apply to your problem, just figure out what the above quote is referring to. The last example on that

Re: [R] cosinor analysis

2009-01-08 Thread Charles C. Berry
On Thu, 8 Jan 2009, Anne Berger wrote: Hallo, I didn?t found any facilities for Halbergs cosinor analysis in R. This analysis is well known in the Chronobiology as the least square approximation of time series using cosine function of known period (in my case of 24hours-period). I

Re: [R] Excluding data with apply

2009-01-08 Thread David Winsemius
Well, testing would have been wise.. The last variable in the sample dataframe had 5 elements and the use negation of a character vector is not proper. df[ , -c(a,d)] Error in -c(a, d) : invalid argument to unary operator So if you limit yourself to negative indexing of numeric references

[R] F-test

2009-01-08 Thread Swanton0822
hi, i was wondering if i have been giving ANOVA table in R: Response:MPG DF Sum SqMean Sq F vaule Pr(F) Model 1 216.750216.7506.1272 0.04811* Model.Mixture4

Re: [R] R in the NY Times

2009-01-08 Thread Louis Bajuk-Yorgan
As the product manager for S+, I'd like to comment as well. I think the burgeoning interest in R demonstrates that there's demand for analytics to solve real, business-critical problems in a broad spectrum of companies and roles, and that some of the incumbent analytics offerings, in particular

Re: [R] R in the NY Times

2009-01-08 Thread Carlos J. Gil Bellosta
On Thu, 2009-01-08 at 10:42 -0600, Stas Kolenikov wrote: A really good measure for R will be the total # of the downloads of r-base for all platforms from all CRAN mirrors (and I would expect that # can be found from the servers' logs). Hello, You obviate here that many of us are downloading

Re: [R] R in the NY Times

2009-01-08 Thread ohri2007
Yes I think R as a package can really learn from SAS and SPSS in making GUI more user friendly , even at the risk of dumbing down some complexity.. also as a consultant I know that selling software requires a lot of marketing follow ups..which is why R has lagged behind in actual implementation

Re: [R] R in the NY Times

2009-01-08 Thread Rahul-A.Agarwal
I believe R as a package has everything people with little knowledge of programming can handle quite easily. Moreover even if someone has no programming knowledge can learn R without much effort. I also believe if people in corporate world start using R instead of other complex software which are

Re: [R] Convert to as.Date

2009-01-08 Thread Rowe, Brian Lee Yung (Portfolio Analytics)
You can use zoo to do this: library(zoo) DA - c(1991q1, 1993q2) as.Date(as.yearqtr(DA)) -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Shruthi Jayaram Sent: Thursday, January 08, 2009 1:46 PM To: r-help@r-project.org Subject: [R]

Re: [R] R in the NY Times

2009-01-08 Thread Marc Schwartz
on 01/08/2009 01:12 PM Andrew Choens wrote: On Thu, 2009-01-08 at 10:42 -0600, Stas Kolenikov wrote: A really good measure for R will be the total # of the downloads of r-base for all platforms from all CRAN mirrors (and I would expect that # can be found from the servers' logs). Given that it

Re: [R] R in the NY Times

2009-01-08 Thread Carlos J. Gil Bellosta
On Thu, 2009-01-08 at 13:52 -0600, Marc Schwartz wrote: Reading the posts on SAS-L since yesterday via Google RSS, where the NYT article was also posted, some have noted that SAS itself offers online support forums (http://support.sas.com/forums/index.jspa). From a quick review, it looks

Re: [R] Another newbie question

2009-01-08 Thread AllenL
Thank you all! In future I will include examples of my code to make things simpler for you. This is what I settled on: Sp.presence-Data[,14:31] ##The subset of my data set I'm interested in (the presence/absence data) Sp.presence$Species-apply(Sp.presence,1,function(x)

[R] help

2009-01-08 Thread Felipe Carrillo
Hi: I am going through some of the xtable examples and I can't make the one below work. I need to create a longtable on the fly keeping the column headers for all the pages and I thought this example could give some ideas on how to do it. I am using Sweave and xtable to create my tables and

Re: [R] Fw: Re: R Stacked Histogram

2009-01-08 Thread Jason Rupert
Thank you again for your response.   This worked great.   Quick question about the legend for qplot.  Instead of being outside the plot, is it possible to move the location of the legend to the upper left or right corner of the plot?  Could you possibly provide an example.   Thank you again

[R] 'equivalent' sign in plot title

2009-01-08 Thread jjh21
Hello, I would like to add the 'equivalent' sign (i.e., three horizontal lines, almost like an equal sign) to a plot. Does R have anything like the LaTeX command '\equiv' that I could put into a plot title? Thanks -- View this message in context:

Re: [R] 'equivalent' sign in plot title

2009-01-08 Thread Peter Alspach
Try plot(1:10, main=expression(x%==%y)) For more details see ?plotmath or demo(plotmath) HTH ... Peter Alspach -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of jjh21 Sent: Friday, 9 January 2009 10:07 a.m. To:

[R] Want to post

2009-01-08 Thread Alon Ben-Ari
Dont be so American Cut the paper work lets get going [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Alignment of image plot overlay

2009-01-08 Thread tsippel
I'm having trouble with alignment of a trend line overlayed onto an image plot. The two should be plotted on the same x-axis (time-series). However, the trend line begins about an inch into the image plot x-axis and ends about an inch off of end of the image plot. Once I have the alignment

Re: [R] Alignment of image plot overlay

2009-01-08 Thread Greg Snow
The updateusr function in the TeachingDemos package can be used to change the user coordinates of the first plot to match what you think they should before adding to the plot. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org

Re: [R] Bar Plot with Connected Points on 1 Y-Axis

2009-01-08 Thread jimdare
Cheers Marc, it works perfectly now. Thanks for your help! Marc Schwartz wrote: You can generalize the approach by using something like: ... ylim = c(0, max(DF$TACC, DF$Catch) * 1.1) ... That would allow you to use the max value of the two columns, multiplied by a fudge

[R] Creating file names.

2009-01-08 Thread greggallen
I have a function that is called from a sourced script, whenever certain conditions are met. I have this inside it: if(P 0.9)save( file = today() , myLMSmodel ); This gives me files that look like this: Thu Jan 8 14:04:43 2009 That works fine on my Mac, but it would be much more convenient

Re: [R] help

2009-01-08 Thread Duncan Mackay
Unfortunately a little long Try this (I have used Miktex 2.7 on win XP and R 2.8 for this) I have not gone into the detail of the first chunk I just created the header strings in R and copied to make the Sweave file lines With a little thought it can be done I suppose ## Sweave script

Re: [R] Creating file names.

2009-01-08 Thread Peter Alspach
Gregg You could either manipulate the character returned by today(), or try using Sys.time() instead. For example, on Windows, paste(strsplit(date(), ' ')[[1]][c(2,3,5,4)], collapse='_') gives Jan_09_2009_13:40:57 which approaches what you want. Alternatively format(Sys.time(),

[R] Graph Script/X-axis label Problem

2009-01-08 Thread jimdare
Hi, I created the script below to plot multiple graphs. It seems to work pretty well (see BYX7) except when the max value for catch or TACC is very small. In this case the x-axis labels are removed from the x-axis (see BYX10). I can't figure out why given that y co-ordinate for the text

Re: [R] Creating file names.

2009-01-08 Thread Don MacQueen
Your easiest starting point is probably file = format(Sys.time()) and then look at the online help for the strptime function to learn about formatting options. For example format( Sys.time() , '%m%d%y' ) Use paste() and some formatting of Sys.time() to construct the form you are

[R] The ATT researchers and the New York Times

2009-01-08 Thread Robert Wilkins
Is anyone in the leadership of the R-project going to contact the New York Times and clarify that the article gave remarkably short shrift to the people who designed the user interface for R, to a large extent ATT researchers from an earlier generation? It would be the appropriate thing to do.

[R] ATT Researchers and the New York Times

2009-01-08 Thread Robert Wilkins
Is anyone in the leadership of the R-project going to contact the New York Times and clarify that the article gave remarkably short shrift to the people who designed the user interface for R, to a large extent ATT researchers from an earlier generation? It would be the appropriate thing to do.

[R] Using apply for two datasets

2009-01-08 Thread Christos Argyropoulos
It depends on how the data are arranged ## x-matrix(c(1,2,3,2,8,2,4,5,6),nrow=3) y-matrix(c(10,2,13,0,8,4,4.2,5.2,6.2),nrow=3) q-mapply(t.test,as.data.frame(x),as.data.frame(y)) q ## The ith column of q contain the results of applying t.test

Re: [R] Creating file names.

2009-01-08 Thread Tony Breyal
Hi, if i understood correctly, you probably want the something like: ## R format(Sys.time(), LMSMOD_%b_%d_%y_%H%M) (works on windows, don't know about Mac) Hope that helps a little, Tony Breyal On 8 Jan, 23:47, greggal...@gmail.com wrote: I have a function that is called from a sourced

[R] survey statistics, rate/proportions with standard errors

2009-01-08 Thread Robert Wilkins
what does R have to compare with , say , proc surveymeans, estimate survey means/proportions with standard errors, using Taylor methods? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] What does R have for age-adjusted survey analysis?

2009-01-08 Thread Robert Wilkins
A procedure that , after adjusting for sampling weights, also explicitly does an age adjustment to conform with an age distribution of an older census? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] Extracting File Basename without Extension

2009-01-08 Thread Gundala Viswanath
Dear all, The basename() function returns the extension also: myfile - path1/path2/myoutput.txt basename(myfile) [1] myoutput.txt Is there any other function where it just returns plain base: myoutput i.e. without 'txt' - Gundala Viswanath Jakarta - Indonesia

Re: [R] ATT Researchers and the New York Times

2009-01-08 Thread stephen sefick
It has worked wonders for me over the last years. It clunkily and reliably delivers messages to my inbox that are either someone elses question about R or an answer to one of my questions. clunkily yours Stephen Sefick On Thu, Jan 8, 2009 at 8:58 PM, Robert Wilkins irishhac...@gmail.com wrote:

Re: [R] Extracting File Basename without Extension

2009-01-08 Thread jim holtman
You can use 'sub' to get rid of the extensions: sub(^([^.]*).*, \\1, 'filename.extension') [1] filename sub(^([^.]*).*, \\1, 'filename.extension.and.more') [1] filename sub(^([^.]*).*, \\1, 'filename without extension') [1] filename without extension On Thu, Jan 8, 2009 at 9:10 PM, Gundala

Re: [R] Problem with aovlmer.fnc in languageR

2009-01-08 Thread Chuck Clifton
I have been corresponding with Harald Baayen about this problem, and it appears that something is broken with the mcmcsamp function that pvals.fnc calls. Running pvals.fnc on an lmer object actually changes the contents of the object (get a summary before and after running pvals.fnc and you'll

Re: [R] survey statistics, rate/proportions with standard errors

2009-01-08 Thread Stas Kolenikov
On 1/8/09, Robert Wilkins iwriteco...@gmail.com wrote: what does R have to compare with , say , proc surveymeans, estimate survey means/proportions with standard errors, using Taylor methods? survey package by Thomas Lumley has pretty much everything you would need. Way more than SAS, at any

[R] Mac woes

2009-01-08 Thread Stuart Vyse
I am new to R but have spent a little time messing around with it. I am running version 2.8.1 under Leopard OS 10.5 on a power PC and also on my laptop, which is Intel. Unless, I am wrong, this version is just a bit buggy. 1. In the R Commander, there is no direct entry of data. The Edit

[R] VCOV Source Code

2009-01-08 Thread Murray Cooper
Mostly I am interested in using R for statistics. I am also interested in being able to look at source code. I hope to be able to write extensions. I tried the suggestion below but was unable to access vcov.lm. methods(vcov) [1] vcov.Arima* vcov.glm* vcov.lm*vcov.mlm* vcov.nls*

Re: [R] VCOV Source Code

2009-01-08 Thread Henrik Bengtsson
stats::vcov.lm Error: 'vcov.lm' is not an exported object from 'namespace:stats' stats:::vcov.lm function (object, ...) { so - summary.lm(object, corr = FALSE) so$sigma^2 * so$cov.unscaled } environment: namespace:stats /Henrik On Thu, Jan 8, 2009 at 1:44 AM, Carlos J. Gil Bellosta

[R] Programming Question (setting ylim generally)

2009-01-08 Thread stephen sefick
library(StreamMetabolism) day - (structure(c(10.08, 10.08, 10.05, 10.03, 10, 9.98, 9.96, 9.95, 9.95, 9.96, 9.96, 9.98, 10.01, 10.05, 10.06, 10.09, 10.11, 10.11, 10.13, 10.13, 10.15, 10.15, 10.13, 10.14, 10.11, 10.13, 10.14, 10.13, 10.12, 10.13, 10.14, 10.16, 10.18, 10.19, 10.23, 10.27, 10.33,

  1   2   >